from lorenz import Lorenz import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D A = Lorenz([-1,1,0]) u1 = A.solve(50,.01) B = Lorenz([-1.001,1.001,0.001]) u2 = B.solve(50,.01) print('---------------------------------------------------------------') print('',u1,'\n',u2) print('---------------------------------------------------------------') print(u1[0][0],u2[0][0]) print(u1[-1][0],u2[-1][0]) fig = plt.figure() ax = fig.gca(projection='3d') ax.set_xlabel("X Axis") ax.set_ylabel("Y Axis") ax.set_zlabel("Z Axis") ax.set_title("Lorenz Attractor") plt.plot(u1[:,0],u1[:,1],u1[:,2]) plt.plot(u2[:,0],u2[:,1],u2[:,2]) plt.show()
from lorenz import Lorenz sigma = 10 rho = 28 beta = 8/3 L1 = Lorenz([-1,1,0],sigma,rho,beta) u1 = L1.solve(50,.01) L2 = Lorenz([-1.001,1.001,.001],sigma,rho,beta) u2 = L2.solve(50,.01) print(u1[0,0],u2[0,0]) print(u1[-1,0],u2[-1,0]) print(L1.df(2)) print(L1.test(u1))