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