{VERSION 5 0 "IBM INTEL LINUX" "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 1 }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 "Heading 1" 0 3 1 {CSTYLE "" -1 -1 
"" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }1 0 0 0 8 4 0 0 0 0 0 0 -1 0 }
{PSTYLE "Heading 2" 3 4 1 {CSTYLE "" -1 -1 "" 1 14 0 0 0 0 0 0 0 0 0 
0 0 0 0 1 }0 0 0 -1 8 2 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 1 }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 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Title" 0 18 1 
{CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 1 0 0 0 0 0 0 1 }3 0 0 -1 12 12 
0 0 0 0 0 0 19 0 }{PSTYLE "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 
0 0 0 0 0 0 0 0 0 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 
257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }3 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#%5Math~696.200~and~700G" }}}{PARA 18 "" 0 "" {TEXT -1 
26 "Euler - First Order Method" }}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 17 
"Set up domain ..." }}{PARA 257 "" 0 "" {TEXT -1 48 "a < x < b, \nleng
th L=b-a, \nN steps of size h=L/N" }}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 48 "N := 20:\na := 0.0:\nb := 1.0:\nL := b-a:\nh := L/N:
" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 32 "Define differential equatio
n ..." }}{PARA 256 "" 0 "" {TEXT -1 14 "dy/dx = f(x,y)" }}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "f := (x,y) -> y:\nexact := x -> exp
(x):" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 33 "Euler Method (1st order
) Code ..." }}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 33 "define arrays for x
 and y values " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "x := array
(0..N):\ny := array(0..N):\nx[0] := a:\ny[0] := exact(a):" }}}}{SECT 
0 {PARA 4 "" 0 "" {TEXT -1 17 "integrate N steps" }}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 83 "for i from 1 to N do\n  x[i] := x[i-1]+h;  \n \+
 y[i] := y[i-1]+h*f(x[i-1],y[i-1]);\nod:" }}}}}{SECT 0 {PARA 3 "" 0 "
" {TEXT -1 34 "Compute error at the last step ..." }}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 112 "# error at last step!\n# WARNING: `error` mi
ght conflict with Maple's meaning of that name\nerr := exact(b)-y[N];
\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$errG$\")CT)\\'!\"*" }}}}
{SECT 0 {PARA 3 "" 0 "" {TEXT -1 54 "Now plot the numerical solution a
nd the exact solution" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "wit
h(plots):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoor
ds has been redefined\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 116 
"plot1 := plot( exact(x), x=a..b, color=blue):\nplot2 := plot( [[x[n],
y[n]] $n=0..N],x=a..b):\ndisplay(\{plot1,plot2\});\n" }}{PARA 13 "" 1 
"" {GLPLOT2D 331 331 331 {PLOTDATA 2 "6&-%'CURVESG6$777$$\"\"!F)$\"\"
\"F)7$$\"3G+++++++]!#>$\"3/++++++]5!#<7$$\"3/+++++++5!#=$\"3/+++++]-6F
27$$\"3'**************\\\"F6$\"3\"*********\\id6F27$$\"37+++++++?F6$\"
3,+++]i]:7F27$$\"3++++++++DF6$\"34+++i:Gw7F27$$\"3))**************HF6$
\"3#*******Rc4S8F27$$\"3y*************\\$F6$\"32+++A/529F27$$\"3C+++++
++SF6$\"3)******HWbuZ\"F27$$\"37+++++++XF6$\"32+++:#G8b\"F27$$\"3+++++
+++]F6$\"3!******fi%*)G;F27$$\"3W+++++++bF6$\"31+++d$R.r\"F27$$\"3y***
***********fF6$\"3#******\\Kcez\"F27$$\"3C+++++++lF6$\"3.+++T\"\\c)=F2
7$$\"3c**************pF6$\"33+++)fJ*z>F27$$\"3++++++++vF6$\"3'******z<
G*y?F27$$\"3W+++++++!)F6$\"39+++(euG=#F27$$\"3y*************\\)F6$\"31
+++;$=?H#F27$$\"3C+++++++!*F6$\"3;+++K#>mS#F27$$\"3c*************\\*F6
$\"3()*****R>]p_#F27$F*$\"3()*****Rq(H`EF2-%'COLOURG6&%$RGBG$\"#5!\"\"
F(F(-F$6$7SF'7$$\"3dmmm;arz@F/$\"3)y_2wWO?-\"F27$$\"3[LL$e9ui2%F/$\"3a
NdbY\\gT5F27$$\"3nmmm\"z_\"4iF/$\"3$=fT9tfS1\"F27$$\"3[mmmT&phN)F/$\"3
<Hr`&G_r3\"F27$$\"3BLLe*=)H\\5F6$\"3L^LEiEj56F27$$\"3gmm\"z/3uC\"F6$\"
3r]yZ%yaG8\"F27$$\"3%)***\\7LRDX\"F6$\"3loPGkJLc6F27$$\"3^mm\"zR'ok;F6
$\"3-Am\"\\\\E6=\"F27$$\"3w***\\i5`h(=F6$\"36)e/'[$pj?\"F27$$\"3YLLL3E
n$4#F6$\"3tRSV5x*GB\"F27$$\"3qmm;/RE&G#F6$\"3*>3j2pYnD\"F27$$\"3\")***
**\\K]4]#F6$\"3Tm9jYu9%G\"F27$$\"3$******\\PAvr#F6$\"3%*p#e<$=E78F27$$
\"3(******\\nHi#HF6$\"3gYw>&\\P*R8F27$$\"3jmm\"z*ev:JF6$\"3AKG01]dl8F2
7$$\"3?LLL347TLF6$\"3A9H.#p*p'R\"F27$$\"3,LLLLY.KNF6$\"3KLbfo2iB9F27$$
\"3v***\\7o7Tv$F6$\"3#ouv#H**eb9F27$$\"3'GLLLQ*o]RF6$\"3s1/UDl[%[\"F27
$$\"3A++D\"=lj;%F6$\"3\\;8&[1^o^\"F27$$\"31++vV&R<P%F6$\"3PKsg\"RD$[:F
27$$\"3XLL$e9Ege%F6$\"3JJ$))*o>'=e\"F27$$\"3HLLeR\"3Gy%F6$\"3+.c]c%)H8
;F27$$\"3emm;/T1&*\\F6$\"3$**eh)zw!zk\"F27$$\"3&em;zRQb@&F6$\"3G,G!GGV
Yo\"F27$$\"3\\***\\(=>Y2aF6$\"3L6+U5yG<<F27$$\"39mm;zXu9cF6$\"3&Q>jrpb
Kv\"F27$$\"3l******\\y))GeF6$\"33!=v7P07z\"F27$$\"3'*)***\\i_QQgF6$\"3
XOlB#\\E\"H=F27$$\"3@***\\7y%3TiF6$\"3neT_<6em=F27$$\"36****\\P![hY'F6
$\"3DV)H1Jn!4>F27$$\"3jKLL$Qx$omF6$\"3vZcmrs1[>F27$$\"3!)*****\\P+V)oF
6$\"39RAPIze!*>F27$$\"3?mm\"zpe*zqF6$\"3hLeRd*=*H?F27$$\"3%)*****\\#\\
'QH(F6$\"3zh'*f?z!Q2#F27$$\"3GKLe9S8&\\(F6$\"3^**33Q,(f6#F27$$\"3R***
\\i?=bq(F6$\"3m>.C'Qe4;#F27$$\"3\"HLL$3s?6zF6$\"3C(p,G?ne?#F27$$\"3a**
*\\7`Wl7)F6$\"3#[n%ezG)QD#F27$$\"3#pmmm'*RRL)F6$\"3^8U**za6,BF27$$\"3Q
mm;a<.Y&)F6$\"3+^n_#[T/N#F27$$\"3<LLe9tOc()F6$\"3wf*=(>KS+CF27$$\"3t**
****\\Qk\\*)F6$\"3u&G]\"H'[sW#F27$$\"3DLL$3dg6<*F6$\"3bh$4Z;k?]#F27$$
\"3HmmmmxGp$*F6$\"38[Xr/78_DF27$$\"3B++D\"oK0e*F6$\"3j4')HYrh1EF27$$\"
3B++v=5s#y*F6$\"3!*=xB3j&)fEF27$F*$\"34X!f%G=G=FF2-F\\r6&F^rF(F($\"*++
++\"!\")-%+AXESLABELSG6%Q\"x6\"Q!F\\bl-%%FONTG6#%(DEFAULTG-%%VIEWG6$;F
(F_rFabl" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Cur
ve 1" "Curve 2" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}}
{MARK "6 3 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 
2 33 1 1 }