{VERSION 6 0 "IBM INTEL LINUX" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 14 "11assgnA09.mws" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 57 "For MATH 333 Assi gnment #11 3.3: 23, 24 -- 11assgnA09.mws" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 81 "\"Find the general real-valued sol ution of each ODE. Plot solution curves in the " }{XPPEDIT 18 0 "ty; " "6#%#tyG" }{TEXT -1 11 "-plane for " }}{PARA 0 "" 0 "" {TEXT -1 6 "- 1 <= " }{XPPEDIT 18 0 "t;" "6#%\"tG" }{TEXT -1 13 " <= 5, where " } {XPPEDIT 18 0 "y(0) = 1;" "6#/-%\"yG6#\"\"!\"\"\"" }{TEXT -1 107 ", y' (0) = -6, -3, 0, 3, 6. Plot the corresponding orbits in a rectangle that shows their main features.\"" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "with(DEtools, DEplot):" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 26 "Common Initial-Value Stuff" }}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 15 "TRG := t=-1..5:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "yZRO := 1: ypZROlist := [-6, -3, 0, 3, 6]:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "PTS := map( x -> [1,x], ypZR Olist);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "IClist := map(x \+ -> [y(0)=1, D(y)(0) = x], ypZROlist);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "ICyvLIST := map(x -> [y(0)=1, v(0) = x], ypZROlist); " }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 7 "3.3: 23" }}{EXCHG {PARA 0 " " 0 "" {TEXT -1 25 "The five solution curves:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "DE23 := (D@@2)(y)(t) + D(y)(t) = 0: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "DE := DE23;" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 40 "IVPlist := map(q -> \{DE,op(q)\}, IClist):" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "solnLIST := map( u -> dsolv e(u, y(t)), IVPlist ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "p lotLIST := map(z -> plot(rhs(z), TRG) ,solnLIST): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "plots[display](plotLIST);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 11 "The orbits:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "RHS := solve( DE, (D@@2)(y)(t)):" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 33 "RHS := subs(D(y)(t) = v(t), RHS):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "DEA := D(y)(t) = v(t);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "DEB := D(v)(t) = RHS;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "SYS := \{DEA, DEB\};" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "DEplot(SYS, [y(t), v(t)], TR G, ICyvLIST, stepsize=0.05);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 7 "3.3: 24" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "The five solution curves:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "DE24 : = (D@@2)(y)(t) + 2*D(y)(t) + 65*y(t)= 0: " }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 11 "DE := DE24;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "IVPlist := map(q -> \{DE,op(q)\}, IClist):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "solnLIST := map( u -> dsolve(u, y(t)), IVPl ist ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "plotLIST := map(z -> plot(rhs(z), TRG) ,solnLIST): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "plots[display](plotLIST);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 11 "The orbits:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "RHS := solve( DE, (D@@2)(y)(t)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "RHS := subs(D(y)(t) = v(t), RHS):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "DEA := D(y)(t) = v(t);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 21 "DEB := D(v)(t) = RHS;" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 18 "SYS := \{DEA, DEB\};" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "DEplot(SYS, [y(t), v(t)], TRG, ICyvLIST, stepsiz e=0.05);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}}{MARK "6 " 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }