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); |
| > | 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); |
| > | display(Nerk, orientation=[0,90]); The "orbit" in xy-space |
| > | display(Nerk, orientation=[90,90]); The ty-space "component curve" |
| > | display(Nerk, orientation=[-90,0]); The tx-space "component curve" |
| > |
| > |
| > |
| > |