Friday, December 13, 2019

Plotting Recurrence Tables 3


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

No comments:

Post a Comment