Solution Curves -- with Animation!

> IC1;

> P1:= DEplot( MODEL, VAR, DOMAIN, [ IC1 ], stepsize=0.01,

> scene=[ t, x ], arrows=NONE, linecolor=BLUE ):

> P2:= DEplot( MODEL, VAR, DOMAIN, [ IC1 ], stepsize=0.01,

> scene=[ t, y ], arrows=NONE, linecolor=GREEN ):

> P3:= DEplot( MODEL, VAR, DOMAIN, [ IC1 ], stepsize=0.01,

> scene=[ t, z ], arrows=NONE, linecolor=RED ):

> LS1 := display( [ P1, P2, P3 ] , title = `Solution Curves for Lorenz System (dy=0)` ): # combined solution curves

> LS1;

>

> IC2;

> P1:= DEplot( MODEL, VAR, DOMAIN, [ IC2 ], stepsize=0.01,

> scene=[ t, x ], arrows=NONE, linecolor=BLUE ):

> P2:= DEplot( MODEL, VAR, DOMAIN, [ IC2 ], stepsize=0.01,

> scene=[ t, y ], arrows=NONE, linecolor=GREEN ):

> P3:= DEplot( MODEL, VAR, DOMAIN, [ IC2 ], stepsize=0.01,

> scene=[ t, z ], arrows=NONE, linecolor=RED ):

> LS2 := display( [ P1, P2, P3 ] , title = `Solution Curves for Lorenz System (dy=0.1)` ): # combined solution curves

> display( LS2 );

>

> IC3;

> P1:= DEplot( MODEL, VAR, DOMAIN, [ IC3 ], stepsize=0.01,

> scene=[ t, x ], arrows=NONE, linecolor=BLUE ):

> P2:= DEplot( MODEL, VAR, DOMAIN, [ IC3 ], stepsize=0.01,

> scene=[ t, y ], arrows=NONE, linecolor=GREEN ):

> P3:= DEplot( MODEL, VAR, DOMAIN, [ IC3 ], stepsize=0.01,

> scene=[ t, z ], arrows=NONE, linecolor=RED ):

> LS3 := display( [ P1, P2, P3 ] , title = `Solution Curves for Lorenz System (dy=0.01)` ): # combined solution curves

> LS3;

>

> IC4;

> P1:= DEplot( MODEL, VAR, DOMAIN, [ IC4 ], stepsize=0.01,

> scene=[ t, x ], arrows=NONE, linecolor=BLUE ):

> P2:= DEplot( MODEL, VAR, DOMAIN, [ IC4 ], stepsize=0.01,

> scene=[ t, y ], arrows=NONE, linecolor=GREEN ):

> P3:= DEplot( MODEL, VAR, DOMAIN, [ IC4 ], stepsize=0.01,

> scene=[ t, z ], arrows=NONE, linecolor=RED ):

> LS4 := display( [ P1, P2, P3 ] , title = `Solution Curves for Lorenz System (dy=0.001)` ): # combined solution curves

> LS4;

>

> display( [ seq( LS.i, i=1..4 ) ], insequence=true );

> display( [ seq( display( LS.i, view=[ 0.. 5,-25..50] ), i=1..4 ) ], insequence=true );

> display( [ seq( display( LS.i, view=[10..15,-25..50] ), i=1..4 ) ], insequence=true );

> display( [ seq( display( LS.i, view=[20..25,-25..50] ), i=1..4 ) ], insequence=true );

>