{VERSION 5 0 "SUN SPARC SOLARIS" "5.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 Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 " " 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 } {PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 75 " \+ Math 2280: Computer Project 1" }}{PARA 0 "" 0 "" {TEXT -1 70 " E uler's Method" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 106 "In this project you will compute some numerical solution s to initial value problems using Euler's method. " }}{PARA 0 "" 0 "" {TEXT -1 60 "First let's restart Maple and load the appropriate packag es." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "r estart: with (DEtools): with (plots): with (linalg):" }}{PARA 7 "" 1 " " {TEXT -1 50 "Warning, the name changecoords has been redefined\n" }} {PARA 7 "" 1 "" {TEXT -1 45 "Warning, the name adjoint has been redefi ned\n" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names nor m and trace have been redefined and unprotected\n" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 116 "We will consider the differential equation y' = 2 y. Here's a plot of the slope field for this differential equation." } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "f := y \+ -> 2*y;\nfieldplot ([1,f(y)], x = -2..2, y = -2..2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fGf*6#%\"yG6\"6$%)operatorG%&arrowGF(,$*&\"\"#\" \"\"9$F/F/F(F(F(" }}{PARA 13 "" 1 "" {GLPLOT2D 381 381 381 {PLOTDATA 2 "6$-%'CURVESG6\\dl7%7$$!3wI4zl+`D?!#<$!3'pFOotzy*=F*7$$!3Cp!4U$*pW(> F*$!3-BP;j-7-@F*7$$!3Tbd=6Z**f>F*$!3(>_XcP=:0#F*7%7$$!33d9+]P+:=F*F+7$ $!3d&f>%=O%Rw\"F*F07$$!3u\"G'R&Ro%\\@MuZ/;F*F+7$$!3!>7 IEI8MF]oF+7$$!3k)zK*ySf-HF]oF07$$!3Uh'*p[=%yv#F]oF57%7$$!3Jxm &oBKzI\"F]oF+7$$!3N;1Q5#4L(z!#>F07$$!3AW#\\!3pyDlFiqF57%7$$\"3:%f!Q5#4 L(zFiqF+7$$\"34vm&oBKzI\"F]oF07$$\"3I7)*3nWo_9F]oF57%7$$\"3U'zK*ySf-HF ]oF+7$$\"3a69v%R&>8MF]oF07$$\"3u[X)\\iZzb$F]oF57%7$$\"39Lv#oBdy+&F]oF+ 7$$\"3E[hk_&e%=bF]oF07$$\"3Z&GzGy5Km&F]oF57%7$$\"3()pAs%R?J6(F]oF+7$$ \"3)\\)3a5ASxSRZoxF]oF57%7$$\"3f1qh_NQ=#*F]oF+7$$\"3s@ cVo[)*G(*F]oF07$$\"3$*e(o')4PP()*F]oF57%7$$\"3Lu60rYOK6F*F+7$$\"3$e.LE ![U$=\"F*F07$$\"3m\\jlD+!z>\"F*F57%7$$\"3+[1%o)4*GM\"F*F+7$$\"3]4DU=6& RR\"F*F07$$\"3NBeWTjU39F*F57%7$$\"3n@,j-tT`:F*F+7$$\"3>$)>@MuZ/;F*F07$ $\"3-(HNsl_*=;F*F57%7$$\"3N&f>%=O%Rw\"F*F+7$$\"3'oX,+v.]\"=F*F07$$\"3p qZ-t*y%H=F*F57%7$$\"3-p!4U$*pW(>F*F+7$$\"3wI4zl+`D?F*F07$$\"3OWU\"))G0 +/#F*F57%7$F($!3=eCVAH5)p\"F*7$F.$!3X%f))fWW3)=F*7$$!3Dp@`'3#oi>F*$!37 !R4\"Hs*[$=F*7%7$F9Fcx7$FA`\\XH;a\"F*F[y7%7$FMFcx7$FPFfx7$$!3A[P;RJ5J8F*F[y7%7$FWFcx7$FZFfx7 $$!3cuUPBod?6F*F[y7%7$F[oFcx7$F_oFfx7$$!3$)3![e2005*F]oF[y7%7$FfoFcx7$ FioFfx7$$!35sK&z\">C&*pF]oF[y7%7$F`pFcx7$FcpFfx7$$!3$f`e+wy**)[F]oF[y7 %7$FjpFcx7$F]qFfx7$$!3?*zj@g:Zy#F]oF[y7%7$FdqFcx7$FgqFfx7$$!3o?1pUW_%z 'FiqF[y7%7$F_rFcx7$FbrFfx7$$\"3![nDOr5eU\"F]oF[y7%7$FirFcx7$F\\sFfx7$$ \"3'4T?:(Q2JNF]oF[y7%7$FcsFcx7$FfsFfx7$$\"3C[^THqLOcF]oF[y7%7$F]tFcx7$ F`tFfx7$$\"3'\\))4t=+;u(F]oF[y7%7$FgtFcx7$FjtFfx7$$\"3q@Y?XL'o%)*F]oF[ y7%7$FauFcx7$FduFfx7$$\"3#e$*4.l7_>\"F*F[y7%7$F[vFcx7$F^vFfx7$$\"3_4%* 4m*QdS\"F*F[y7%7$FevFcx7$FhvFfx7$$\"3>$))))=Glih\"F*F[y7%7$F_wFcx7$Fbw Ffx7$$\"3'oNywf\"zE=F*F[y7%7$FiwFcx7$F\\xFfx7$$\"3`IyY8zJP?F*F[y7%7$F( $!3kR'G!3hK)\\\"F*7$F.$!3nlM\")G'o&f;F*7$$!33$ey=Yp`'>F*$!31eKd#3w#=;F *7%7$F9Fa`l7$FkK7\"F*Fi`l7%7$F[oFa`l7$F_oFd`l7$$!3Lx>EjqFiqFi`l7%7$F_rFa`l7$FbrFd`l7$$\"3-P:;gp$*)R\"F]oFi`l7 %7$FirFa`l7$F\\sFd`l7$$\"3=ti0=,?/NF]oFi`l7%7$FcsFa`l7$FfsFd`l7$$\"3\" *45&fFj%4cF]oFi`l7%7$F]tFa`l7$F`tFd`l7$$\"3jYd%QVEZr(F]oFi`l7%7$FgtFa` l7$FjtFd`l7$$\"3O$[S)*F]oFi`l7%7$FauFa`l7$FduFd`l7$$\"3+AN'\\FDD >\"F*Fi`l7%7$F[vFa`l7$F^vFd`l7$$\"3o&*Hv!f^IS\"F*Fi`l7%7$FevFa`l7$FhvF d`l7$$\"3NpCa1zd8;F*Fi`l7%7$F_wFa`l7$FbwFd`l7$$\"3-V>LAU5C=F*Fi`l7%7$F iwFa`l7$F\\xFd`l7$$\"39<97Q0jM?F*Fi`l7%7$F($!33@[i$H\\&)H\"F*7$F.$!3)o LQ;\"GHQ9F*7$$!3\"p*\\APo0o>F*$!3*f7Pg$\\l,9F*7%7$F9F_hl7$F^**>L(FiqFghl7%7 $F_rF_hl7$FbrFbhl7$$\"3A*R(p1K1s8F]oFghl7%7$FirF_hl7$F\\sFbhl7$$\"3&f8 #fkjKxMF]oFghl7%7$FcsF_hl7$FfsFbhl7$$\"3oso[A&*e#e&F]oFghl7%7$F]tF_hl7 $F`tFbhl7$$\"3S4;Q!o_yo(F]oFghl7%7$FgtF_hl7$FjtFbhl7$$\"37YjFQe6$z*F]o Fghl7%7$FauF_hl7$FduFbhl7$$\"3;3rh**y$)*=\"F*Fghl7%7$F[vF_hl7$F^vFbhl7 $$\"3&=e1a@k.S\"F*Fghl7%7$FevF_hl7$FhvFbhl7$$\"3_bg>J0*3h\"F*Fghl7%7$F _wF_hl7$FbwFbhl7$$\"3>Hb)p%oT@=F*Fghl7%7$FiwF_hl7$F\\xFbhl7$$\"3I.]xiJ %>.#F*Fghl7%7$F($!3_-5AzCx)4\"F*7$F.$!353KY%*p,<7F*7$$!3`59d7Uuq>F*$!3 o$*4]*yL]=\"F*7%7$F9F]`m7$Fy'*y@gk\\e:.#3;F*Fe`m7%7$F_wF]`m7$FbwF``m7$$\"3e:\"R;ZH(== F*Fe`m7%7$FiwF]`m7$F\\xF``m7$$\"3[*eGuyb#H?F*Fe`m7%7$F($!3[Q=<[m&**)*) F]o7$F.$!3E%z!)Gx6u&**F]o7$$!3OCy\"zeJM(>F*$!3O<'['Hk7%o*F]o7%7$F9F[hm 7$FQUEj*=h(>F*$!3x(H(Gk\\\"z^(F]o7%7$F9Fi_n7$F(poJ%Q([B*F]oFa`n7%7$FfoFi_n7$FioF\\`n7$$!3YgRF&o5'H rF]oFa`n7%7$F`pFi_n7$FcpF\\`n7$$!3uB#zt_ZV-&F]oFa`n7%7$FjpFi_n7$F]qF\\ `n7$$!3d([%[pV3>HF]oFa`n7%7$FdqFi_n7$FgqF\\`n7$$!3$H](*e67#Q\")FiqFa`n 7%7$F_rFi_n7$FbrF\\`n7$$\"3U')\\IY>W\"H\"F]oFa`n7%7$FirFi_n7$F\\sF\\`n 7$$\"3gA(*>/^q'R$F]oFa`n7%7$FcsFi_n7$FfsF\\`n7$$\"3))fW4i#o>]&F]oFa`n7 %7$F]tFi_n7$F`tF\\`n7$$\"3h'>*)*>9B2wF]oFa`n7%7$FgtFi_n7$FjtF\\`n7$$\" 3LLR)yd%\\7(*F]oFa`n7%7$FauFi_n7$FduF\\`n7$$\"3)o'ydtdx\"=\"F*Fa`n7%7$ F[vFi_n7$F^vF\\`n7$$\"3cStO*3-BR\"F*Fa`n7%7$FevFi_n7$FhvF\\`n7$$\"3C9o :0%GGg\"F*Fa`n7%7$F_wFi_n7$FbwF\\`n7$$\"3\"zGY4saL\"=F*Fa`n7%7$FiwFi_n 7$F\\xF\\`n7$$\"3!=wNn.\")Q-#F*Fa`n7%7$F($!3=ma4g.U%*\\F]o7$F.$!3l>#y$ Ha*=`&F]o7$$!3!=l5'Qj!)y>F*$!3'f(f#*)\\.oXE\"F]oF_hn7%7$FirFggn7$F\\sFjgn7$$\"3P&eN 2NJ)pLF]oF_hn7%7$FcsFggn7$FfsFjgn7$$\"3a@.j3X4vaF]oF_hn7%7$F]tFggn7$F` tFjgn7$$\"3Ee]_mwN!e(F]oF_hn7%7$FgtFggn7$FjtFjgn7$$\"3+&z>W#3i&o*F]oF_ hn7%7$FauFggn7$FduFjgn7$$\"3F`9B)R)3z6F*F_hn7%7$F[vFggn7$F^vFjgn7$$\"3 &p#4-9Zh*Q\"F*F_hn7%7$FevFggn7$FhvFjgn7$$\"3j+/\")H59+;F*F_hn7%7$F_wFg gn7$FbwFjgn7$$\"3Iu)*fXtm5=F*F_hn7%7$FiwFggn7$F\\xFjgn7$$\"3(zM*QhO>@? F*F_hn7%7$F($!3/!Gdg@_m*HF]o7$F.$!3MKpids8>LF]o7$$!3klq&Rr$\\\")>F*$!3 #elkN.#\\&=$F]o7%7$F9Fe_o7$Ffn\")Rn4x\"F*F]`o7%7$FCFe_o7$F FFh_o7$$!3H=\"yB3T/c\"F*F]`o7%7$FMFe_o7$FPFh_o7$$!3hW')emZ\"*\\8F*F]`o 7%7$FWFe_o7$FZFh_o7$$!3%4<*z]%)QR6F*F]`o7%7$F[oFe_o7$F_oFh_o7$$!3vsp4] 8i)G*F]oF]`o7%7$FfoFe_o7$FioFh_o7$$!3/OA?#>eL=(F]oF]`o7%7$F`pFe_o7$Fcp Fh_o7$$!3J*\\2V.&4y]F]oF]`o7%7$FjpFe_o7$F]qFh_o7$$!3eiFTw=$G(HF]oF]`o7 %7$FdqFe_o7$FgqFh_o7$$!3Bd-=&=(ov')FiqF]`o7%7$F_rFe_o7$FbrFh_o7$$\"3+6 nPRWpP7F]oF]`o7%7$FirFe_o7$F\\sFh_o7$$\"3eZ9F(fdHM$F]oF]`o7%7$FcsFe_o7 $FfsFh_o7$$\"3J%=m^v?#[aF]oF]`o7%7$F]tFe_o7$F`tFh_o7$$\"3/@418R[`vF]oF ]`o7%7$FgtFe_o7$FjtFh_o7$$\"3wdc&42Z(e'*F]oF]`o7%7$FauFe_o7$FduFh_o7$$ \"3WR])G-,k<\"F*F]`o7%7$F[vFe_o7$F^vFh_o7$$\"368XnQt#pQ\"F*F]`o7%7$Fev Fe_o7$FhvFh_o7$$\"3y')RYaOX(f\"F*F]`o7%7$F_wFe_o7$FbwFh_o7$$\"3ZgMDq*z z!=F*F]`o7%7$FiwFe_o7$F\\xFh_o7$$\"39MH/'G1&=?F*F]`o7%7$F($!3WS4>?2%)) )**Fiq7$F.$!3!\\kve3zj5\"F]o7$$!3ZzMI*3\"=%)>F*$!3rML?o0G>5F]o7%7$F9Fc go7$F9-EN\"F*F[ho7%7$FWFcgo7$FZFfgo7$$!3y%eXh #e2U6F*F[ho7%7$F[oFcgo7$F_oFfgo7$$!3**46c.^\\:$*F]oF[ho7%7$FfoFcgo7$Fi oFfgo7$$!3EtjmX>B5sF]oF[ho7%7$F`pFcgo7$FcpFfgo7$$!3aO;x(yo\\5&F]oF[ho7 %7$FjpFcgo7$F]qFfgo7$$!3O+p()Hcq**HF]oF[ho7%7$FdqFcgo7$FgqFfgo7$$!3oL; #)>ZUW*)FiqF[ho7%7$F_rFcgo7$FbrFfgo7$$\"3OtD\"fo?3@\"F]oF[ho7%7$FirFcg o7$F\\sFfgo7$$\"3!)4t!Q%Q3;LF]oF[ho7%7$FcsFcgo7$FfsFfgo7$$\"33Z?q,qM@a F]oF[ho7%7$F]tFcgo7$F`tFfgo7$$\"3\"Qy'ff,hEvF]oF[ho7%7$FgtFcgo7$FjtFfg o7$$\"3a?:\\4>?2%)))**Fiq7$F.$\"3aUc( e3zj5\"F]o7$$!33$*)\\YYoo)>F*$\"3F')z:(*3$p9\"F]o7%7$F9Fa_p7$F/')[@Mw\\,KwW6F*Fi_p7% 7$F[oFa_p7$F_oFd_p7$$!3K[_-d)oBM*F]oFi_p7%7$FfoFa_p7$FioFd_p7$$!3g608* p0rB(F]oFi_p7%7$F`pFa_p7$FcpFd_p7$$!3)[xN7aU=8&F]oFi_p7%7$FjpFa_p7$F]q Fd_p7$$!3:Q5M$Qzl-$F]oFi_p7%7$FdqFa_p7$FgqFd_p7$$!385IYaA;8#*FiqFi_p7% 7$F_rFa_p7$FbrFd_p7$$\"3qN%[C$p%R=\"F]oFi_p7%7$FirFa_p7$F\\sFd_p7$$\"3 -sJM!45#*G$F]oFi_p7%7$FcsFa_p7$FfsFd_p7$$\"3u3zB[KZ%R&F]oFi_p7%7$F]tFa _p7$F`tFd_p7$$\"3YXE81kt*\\(F]oFi_p7%7$FgtFa_p7$FjtFd_p7$$\"3?#QFSc** \\g*F]oFi_p7%7$FauFa_p7$FduFd_p7$$\"3*>@#>si-r6F*Fi_p7%7$F[vFa_p7$F^vF d_p7$$\"3n&o\")ze_:Q\"F*Fi_p7%7$FevFa_p7$FhvFd_p7$$\"3Mf6x.*y?f\"F*Fi_ p7%7$F_wFa_p7$FbwFd_p7$$\"3-L1c>_g-=F*Fi_p7%7$FiwFa_p7$F\\xFd_p7$$\"3# p5]``JJ,#F*Fi_p7%7$F($\"3#yFdg@_m*HF]o7$F.$\"37Ipids8>LF]o7$$!3#pI'**R eb*)>F*$\"3Q2$>DOUJJ$F]o7%7$F9F_gp7$FU@Wy\"**z\"F*Fggp7%7$FiwF_gp7$F\\xFbgp7$$\"33$p.+;W/,#F*Fggp7%7$ F($\"3'RY&4g.U%*\\F]o7$F.$\"3V<#y$Ha*=`&F]o7$$!3w?FM:KC#*>F*$\"3^F1)y# QNzaF]o7%7$F9F]_q7$Fr:F*Fe_q7%7$FMF]_q7$FPF`_q7$$!3t*HuzEk1O\"F*Fe_q7%7$FWF]_q7$ FZF`_q7$$!31E[=_z8]6F*Fe_q7%7$F[oF]_q7$F_oF`_q7$$!3)Q_`RO;hR*F]oFe_q7% 7$FfoF]_q7$FioF`_q7$$!3=(yeg?`3H(F]oFe_q7%7$F`pF]_q7$FcpF`_q7$$!3W]S;[ +f&=&F]oFe_q7%7$FjpF]_q7$F]qF`_q7$$!3<8$p-*oK!3$F]oFe_q7%7$FdqF]_q7$Fg qF`_q7$$!3/jduBtj](*FiqFe_q7%7$F_rF]_q7$FbrF`_q7$$\"3Ug,_D%*>I6F]oFe_q 7%7$FirF]_q7$F\\sF`_q7$$\"3+(*[T$eiaB$F]oFe_q7%7$FcsF]_q7$FfsF`_q7$$\" 3=L'48uD2M&F]oFe_q7%7$F]tF]_q7$F`tF`_q7$$\"3!*pV?**)))fW(F]oFe_q7%7$Fg tF]_q7$FjtF`_q7$$\"3j1\"*4d?D^&*F]oFe_q7%7$FauF]_q7$FduF`_q7$$\"3L%Q* \\@:ll6F*Fe_q7%7$F[vF]_q7$F^vF`_q7$$\"3+e))GPyF*$\"3A[>C$Hlbk(F]o7%7$F9F[gq7$F +\"F]oFcgq7%7$F_rF[gq7$FbrF^gq7$$\"3yAg0scK.6F]oFcgq7%7$FirF[gq7$F\\sF ^gq7$$\"3Af2&*H))e3KF]oFcgq7%7$FcsF[gq7$FfsF^gq7$$\"3&f\\Xy)>&QJ&F]oFc gq7%7$F]tF[gq7$F`tF^gq7$$\"3oK-uX^6>uF]oFcgq7%7$FgtF[gq7$FjtF^gq7$$\"3 Sp\\j.$yV_*F]oFcgq7%7$FauF[gq7$FduF^gq7$$\"3\\qH:YT'H;\"F*Fcgq7%7$F[vF [gq7$F^vF^gq7$$\"3Y!\\t8F*Fcgq7%7$FevF[gq7$FhvF^gq7$$\"3%y\">txn, %e\"F*Fcgq7%7$F_wF[gq7$FbwF^gq7$$\"3_\"R@N4VXz\"F*Fcgq7%7$FiwF[gq7$F\\ xF^gq7$$\"3Tl3J4%p]+#F*Fcgq7%7$F($\"3EO=<[m&**)*)F]o7$F.$\"3/#z!)Gx6u& **F]o7$$!3?[b.mzh(*>F*$\"3.qKgenx6)*F]o7%7$F9Fi^r7$F(RuS8$F]oFa_r7%7$FdqFi^r7$FgqF\\_r7$$!3t^GIR7\")G5F ]oFa_r7%7$F_rFi^r7$FbrF\\_r7$$\"3)\\)=f=>Xw5F]oFa_r7%7$FirFi^r7$F\\sF \\_r7$$\"3W@m[w]r\"=$F]oFa_r7%7$FcsFi^r7$FfsF\\_r7$$\"3re8QM#ypG&F]oFa _r7%7$F]tFi^r7$F`tF\\_r7$$\"3W&4wARTAR(F]oFa_r7%7$FgtFi^r7$FjtF\\_r7$$ \"3;K3<]X](\\*F]oFa_r7%7$FauFi^r7$FduF\\_r7$$\"3*ob13xw-;\"F*Fa_r7%7$F [vFi^r7$F^vF\\_r7$$\"3bIgf'3.3P\"F*Fa_r7%7$FevFi^r7$FhvF\\_r7$$\"3B/bQ -%H8e\"F*Fa_r7%7$F_wFi^r7$FbwF\\_r7$$\"3!z(\\<=d&=z\"F*Fa_r7%7$FiwFi^r 7$F\\xF\\_r7$$\"3e^W'R.#Q-?F*Fa_r7%7$F($\"3G-5AzCx)4\"F*7$F.$\"3)y?jW* p,<7F*7$$!3\"='>QT`I+?F*$\"31fkRA))z(>\"F*7%7$F9Fgfr7$F^kYu9P(F] oF_gr7%7$F`pFgfr7$FcpFjfr7$$!3Cjkb38@m_F]oF_gr7%7$FjpFgfr7$F]qFjfr7$$! 3^E\")QwOltF]oF_gr7%7$FgtFgfr7$FjtFjfr7$$\"3$Qp1nzI1Z*F]oF _gr7%7$FauFgfr7$FduFjfr7$$\"30V,Y&R*ed6F*F_gr7%7$F[vFgfr7$F^vFjfr7$$\" 3s;'\\7r:\"o8F*F_gr7%7$FevFgfr7$FhvFjfr7$$\"3S!4Rq-U'y:F*F_gr7%7$F_wFg fr7$FbwFjfr7$$\"32k&GGMo\"*y\"F*F_gr7%7$FiwFgfr7$F\\xFjfr7$$\"3'z.='eY p**>F*F_gr7%7$F($\"3'3#[i$H\\&)H\"F*7$F.$\"3mO$Q;\"GHQ9F*7$$!3kv$Gnr#* H+#F*$\"39\"fK*o*>WT\"F*7%7$F9Fe^s7$Fe\"F*F]_s7%7$FMFe^s7$FPFh^s7$$!3%[&*f$pPTr8F *F]_s7%7$FWFe^s7$FZFh^s7$$!3<\"[qNX()3;\"F*F]_s7%7$F[oFe^s7$F_oFh^s7$$ !3#R25yP6O]*F]oF]_s7%7$FfoFe^s7$FioFh^s7$$!3?P`\"*>#[$)R(F]oF]_s7%7$F` pFe^s7$FcpFh^s7$$!3Y+1-i]3$H&F]oF]_s7%7$FjpFe^s7$F]qFh^s7$$!3vje7/>#y= $F]oF]_s7%7$FdqFe^s7$FgqFh^s7$$!3,F6BY(eD3\"F]oF]_s7%7$F_rFe^s7$FbrFh^ s7$$\"3q4Om6WqA5F]oF]_s7%7$FirFe^s7$F\\sFh^s7$$\"3UY$e&pv'z7$F]oF]_s7% 7$FcsFe^s7$FfsFh^s7$$\"3:$3`usIKB&F]oF]_s7%7$F]tFe^s7$F`tFh^s7$$\"3))> yM&)Q\\QtF]oF]_s7%7$FgtFe^s7$FjtFh^s7$$\"3gcDCVqvV%*F]oF]_s7%7$FauFe^s 7$FduFh^s7$$\"3AHP6??!\\:\"F*F]_s7%7$F[vFe^s7$F^vFh^s7$$\"3*G?.fLGaO\" F*F]_s7%7$FevFe^s7$FhvFh^s7$$\"3cwEp^Y&fd\"F*F]_s7%7$F_wFe^s7$FbwFh^s7 $$\"3C]@[n4['y\"F*F]_s7%7$FiwFe^s7$F\\xFh^s7$$\"39C;F$G2q*>F*F]_s7%7$F ($\"3UR'G!3hK)\\\"F*7$F.$\"3XlM\")G'o&f;F*7$$!3[*yu?4!o0?F*$\"3@B(oa6T 5j\"F*7%7$F9Fcfs7$FADuF]oF[gs7%7$F`pFcfs7$FcpFffs7$$!3\")Q Z[:)e*>`F]oF[gs7%7$FjpFcfs7$F]qFffs7$$!3`,+fdcp9KF]oF[gs7%7$FdqFcfs7$F gqFffs7$$!3nk_p*\\K%46F]oF[gs7%7$F_rFcfs7$FbrFffs7$$\"3e?Z*>e1$e**FiqF [gs7%7$FirFcfs7$F\\sFffs7$$\"3k3U4;Q4,JF]oF[gs7%7$FcsFcfs7$FfsFffs7$$ \"3\"[%*))R(pN1_F]oF[gs7%7$F]tFcfs7$F`tFffs7$$\"3a\"o$)=8?;J(F]oF[gs7% 7$FgtFcfs7$FjtFffs7$$\"3E=%y(*G$)oT*F]oF[gs7%7$FauFcfs7$FduFffs7$$\"3h :twWY@_6F*F[gs7%7$F[vFcfs7$F^vFffs7$$\"3G*yc0'4ui8F*F[gs7%7$FevFcfs7$F hvFffs7$$\"3&HEYjFnKd\"F*F[gs7%7$F_wFcfs7$FbwFffs7$$\"3iOd8#f$z$y\"F*F [gs7%7$FiwFcfs7$F\\xFffs7$$\"3_5_#z!*>V*>F*F[gs7%7$F($\"3'zXKC#H5)p\"F *7$F.$\"3B%f))fWW3)=F*7$$!3J.7UnuO3?F*$\"3Gb[+iAmZ=F*7%7$F9Fa^t7$F%yju%G(F]oFi^t7%7$FgtFa^t7$FjtFd^t7$$\"3.\"G 9j`4+R*F]oFi^t7%7$FauFa^t7$FduFd^t7$$\"3x,4Ups_\\6F*Fi^t7%7$F[vFa^t7$F ^vFd^t7$$\"3Wv.@&e`+O\"F*Fi^t7%7$FevFa^t7$FhvFd^t7$$\"37\\)**4!*z0d\"F *Fi^t7%7$F_wFa^t7$FbwFd^t7$$\"3zA$*y;i5\"y\"F*Fi^t7%7$FiwFa^t7$F\\xFd^ t7$$\"3o'zyD`K;*>F*Fi^t7%7$F($\"3uwi$otzy*=F*7$F.$\"3-BP;j-7-@F*7$$!39 *R&*eZz 8F*Fgft7%7$FWF_ft7$FZFbft7$$!3XA(4'z&\\*o6F*Fgft7%7$F[oF_ft7$F_oFbft7$ $!3r'[-#QEB%e*F]oFgft7%7$FfoF_ft7$FioFbft7$$!3**\\xI![p*yuF]oFgft7%7$F `pF_ft7$FcpFbft7$$!3F8ITAjqt`F]oFgft7%7$FjpF_ft7$F]qFbft7$$!35x#=X;V%o KF]oFgft7%7$FdqF_ft7$FgqFbft7$$!3#)RNi1+=j6F]oFgft7%7$F_rF_ft7$FbrFbft 7$$\"31p>r7:$3U*FiqFgft7%7$FirF_ft7$F\\sFbft7$$\"32Lf;4jMZIF]oFgft7%7$ FcsF_ft7$FfsFbft7$$\"3Nq11n%4E:&F]oFgft7%7$F]tF_ft7$F`tFbft7$$\"322a& \\isyD(F]oFgft7%7$FgtF_ft7$FjtFbft7$$\"3!Q9]GyNJO*F]oFgft7%7$FauF_ft7$ FduFbft7$$\"3%z[uS*)Ro9\"F*Fgft7%7$F[vF_ft7$F^vFbft7$$\"3hhR')4iOd8F*F gft7%7$FevF_ft7$FhvFbft7$$\"3GNMlDD*yc\"F*Fgft7%7$F_wF_ft7$FbwFbft7$$ \"3'*3HWT)=%yF*Fgft-%+AXE SLABELSG6$%\"xG%\"yG" 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 110 "We're going to solve this differential equation for 0 < x < 1, with a steps ize of 1/100 = .01 (or, 100 steps)." }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "n:= 100; h := .01;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"$+\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% \"hG$\"\"\"!\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 65 "Next we set up the data array for our solution and initialize it." }{MPLTEXT 1 0 0 " " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "y:= \+ vector[n+1];\nx := vector[n+1];\nx[1] := 0;\ny[1] := 1;" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"yG&%'vectorG6#\"$,\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG&%'vectorG6#\"$,\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"xG6#\"\"\"\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"yG6 #\"\"\"F'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 58 "Now we fill in the d ata array with the iteration procedure" }}{PARA 0 "" 0 "" {TEXT -1 59 " x[i+1] = x[i] + h, y[i+1] = y[i] + hf(y[i])." } {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 57 "The basic for mat of a do loop in Maple is the following: " }}{PARA 0 "" 0 "" {TEXT -1 52 " for i from 1 to n do od:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 73 "Next we will plot the numerical solution and the real sol ution together. " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "numerical := pointplot(\{seq([x[i],y[i]], i = 1..n+1) \}):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "actual := plot(exp( 2*x), x = 0..1):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 80 "display (\{numerical,actual\}, title = \"numerical solutiuon vs. actual solut ion\");" }}{PARA 13 "" 1 "" {GLPLOT2D 381 381 381 {PLOTDATA 2 "6'-%'PO INTSG6aq7$$\"#&)!\"#$\"+w(yGQ&!\"*7$$\"#%)F)$\"+L@Lx_F,7$$\"#dF)$\"+!f )y\"4$F,7$$\"#cF)$\"+%Gl6.$F,7$$\"#$)F)$\"+K]&Q<&F,7$$\"##)F)$\"+%*oSs ]F,7$$\"#`F)$\"+ZZLcGF,7$$\"#_F)$\"+$=G.!GF,7$$\"#BF)$\"+i#**od\"F,7$$ \"#AF)$\"+p'zfa\"F,7$$\"#\")F)$\"+Nz%H(\\F,7$$\"#!)F)$\"+_\"Ra([F,7$$ \"#JF)$\"+8))eZ=F,7$$\"#IF)$\"+\"eh8\"=F,7$$\"\"\"F)$\"$-\"F)7$$\"\"#F )$\"&//\"!\"%7$$\"\"$F)$\"(371\"!\"'7$$\"\"%F)$\"*;KC3\"!\")7$$\"\"&F) $\"+.33/6F,7$$\"\"'F)$\"+>C;E6F,7$$\"\"(F)$\"+nco[6F,7$$\"\")F)$\"+!Qf ;<\"F,7$$\"\"*F)$\"+oD4&>\"F,7$$\"#5F)$\"+>W**=7F,7$$\"#6F)$\"+2VPV7F, 7$$\"#7F)$\"+$zT#o7F,7$$\"#8F)$\"+Hmg$H\"F,7$$\"#9F)$\"+i(y%>8F,7$$\"# :F)$\"+P$oeM\"F,7$$\"#;F)$\"+/dys8F,7$$\"#F)$\"+r6\"oX\"F,7$$\"#?F)$\"+%RZf[\"F,7$$\"#@F)$\"+Ujm::F ,7$$\"#'*F)$\"+wJ$Hp'F,7$$\"#(*F)$\"+S=zEoF,7$$\"#)*F)$\"+xwKjpF,7$$\" #**F)$\"+IUf-rF,7$$\"$+\"F)$\"+:hkWsF,7$$\"#pF)$\"+5b8@RF,7$$\"#oF)$\" ++0DWQF,7$$\"#$*F)$\"+;/!pI'F,7$$\"##*F)$\"+-dB$='F,7$$\"#vF)$\"+_a$eT %F,7$$\"#uF)$\"+V/DHVF,7$$\"#lF)$\"+1J_AOF,7$$\"#kF)$\"+TK\\^NF,7$$\"# RF)$\"+mZuk@F,7$$\"#QF)$\"+!z)HA@F,7$$\"#tF)$\"+zJOWUF,7$$\"#sF)$\"+s. 9hTF,7$$\"#VF)$\"+_$*=VBF,7$$\"#UF)$\"+jWC(H#F,7$$\"#\"*F)$\"+'y&*>1'F ,7$$\"#!*F)$\"+BJ8VfF,7$$\"#zF)$\"+1B%)zZF,7$$\"#yF)$\"+3*>ho%F,7$$\"# \\F)$\"+!z6)QEF,7$$\"#[F)$\"+#Qqqe#F,7$$\"#NF)$\"+]&*))**>F,7$$\"#MF)$ \"+Hgng>F,7$$\"#DF)$\"+#*fgS;F,7$$\"#CF)$\"+ZsV3;F,7$$\"#xF)$\"+/_B%f% F,7$$\"#wF)$\"+h@:/XF,7$$\"#^F)$\"+(y>au#F,7$$\"#]F)$\"+E!)e\"p#F,7$$ \"#XF)$\"+-U&yV#F,7$$\"#WF)$\"+RJ0!R#F,7$$\"#FF)$\"+uk)oq\"F,7$$\"#EF) $\"+7\"=Mn\"F,7$$\"#TF)$\"+a/?_AF,7$$\"#SF)$\"+h'R!3AF,7$$\"#nF)$\"+RI ()oPF,7$$\"#mF)$\"+oN(\\p$F,7$$\"#fF)$\"+\\op;KF,7$$\"#eF)$\"+iVi`JF,7 $$\"#LF)$\"+,9BA>F,7$$\"#KF)$\"+*eSX)=F,7$$\"#&*F)$\"+#>*phlF,7$$\"#%* F)$\"+C%QIV'F,7$$\"#()F)$\"+zqM+cF,7$$\"#')F)$\"+_j`!\\&F,7$$\"#bF)$\" +r1trHF,7$$\"#aF)$\"+U9Y8HF,7$$\"#PF)$\"+)3&o!3#F,7$$\"#OF)$\"+Tt))R?F ,7$$\"#hF)$\"+-9lYLF,7$$\"#gF)$\"+'yI5G$F,7$$\"#HF)$\"+(oWex\"F,7$$\"# GF)$\"+.U-Th&=[$F,7$$\"#iF)$\"+IWe8MF,-%'CURVE SG6$7SF_jl7$$\"3dmmm;arz@!#>$\"3%\\@Q-]eX/\"!#<7$$\"3[LL$e9ui2%Fd[m$\" 3yw$*ok3%\\3\"Fg[m7$$\"3nmmm\"z_\"4iFd[m$\"3#Gbv?6BA8\"Fg[m7$$\"3[mmmT &phN)Fd[m$\"3c.'[>4+>=\"Fg[m7$$\"3BLLe*=)H\\5!#=$\"3]zh[5\\]L7Fg[m7$$ \"3fmm\"z/3uC\"Fj\\m$\"3[dVri*fLG\"Fg[m7$$\"3%)***\\7LRDX\"Fj\\m$\"3E8 D#oQ1rL\"Fg[m7$$\"3]mm\"zR'ok;Fj\\m$\"3EDS*pzf]R\"Fg[m7$$\"3v***\\i5`h (=Fj\\m$\"37eED0qKb9Fg[m7$$\"3WLLL3En$4#Fj\\m$\"36(\\$Rkn.?:Fg[m7$$\"3 pmm;/RE&G#Fj\\m$\"3!y&RuWATz:Fg[m7$$\"3\")*****\\K]4]#Fj\\m$\"3h9,pkY. \\;Fg[m7$$\"3$******\\PAvr#Fj\\m$\"3=q&*3:6.Acr4&zk=Fg[m7$$\" 3?LLL347TLFj\\m$\"3PP[sH+x]>Fg[m7$$\"3+LLLLY.KNFj\\m$\"3&p`wF4'pE?Fg[m 7$$\"3v***\\7o7Tv$Fj\\m$\"3Lt*3A/U(=@Fg[m7$$\"3&GLLLQ*o]RFj\\m$\"3'>va TC+P?#Fg[m7$$\"3@++D\"=lj;%Fj\\m$\"3ns0%H:P3I#Fg[m7$$\"31++vV&Rt-EFg[m7$$\"3cmm;/T1&*\\Fj\\m$\"3;3oM@(*f:FFg[m7 $$\"3%em;zRQb@&Fj\\m$\"33*f#H!*H-QGFg[m7$$\"3[***\\(=>Y2aFj\\m$\"3(o:< QUx!\\HFg[m7$$\"39mm;zXu9cFj\\m$\"3I(>N'Rb!R2$Fg[m7$$\"3k******\\y))Ge Fj\\m$\"3qb]3#o;%3KFg[m7$$\"3'*)***\\i_QQgFj\\m$\"3UO:gCPqXLFg[m7$$\"3 @***\\7y%3TiFj\\m$\"3QPxHo]7%[$Fg[m7$$\"35****\\P![hY'Fj\\m$\"3Gia^'*z `WOFg[m7$$\"3jKLL$Qx$omFj\\m$\"3]S[$\\4m\\z$Fg[m7$$\"3!)*****\\P+V)oFj \\m$\"3gMOa3.WiRFg[m7$$\"3?mm\"zpe*zqFj\\m$\"3Z))\\ft4d?TFg[m7$$\"3%)* ****\\#\\'QH(Fj\\m$\"3b#*=a\"Hz1I%Fg[m7$$\"3GKLe9S8&\\(Fj\\m$\"3\")*4] _iHtZ%Fg[m7$$\"3R***\\i?=bq(Fj\\m$\"3'oOiq9T(pYFg[m7$$\"3\"HLL$3s?6zFj \\m$\"39&='R;,&e'[Fg[m7$$\"3a***\\7`Wl7)Fj\\m$\"3o1\\)[.))*z]Fg[m7$$\" 3#pmmm'*RRL)Fj\\m$\"3w***oAXK^H&Fg[m7$$\"3Qmm;a<.Y&)Fj\\m$\"3M_Dyi^dCb Fg[m7$$\"3w&Fg[m7$$\"3t******\\Qk\\*)Fj\\m$ \"3(H&z(G&e-*)fFg[m7$$\"3CLL$3dg6<*Fj\\m$\"3!evAV3D.E'Fg[m7$$\"3Hmmmmx Gp$*Fj\\m$\"3/Kx2'otL^'Fg[m7$$\"3A++D\"oK0e*Fj\\m$\"3Q2zPZHX%z'Fg[m7$$ \"3A++v=5s#y*Fj\\m$\"3,))>Z!eN[2(Fg[m7$Fbjl$\"3S]1$*)4c!*Q(Fg[m-%'COLO URG6&%$RGBG$Fjr!\"\"F`jlF`jl-%&TITLEG6#QHnumerical~solutiuon~vs.~actua l~solution6\"-%+AXESLABELSG6$Q\"xF\\[nQ!F\\[n-%%VIEWG6$;F`jlFbjl%(DEFA ULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1 " "Curve 2" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 240 "You have to look \+ fairly closely to tell that the numerical solution lies below the actu al solution, but it does. This is easier to see if you look on a large r interval, or at a differential equation with a steeper slope field. \+ Play around!." }}}}{MARK "15 0 0" 240 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }