Phase Portrait -- with Animation!
> IC1 := [ x(0) = 0, y(0) = 1, z(0)=0 ] ;
> PHASEplot1 := DEplot3d( MODEL, VAR, DOMAIN, [ IC1 ], WINDOW,
> arrows=NONE, stepsize=0.01,
> scene=[ x, y, z ], linecolor=BLUE,
> title = `Phase Portrait for Lorenz System (dy=0)` ):
>
> IC2 := [ x(0) = 0, y(0) = 1.1, z(0)=0 ] ;
> PHASEplot2 := DEplot3d( MODEL, VAR, DOMAIN, [ IC2 ], WINDOW,
> arrows=NONE, stepsize=0.01,
> scene=[ x, y, z ], linecolor=BLUE,
> title = `Phase Portrait for Lorenz System (dy=0.1)` ):
>
> IC3 := [ x(0) = 0, y(0) = 1.01, z(0)=0 ] ;
> PHASEplot3 := DEplot3d( MODEL, VAR, DOMAIN, [ IC3 ], WINDOW,
> arrows=NONE, stepsize=0.01,
> scene=[ x, y, z ], linecolor=BLUE,
> title = `Phase Portrait for Lorenz System (dy=0.01)` ):
>
> IC4 := [ x(0) = 0, y(0) = 1.001, z(0)=0 ] ;
> PHASEplot4 := DEplot3d( MODEL, VAR, DOMAIN, [ IC4 ], WINDOW,
> arrows=NONE, stepsize=0.01,
> scene=[ x, y, z ], linecolor=BLUE,
> title = `Phase Portrait for Lorenz System (dy=0.001)` ):
>
> display( [ seq( PHASEplot.i, i=1..4 ) ], insequence=true ) ;
>