{VERSION 4 0 "IBM INTEL NT" "4.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 256 "" 0 1 0 0 0 0 0 2 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 "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 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 
2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 2 1 0 1 0 2 2 0 1 }{PSTYLE "Warning" 
-1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 1 3 1 }
1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "
" -1 -1 "Times" 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 }{PSTYLE "Title
" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 1 1 1 1 }
3 1 0 0 12 12 1 0 1 0 2 2 19 1 }{PSTYLE "Author" -1 19 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 8 8 1 0 1 0 2 
2 0 1 }}
{SECT 0 {PARA 18 "" 0 "" {TEXT -1 47 "Euler's Method for \nNumerical I
ntegrating ODE's" }}{PARA 19 "" 0 "" {TEXT -1 0 "" }}{SECT 0 {PARA 3 "
" 0 "" {TEXT -1 14 "Initialization" }}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 36 "restart:\nwith(plots):\nwith(DEtools):" }}{PARA 7 "" 
1 "" {TEXT -1 50 "Warning, the name changecoords has been redefined\n
" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 22 "Derive Euler's Method:" }}
{SECT 0 {PARA 4 "" 0 "" {TEXT -1 33 "First order differential equation
" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "diff(y(x),x) = f(x,y(x))
;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%%diffG6$-%\"yG6#%\"xGF*-%\"fG6
$F*F'" }}}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 31 "Finite difference appr
oximation" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "finite_diff_eq \+
:= (y[i+1]-y[i])/h = f(x[i],y[i]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#
>%/finite_diff_eqG/*&,&&%\"yG6#,&%\"iG\"\"\"F-F-F-&F)6#F,!\"\"F-%\"hGF
0-%\"fG6$&%\"xGF/F." }}}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 9 "Solve for
" }{TEXT 256 1 " " }{XPPEDIT 18 0 "y[i+1];" "6#&%\"yG6#,&%\"iG\"\"\"F(
F(" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "y[i+1] := solve(finite
_diff_eq,y[i+1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"yG6#,&%\"iG
\"\"\"F)F),&&F%6#F(F)*&-%\"fG6$&%\"xGF,F+F)%\"hGF)F)" }}}}}{SECT 0 
{PARA 3 "" 0 "" {TEXT -1 26 "Set up initial conditions:" }}{SECT 0 
{PARA 4 "" 0 "" {TEXT -1 58 "Step size, number of intervals, initial c
onditions, f(x,y)" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "h := 0.
1;\nN := 20;\nx0 := 0.0;\ny0 := 1.0;\nf := (x,y) -> cos(x+y);" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"hG$\"\"\"!\"\"" }}{PARA 11 "" 1 "
" {XPPMATH 20 "6#>%\"NG\"#?" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x0G$
\"\"!F&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y0G$\"#5!\"\"" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#>%\"fGR6$%\"xG%\"yG6\"6$%)operatorG%&arrowGF
)-%$cosG6#,&9$\"\"\"9%F2F)F)F)" }}}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 
17 "Initialize arrays" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "x[0
] := x0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"xG6#\"\"!$F'F'" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "y[0] := y0;" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#>&%\"yG6#\"\"!$\"#5!\"\"" }}}}}{SECT 0 {PARA 3 "" 
0 "" {TEXT -1 20 "Apply Euler's method" }}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 81 "for i from 0 to N do\n  x[i+1] := x[i] + h:\n  y[i+1]
 := y[i] + f(x[i],y[i])*h:\nod:" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 
20 "Plot the data points" }}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 27 "Creat
e array of data points" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "l \+
:= [[x[n],y[n]] $n=0..N]:" }}}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 13 "No
w, plot it!" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "plot(l,style=
point);\nplot1 := plot(l,style=point,color=red):" }}{PARA 13 "" 1 "" 
{GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$777$$\"\"!F)$\"\"\"F)
7$$\"3/+++++++5!#=$\"3/+++J-.a5!#<7$$\"35+++++++?F/$\"3#*******y2^%4\"
F27$$\"3))**************HF/$\"3(******\\<*y@6F27$$\"3A+++++++SF/$\"3(*
*****4\"[jO6F27$$\"3++++++++]F/$\"3-+++)H]+9\"F27$$\"3w**************f
F/$\"3*******\\VIJ8\"F27$$\"3a**************pF/$\"3&******zAop6\"F27$$
\"3U+++++++!)F/$\"33+++A*)f#4\"F27$$\"3A+++++++!*F/$\"3++++37(41\"F27$
F*$\"3(******>>OH-\"F27$$\"33+++++++6F2$\"3g*****zaqCz*F/7$$\"3%******
********>\"F2$\"30+++0gk0$*F/7$$\"3/+++++++8F2$\"3-+++`hlu()F/7$$\"3!*
*************R\"F2$\"3Y*****p(=_/#)F/7$$\"3++++++++:F2$\"3_+++q%4'*f(F
/7$$\"33+++++++;F2$\"3!******Hz;P'pF/7$$\"3%**************p\"F2$\"3C++
+43:+jF/7$$\"3/+++++++=F2$\"3t******[ez6cF/7$$\"3!***************=F2$
\"3!******R(G<,\\F/7$$\"\"#F)$\"3!******zY!\\qTF/-%'COLOURG6&%$RGBG$\"
#5!\"\"F(F(-%+AXESLABELSG6$Q!6\"Fdr-%&STYLEG6#%&POINTG-%%VIEWG6$%(DEFA
ULTGF]s" 1 5 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curv
e 1" }}}}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 27 "Maple (Analytical) Sol
ution" }}{PARA 0 "" 0 "" {TEXT -1 48 "Use the same f, but different de
pendent variable" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "unassign
(x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "ode:=diff(u(x),x) =
 f(x,u(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$odeG/-%%diffG6$-%\"u
G6#%\"xGF,-%$cosG6#,&F,\"\"\"F)F1" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 9 "u0 := y0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u0G$\"
#5!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "dsolve( \{ode,u(
0)=u0\}, u(x) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"xG,&F'
!\"\"*&\"\"#\"\"\"-%'arctanG6#,&F'F,-%$tanG6##F,F+F,F,F," }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "ans := unapply(rhs(%),x);" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#>%$ansGR6#%\"xG6\"6$%)operatorG%&arrowGF(,&9
$!\"\"*&\"\"#\"\"\"-%'arctanG6#,&F-F1-%$tanG6##F1F0F1F1F1F(F(F(" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 70 "plot( ans(x), x=0..N*h);\npl
ot2 := plot( ans(x), x=0..N*h, color=blue):" }}{PARA 13 "" 1 "" 
{GLPLOT2D 400 300 300 {PLOTDATA 2 "6%-%'CURVESG6$7S7$$\"\"!F)$\"\"\"F)
7$$\"39LLLL3VfV!#>$\"3:MG.)z?B-\"!#<7$$\"3'pmm;H[D:)F/$\"3dsJ7w^tR5F27
$$\"3LLLLe0$=C\"!#=$\"3;'*f=905d5F27$$\"3ILLL3RBr;F;$\"3i[1U>kBs5F27$$
\"3Ymm;zjf)4#F;$\"3cv*)e,o+&3\"F27$$\"3=LL$e4;[\\#F;$\"3Y+\"[*[$[[4\"F
27$$\"3p****\\i'y]!HF;$\"3P%pD8qlI5\"F27$$\"3,LL$ezs$HLF;$\"3O=M_6_^46
F27$$\"3_****\\7iI_PF;$\"3<y(*HG*RR6\"F27$$\"3#pmmm@Xt=%F;$\"3ln$e3\\!
[;6F27$$\"3QLLL3y_qXF;$\"3#\\_8S>$4<6F27$$\"3i******\\1!>+&F;$\"3T$Q^n
xMg6\"F27$$\"3()******\\Z/NaF;$\"3t4M)*>m=86F27$$\"3'*******\\$fC&eF;$
\"3D+C0kR#)36F27$$\"3ELL$ez6:B'F;$\"3v6jRPka.6F27$$\"3Smmm;=C#o'F;$\"3
_9K_DWr&4\"F27$$\"3-mmmm#pS1(F;$\"3=f2/<d\"y3\"F27$$\"3]****\\i`A3vF;$
\"3y\\)4YOQs2\"F27$$\"3slmmm(y8!zF;$\"3%pNg4H'om5F27$$\"3V++]i.tK$)F;$
\"3mL.1Of)Q0\"F27$$\"39++](3zMu)F;$\"3))y)>Mcj0/\"F27$$\"3#pmm;H_?<*F;
$\"3wl(G?rTb-\"F27$$\"3emm;zihl&*F;$\"3y\"=6`-'y55F27$$\"39LLL3#G,***F
;$\"3Y<U4Z)z)Q**F;7$$\"3<LLezw5V5F2$\"3gY@Vpc3`(*F;7$$\"3!****\\PQ#\\
\"3\"F2$\"3Wyg'osPKe*F;7$$\"3BLL$e\"*[H7\"F2$\"3G._3'p:<R*F;7$$\"3#***
****pvxl6F2$\"3oWsm59X&=*F;7$$\"3z****\\_qn27F2$\"39o]6&y\"zv*)F;7$$\"
3%)***\\i&p@[7F2$\"3!ye=E%\\(ew)F;7$$\"3#)****\\2'HKH\"F2$\"3!*pt%fd\\
]_)F;7$$\"3_mmmwanL8F2$\"3&Qu$3)\\F?I)F;7$$\"3'******\\2goP\"F2$\"3oT?
,/XId!)F;7$$\"3CLLeR<*fT\"F2$\"3%pIVeFF*HyF;7$$\"3'******\\)Hxe9F2$\"3
'=_A^E0bd(F;7$$\"3Ymm\"H!o-*\\\"F2$\"3UjhlZ(*yItF;7$$\"3))***\\7k.6a\"
F2$\"3a<II^&>(pqF;7$$\"3emmmT9C#e\"F2$\"3Jx!GfX$\\4oF;7$$\"3\"****\\i!
*3`i\"F2$\"39BtT*>;?`'F;7$$\"3QLLL$*zym;F2$\"3%eP*Q)42,E'F;7$$\"3GLL$3
N1#4<F2$\"3+ygzXG^xfF;7$$\"3kmm\"HYt7v\"F2$\"3ihQ)fjIHp&F;7$$\"3%*****
**p(G**y\"F2$\"3\"4)Q-[syFaF;7$$\"3lmm;9@BM=F2$\"3DGt@sCz>^F;7$$\"3ELL
L`v&Q(=F2$\"3.S\\$Hm\\2%[F;7$$\"30++DOl5;>F2$\"3H)Hj[^t'RXF;7$$\"3/++v
.Uac>F2$\"3WVRb-A?[UF;7$$\"\"#F)$\"3gA0*H*4_JRF;-%'COLOURG6&%$RGBG$\"#
5!\"\"F(F(-%+AXESLABELSG6$Q\"x6\"Q!6\"-%%VIEWG6$;F($\"#?F_[l%(DEFAULTG
" 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 "" {MPLTEXT 1 0 23 "display(\{plot1,plot2\});
" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$
7S7$$\"\"!F)$\"\"\"F)7$$\"39LLLL3VfV!#>$\"3:MG.)z?B-\"!#<7$$\"3'pmm;H[
D:)F/$\"3dsJ7w^tR5F27$$\"3LLLLe0$=C\"!#=$\"3;'*f=905d5F27$$\"3ILLL3RBr
;F;$\"3i[1U>kBs5F27$$\"3Ymm;zjf)4#F;$\"3cv*)e,o+&3\"F27$$\"3=LL$e4;[\\
#F;$\"3Y+\"[*[$[[4\"F27$$\"3p****\\i'y]!HF;$\"3P%pD8qlI5\"F27$$\"3,LL$
ezs$HLF;$\"3O=M_6_^46F27$$\"3_****\\7iI_PF;$\"3<y(*HG*RR6\"F27$$\"3#pm
mm@Xt=%F;$\"3ln$e3\\![;6F27$$\"3QLLL3y_qXF;$\"3#\\_8S>$4<6F27$$\"3i***
***\\1!>+&F;$\"3T$Q^nxMg6\"F27$$\"3()******\\Z/NaF;$\"3t4M)*>m=86F27$$
\"3'*******\\$fC&eF;$\"3D+C0kR#)36F27$$\"3ELL$ez6:B'F;$\"3v6jRPka.6F27
$$\"3Smmm;=C#o'F;$\"3_9K_DWr&4\"F27$$\"3-mmmm#pS1(F;$\"3=f2/<d\"y3\"F2
7$$\"3]****\\i`A3vF;$\"3y\\)4YOQs2\"F27$$\"3slmmm(y8!zF;$\"3%pNg4H'om5
F27$$\"3V++]i.tK$)F;$\"3mL.1Of)Q0\"F27$$\"39++](3zMu)F;$\"3))y)>Mcj0/
\"F27$$\"3#pmm;H_?<*F;$\"3wl(G?rTb-\"F27$$\"3emm;zihl&*F;$\"3y\"=6`-'y
55F27$$\"39LLL3#G,***F;$\"3Y<U4Z)z)Q**F;7$$\"3<LLezw5V5F2$\"3gY@Vpc3`(
*F;7$$\"3!****\\PQ#\\\"3\"F2$\"3Wyg'osPKe*F;7$$\"3BLL$e\"*[H7\"F2$\"3G
._3'p:<R*F;7$$\"3#*******pvxl6F2$\"3oWsm59X&=*F;7$$\"3z****\\_qn27F2$
\"39o]6&y\"zv*)F;7$$\"3%)***\\i&p@[7F2$\"3!ye=E%\\(ew)F;7$$\"3#)****\\
2'HKH\"F2$\"3!*pt%fd\\]_)F;7$$\"3_mmmwanL8F2$\"3&Qu$3)\\F?I)F;7$$\"3'*
*****\\2goP\"F2$\"3oT?,/XId!)F;7$$\"3CLLeR<*fT\"F2$\"3%pIVeFF*HyF;7$$
\"3'******\\)Hxe9F2$\"3'=_A^E0bd(F;7$$\"3Ymm\"H!o-*\\\"F2$\"3UjhlZ(*yI
tF;7$$\"3))***\\7k.6a\"F2$\"3a<II^&>(pqF;7$$\"3emmmT9C#e\"F2$\"3Jx!GfX
$\\4oF;7$$\"3\"****\\i!*3`i\"F2$\"39BtT*>;?`'F;7$$\"3QLLL$*zym;F2$\"3%
eP*Q)42,E'F;7$$\"3GLL$3N1#4<F2$\"3+ygzXG^xfF;7$$\"3kmm\"HYt7v\"F2$\"3i
hQ)fjIHp&F;7$$\"3%*******p(G**y\"F2$\"3\"4)Q-[syFaF;7$$\"3lmm;9@BM=F2$
\"3DGt@sCz>^F;7$$\"3ELLL`v&Q(=F2$\"3.S\\$Hm\\2%[F;7$$\"30++DOl5;>F2$\"
3H)Hj[^t'RXF;7$$\"3/++v.Uac>F2$\"3WVRb-A?[UF;7$$\"\"#F)$\"3gA0*H*4_JRF
;-%'COLOURG6&%$RGBGF(F($\"*++++\"!\")-F$6%77F'7$$\"3/+++++++5F;$\"3/++
+J-.a5F27$$\"35+++++++?F;$\"3#*******y2^%4\"F27$$\"3))**************HF
;$\"3(******\\<*y@6F27$$\"3A+++++++SF;$\"3(******4\"[jO6F27$$\"3++++++
++]F;$\"3-+++)H]+9\"F27$$\"3w**************fF;$\"3*******\\VIJ8\"F27$$
\"3a**************pF;$\"3&******zAop6\"F27$$\"3U+++++++!)F;$\"33+++A*)
f#4\"F27$$\"3A+++++++!*F;$\"3++++37(41\"F27$F*$\"3(******>>OH-\"F27$$
\"33+++++++6F2$\"3g*****zaqCz*F;7$$\"3%**************>\"F2$\"30+++0gk0
$*F;7$$\"3/+++++++8F2$\"3-+++`hlu()F;7$$\"3!**************R\"F2$\"3Y**
***p(=_/#)F;7$$\"3++++++++:F2$\"3_+++q%4'*f(F;7$$\"33+++++++;F2$\"3!**
****Hz;P'pF;7$$\"3%**************p\"F2$\"3C+++43:+jF;7$$\"3/+++++++=F2
$\"3t******[ez6cF;7$$\"3!***************=F2$\"3!******R(G<,\\F;7$Fez$
\"3!******zY!\\qTF;-Fjz6&F\\[lF][lF(F(-%&STYLEG6#%&POINTG-%+AXESLABELS
G6%Q\"x6\"Q!6\"%(DEFAULTG-%%VIEWG6$;F($\"#?!\"\"F`bl" 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 "" {MPLTEXT 1 0 0 "" }}}}}{MARK "8" 0 }
{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }