{VERSION 6 0 "IBM INTEL NT" "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 } {CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 260 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 270 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 271 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{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 "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Ti mes" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 } {PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 55 "ENGI 3424 Example 2.5.2 N umeric solution of the ODE " }{XPPEDIT 18 0 "x^2;" "6#*$)%\"xG\"\"# \"\"\"" }{TEXT -1 1 " " }{XPPEDIT 18 0 "d^2*y/(dx^2);" "6#*(%\"dG\"\"# %\"yG\"\"\"*$)%#dxGF%F'!\"\"" }{TEXT -1 2 " " }{XPPEDIT 18 0 "-6*x;" "6#,$*&\"\"'\"\"\"%\"xGF&!\"\"" }{TEXT -1 1 " " }{XPPEDIT 18 0 "dy/dx+ 12*y;" "6#,&*&%#dyG\"\"\"%#dxG!\"\"F&*&\"#7F&%\"yGF&F&" }{TEXT -1 14 " = 24 with " }{TEXT 260 2 " y" }{TEXT -1 10 "(1) = 2, " }{TEXT 261 2 "y'" }{TEXT -1 7 "(1) = 1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(DEtools):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 140 " DEplot(x^2*diff(y(x),x$2) - 6*x*diff(y(x),x) + 12*y(x) - 24, y(x),\nx= 0.001..2,[[y(1)=2,D(y)(1)=1]],\ny=0..5,stepsize=.05, linecolour=[blue] );" }}{PARA 13 "" 1 "" {GLPLOT2D 404 404 404 {PLOTDATA 2 "6%-%'CURVESG 6M7$7$$\"\"\"!\"$$\"3p(p!4)z_++#!#<7$$\"#^F*$\"3iBYrdk-+?F-7$F.7$$\"$, \"F*$\"3'\\=0+UF\"**>F-7$F47$$\"$^\"F*$\"3Y3F!4$\\5(*>F-7$F:7$$\"$,#F* $\"3I2myu&HN*>F-7$F@7$$\"$^#F*$\"3MNE7[(o\"))>F-7$FF7$$\"$,$F*$\"3+(=` ._Z4)>F-7$FL7$$\"$^$F*$\"3if]\\HG%>(>F-7$FR7$$\"$,%F*$\"33A2/:AQh>F-7$ FX7$$\"$^%F*$\"3uS_@^Mk\\>F-7$Fhn7$$\"$,&F*$\"3;(QixTas$>F-7$F^o7$$\"$ ^&F*$\"3^SRQZI*[#>F-7$Fdo7$$\"$,'F*$\"3cIu)=I(Q8>F-7$Fjo7$$\"$^'F*$\"3 @>$*Qg^r.>F-7$F`p7$$\"$,(F*$\"3Sy618Y+(*=F-7$Ffp7$$\"$^(F*$\"3+J$>wlLX *=F-7$F\\q7$$\"$,)F*$\"3])zdsHIx*=F-7$Fbq7$$\"$^)F*$\"3<'f!QRD<3>F-7$F hq7$$\"$,*F*$\"3s&e'z%R)eF>F-7$F^r7$$\"$^*F*$\"37HB>xe&y&>F-7$Fdr7$$\" %,5F*$\"3w**4+.I+,?F-7$Fjr7$$\"%^5F*$\"3I55w>w?f?F-7$F`s7$$\"%,6F*$\"3 _s,YzyzM@F-7$Ffs7$$\"%^6F*$\"3E^t>[ BF-7$Fbt7$$\"%^7F*$\"3Qh:9,ET\"\\#F-7$Fht7$$\"%,8F*$\"3A&>uWbDGm#F-7$F ^u7$$\"%^8F*$\"3Ms)QKB9b'GF-7$Fdu7$$\"%,9F*$\"3W7+rgmq-JF-7$Fju7$$\"%^ 9F*$\"35)[)Qh3yxLF-7$F`v7$$\"%,:F*$\"3+S&)*4'[E%p$F-7$Ffv7$$\"%^:F*$\" 3%))*\\['oOe0%F-7$F\\w7$$\"%,;F*$\"3aw2,mVKmWF-7$Fbw7$$\"%^;F*$\"3UdY& *GfqH\\F-7$Fhw7$$\"%, " 0 "" {MPLTEXT 1 0 198 "DEplot(x^2*diff(y(x),x$2) - 6*x*diff(y(x),x) + 12*y(x) - 24, y(x),\nx=0.001..3,[[y(1)=2,D(y)(1)=1 ], [y(1)=0.1,D(y)(1)=1], [y(1)=2,D(y)(1)=0.01]],\ny=0..10,stepsize=.05 , linecolour=[blue, green, red]);" }}{PARA 13 "" 1 "" {GLPLOT2D 404 404 404 {PLOTDATA 2 "6'-%'CURVESG6[o7$7$$\"\"\"!\"$$\"3p(p!4)z_++#!#<7 $$\"#^F*$\"3iBYrdk-+?F-7$F.7$$\"$,\"F*$\"3'\\=0+UF\"**>F-7$F47$$\"$^\" F*$\"3Y3F!4$\\5(*>F-7$F:7$$\"$,#F*$\"3I2myu&HN*>F-7$F@7$$\"$^#F*$\"3MN E7[(o\"))>F-7$FF7$$\"$,$F*$\"3+(=`._Z4)>F-7$FL7$$\"$^$F*$\"3if]\\HG%>( >F-7$FR7$$\"$,%F*$\"33A2/:AQh>F-7$FX7$$\"$^%F*$\"3uS_@^Mk\\>F-7$Fhn7$$ \"$,&F*$\"3;(QixTas$>F-7$F^o7$$\"$^&F*$\"3^SRQZI*[#>F-7$Fdo7$$\"$,'F*$ \"3cIu)=I(Q8>F-7$Fjo7$$\"$^'F*$\"3@>$*Qg^r.>F-7$F`p7$$\"$,(F*$\"3Sy618 Y+(*=F-7$Ffp7$$\"$^(F*$\"3+J$>wlLX*=F-7$F\\q7$$\"$,)F*$\"3])zdsHIx*=F- 7$Fbq7$$\"$^)F*$\"3<'f!QRD<3>F-7$Fhq7$$\"$,*F*$\"3s&e'z%R)eF>F-7$F^r7$ $\"$^*F*$\"37HB>xe&y&>F-7$Fdr7$$\"%,5F*$\"3w**4+.I+,?F-7$Fjr7$$\"%^5F* $\"3I55w>w?f?F-7$F`s7$$\"%,6F*$\"3_s,YzyzM@F-7$Ffs7$$\"%^6F*$\"3E^t>[BF-7$Fbt7$$\"%^7F*$\"3Qh:9,ET\"\\#F- 7$Fht7$$\"%,8F*$\"3A&>uWbDGm#F-7$F^u7$$\"%^8F*$\"3Ms)QKB9b'GF-7$Fdu7$$ \"%,9F*$\"3W7+rgmq-JF-7$Fju7$$\"%^9F*$\"35)[)Qh3yxLF-7$F`v7$$\"%,:F*$ \"3+S&)*4'[E%p$F-7$Ffv7$$\"%^:F*$\"3%))*\\['oOe0%F-7$F\\w7$$\"%,;F*$\" 3aw2,mVKmWF-7$Fbw7$$\"%^;F*$\"3UdY&*GfqH\\F-7$Fhw7$$\"%,Q'z'GG\"=.'F-7$Fdx7$$\"%,=F*$\"3Gr5GIUCzmF-7$ Fjx7$$\"%^=F*$\"3#QQQ\\k\")pR(F-7$F`y7$$\"%,>F*$\"3-@W33Jv*=)F-7$Ffy7$ $\"%^>F*$\"3qd!fjlOC1*F-7$F\\z7$$\"%,?F*$\"3KE3wKg+-5!#;7$Fbz7$%*undef inedGFjz7$FizFizF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[ [lF[[lF[[lF[[lF[[l-%'COLOURG6&%$RGBG$\"\"!Fa[lF`[l$\"*++++\"!\")-%&STY LEG6#%%LINEG-%*THICKNESSG6#\"\"$-F$6[o7$7$F($\"3A*oBYTx/+#F-7$F/$\"3P& *e\"RMX-+#F-7$Fc\\l7$F5$\"3))RfD&e6B*>F-7$Fg\\l7$F;$\"3cV[)y3AT(>F-7$F []l7$FA$\"3)*)yJ&\\zDT>F-7$F_]l7$FG$\"3)R,Z=Z82*=F-7$Fc]l7$FM$\"3$R6<, Qb0#=F-7$Fg]l7$FS$\"390+kq=()H,M+Xh(=;F-7$F_^l7$Fin $\"3a/'*\\(pH$)[\"F-7$Fc^l7$F_o$\"3g)Q.?7)oS8F-7$Fg^l7$Feo$\"359Fq'p`* y6F-7$F[_l7$F[p$\"3uCaQZ'[s+\"F-7$F__l7$Fap$\"3;wS8xM+2$)!#=7$Fc_l7$Fg p$\"3%Qs$HOETalFf_l7$Fh_l7$F]q$\"39\"o97A*3')[Ff_l7$F\\`l7$Fcq$\"3qfE9 G*eMQ$Ff_l7$F``l7$Fiq$\"3heZpmx*z8#Ff_l7$Fd`l7$F_r$\"3W%*4/))>B^7Ff_l7 $Fh`l7$Fer$\"3WN>]y.QZ$)!#>7$F\\al7$F[s$\"3\"*Hp1#=W,,\"Ff_l7$Faal7$Fa s$\"3_(=AyZ;\"4>Ff_l7$Feal7$Fgs$\"3?0mi$Q!RtOFf_l7$Fial7$F]t$\"3;P#o_* *RZX'Ff_l7$F]bl7$Fct$\"3HlVM6#*\\T5F-7$Fabl7$Fit$\"3G59&*fBgs:F-7$Febl 7$F_u$\"3M2*Hv8wpD#F-7$Fibl7$Feu$\"3x$Q$GU#=Q6$F-7$F]cl7$F[v$\"3eA>Ty8 LjTF-7$Facl7$Fav$\"36`y3cKsEaF-7$Fecl7$Fgv$\"3-gHV\"f1i#pF-7$Ficl7$F]w $\"3U3w\\1\"**\\o)F-7$F]dl7$Fcw$\"3yxT#HNKF2\"Fgz7$FadlFizF[[lF[[lF[[l F[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[[lF[ [lF[[lF[[lF[[lF[[lF[[lF[[l-F][l6&F_[lF`[lFb[lF`[lFe[lFi[l-F$6[o7$7$F($ \"3W24)z_+++#F-7$F/$\"3TYrdk-++?F-7$F]el7$F5$\"3%=0+UF\"****>F-7$Fael7 $F;$\"3-F!4$\\5(***>F-7$Feel7$FA$\"3Amyu&HN***>F-7$Fiel7$FG$\"3SE7[(o \"))**>F-7$F]fl7$FM$\"3-KN?v%4)**>F-7$Fafl7$FS$\"3k]\\HG%>(**>F-7$Fefl 7$FY$\"3=2/:AQh**>F-7$Fifl7$Fin$\"3U_@^Mk\\**>F-7$F]gl7$F_o$\"3yBwF-7$Fagl7$Feo$\"3ARQZI*[#**>F-7$Fegl7$F[p$\"37u)=I(Q8**>F-7$Figl7$F ap$\"3;$*Qg^r.**>F-7$F]hl7$Fgp$\"3o618Y+(*)*>F-7$Fahl7$F]q$\"35$>wlLX* )*>F-7$Fehl7$Fcq$\"3#zdsHIx*)*>F-7$Fihl7$Fiq$\"3-1QRD<3**>F-7$F]il7$F_ r$\"3%e'z%R)eF**>F-7$Fail7$Fer$\"3EB>xe&y&**>F-7$Feil7$F[s$\"355+.I+,+ ?F-7$Fiil7$Fas$\"3:5w>w?f+?F-7$F]jl7$Fgs$\"3!=g%zyzM,?F-7$Fajl7$F]t$\" 3&HFQy^-B+#F-7$Fejl7$Fct$\"3Lj>^t>[.?F-7$Fijl7$Fit$\"3R:9,ET\"\\+#F-7$ F][m7$F_u$\"3*=uWbDGm+#F-7$Fa[m7$Feu$\"3i)QKB9b'3?F-7$Fe[m7$F[v$\"3=+r gmq-6?F-7$Fi[m7$Fav$\"3*[)Qh3yx8?F-7$F]\\m7$Fgv$\"3^&)*4'[E%p,#F-7$Fa \\m7$F]w$\"3N]['oOe0-#F-7$Fe\\m7$Fcw$\"3(z5gOCjY-#F-7$Fi\\m7$Fiw$\"3=Z &*GfqHH?F-7$F]]m7$F_x$\"3'R4fS4,X.#F-7$Fa]m7$Fex$\"3E(z'GG\"=./#F-7$Fe ]m7$F[y$\"3m6GIUCzY?F-7$Fi]m7$Fay$\"3O&Q\\k\")pR0#F-7$F]^m7$Fgy$\"3iX3 3Jv*=1#F-7$Fa^m7$F]z$\"3%GfjlOC12#F-7$Fe^m7$Fcz$\"39&3wKg+-3#F-7$Fi^m7 $$\"%^?F*$\"33O)3;-y14#F-7$F]_m7$$\"%,@F*$\"3t-W'>!*4@5#F-7$Fc_m7$$\"% ^@F*$\"3)GQdZ-^X6#F-7$Fi_m7$$\"%,AF*$\"3v1WTqw0G@F-7$F_`m7$$\"%^AF*$\" 3*>8u$>woU@F-7$Fe`m7$$\"%,BF*$\"37Ts3_,]e@F-7$F[am7$$\"%^BF*$\"31Ok,\\ gbv@F-7$Faam7$$\"%,CF*$\"3ULkj!f#F-7$Fgam7$$\"%^CF*$\"3\\jRVd&[O@#F -7$F]bm7$$\"%,DF*$\"3!ew2*HU\"[B#F-7$Fcbm7$$\"%^DF*$\"3#yeo&)Q\"[dAF-7 $Fibm7$$\"%,EF*$\"37!=RRJ= " 0 "" {MPLTEXT 1 0 61 "ode1 := x^2*diff(y(x),x$2) - 6*x*diff(y(x),x) + 12*y(x) - 24;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode1G,**&)%\"xG\"\"#\"\"\"-%%diffG6$-%\"yG6# F(-%\"$G6$F(F)F*F**(\"\"'F*F(F*-F,6$F.F(F*!\"\"*&\"#7F*F.F*F*\"#CF8" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "ic1 := y(1) = 2, D(y)(1) = 1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ic1G6$/-%\"yG6#\"\"\"\"\"#/- -%\"DG6#F(F)F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "dsol1 := \+ dsolve(\{ode1, ic1\}, numeric, range=0.01..3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&dsol1Gf*6#%(x_rkf45G6)%$resG%%dataG%%varsG%)solnproc G%)outpointG%&ndsolG%\"iG6#%inCopyright~(c)~2000~by~Waterloo~Maple~Inc .~All~rights~reserved.G6\"C*>%8_EnvDSNumericSaveDigitsG%'DigitsG>F6\"# 9@%/%-_EnvInFsolveG%%trueG>8(-&%&evalfG6#F56#9$>F>-FAFC>8%`6$%$GetG%$S etGb6#%+thismoduleG6#%%DataGF26$%$GetG%$SetGF2F2F2F2%,_m146497388G6#%% DataG>8&-_FHFS6#Q(sysvarsF2>8'-Fen6#Q/soln_procedureF2@$4-%%typeG6$F>. %(numericGC$@,-%'memberG6$FD7+Q&startF2.%&startGQ'methodF2.%'methodGQ% leftF2Q&rightF2Q)leftdataF2Q*rightdataF2Q+enginedataF2C$>8$-Fin6#-%(co nvertG6$FD.%'stringG@&-F`o6$Fgp.%&arrayGO-%%evalG6$Fgp\"\"\"0FgpQ)proc nameF2OFgp-Fgo6$FD7&Q%lastF2.%%lastGQ(initialF2.%(initialGC$>FgpFhp@$- F`o6$Fgp.%%listGO7#-%$seqG6$/&FY6#,&8*FhqFhqFhq&FgpFcs/Fes;\"\"!\"\"#3 -F`o6$F>%\"=G-Fgo6$-%$lhsG6#F>7$FbrFcrC$@%-F`o6$-%$rhsGFctFjr>Fgp-Fin6 #/Fbr7$Fis-%#opG6#FitYQKinitial~values~must~be~specified~in~a~listF2@$ FhrOF]s/F>Q.solnprocedureF2O-Ffq6#Fin/F>FgnOFY@%09!%(unknownGO-.F`vFCC $>8)-%(pointtoG6#&-Fen6#Q0soln_proceduresF26#FisO-.FgvFCZ%C$>Fgp-FinFc tF]sF2YF2F2F2F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "dsol1(2); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%/%\"xG$\"\"#\"\"!/-%\"yG6#F%$\"3 &REUN(3++5!#;/-%%diffG6$F*F%$\"3\">rauU-++#F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "dsol1(0.5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7 %/%\"xG$\"\"&!\"\"/-%\"yG6#F%$\"3*Qs*)\\%**\\P>!# " 0 " " {MPLTEXT 1 0 199 "phaseportrait([D(x)(t) = y(t), \nD(y)(t) = (24-12* x(t)+6*t*y(t))/t^2], [x(t),y(t)], t=0.01..2, \n[[x(1)=2, y(1)=1]], x=1 ..3, y=0..6, \nstepsize=.05, linecolour=[blue], \ntitle=`ENGI 3424 Exa mple 2.5.2`);" }}{PARA 13 "" 1 "" {GLPLOT2D 404 404 404 {PLOTDATA 2 "6 &-%'CURVESG6M7$7$%*undefinedGF(F'F&F&F&F&F&F&F&F&F&F&F&F&7$F'7$$\"37Du )[H2i*=!#<$!3_w2/$GGb1)!#>7$F*7$$\"3=F%*p=ok%*=F-$\"3z&R-J6:.J#F07$F27 $$\"3A)\\bs*o-**=F-$\"3?G8$o9BYd\"!#=7$F87$$\"3)*4-UTD&4\">F-$\"3CSWp0 8AaKF=7$F?7$$\"3o#*[ck(y@$>F-$\"3cj!>'>b#)*H&F=7$FE7$$\"3m\"RDBe5Y'>F- $\"3!\\?3&HbVTxF=7$FK7$$\"3\"=Eu*4II5?F-$\"3]:5O+S!41\"F-7$FQ7$$\"3u@T /\")3Yr?F-$\"3)4Ka_.lKR\"F-7$FW7$$\"3')>8fa$R/:#F-$\"3'GjI1`EUx\"F-7$F gn7$$\"35-@fZMu\\AF-$\"36E'4#*[)y1AF-7$F]o7$$\"3CjTLy\"G?P#F-$\"3]yp6= 4&Rp#F-7$Fco7$$\"3ZwzRm&)4?DF-$\"3E?R&y#QrQKF-7$Fio7$$\"3!=bTEj4pp#F-$ \"3i@T`Js2WQF-7$F_p7$$\"3z[z=*Rmb!HF-$\"3he?JW6/8XF-7$Fep7$$\"3_,]T*)Q K\\JF-$\"3!eG]He0'[_F-7$F[q7$$\"3iSo%z7(eJMF-$\"35j')[l0x`gF-7$FaqF'F& F&F&F&F&F&F&F&F&F&-%'COLOURG6&%$RGBG$\"\"!F\\rF[r$\"*++++\"!\")-%&STYL EG6#%%LINEG-%*THICKNESSG6#\"\"$-%%VIEWG6$;$\"+++++!*!#5$\"+++++J!\"*;$ !+++++IF^s$\"+++++jFas-%&TITLEG6#%8ENGI~3424~Example~2.5.2G-%+AXESLABE LSG6$Q\"x6\"Q\"yF_t" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 40 "Exact \+ solution of this Cauchy-Euler ODE:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "dsolve(ode1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-% \"yG6#%\"xG,(*&)F'\"\"$\"\"\"%$_C2GF,F,*&)F'\"\"%F,%$_C1GF,F,\"\"#F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "dsolve(\{ode1, ic1\});" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG,(*$)F'\"\"$\"\"\"!\"\" *$)F'\"\"%F,F,\"\"#F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}}{MARK "16" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }