Visualization of Recurrence Tables 2

xxxxxxxxxx

 

1

def gen(max=Infinity):

2

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

3

    while n<=max:

4

        yield xyn; n+=1

5

        xn,yn=xyn[0],xyn[1]

6

        xyn=[.7*xn+yn,-.8+xn^2]

7

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

8

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

9

      for i in [1..N]]

10

p=polygon([[-1.5,-1],[-1.5,1.5],[1,1.5],

11

           [1,-1],[-1.5,-1]],

12

          color='silver',alpha=.05)

13

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

14

       for i in [1..N-1]])

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

---

Messages