Art of Recurrence Tables 4

xxxxxxxxxx

 

1

def gen(a,max=Infinity):

2

    n=1; z0=1+1*I; xyn=[(a*z0).real(),(a*z0).imag()]

3

    while n<=max:

4

        yield xyn; n+=1; xn,yn=xyn[0],xyn[1]; zn=xn+yn*I

5

        zn=a+(1-a)*zn; xyn=[zn.real(),zn.imag()]

6

N,A=3000,.01+.15*I; XY=gen(A,N); XYN=[el for el in XY]

7

sum([point([XYN[i]],size=1,color=hue(1-i/N)) 

8

     for i in [1..N-1]]).show(axes=False)

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages