DEplot3d_928.mws

DEplot3d_928.mws
Wed Sep 28 14:14:56 MDT 2005

>    restart;

####################################################

>    with(plots, display):

>    with(DEtools, DEplot3d):

###############################################################

 

Here the DE is 3x'' + 5x' - 2x = 0 and y = x'

:

>    DE1 := D(x)(t) = y(t);

DE1 := D(x)(t) = y(t)

>    DE2 := D(y)(t) = 2*x(t) - 5*y(t);

DE2 := D(y)(t) = 2*x(t)-5*y(t)

>    SYS := {DE1, DE2}:  VARS := {x(t), y(t)}:  

>    IC := [[ x(0) = -1, y(0)= 1/2]]:

>    TRG := -5..5:  XRG := -2..2:  YRG := TRG:

>    OPT := scene=[t, x(t), y(t)]:

>    Nerk := DEplot3d(SYS, VARS, t=TRG, IC, x=XRG,  y=YRG, OPT):

>    display(Nerk);

[Maple Plot]

>    display(Nerk, orientation=[0,90]);   The "orbit" in xy-space

[Maple Plot]

>    display(Nerk, orientation=[90,90]);   The ty-space "component curve"

[Maple Plot]

>    display(Nerk, orientation=[-90,0]);   The tx-space "component curve"

[Maple Plot]

>   

>   

>   

>