fabiohono@gmail.com

Main

Curriculum

Hobbies

Photos

Links

Courses

Files

Simulations

Versions

Lorenz Attractor



The "Butterfly Effect" is often ascribed to the American meteorologist, Edward Lorenz, who developed in 1963 the system of 3 differential equations below:

diff(X(t),t) = sigma*(Y-X)

diff(Y(t),t) = -X*Z+r*X-Y

diff(Z(t),t) = X*Y-b*Z

For further details on the meaning of each equation please check Weisstein (refered below).

In this article Lorenz wrote that "if the theory were correct, one flap of a seagull's wings would be enough to alter the course of the weather forever". In 1972, in his speech delivered at the American Association for the Advancement of Science in Washington, D.C the seagull "became" a butterfly an left a question:

"Does the Flap of a Butterfly’s Wings in Brazil set off a Tornado in Texas?"

Initialization

> restart;

> with( DEtools ):

> with( plots ):

> with( linalg ):

> with( PDEtools ):
Warning, the protected names norm and trace have been redefined and unprotected

Warning, the name changecoords has been redefined

Warning, the name adjoint has been redefined

Warning, the protected names norm and trace have been redefined and unprotected

Parameters

> sigma:=5:

> r:=15:

> b:=1:

Model Dynamics

> eqs:={diff(X(t),t)=-sigma*X(t)+sigma*Y(t),
diff(Y(t),t)=-X(t)*Z(t)+r*X(t)-Y(t), diff(Z(t),t)=X(t)*Y(t)-b*Z(t)}:

> iniset:=[[X(0)=-4,Y(0)=-2.5,Z(0)=15]]:

> DEplot3d(eqs,{X(t),Y(t),Z(t)},t=0..180,X=-8..8,Y=-11..11,Z=4..24,iniset,stepsize=0.1, method=rkf45,linecolor=navy,thickness=1,axes=FRAME);

Powered by: JavaView

obs: If the graph doesn't come out press "reload" in your web browser. You must have enabled/instaled the Java Virtual Machine.
To rotate, rescale or pick a vertex click the right-button and move your mouse.
The picture may take a few seconds to come out.

References:

- Lorenz, E. N. "Deterministic Nonperiodic Flow." J. Atmos. Sci. 20, 130-141, 1963.

- Weisstein, E. W. "Lorenz Attractor." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/LorenzAttractor.html