Friday, December 13, 2019

Plotting Recurrence Tables 5


import numpy,pylab
def gen(a,b,n):
    i=1; xyi=[0.,0.]
    while i<=n:
        yield xyi; i+=1; xi,yi=xyi[0],xyi[1]
        xyi=[(1+.01*a)*abs(xi)-(1+.00001*b)*yi+1,xi]
for n in [50000,100000]: 
    a,b=7,5; xy=gen(a,b,n); xyn=numpy.array([el for el in xy])      
    pylab.figure(figsize=(10,10)); ax=pylab.gca(); ax.set_facecolor('ghostwhite')
    pylab.scatter(xyn[:,0],xyn[:,1],s=.05,c=[numpy.random.random(3)])
    pylab.grid(c='slategray',alpha=.4); pylab.show()

No comments:

Post a Comment