A complete Lab manual for M.Sc physics Computational Physics Practical

Prepared by Prof. Purushothaman sir. Excellent work

http://www.scribd.com/doc/79536004/mscpractical2011-12

Logistic map- Bifurcation diagram

from math import*
from pylab import*
def logistic(x,r):
return r*x*(1-x)
#x0=input("Enter intial value x0= ")
rp=[]
x=[]
r=2
while r<=6:
x0=0.35
for i in range(250):
x0=logistic(x0,r)
for i in range(251,600):
x0=logistic(x0,r)
x.append(x0)
rp.append(r)
r+=0.01
clf()
scatter(rp,x,marker='o',s=0.005)
show()