漸化式となる関数からnumpy配列を高速に生成する方法を探しています。

ある実数値 a0 と 関数 f(x) が定まっているとして、以下のコードをなるたけ実行効率をよくするにはどうすればいいでしょうか。

a = a0
L = [a]
for i in range(100):
    a = f(a)
    L += [a]

何卒、ご教示ください。