Art of Recurrence Tables 2

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-(.2+.01*a)*xn^2+yn,(.9991+.0001*b)*xn]

6

N=3000; XY=gen(-11,-19987,N); XYN=[el for el in XY]

7

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

8

(polygon([[-.25,-1.25],[-.25,.25],[1.25,.25],

9

          [1.25,-1.25],[-.25,-1.25]],

10

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

11

sum([point([XYN[i]],size=2,color=cools[i]) 

12

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

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

---

Messages