Art of Recurrence Tables 3

xxxxxxxxxx

 

1

def gen(a,b,max=Infinity):

2

    n=1; xyn=[0.,0.]

3

    while n<=max:

4

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

5

        xyn=[(1+.01*a)*abs(xn)-(1+.00001*b)*yn+1,xn]

6

N,A,B=3000,7,5; XY=gen(A,B,N); XYN=[el for el in XY]

7

cols=[colormaps.autumn(1.*i/N)[:3] for i in [1..N]]

8

(polygon([[-5,-5],[-5,15],[15,15],[15,-5],[-5,-5]],

9

         color='silver',alpha=.1)+\

10

sum([point([XYN[i]],size=1,color=cols[i]) 

11

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

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

---

Messages