《数学物理趣谈》一书中提到的混沌效应,参见书中P115-P117.
# -*- coding: utf-8 -*-
"""
Created on Tue Nov 20 12:41:40 2018
@author: yujin.wang
"""
import numpy as np
import matplotlib.pyplot as plt
y=[]
xx=[]
t = np.linspace(0.01,4,num=100)
for k in t:
x = [0.0000001]
z=[]
for iter in range(300):
x.append(k*x[-1]-k*x[-1]*x[-1])
z.append(k)
y.extend(x[-100:])
xx.extend(z[-100:])
plt.scatter(xx,y)
plt.show()