LORENZ - Random Differential Equation Behaviour

Walas (1991) has shown that the following set of ordinary differential equations gives random chaotic behaviour. 

The values of the parameters b, r and s and the initial values for xi, X2 and X3 in LORENZ are identical to those employed by Walas. 

As noted by Walas, the solution exhibits random behaviour and the solution jumps erratically between two critical point regions, as shown very clearly by plotting the variables as a phase-plane plot. 

CONSTANT TPIN=30, CINT=0.5 1 SIM INTERACT RESET GOTO 1 INITIAL 

Enjoy playing with this system using different parameter values, including CINT, and the differing numerical integration routines. For further study of the phenomenon read the references of Walas (1991) and Denn (1987). 

---