{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 } {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 -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 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 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 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 28 "The maple process of solving" } }{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 45 " y\" \+ = y' + 2y + cos(x), x in [0, Pi/2]" }}{PARA 0 "" 0 "" {TEXT -1 35 " y(0) = -0.3, y(Pi/2)=-0.1" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "a:=0.0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aG\"\"!" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "b:=evalf(Pi/2); # Ca pital P" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG$\"+Fjzq:!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "n:=5; # can be changed" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "h:=(b-a)/n;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"hG$\"+aEfTJ!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "x:=array(0..n,[ seq(a+i*h,i=0..n) ]);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"xG-%&ARRAYG6$7#;\"\"!\"\"&7(/6#F*F*/6#\"\"\"$ \"+aEfTJ!#5/6#\"\"#$\"+3`=$G'F4/6#\"\"$$\"+izxC%*F4/6#\"\"%$\"+iqjc7! \"*/6#F+$\"+Fjzq:FD" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "y:=a rray(0..n); # empty now" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"yG-% &arrayG6$;\"\"!\"\"&7\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 150 "solve( \{ y[0]=-0.3,\nseq( (2-h)*y[i+1]-4*(1+h^2)*y[i]+(2+h)*y[i-1 ] =2*h^2*cos(x[i]) , i=1..n-1 ),\ny[n]=-0.1\}, \{ seq(y[ i],i=0..n) \} );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#<(/&%\"yG6#\"\"!$ !+++++I!#5/&F&6#\"\"&$!+++++5F+/&F&6#\"\"\"$!+!4/Q;$F+/&F&6#\"\"%$!+^u vw=F+/&F&6#\"\"$$!+Wx/sDF+/&F&6#\"\"#$!+b<(f,$F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "assign(%); # assigning above values to y's " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "eval(y); # the y val ues (a table)" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%&ARRAYG6$7#;\"\"! \"\"&7(/6#F($!+++++I!#5/6#\"\"\"$!+!4/Q;$F//6#\"\"#$!+b<(f,$F//6#\"\"$ $!+Wx/sDF//6#\"\"%$!+^uvw=F//6#F)$!+++++5F/" }}}{EXCHG {PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 25 "We want to see the graph:" }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "points:=[ seq([x[i],y[i]],i=0..n ) ];" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%'pointsG7(7$\"\"!$!+++++I!#57$$\"+aEfTJF*$ !+!4/Q;$F*7$$\"+3`=$G'F*$!+b<(f,$F*7$$\"+izxC%*F*$!+Wx/sDF*7$$\"+iqjc7 !\"*$!+^uvw=F*7$$\"+Fjzq:F=$!+++++5F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "plot(points);" }}{PARA 13 "" 1 "" {GLPLOT2D 370 278 278 {PLOTDATA 2 "6%-%'CURVESG6$7(7$\"\"!$!1+++++++I!#;7$$\"1+++aEfTJF+ $!1+++!4/Q;$F+7$$\"1+++3`=$G'F+$!1+++b<(f,$F+7$$\"1,++izxC%*F+$!1+++Wx /sDF+7$$\"1+++iqjc7!#:$!1+++^uvw=F+7$$\"1+++Fjzq:F>$!1+++++++5F+-%'COL OURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$%!GFP-%%VIEWG6$%(DEFAULTGFT " 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "plot(points,style=point);" }}{PARA 13 "" 1 "" {GLPLOT2D 370 278 278 {PLOTDATA 2 "6&-%'CURVESG6$7(7$\"\"!$!1++++ +++I!#;7$$\"1+++aEfTJF+$!1+++!4/Q;$F+7$$\"1+++3`=$G'F+$!1+++b<(f,$F+7$ $\"1,++izxC%*F+$!1+++Wx/sDF+7$$\"1+++iqjc7!#:$!1+++^uvw=F+7$$\"1+++Fjz q:F>$!1+++++++5F+-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$%!GFP -%&STYLEG6#%&POINTG-%%VIEWG6$%(DEFAULTGFX" 1 5 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 49 "We happen to know the exact solution in a formula" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 31 " y(x) = -(sin(x)+3cos(x))/10" }}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{PARA 0 "" 0 "" {TEXT -1 27 "Let's compare the solution." }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "g:=t -> -(sin(t)+3*cos(t))/10; # d efine the exact solution" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"gGR6#% \"tG6\"6$%)operatorG%&arrowGF(,&-%$sinG6#9$#!\"\"\"#5-%$cosGF/#!\"$F3F (F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "plot(g,0..Pi/2);" }}{PARA 13 "" 1 "" {GLPLOT2D 370 278 278 {PLOTDATA 2 "6%-%'CURVESG6$7S 7$\"\"!$!1+++++++I!#;7$$\"1GK5j*))QU$!#<$!1T2u>RZKIF+7$$\"1XXYUk*HS'F/ $!1\\)3wcQy0$F+7$$\"1N68yVJ`(*F/$!1;z$)>37$3$F+7$$\"16YeRSe78F+$!1a9la d20JF+7$$\"1,,hQPB[;F+$!1my![:?M7$F+7$$\"1B*ed(RUf>F+$!1Gs>s[GPJF+7$$ \"1IX[TMk\"G#F+$!1iV\"=bR%[JF+7$$\"1hw(QF$)[h#F+$!1>y/Nu`cJF+7$$\"1G?! >Saq%HF+$!1fcJD77hJF+7$$\"12n4OKt)G$F+$!1/6EVu>iJF+7$$\"1Y]BRTo*e$F+$! 1%[n$ox3gJF+7$$\"1t79wN[GRF+$!1M'e9`)GaJF+7$$\"1](GF+7$$\"1\"3@7LGi%yF+$!1K44dF_HGF+7$$\"16 HIT&[D>)F+$!11v(ovM*yFF+7$$\"1w)\\AI@S\\)F+$!1Js;uA=KFF+7$$\"1o()=V,i> ))F+$!1!pt&[1!*yEF+7$$\"1`\"[dg&*f:*F+$!17g_!)Q(3i#F+7$$\"1Tt6rL2&[*F+ $!1*4<18O7c#F+7$$\"1la(*HIZ.)*F+$!1Wq&G,$*3]#F+7$$\"1m\"[l:+d,\"!#:$!1 'yw(e[#4V#F+7$$\"1\"3XyEmu/\"F`u$!18(f.`haO#F+7$$\"1K_*>P$Q\"3\"F`u$!1 t(Q()\\KHH#F+7$$\"1M\"G%4t676F`u$!1$GEQJI\\A#F+7$$\"1p*Q4iq'*z.@\"F`u$!1jBI(oTP *>F+7$$\"1Q$[)>&*oU7F`u$!1mFlZ$3M\">F+7$$\"1^lOFY^w7F`u$!1)QF]Cor#=F+7 $$\"1\")f6EA448F`u$!1Tt$\\\\J@u\"F+7$$\"1;T7EvSU8F`u$!1V6E@yD`;F+7$$\" 1`wiepWv8F`u$!1(R)oBxIj:F+7$$\"1$obdw1eS\"F`u$!1ZKJrh9z9F+7$$\"1%)o#3` -1W\"F`u$!17W;bc,\"Q\"F+7$$\"1#*)4zFC\"F+7$$\"13Y:@imO:F`u$!1\"=\"4k!)y,6F+7$$\"1++lBjz q:F`u$!1!pM%4+++5F+-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$%!G Fb[l-%%VIEWG6$;F($\"+Fjzq:!\"*%(DEFAULTG" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "plot1:=plot(points):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot2:=plot(g,0..Pi/2):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "plots[display](\{plot1,plot2\});" }}{PARA 13 "" 1 "" {GLPLOT2D 370 278 278 {PLOTDATA 2 "6&-%'CURVESG6$7S7$\"\"!$!1+++++++I!#;7$$\"1GK 5j*))QU$!#<$!1T2u>RZKIF+7$$\"1XXYUk*HS'F/$!1\\)3wcQy0$F+7$$\"1N68yVJ`( *F/$!1;z$)>37$3$F+7$$\"16YeRSe78F+$!1a9lad20JF+7$$\"1,,hQPB[;F+$!1my![ :?M7$F+7$$\"1B*ed(RUf>F+$!1Gs>s[GPJF+7$$\"1IX[TMk\"G#F+$!1iV\"=bR%[JF+ 7$$\"1hw(QF$)[h#F+$!1>y/Nu`cJF+7$$\"1G?!>Saq%HF+$!1fcJD77hJF+7$$\"12n4 OKt)G$F+$!1/6EVu>iJF+7$$\"1Y]BRTo*e$F+$!1%[n$ox3gJF+7$$\"1t79wN[GRF+$! 1M'e9`)GaJF+7$$\"1]( GF+7$$\"1\"3@7LGi%yF+$!1K44dF_HGF+7$$\"16HIT&[D>)F+$!11v(ovM*yFF+7$$\" 1w)\\AI@S\\)F+$!1Js;uA=KFF+7$$\"1o()=V,i>))F+$!1!pt&[1!*yEF+7$$\"1`\"[ dg&*f:*F+$!17g_!)Q(3i#F+7$$\"1Tt6rL2&[*F+$!1*4<18O7c#F+7$$\"1la(*HIZ.) *F+$!1Wq&G,$*3]#F+7$$\"1m\"[l:+d,\"!#:$!1'yw(e[#4V#F+7$$\"1\"3XyEmu/\" F`u$!18(f.`haO#F+7$$\"1K_*>P$Q\"3\"F`u$!1t(Q()\\KHH#F+7$$\"1M\"G%4t676 F`u$!1$GEQJI\\A#F+7$$\"1p*Q4iq'*z.@\"F`u$!1jBI(oTP*>F+7$$\"1Q$[)>&*oU7F`u$!1mFl Z$3M\">F+7$$\"1^lOFY^w7F`u$!1)QF]Cor#=F+7$$\"1\")f6EA448F`u$!1Tt$\\\\J @u\"F+7$$\"1;T7EvSU8F`u$!1V6E@yD`;F+7$$\"1`wiepWv8F`u$!1(R)oBxIj:F+7$$ \"1$obdw1eS\"F`u$!1ZKJrh9z9F+7$$\"1%)o#3`-1W\"F`u$!17W;bc,\"Q\"F+7$$\" 1#*)4zFC\"F+7$$\" 13Y:@imO:F`u$!1\"=\"4k!)y,6F+7$$\"1++lBjzq:F`u$!1!pM%4+++5F+-%'COLOURG 6&%$RGBG$\"#5!\"\"F(F(-F$6$7(F'7$$\"1+++aEfTJF+$!1+++!4/Q;$F+7$$\"1+++ 3`=$G'F+$!1+++b<(f,$F+7$$\"1,++izxC%*F+$!1+++Wx/sDF+7$$\"1+++iqjc7F`u$ !1+++^uvw=F+7$$\"1+++Fjzq:F`u$!1+++++++5F+Fhz-%+AXESLABELSG6$%!GF^]l-% %VIEWG6$;F($\"+Fjzq:!\"*%(DEFAULTG" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}}{MARK "21 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }