ロジスティック写像.py

import matplotlib.pyplot as plt
import random

a1 = 1.0
a2 = 4.0

xs = []
a = []
for i in range(100000):
	an = a1 + (a2 - a1) * float(i) / 100000
	print i
	x = random.random()
	for j in range(1000):
        	x = an * x * (1 - x)
	a.append(an)
	xs.append(x)

plt.plot(a,xs,",")
plt.show()