Art of Recurrence Tables

xxxxxxxxxx

 

1

def gen(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*xn^2+yn,.9991*xn]

6

N=3000; XY=gen(N); XYN=[el for el in XY]

7

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

8

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

9

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

10

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

11

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

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

---

Messages