December 2017
Intermediate to advanced
386 pages
10h 42m
English
This is the main code of the solution:
import numpyfrom numpy import linspaceimport scipyfrom scipy.integrate import odeintimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dsigma=10.0b=8/3.0r=28.0f = lambda x,t: [sigma*(x[1]-x[0]), r*x[0]-x[1]-x[0]*x[2], x[0]*x[1]-b*x[2]]
This is the final result:

This precisely shows the chaotic behavior of the solutions. Let's observe the fluctuations of the vertical temperature in detail, along with the fluctuation of the horizontal temperature against the vertical. Issue the following commands:
plt.rcParams['figure.figsize'] = (10.0, 5.0)plt.subplot(121); plt.plot(t,Z) ...
Read now
Unlock full access