{VERSION 2 3 "IBM INTEL NT" "2.3" } {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 "Times" 1 12 0 0 0 1 2 1 2 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 P lot" 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 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 12 128 0 128 1 2 1 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 11 0 128 128 1 2 1 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 } } {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "Euler's method - Example 1 " }}{PARA 0 "" 0 "" {TEXT -1 36 "------------------------------------ " }}{PARA 0 "" 0 "" {TEXT -1 32 "Solve the initial value problem:" }} {PARA 0 "" 0 "" {TEXT -1 52 " dy/dx = 2y y( 0) = 1 " }}{PARA 0 "" 0 "" {TEXT -1 125 "in the domain [ 0 , 2 ] u sing various values of step size h. Compare the approximations with th e true solution y = e ^ (2x)." }}{PARA 0 "" 0 "" {TEXT -1 19 " \+ " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "ExactSolution := (x) -> exp( 2*x);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%.ExactSolutionG:6#%\"xG6\"6$%)operatorG%&arrowGF(-%$e xpG6#,$9$\"\"#F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 63 "The exact s olution is indicated by the blue curve in the graph." }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "picture1 \+ := plot(ExactSolution(x), x=0..2, color=COLOR(RGB,0,0,1)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "h:=0.5; # the step size" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"hG$\"\"&!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "oldx:=0; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%o ldxG\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "oldy:=1; # th e initial condition" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG\"\"\" " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "lastx:=2;" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%&lastxG\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 6 "L:=[];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% \"LG7\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "L:=[op(L),[oldx, oldy]];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7#7$\"\"!\"\"\"" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "slope:=(x,y) -> 2*y;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%&slopeG:6$%\"xG%\"yG6\"6$%)operatorG %&arrowGF),$9%\"\"#F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 " while oldx < lastx do oldy:=oldy+h*slope(oldx,oldy);oldx:=oldx+h;L:=[o p(L),[oldx,oldy]];od;" }}{PARA 11 "" 0 "" {TEXT -1 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"#?!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"&!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7$7 $\"\"!\"\"\"7$$\"\"&!\"\"$\"#?F," }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% %oldyG$\"$+%!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#5!\"\" " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7%7$\"\"!\"\"\"7$$\"\"&!\"\" $\"#?F,7$$\"#5F,$\"$+%!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG $\"%+!)!\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#:!\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7&7$\"\"!\"\"\"7$$\"\"&!\"\"$\"# ?F,7$$\"#5F,$\"$+%!\"#7$$\"#:F,$\"%+!)!\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"'++;!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%%oldxG$\"#?!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7'7$\"\"! \"\"\"7$$\"\"&!\"\"$\"#?F,7$$\"#5F,$\"$+%!\"#7$$\"#:F,$\"%+!)!\"$7$F-$ \"'++;!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "picture2:=plo t(L,style=point,symbol=diamond );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#> %)picture2G-%/INTERFACE_PLOTG6%-%'CURVESG6$7'7$\"\"!$\"\"\"F-7$$\"1+++ ++++]!#;$\"\"#F-7$F.$\"\"%F-7$$\"1+++++++:!#:$\"\")F-7$F4$\"#;F--%'COL OURG6&%$RGBG$\"#5!\"\"F-F--%'SYMBOLG6#%(DIAMONDG-%&STYLEG6#%&POINTG" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "display(\{picture1,picture 2\});" }}{PARA 13 "" 1 "" {INLPLOT "6&-%'CURVESG6&7'7$\"\"!$\"\"\"F(7$ $\"1+++++++]!#;$\"\"#F(7$F)$\"\"%F(7$$\"1+++++++:!#:$\"\")F(7$F/$\"#;F (-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%'SYMBOLG6#%(DIAMONDG-%&STYLEG6#%&PO INTG-F$6$7SF'7$$\"1LLLL3VfV!#<$\"14*>/Y-64\"F77$$\"1nmm\"H[D:)FR$\"1LH ()zm4x6F77$$\"1LLLe0$=C\"F.$\"1]!=bG\"F77$$\"1LLL3RBr;F.$\"1<\"[t#y )oR\"F77$$\"1mm;zjf)4#F.$\"1=N>kV`@:F77$$\"1LL$e4;[\\#F.$\"1(>)*Qz7qk \"F77$$\"1++]i'y]!HF.$\"1d:o*[`yy\"F77$$\"1LL$ezs$HLF.$\"16!*=P=>Y>F77 $$\"1++]7iI_PF.$\"1rw@#pwz6#F77$$\"1nmm;_M(=%F.$\"1_y5l<^5BF77$$\"1LLL 3y_qXF.$\"1[r%)QIa%\\#F77$$\"1+++]1!>+&F.$\"1R'\\`K:$>FF77$$\"1+++]Z/N aF.$\"1GmLh6RlHF77$$\"1+++]$fC&eF.$\"1'op)HydBKF77$$\"1LL$ez6:B'F.$\"1 uxPa2YxMF77$$\"1mmm;=C#o'F.$\"1VE()3P]0QF77$$\"1mmmm#pS1(F.$\"1\"ovC0( \\2TF77$$\"1++]i`A3vF.$\"1_QUTy1*[%F77$$\"1mmmm(y8!zF.$\"19x3mZHc[F77$ $\"1++]i.tK$)F.$\"1\"o^Vg^QH&F77$$\"1++](3zMu)F.$\"1&=`d^-ru&F77$$\"1n mm\"H_?<*F.$\"1h)[jnT9E'F77$$\"1nm;zihl&*F.$\"1M/5%R8Ux'F77$$\"1LLL3#G ,***F.$\"1j)yl%=[utF77$$\"1LLezw5V5F7$\"1X$HM\\uV0)F77$$\"1++v$Q#\\\"3 \"F7$\"1S!o;ldqp)F77$$\"1LL$e\"*[H7\"F7$\"1`**fm_*)[%*F77$$\"1+++qvxl6 F7$\"1o\")egPRH5!#97$$\"1++]_qn27F7$\"1f7GTLP>6Ffv7$$\"1++Dcp@[7F7$\"1 rS:\\F\"R@\"Ffv7$$\"1++]2'HKH\"F7$\"1%\\/3sl#G8Ffv7$$\"1nmmwanL8F7$\"1 *\\qhww,W\"Ffv7$$\"1+++v+'oP\"F7$\"1Qr()>L4q:Ffv7$$\"1LLeR<*fT\"F7$\"1 bL%[^5zp\"Ffv7$$\"1+++&)Hxe9F7$\"1G(poB%e\\=Ffv7$$\"1mm\"H!o-*\\\"F7$ \"1?UJd![Y+#Ffv7$$\"1++DTO5T:F7$\"1Nn0Q#[1=#Ffv7$$\"1nmmT9C#e\"F7$\"1x K_N(\\wO#Ffv7$$\"1++D1*3`i\"F7$\"1#QdM%yi!e#Ffv7$$\"1LLL$*zym;F7$\"1Js qoF%Q!GFfv7$$\"1LL$3N1#4LFfv7$$\"1+++q(G**y\"F7$\"1\\SkmI%oe$Ffv7$$\"1nm;9@BM=F7$\"1s+9;q ;>RFfv7$$\"1LLL`v&Q(=F7$\"1\")*[un.CC%Ffv7$$\"1++DOl5;>F7$\"1LKT&3fkh% Ffv7$$\"1++v.Uac>F7$\"1T72\\)H`+&Ffv7$F/$\"1CWJ.]\")faFfv-%&COLORG6&F@ F(F(F*-%+AXESLABELSG6$%\"xG%!G-%%VIEWG6$;F(F/%(DEFAULTG" 2 304 304 304 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 108 "You can copy from the above all the text in red and chan ge the value of step size h for your own experiment." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 136 "We now decrese the step size h and re-run the programme. It is necessary to \+ reset the values of oldx and oldy to the initial condition: " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "h:=0.1; # the step size" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"hG$\"\"\"!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "oldx:=0; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "oldy:=1 ; # the initial condition" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG \"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "L:=[];" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7\"" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 23 "L:=[op(L),[oldx,oldy]];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7#7$\"\"!\"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "slope:=(x,y) -> 2*y;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&slopeG:6$%\"xG%\"yG6\"6$%)operatorG%&arrowGF),$9%\"\"#F)F)" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 "while oldx < lastx do oldy:= oldy+h*slope(oldx,oldy);oldx:=oldx+h;L:=[op(L),[oldx,oldy]];od;" }} {PARA 11 "" 0 "" {TEXT -1 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ol dyG$\"#7!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"\"!\"\" " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7$7$\"\"!\"\"\"7$$F(!\"\"$\" #7F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"$W\"!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"#!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7%7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\" !\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"%G%%oldxG$\"\"$!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7&7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\" !\"#7$$\"\"$F+$\"%G%%oldyG$ \"&O2#!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"%!\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7'7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+ 7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G%%oldyG$\"'K)[#!\"&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"&!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG7(7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$ F+$\"%G%%oldyG$\"(%)f)H!\"'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"'!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"LG7)7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$ \"%G%%oldyG$\")3=$e$!\"(" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\"(!\"\"" }}{PARA 12 "" 1 " " {XPPMATH 20 "6#>%\"LG7*7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W \"!\"#7$$\"\"$F+$\"%G%%oldyG$\"*'p\")*H%!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"\")!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG7+7 $\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G%%oldyG$\"+_.yf^!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%%oldxG$\"\"*!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG7,7$\" \"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G%%oldyG$\"+Akt\">'!\"*" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%%oldxG$\"#5!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG7-7$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\" \"$F+$\"%G'Fgn" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"+1P3Iu!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$ \"#6!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG7.7$\"\"!\"\"\"7$$F (!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG $\"+Z/5;*)!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#7!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG7/7$\"\"!\"\"\"7$$F(!\"\"$\"#7 F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\" #6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %%oldyG$\"+a?$*p5!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#8 !\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG707$\"\"!\"\"\"7$$F(!\" \"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fg n7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$$\"#8F+$\"+a?$*p5FW" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"+l%=RG\"!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#9!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#> %\"LG717$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$ \"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$ $\"#8F+$\"+a?$*p5FW7$$\"#9F+$\"+l%=RG\"FW" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"+e@qS:!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%%oldxG$\"#:!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG727$\"\" !\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$$\"#8F+$\"+a?$ *p5FW7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e@qS:FW" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"+!fU)[=!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#;!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG737$\" \"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$$\"#8F+$\"+a ?$*p5FW7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e@qS:FW7$$\"#;F+$\"+!fU)[=FW " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldyG$\"+36h=A!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#%\"LG747$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\" !\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F, $\"+Z/5;*)Fgn7$$\"#8F+$\"+a?$*p5FW7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e @qS:FW7$$\"#;F+$\"+!fU)[=FW7$$\"#%%oldyG$\"+ILLiE!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%%oldxG$\"#=!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG757$\"\" !\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$$\"#8F+$\"+a?$ *p5FW7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e@qS:FW7$$\"#;F+$\"+!fU)[=FW7$ $\"#%%oldyG$\"+'***z%>$!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%o ldxG$\"#>!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"LG767$\"\"!\"\" \"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$\"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$$\"#8F+$\"+a?$*p5FW 7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e@qS:FW7$$\"#;F+$\"+!fU)[=FW7$$\"#< F+$\"+36h=AFW7$$\"#=F+$\"+ILLiEFW7$$\"#>F+$\"+'***z%>$FW" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%%oldyG$\"+&**fP$Q!\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%oldxG$\"#?!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#> %\"LG777$\"\"!\"\"\"7$$F(!\"\"$\"#7F+7$$\"\"#F+$\"$W\"!\"#7$$\"\"$F+$ \"%G'Fgn7$$\"#6F+$\"+1P3IuFgn7$F,$\"+Z/5;*)Fgn7$ $\"#8F+$\"+a?$*p5FW7$$\"#9F+$\"+l%=RG\"FW7$$\"#:F+$\"+e@qS:FW7$$\"#;F+ $\"+!fU)[=FW7$$\"#F+$\"+'***z %>$FW7$$\"#?F+$\"+&**fP$QFW" }}}{EXCHG {PARA 12 "" 1 "" {TEXT -1 70 "T he array L gives the solution obtained by Euler's method in the form: " }}{PARA 0 "" 0 "" {TEXT -1 54 "L=[ [x_0, y_0], [x_1, y_1], [x_2, y _2], .... ]." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 12 "" 1 "" {TEXT -1 0 "" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "picture3:=plot(L,style=poi nt,symbol=circle );" }{XPPMATH 20 "6#>%)picture3G-%/INTERFACE_PLOTG6%- %'CURVESG6$777$\"\"!$\"\"\"F-7$$\"1+++++++5!#;$\"1+++++++7!#:7$$\"1+++ ++++?F3$\"1++++++S9F67$$\"1+++++++IF3$\"1++++++G'F67$$\"1+++++++6F6$\"1+ ++1P3IuF67$F4$\"1+++Z/5;*)F67$$\"1+++++++8F6$\"1+++a?$*p5!#97$$\"1++++ +++9F6$\"1+++l%=RG\"Fio7$$\"1+++++++:F6$\"1+++e@qS:Fio7$$\"1+++++++;F6 $\"1+++!fU)[=Fio7$$\"1+++++++F6$\"1+++'***z%>$Fio7$$\"\"#F-$\"1+++&**fP$QFio- %'COLOURG6&%$RGBG$\"#5!\"\"F-F--%'SYMBOLG6#%'CIRCLEG-%&STYLEG6#%&POINT G" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "display(\{picture1,pic ture2,picture3\});" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 13 "" 1 " " {INLPLOT "6'-%'CURVESG6&7'7$\"\"!$\"\"\"F(7$$\"1+++++++]!#;$\"\"#F(7 $F)$\"\"%F(7$$\"1+++++++:!#:$\"\")F(7$F/$\"#;F(-%'COLOURG6&%$RGBG$\"#5 !\"\"F(F(-%'SYMBOLG6#%(DIAMONDG-%&STYLEG6#%&POINTG-F$6&77F'7$$\"1+++++ ++5F.$\"1+++++++7F77$$\"1+++++++?F.$\"1++++++S9F77$$\"1+++++++IF.$\"1+ +++++G'F77$$\" 1+++++++6F7$\"1+++1P3IuF77$FR$\"1+++Z/5;*)F77$$\"1+++++++8F7$\"1+++a?$ *p5!#97$$\"1+++++++9F7$\"1+++l%=RG\"Fdq7$F5$\"1+++e@qS:Fdq7$$\"1++++++ +;F7$\"1+++!fU)[=Fdq7$$\"1+++++++F7$\"1+++'***z%>$Fdq7$F/$\"1+++&**fP$QFdqF=- FE6#%'CIRCLEGFH-F$6$7SF'7$$\"1LLLL3VfV!#<$\"14*>/Y-64\"F77$$\"1nmm\"H[ D:)F]t$\"1LH()zm4x6F77$$\"1LLLe0$=C\"F.$\"1]!=bG\"F77$$\"1LLL3RBr;F .$\"1<\"[t#y)oR\"F77$$\"1mm;zjf)4#F.$\"1=N>kV`@:F77$$\"1LL$e4;[\\#F.$ \"1(>)*Qz7qk\"F77$$\"1++]i'y]!HF.$\"1d:o*[`yy\"F77$$\"1LL$ezs$HLF.$\"1 6!*=P=>Y>F77$$\"1++]7iI_PF.$\"1rw@#pwz6#F77$$\"1nmm;_M(=%F.$\"1_y5l<^5 BF77$$\"1LLL3y_qXF.$\"1[r%)QIa%\\#F77$$\"1+++]1!>+&F.$\"1R'\\`K:$>FF77 $$\"1+++]Z/NaF.$\"1GmLh6RlHF77$$\"1+++]$fC&eF.$\"1'op)HydBKF77$$\"1LL$ ez6:B'F.$\"1uxPa2YxMF77$$\"1mmm;=C#o'F.$\"1VE()3P]0QF77$$\"1mmmm#pS1(F .$\"1\"ovC0(\\2TF77$$\"1++]i`A3vF.$\"1_QUTy1*[%F77$$\"1mmmm(y8!zF.$\"1 9x3mZHc[F77$$\"1++]i.tK$)F.$\"1\"o^Vg^QH&F77$$\"1++](3zMu)F.$\"1&=`d^- ru&F77$$\"1nmm\"H_?<*F.$\"1h)[jnT9E'F77$$\"1nm;zihl&*F.$\"1M/5%R8Ux'F7 7$$\"1LLL3#G,***F.$\"1j)yl%=[utF77$$\"1LLezw5V5F7$\"1X$HM\\uV0)F77$$\" 1++v$Q#\\\"3\"F7$\"1S!o;ldqp)F77$$\"1LL$e\"*[H7\"F7$\"1`**fm_*)[%*F77$ $\"1+++qvxl6F7$\"1o\")egPRH5Fdq7$$\"1++]_qn27F7$\"1f7GTLP>6Fdq7$$\"1++ Dcp@[7F7$\"1rS:\\F\"R@\"Fdq7$$\"1++]2'HKH\"F7$\"1%\\/3sl#G8Fdq7$$\"1nm mwanL8F7$\"1*\\qhww,W\"Fdq7$$\"1+++v+'oP\"F7$\"1Qr()>L4q:Fdq7$$\"1LLeR <*fT\"F7$\"1bL%[^5zp\"Fdq7$$\"1+++&)Hxe9F7$\"1G(poB%e\\=Fdq7$$\"1mm\"H !o-*\\\"F7$\"1?UJd![Y+#Fdq7$$\"1++DTO5T:F7$\"1Nn0Q#[1=#Fdq7$$\"1nmmT9C #e\"F7$\"1xK_N(\\wO#Fdq7$$\"1++D1*3`i\"F7$\"1#QdM%yi!e#Fdq7$$\"1LLL$*z ym;F7$\"1JsqoF%Q!GFdq7$$\"1LL$3N1#4LFdq7$$\"1+++q(G**y\"F7$\"1\\SkmI%oe$Fdq7$$\"1nm;9@BM =F7$\"1s+9;q;>RFdq7$$\"1LLL`v&Q(=F7$\"1\")*[un.CC%Fdq7$$\"1++DOl5;>F7$ \"1LKT&3fkh%Fdq7$$\"1++v.Uac>F7$\"1T72\\)H`+&Fdq7$F/$\"1CWJ.]\")faFdq- %&COLORG6&F@F(F(F*-%+AXESLABELSG6$%\"xG%!G-%%VIEWG6$;F(F/%(DEFAULTG" 2 304 304 304 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 250 76 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 191 "The circlces give the approximation with step size h=0.1. It moves closer to the true solution (the blue curve ) compared to the the diamonds, which is the approximation with step s ize h=0.25." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "0 4 0" 111 }{VIEWOPTS 1 1 0 1 1 1803 }