Visualization of Recurrence Tables

xxxxxxxxxx

 

1

def gen(max=Infinity):

2

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

3

    while n<=max:

4

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

5

        xyn=[xn+.684*(xn-xn^2+yn),

6

             yn+.684*(yn-yn^2+xn)]

7

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

8

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

9

       for i in [1..N]]

10

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

11

          color='steelblue',alpha=.1)

12

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

13

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

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

---

Messages