{VERSION 5 0 "IBM INTEL NT" "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 "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 "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 } } {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 16 "Class 12/15/2004" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 36 "Predator- Prey problem in vector form" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "f := Vector(2): #RHS vecto r function" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "f[1] := (t,x, y)-> x - 0.1*x*y;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"fG6#\"\"\"f* 6%%\"tG%\"xG%\"yG6\"6$%)operatorG%&arrowGF-,&9%F'*($F'!\"\"F'F2F'9&F'F 5F-F-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "f[2] := (t,x,y)- > 0.02*x*y - 0.5*y;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"fG6#\"\"#f *6%%\"tG%\"xG%\"yG6\"6$%)operatorG%&arrowGF-,&*($F'!\"#\"\"\"9%F59&F5F 5*&$\"\"&!\"\"F5F7F5F;F-F-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "x0 := Vector([20,4]); # initial values" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x0G-%'RTABLEG6%\")S)ev\"-%'MATRIXG6#7$7#\"#?7#\"\"%& %'VectorG6#%'columnG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "a, \+ b := 0, 10; # interval domain" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>6 $%\"aG%\"bG6$\"\"!\"#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "n \+ := 50:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "read(\"f:/340/rkm xm.txt\"):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "t,x := rkmxm( f,a,b,x0,n);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%L~End~of~Runge-Kutta~ method~for~mxm~IVP-ODE~G" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>6$%\"tG% \"xG6$-%'RTABLEG6*\")!)=c<%)anythingG&%'VectorG6#%'columnG%,rectangula rG%.Fortran_orderG7\"\"\"\";F4\"#^-F)6*\")g@c " 0 "" {MPLTEXT 1 0 39 "pts1 := [seq([t[j],x[1][j] ],j=1..n+1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "plot(pts1) ;" }}{PARA 13 "" 1 "" {GLPLOT2D 307 156 156 {PLOTDATA 2 "6%-%'CURVESG6 $7U7$$\"\"!F)$\"#?F)7$$\"35+++++++?!#=$\"39+++^9\\cA!#;7$$\"3A+++++++S F/$\"39+++tQ\"ya#F27$$\"3w**************fF/$\"37+++!)fGwGF27$$\"3U++++ +++!)F/$\"3.+++;I>VKF27$$\"\"\"F)$\"3#)*****RqV![OF27$$\"3%*********** ***>\"!#<$\"3!******zH;v3%F27$$\"3!**************R\"FJ$\"3M+++\"RFSb%F 27$$\"33+++++++;FJ$\"3#)*****H)4\"R.&F27$$\"3/+++++++=FJ$\"3))*****R() [`]&F27$$\"\"#F)$\"3s*****zFHl$fF27$$\"3;+++++++AFJ$\"3%)*****41/^G'F2 7$$\"3!**************R#FJ$\"3J******\\ry+lF27$$\"33+++++++EFJ$\"3u**** *p7FL`'F27$$\"3#)*************z#FJ$\"3')*****\\>qnM'F27$$\"\"$F)$\"3F+ ++2pcOfF27$$\"3;+++++++KFJ$\"3M+++*R$>S`F27$$\"3!**************R$FJ$\" 3N+++h\\;JYF27$$\"33+++++++OFJ$\"3K+++S\"fh*QF27$$\"3#)*************z$ FJ$\"3o*******3K$3KF27$$\"\"%F)$\"3\"*******p-\">h#F27$$\"3;+++++++UFJ $\"3%)******[P>A@F27$$\"3M+++++++WFJ$\"3++++'*ohMj?\"F27$$ \"\"&F)$\"30+++ibML5F27$$\"3;+++++++_FJ$\"3U*****>Vv>.*FJ7$$\"3M++++++ +aFJ$\"3!*******zyAf!)FJ7$$\"3k*************f&FJ$\"3')*****z^v/M(FJ7$$ \"3#)*************z&FJ$\"3d*****4573#oFJ7$$\"\"'F)$\"3')******G$[0Y'FJ 7$$\"3;+++++++iFJ$\"3w******fWYJiFJ7$$\"3M+++++++kFJ$\"39+++'GmR6'FJ7$ $\"3k*************f'FJ$\"3k*****zx4]4'FJ7$$\"3#)*************z'FJ$\"35 +++KnimhFJ7$$\"\"(F)$\"3;+++:?!\\K'FJ7$$\"3;+++++++sFJ$\"3g*****fC%Gpl FJ7$$\"3M+++++++uFJ$\"3t******=:9-pFJ7$$\"3k*************f(FJ$\"34+++K v\\GtFJ7$$\"3#)*************z(FJ$\"3#******\\$3!f&yFJ7$$\"\")F)$\"3O** ***H[\"R%\\)FJ7$$\"3G*************>)FJ$\"3C+++6i`c#*FJ7$$\"3M+++++++%) FJ$\"3%******fd]d,\"F27$$\"3k*************f)FJ$\"31+++[m^@6F27$$\"3q++ +++++))FJ$\"33+++Z_,X7F27$$\"\"*F)$\"3)*******[te)Q\"F27$$\"3G******** *****>*FJ$\"3++++y&R[b\"F27$$\"3M+++++++%*FJ$\"3')*****pm5mu\"F27$$\"3 k*************f*FJ$\"3=+++:M!p'>F27$$\"3q+++++++)*FJ$\"3(******pyj(=AF 27$$\"#5F)$\"35+++>630DF2-%'COLOURG6&%$RGBG$F`[l!\"\"F(F(-%+AXESLABELS G6$Q!6\"F\\\\l-%%VIEWG6$%(DEFAULTGFa\\l" 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 39 "pts2 := [seq([t[j],x[2][j]],j=1..n+1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "plot(pts2);" }}{PARA 13 "" 1 "" {GLPLOT2D 307 184 184 {PLOTDATA 2 "6%-%'CURVESG6$7U7$$\"\"!F)$\"\"%F)7 $$\"35+++++++?!#=$\"3$)*****f)3_SR!#<7$$\"3A+++++++SF/$\"3*)*****z5dY# RF27$$\"3w**************fF/$\"35+++-xedRF27$$\"3U+++++++!)F/$\"3I+++Gw nYSF27$$\"\"\"F)$\"3q*****4zf@?%F27$$\"3%**************>\"F2$\"3\")*** **HT/!QWF27$$\"3!**************R\"F2$\"3I+++j=&Hx%F27$$\"33+++++++;F2$ \"30+++UncJ_F27$$\"3/+++++++=F2$\"3i*****zp5\\%eF27$$\"\"#F)$\"3))**** *4UQ*\\mF27$$\"3;+++++++AF2$\"35+++N2'fo(F27$$\"3!**************R#F2$ \"3%******\\Hjd)*)F27$$\"33+++++++EF2$\"33+++*pFf0\"!#;7$$\"3#)******* ******z#F2$\"32+++2>6P7Fio7$$\"\"$F)$\"3********z]6K9Fio7$$\"3;+++++++ KF2$\"3&)*****fw,Xi\"Fio7$$\"3!**************R$F2$\"3#******HZRZz\"Fio 7$$\"33+++++++OF2$\"3/+++'oKe#>Fio7$$\"3#)*************z$F2$\"35+++X@7 3?Fio7$F*$\"35+++Fio7$$\"3k*************f%F2$\"37+++C&33\">Fio7$ $\"3#)*************z%F2$\"33+++`aPA=Fio7$$\"\"&F)$\"3(******RG9Us\"Fio 7$$\"3;+++++++_F2$\"3))*****\\(>a@;Fio7$$\"3M+++++++aF2$\"3$******\\t1 \"=:Fio7$$\"3k*************f&F2$\"3,+++x'*\\;9Fio7$$\"3#)************* z&F2$\"3#******f8W%=8Fio7$$\"\"'F)$\"3,+++\"GK]A\"Fio7$$\"3;+++++++iF2 $\"34+++Mn!p8\"Fio7$$\"3M+++++++kF2$\"3%******RA&Ra5Fio7$$\"3k******** *****f'F2$\"3;*****\\_&3w(*F27$$\"3#)*************z'F2$\"3-+++;76l!*F2 7$$\"\"(F)$\"31+++L\\r4%)F27$$\"3;+++++++sF2$\"3')*****H'e)z!yF27$$\"3 M+++++++uF2$\"35+++FlndsF27$$\"3k*************f(F2$\"33+++F.RcnF27$$\" 3#)*************z(F2$\"3d*****pj=)F2$\"3W+++3p6BbF27$$\"3M+++++++%)F2$\"3S+++I(G^>&F27$ $\"3k*************f)F2$\"3M+++&Gmd!\\F27$$\"3q+++++++))F2$\"3u*****p6x Pl%F27$$\"\"*F)$\"3r******Q!R$QWF27$$\"3G*************>*F2$\"3++++'*f9 fUF27$$\"3M+++++++%*F2$\"3@+++J(=l6%F27$$\"3k*************f*F2$\"3M+++ V=b6SF27$$\"3q+++++++)*F2$\"35+++zBJYRF27$$\"#5F)$\"3&******>&)4T#RF2- %'COLOURG6&%$RGBG$F^[l!\"\"F(F(-%+AXESLABELSG6$Q!6\"Fj[l-%%VIEWG6$%(DE FAULTGF_\\l" 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 0 "" }}{PARA 0 "" 0 "" {TEXT -1 62 "One way to see the solution: plot two function as two cu rves." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "plot(\{pts1,pts2\});" }}{PARA 13 "" 1 "" {GLPLOT2D 307 182 182 {PLOTDATA 2 "6&-%'CURVESG6$7U7$$\"\"!F)$\"#?F)7$$\"35+++++ ++?!#=$\"39+++^9\\cA!#;7$$\"3A+++++++SF/$\"39+++tQ\"ya#F27$$\"3w****** ********fF/$\"37+++!)fGwGF27$$\"3U+++++++!)F/$\"3.+++;I>VKF27$$\"\"\"F )$\"3#)*****RqV![OF27$$\"3%**************>\"!#<$\"3!******zH;v3%F27$$ \"3!**************R\"FJ$\"3M+++\"RFSb%F27$$\"33+++++++;FJ$\"3#)*****H) 4\"R.&F27$$\"3/+++++++=FJ$\"3))*****R()[`]&F27$$\"\"#F)$\"3s*****zFHl$ fF27$$\"3;+++++++AFJ$\"3%)*****41/^G'F27$$\"3!**************R#FJ$\"3J* *****\\ry+lF27$$\"33+++++++EFJ$\"3u*****p7FL`'F27$$\"3#)*************z #FJ$\"3')*****\\>qnM'F27$$\"\"$F)$\"3F+++2pcOfF27$$\"3;+++++++KFJ$\"3M +++*R$>S`F27$$\"3!**************R$FJ$\"3N+++h\\;JYF27$$\"33+++++++OFJ$ \"3K+++S\"fh*QF27$$\"3#)*************z$FJ$\"3o*******3K$3KF27$$\"\"%F) $\"3\"*******p-\">h#F27$$\"3;+++++++UFJ$\"3%)******[P>A@F27$$\"3M+++++ ++WFJ$\"3++++'*ohMj?\"F27$$\"\"&F)$\"30+++ibML5F27$$\"3;++ +++++_FJ$\"3U*****>Vv>.*FJ7$$\"3M+++++++aFJ$\"3!*******zyAf!)FJ7$$\"3k *************f&FJ$\"3')*****z^v/M(FJ7$$\"3#)*************z&FJ$\"3d**** *4573#oFJ7$$\"\"'F)$\"3')******G$[0Y'FJ7$$\"3;+++++++iFJ$\"3w******fWY JiFJ7$$\"3M+++++++kFJ$\"39+++'GmR6'FJ7$$\"3k*************f'FJ$\"3k**** *zx4]4'FJ7$$\"3#)*************z'FJ$\"35+++KnimhFJ7$$\"\"(F)$\"3;+++:?! \\K'FJ7$$\"3;+++++++sFJ$\"3g*****fC%GplFJ7$$\"3M+++++++uFJ$\"3t******= :9-pFJ7$$\"3k*************f(FJ$\"34+++Kv\\GtFJ7$$\"3#)*************z(F J$\"3#******\\$3!f&yFJ7$$\"\")F)$\"3O*****H[\"R%\\)FJ7$$\"3G********** ***>)FJ$\"3C+++6i`c#*FJ7$$\"3M+++++++%)FJ$\"3%******fd]d,\"F27$$\"3k** ***********f)FJ$\"31+++[m^@6F27$$\"3q+++++++))FJ$\"33+++Z_,X7F27$$\"\" *F)$\"3)*******[te)Q\"F27$$\"3G*************>*FJ$\"3++++y&R[b\"F27$$\" 3M+++++++%*FJ$\"3')*****pm5mu\"F27$$\"3k*************f*FJ$\"3=+++:M!p' >F27$$\"3q+++++++)*FJ$\"3(******pyj(=AF27$$\"#5F)$\"35+++>630DF2-%'COL OURG6&%$RGBG$F`[l!\"\"F(F(-F$6$7U7$F(Fiq7$F-$\"3$)*****f)3_SRFJ7$F4$\" 3*)*****z5dY#RFJ7$F9$\"35+++-xedRFJ7$F>$\"3I+++GwnYSFJ7$FC$\"3q*****4z f@?%FJ7$FH$\"3\")*****HT/!QWFJ7$FN$\"3I+++j=&Hx%FJ7$FS$\"30+++UncJ_FJ7 $FX$\"3i*****zp5\\%eFJ7$Fgn$\"3))*****4UQ*\\mFJ7$F\\o$\"35+++N2'fo(FJ7 $Fao$\"3%******\\Hjd)*)FJ7$Ffo$\"33+++*pFf0\"F27$F[p$\"32+++2>6P7F27$F `p$\"3********z]6K9F27$Fep$\"3&)*****fw,Xi\"F27$Fjp$\"3#******HZRZz\"F 27$F_q$\"3/+++'oKe#>F27$Fdq$\"35+++X@73?F27$Fiq$\"35+++F27$Fhr$\"37+++C&33\">F27$F]s$\"33+ ++`aPA=F27$Fbs$\"3(******RG9Us\"F27$Fgs$\"3))*****\\(>a@;F27$F\\t$\"3$ ******\\t1\"=:F27$Fat$\"3,+++x'*\\;9F27$Fft$\"3#******f8W%=8F27$F[u$\" 3,+++\"GK]A\"F27$F`u$\"34+++Mn!p8\"F27$Feu$\"3%******RA&Ra5F27$Fju$\"3 ;*****\\_&3w(*FJ7$F_v$\"3-+++;76l!*FJ7$Fdv$\"31+++L\\r4%)FJ7$Fiv$\"3') *****H'e)z!yFJ7$F^w$\"35+++FlndsFJ7$Fcw$\"33+++F.RcnFJ7$Fhw$\"3d*****p j=&FJ 7$F\\y$\"3M+++&Gmd!\\FJ7$Fay$\"3u*****p6xPl%FJ7$Ffy$\"3r******Q!R$QWFJ 7$F[z$\"3++++'*f9fUFJ7$F`z$\"3@+++J(=l6%FJ7$Fez$\"3M+++V=b6SFJ7$Fjz$\" 35+++zBJYRFJ7$F_[l$\"3&******>&)4T#RFJ-Fd[l6&Ff[lF(Fg[lF(-%+AXESLABELS G6$Q!6\"Fhel-%%VIEWG6$%(DEFAULTGF]fl" 1 2 0 1 10 0 2 9 1 4 2 1.000000 46.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 43 "Another way to see the solution: plot (x,y)" }} {PARA 0 "" 0 "" {TEXT -1 30 "(also called parametric curve)" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "p ts:= [seq([x[1][j],x[2][j]],j=1..n+1)]:\n" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 10 "plot(pts);" }}{PARA 13 "" 1 "" {GLPLOT2D 347 220 220 {PLOTDATA 2 "6%-%'CURVESG6$7U7$$\"#?\"\"!$\"\"%F*7$$\"39+++^9\\cA! #;$\"3$)*****f)3_SR!#<7$$\"39+++tQ\"ya#F0$\"3*)*****z5dY#RF37$$\"37+++ !)fGwGF0$\"35+++-xedRF37$$\"3.+++;I>VKF0$\"3I+++GwnYSF37$$\"3#)*****Rq V![OF0$\"3q*****4zf@?%F37$$\"3!******zH;v3%F0$\"3\")*****HT/!QWF37$$\" 3M+++\"RFSb%F0$\"3I+++j=&Hx%F37$$\"3#)*****H)4\"R.&F0$\"30+++UncJ_F37$ $\"3))*****R()[`]&F0$\"3i*****zp5\\%eF37$$\"3s*****zFHl$fF0$\"3))***** 4UQ*\\mF37$$\"3%)*****41/^G'F0$\"35+++N2'fo(F37$$\"3J******\\ry+lF0$\" 3%******\\Hjd)*)F37$$\"3u*****p7FL`'F0$\"33+++*pFf0\"F07$$\"3')*****\\ >qnM'F0$\"32+++2>6P7F07$$\"3F+++2pcOfF0$\"3********z]6K9F07$$\"3M+++*R $>S`F0$\"3&)*****fw,Xi\"F07$$\"3N+++h\\;JYF0$\"3#******HZRZz\"F07$$\"3 K+++S\"fh*QF0$\"3/+++'oKe#>F07$$\"3o*******3K$3KF0$\"35+++X@73?F07$$\" 3\"*******p-\">h#F0$\"35+++A@F0$\"3'******H?'3 H?F07$$\"3++++'*ohMF07$$\"32+++sV\"\\V\"F0$\"37+++C&3 3\">F07$$\"32+++!H>j?\"F0$\"33+++`aPA=F07$$\"30+++ibML5F0$\"3(******RG 9Us\"F07$$\"3U*****>Vv>.*F3$\"3))*****\\(>a@;F07$$\"3!*******zyAf!)F3$ \"3$******\\t1\"=:F07$$\"3')*****z^v/M(F3$\"3,+++x'*\\;9F07$$\"3d***** 4573#oF3$\"3#******f8W%=8F07$$\"3')******G$[0Y'F3$\"3,+++\"GK]A\"F07$$ \"3w******fWYJiF3$\"34+++Mn!p8\"F07$$\"39+++'GmR6'F3$\"3%******RA&Ra5F 07$$\"3k*****zx4]4'F3$\"3;*****\\_&3w(*F37$$\"35+++KnimhF3$\"3-+++;76l !*F37$$\"3;+++:?!\\K'F3$\"31+++L\\r4%)F37$$\"3g*****fC%GplF3$\"3')**** *H'e)z!yF37$$\"3t******=:9-pF3$\"35+++FlndsF37$$\"34+++Kv\\GtF3$\"33++ +F.RcnF37$$\"3#******\\$3!f&yF3$\"3d*****pj=&F37$$\"31+++[m^@6F0$\"3M+++&Gmd!\\F37$$\"33+++Z_,X7F 0$\"3u*****p6xPl%F37$$\"3)*******[te)Q\"F0$\"3r******Q!R$QWF37$$\"3+++ +y&R[b\"F0$\"3++++'*f9fUF37$$\"3')*****pm5mu\"F0$\"3@+++J(=l6%F37$$\"3 =+++:M!p'>F0$\"3M+++V=b6SF37$$\"3(******pyj(=AF0$\"35+++zBJYRF37$$\"35 +++>630DF0$\"3&******>&)4T#RF3-%'COLOURG6&%$RGBG$\"#5!\"\"$F*F*Fj[l-%+ AXESLABELSG6$Q!6\"F^\\l-%%VIEWG6$%(DEFAULTGFc\\l" 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 0 "" }}{PARA 0 "" 0 "" {TEXT -1 8 "Page 186" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "f := Vector(2): # RHS" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "f[1] \+ := (t,x,y) -> -y + cos(t)*sin(t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> &%\"fG6#\"\"\"f*6%%\"tG%\"xG%\"yG6\"6$%)operatorG%&arrowGF-,&9&!\"\"*& -%$cosG6#9$F'-%$sinGF7F'F'F-F-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "f[2] := (t,x,y) -> x + (sin(t))^2;" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>&%\"fG6#\"\"#f*6%%\"tG%\"xG%\"yG6\"6$%)operatorG%&ar rowGF-,&9%\"\"\"*$)-%$sinG6#9$F'F3F3F-F-F-" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 38 "x0 := Vector([0,0]); # initial values" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x0G-%'RTABLEG6%\")_&H#>-%'MATRIXG6#7$7#\"\"! F-&%'VectorG6#%'columnG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 " a, b := 0, evalf(2*Pi);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>6$%\"aG%\" bG6$\"\"!$\"+3`=$G'!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "n := 100:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "t, x := rkmxm(f ,a,b,x0,n);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%L~End~of~Runge-Kutta~m ethod~for~mxm~IVP-ODE~G" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>6$%\"tG%\" xG6$-%'RTABLEG6*\")#fH#>%)anythingG&%'VectorG6#%'columnG%,rectangularG %.Fortran_orderG7\"\"\"\";F4\"$,\"-F)6*\")#zH#>F,F-F1F2F3F4;F4\"$-#" } }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 35 " First way to draw graph: two curves" }}{PARA 0 "" 0 "" {TEXT -1 0 "" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "curve1 := [seq([t[j],x[1][ j]],j=1..n+1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "curve2 : = [seq([t[j],x[2][j]],j=1..n+1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "plot(\{curve1,curve2\});" }}{PARA 13 "" 1 "" {GLPLOT2D 347 256 256 {PLOTDATA 2 "6&-%'CURVESG6$7aq7$$\"\"!F)F(7$$\"3 )******zI&=$G'!#>$\"39+++;vPp>!#?7$$\"3-+++iqjc7!#=$\"3N+++c%=J#yF07$$ \"3'******>fb\\)=F4$\"31+++r,!*RF47$$\"3[++++49o\")F4$\"35+++iEUf@F47$$\"3#******4Vfk z)F4$\"3\"******zQX6J#F47$$\"3^+++izxC%*F4$\"33+++Sj$HU#F47$$\"3(***** **['4`+\"!#<$\"35+++*>jr[#F47$$\"3-+++-:9o5Fep$\"3'******HZpm\\#F47$$ \"32+++bL(48\"Fep$\"3$******4h6\\W#F47$$\"3'*******3_!Q>\"Fep$\"3$**** **\\U(3EBF47$$\"3/+++iqjc7Fep$\"3\"*******>PDN@F47$$\"33+++:*o%>8Fep$ \"3'*******[6Vo=F47$$\"3!******zw+BQ\"Fep$\"3++++dXpA:F47$$\"3'******4 iK^W\"Fep$\"3/+++gsC'4\"F47$$\"3++++uW'z]\"Fep$\"3++++%))zZ)eF-7$$\"3/ +++Fjzq:Fep$!3!)***********f8&!#D7$$\"35+++!=GOj\"Fep$!3O*****fr>Ln'F- 7$$\"3#******H.gkp\"Fep$!38+++`lT59F47$$\"3)******f)=Hf XI\\A#F47$$\"3-+++RP7A=Fep$!3@+++'*[O0JF47$$\"31+++#fb\\)=Fep$!3?+++kS 3XSF47$$\"3!******\\W(yZ>Fep$!34+++J6SO]F47$$\"3.+++*H>1,#Fep$!3S+++W8 nqgF47$$\"33+++_6Xt?Fep$!3q*****4x+%QrF47$$\"37+++0IGO@Fep$!3l******> \"p$H#)F47$$\"3;+++e[6*>#Fep$!3*)*****H![wK$*F47$$\"3A+++6n%>E#Fep$!3) ******4-LP/\"Fep7$$\"3$)*****RcyZK#Fep$!3%******z\\^J:\"Fep7$$\"3')*** **pT5wQ#Fep$!3%******\\]j.E\"Fep7$$\"3#*******pAW]CFep$!33+++'[.UO\"Fe p7$$F>Fep$!3%******z/DNY\"Fep7$$\"3-+++wf5wDFep$!33+++-r@d:Fep7$$\"32+ ++Hy$*QEFep$!3'******fe>Uk\"Fep7$$\"37+++#op[kb&=Fep7$$\"3!******>Cl-*GFep$!31+++mft1>Fep7$$\"3%******\\4(4`HF ep$!3!*******f\\Fep7$$\"3++++[*Gf,$Fep$!3/+++bdSw>Fep7$$\"3/+++,3wy IFep$!3-+++mN3%*>Fep7$$\"35+++aEfTJFep$!3)******f]*****>Fep7$$\"39+++2 XU/KFep$!3!*******4O3%*>Fep7$$\"3?+++gjDnKFep$!3++++VeSw>Fep7$$\"3!)** ***H@)3ILFep$!33+++\"4vr%>Fep7$$\"3%)*****f1?HR$Fep$!3-+++Sht1>Fep7$$ \"3!*******=>vbMFep$!35+++*pkb&=Fep7$$\"3%******>x$e=NFep$!3!******f`g Uz\"Fep7$$\"33+++EcT\"e$Fep$!3$*******Q)QNs\"Fep7$$\"37+++zuCWOFep$!3( ******p#*>Uk\"Fep7$$\"3<+++K$zqq$Fep$!3)******>[n,Z)eF-7$$\"3C+++(o_!Q[Fep$\"3#*******G%Ri4\"F47$$\"3%)******R X)3!\\Fep$\"3.+++[noA:F47$$\"3M+++$R;P'\\Fep$\"37+++#RB%o=F47$$FRFep$ \"3,+++YgCN@F47$$\"3W+++*4!Q*3&Fep$\"3\"******p')zgK#F47$$\"37+++`>@_^ Fep$\"3')*******>/\\W#F47$$\"3t*****f!Q/:_Fep$\"3#*******RAm'\\#F47$$ \"3A+++fc(yF&Fep$\"3-+++xh:([#F47$$\"3#)*****>^22M&Fep$\"3'*******e&HH U#F47$$\"3K+++l$RNS&Fep$\"3#******p()Q6J#F47$$\"3#******z@rjY&Fep$\"3) )******\\kTf@F47$$\"3T+++rI?HbFep$\"35+++&*osv>F47$$\"3-+++C\\.#f&Fep$ \"3#******>`>#o1fFep$\"3z******p .=HlF-7$$\"3S+++Vg-pfFep$\"3!******4oqZl%F-7$$\"3++++'*y&=.'Fep$\"3/++ +-l&H/$F-7$$\"3g******[(*o%4'Fep$\"3%)*****Rv!))RF07$$\"3?+++3`= $G'Fep$!3#*******[!pN&QFes-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-F$6$7aqF'7$ F+$\"3.+++=a4R7!#@7$F2$\"33+++GS.$))*Fi[m7$F8$\"31+++6y0>LF07$F=$\"3w+ ++S]28yF07$FB$\"3&******RuRC^\"F-7$FG$\"3++++)=.^e#F-7$FL$\"3u*****>Pq A0%F-7$FQ$\"3)*******\\W(*efF-7$FV$\"3-+++b*H8M)F-7$Fen$\"3*******\\?q D7\"F47$Fjn$\"3-+++UQ!GY\"F47$F_o$\"31+++GvLb=F47$Fdo$\"3)******Rj_&*H #F47$Fio$\"3%******p/'p$z#F47$F^p$\"3++++%3()[L$F47$Fcp$\"39+++YN9>RF4 7$Fip$\"3()*****4aE9a%F47$F^q$\"3I+++#[.d>&F47$Fcq$\"3m*****>HG](eF47$ Fhq$\"3]+++[wjrlF47$F]r$\"3S+++n<1xsF47$Fbr$\"3O+++]![A)zF47$Fgr$\"3[* ****zm&px')F47$F\\s$\"3I+++,!*f`$*F47$Fas$\"2)*****4p(******Fep7$Fgs$ \"31+++'4$pg5Fep7$F\\t$\"34+++\"Qfk6\"Fep7$Fat$\"35+++]#\\j;\"Fep7$Fft $\"30+++6(f%47Fep7$F[u$\"3.+++\\)[\\C\"Fep7$F`u$\"3'******Hu\\?F\"Fep7 $Feu$\"3#*******yL3!H\"Fep7$Fju$\"3/+++_.Z)H\"Fep7$F_v$\"3++++q6u'H\"F ep7$Fdv$\"3#******\\)\\a%G\"Fep7$Fiv$\"3!*******>mlh7Fep7$F^w$\"3++++l <)zA\"Fep7$Fcw$\"3%******4FgN=\"Fep7$Fhw$\"34+++[ucG6Fep7$F]x$\"3++++z LJj5Fep7$Fax$\"3G+++(H-C))*F47$Ffx$\"3O+++&)e(4HgF47$Fjy$\"3w******e(oc* [F47$F_z$\"3t*****4(zV9PF47$Fdz$\"33+++%Q%y'\\#F47$Fiz$\"33+++GXda7F47 $F^[l$\"33+++++I5PF_cl7$Fc[l$!3!*******=rca7F47$Fh[l$!3!******>,xn\\#F 47$F]\\l$!3!)******e1V9PF47$Fb\\l$!3?+++K:m&*[F47$Fg\\l$!3e*****>2!4Hg F47$F\\]l$!3++++UX(R5(F47$Fa]l$!3_+++!G_.6)F47$Ff]l$!3y*****fBp\"R!*F4 7$F[^l$!3w******\\eR#))*F47$F`^l$!31+++dFJj5Fep7$Fe^l$!3)******4&ocG6F ep7$Fj^l$!3#******>qfN=\"Fep7$F__l$!3/+++E7)zA\"Fep7$Fd_l$!3)******46c ;E\"Fep7$Fi_l$!3#*******4Xa%G\"Fep7$F^`l$!33+++J2u'H\"Fep7$Fc`l$!35+++ ]*p%)H\"Fep7$Fh`l$!3++++:I3!H\"Fep7$F]al$!35+++?%\\?F\"Fep7$Fbal$!3&** *****o&[\\C\"Fep7$Fgal$!3)******RZf%47Fep7$F\\bl$!3!******p0\\j;\"Fep7 $Fabl$!3++++N#fk6\"Fep7$Ffbl$!33+++(*Hpg5Fep7$F[cl$!2-++?<(******Fep7$ Facl$!3a*****\\'*)f`$*F47$Ffcl$!3;+++>hpx')F47$F[dl$!3U+++!**[A)zF47$F `dl$!3#******p>jqF(F47$Fedl$!3W+++m&R;d'F47$Fidl$!3m*****RpI](eF47$F^e l$!3)******HO1d>&F47$Fcel$!3v*****R*)H9a%F47$Fhel$!3D+++et9>RF47$F]fl$ !3!******HM\"*[L$F47$Fbfl$!3)******pt+Pz#F47$Fgfl$!3-+++Sxb*H#F47$F\\g l$!37+++IIMb=F47$Fagl$!33+++?(4GY\"F47$Ffgl$!3)******fVwD7\"F47$F[hl$! 3w*****H`&RT$)F-7$F`hl$!3I+++HI/ffF-7$Fehl$!3'******pmTB0%F-7$Fjhl$!3% )*****>!p<&e#F-7$F_il$!3%*******ob^7:F-7$Fdil$!3\"******\\;^Q\"yF07$Fi il$!31+++#[[)>LF07$F^jl$!3I*****z<_5*)*Fi[m7$Fcjl$!3))*****fP!>Z7Fi[m7 $Fhjl$!3))*****R)HvO\")F_cl-F][m6&F_[mF(F`[mF(-%+AXESLABELSG6$Q!6\"Fh^ n-%%VIEWG6$%(DEFAULTGF]_n" 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 0 "" }}{PARA 0 "" 0 "" {TEXT -1 52 "Second way to see the solution: th e parametric curve" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "pcurve := [seq([x[1][j],x[2][j]],j=1..n+1)] :" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "plot(pcurve);" }} {PARA 13 "" 1 "" {GLPLOT2D 347 272 272 {PLOTDATA 2 "6%-%'CURVESG6$7aq7 $$\"\"!F)F(7$$\"39+++;vPp>!#?$\"3.+++=a4R7!#@7$$\"3N+++c%=J#yF-$\"33++ +GS.$))*F07$$\"31+++r,!*R$\"31+++6y0>LF-7$$\"31+++h4)H/$F9$\"3w+++ S]28yF-7$$\"3++++o+![l%F9$\"3&******RuRC^\"F97$$\"3d******\\X@HlF9$\"3 ++++)=.^e#F97$$\"3Q+++b@]6')F9$\"3u*****>PqA0%F97$$\"3++++jB$R3\"!#=$ \"3)*******\\W(*efF97$$\"3'******pA#Q98FS$\"3-+++b*H8M)F97$$\"37+++\"H %3X:FS$\"3*******\\?qD7\"FS7$$\"3'******R1D#oFS$\"31+++GvLb=FS7$$\"35+++iEUf@FS$\"3)******Rj_&*H#FS7$ $\"3\"******zQX6J#FS$\"3%******p/'p$z#FS7$$\"33+++Sj$HU#FS$\"3++++%3() [L$FS7$$\"35+++*>jr[#FS$\"39+++YN9>RFS7$$\"3'******HZpm\\#FS$\"3()**** *4aE9a%FS7$$\"3$******4h6\\W#FS$\"3I+++#[.d>&FS7$$\"3$******\\U(3EBFS$ \"3m*****>HG](eFS7$$\"3\"*******>PDN@FS$\"3]+++[wjrlFS7$$\"3'*******[6 Vo=FS$\"3S+++n<1xsFS7$$\"3++++dXpA:FS$\"3O+++]![A)zFS7$$\"3/+++gsC'4\" FS$\"3[*****zm&px')FS7$$\"3++++%))zZ)eF9$\"3I+++,!*f`$*FS7$$!3!)****** *****f8&!#D$\"2)*****4p(******!#<7$$!3O*****fr>Ln'F9$\"31+++'4$pg5Ffs7 $$!38+++`lT59FS$\"34+++\"Qfk6\"Ffs7$$!3*)*****>XI\\A#FS$\"35+++]#\\j; \"Ffs7$$!3@+++'*[O0JFS$\"30+++6(f%47Ffs7$$!3?+++kS3XSFS$\"3.+++\\)[\\C \"Ffs7$$!34+++J6SO]FS$\"3'******Hu\\?F\"Ffs7$$!3S+++W8nqgFS$\"3#****** *yL3!H\"Ffs7$$!3q*****4x+%QrFS$\"3/+++_.Z)H\"Ffs7$$!3l******>\"p$H#)FS $\"3++++q6u'H\"Ffs7$$!3*)*****H![wK$*FS$\"3#******\\)\\a%G\"Ffs7$$!3)* *****4-LP/\"Ffs$\"3!*******>mlh7Ffs7$$!3%******z\\^J:\"Ffs$\"3++++l<)z A\"Ffs7$$!3%******\\]j.E\"Ffs$\"3%******4FgN=\"Ffs7$$!33+++'[.UO\"Ffs$ \"34+++[ucG6Ffs7$$!3%******z/DNY\"Ffs$\"3++++zLJj5Ffs7$$!33+++-r@d:Ffs $\"3G+++(H-C))*FS7$$!3'******fe>Uk\"Ffs$\"3O+++&)e[kb&=Ffs$\"36+++!>(4HgFS7$$!31+++mft1>Ffs$\"3w******e(oc*[FS7$ $!3!*******f\\Ffs$\"3t*****4(zV9PFS7$$!3/+++bdSw>Ffs$\"33+++%Q%y'\\ #FS7$$!3-+++mN3%*>Ffs$\"33+++GXda7FS7$$!3)******f]*****>Ffs$\"33+++++I 5P!#C7$$!3!*******4O3%*>Ffs$!3!*******=rca7FS7$$!3++++VeSw>Ffs$!3!**** **>,xn\\#FS7$$!33+++\"4vr%>Ffs$!3!)******e1V9PFS7$$!3-+++Sht1>Ffs$!3?+ ++K:m&*[FS7$$!35+++*pkb&=Ffs$!3e*****>2!4HgFS7$$!3!******f`gUz\"Ffs$!3 ++++UX(R5(FS7$$!3$*******Q)QNs\"Ffs$!3_+++!G_.6)FS7$$!3(******p#*>Uk\" Ffs$!3y*****fBp\"R!*FS7$$!3)******>[qfN=\"Ffs7$$!3))*****f-_J:\"Ffs$!3/+ ++E7)zA\"Ffs7$$!3\"******4eLP/\"Ffs$!3)******46c;E\"Ffs7$$!3_*****frqF L*FS$!3#*******4Xa%G\"Ffs7$$!30+++C`PH#)FS$!33+++J2u'H\"Ffs7$$!3K+++Ys SQrFS$!35+++]*p%)H\"Ffs7$$!3]******o!y12'FS$!3++++:I3!H\"Ffs7$$!37+++# 33k.&FS$!35+++?%\\?F\"Ffs7$$!3#)*****z@\"4XSFS$!3&*******o&[\\C\"Ffs7$ $!3@+++HAP0JFS$!3)******RZf%47Ffs7$$!3()******Qz$\\A#FS$!3!******p0\\j ;\"Ffs7$$!31+++nTU59FS$!3++++N#fk6\"Ffs7$$!3o+++OoRtmF9$!33+++(*Hpg5Ff s7$$!3c***********\\H)Fd[l$!2-++?<(******Ffs7$$\"3#******>n,Z)eF9$!3a* ****\\'*)f`$*FS7$$\"3#*******G%Ri4\"FS$!3;+++>hpx')FS7$$\"3.+++[noA:FS $!3U+++!**[A)zFS7$$\"37+++#RB%o=FS$!3#******p>jqF(FS7$$\"3,+++YgCN@FS$ !3W+++m&R;d'FS7$$\"3\"******p')zgK#FS$!3m*****RpI](eFS7$$\"3')*******> /\\W#FS$!3)******HO1d>&FS7$$\"3#*******RAm'\\#FS$!3v*****R*)H9a%FS7$$ \"3-+++xh:([#FS$!3D+++et9>RFS7$$\"3'*******e&HHU#FS$!3!******HM\"*[L$F S7$$\"3#******p()Q6J#FS$!3)******pt+Pz#FS7$$\"3))******\\kTf@FS$!3-+++ Sxb*H#FS7$$\"35+++&*osv>FS$!37+++IIMb=FS7$$\"3#******>`>#o!p<&e#F97$$\"3!** ****4oqZl%F9$!3%*******ob^7:F97$$\"3/+++-l&H/$F9$!3\"******\\;^Q\"yF-7 $$\"3%)*****Rv!))RLF-7$$\"3.+++&fvH#yF-$!3I*****z<_5*)* F07$$\"3#******zh'Gp>F-$!3))*****fP!>Z7F07$$!3#*******[!pN&QFcs$!3))** ***R)HvO\")Fd[l-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$Q!6\"Fi [m-%%VIEWG6$%(DEFAULTGF^\\m" 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 0 "" }}}}{MARK "0 1 0" 16 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }{RTABLE_HANDLES 17558840 17561880 17562160 19229552 19229592 19229792 }{RTABLE M7R0 I5RTABLE_SAVE/17558840X*%)anythingG6"6"[gl!#%!!!"#"#"#?""%F& } {RTABLE M7R0 I5RTABLE_SAVE/17561880X*%)anythingG6"6"[gl!#%!!!"T"T$""!F($"+++++?!#5$"+++++SF+ $"+++++gF+$"+++++!)F+$"+++++5!"*$"+++++7F4$"+++++9F4$"+++++;F4$"+++++=F4$F*F4$" +++++AF4$"+++++CF4$"+++++EF4$"+++++GF4$"+++++IF4$"+++++KF4$"+++++MF4$"+++++OF4$ "+++++QF4$F-F4$"+++++UF4$"+++++WF4$"+++++YF4$"+++++[F4$"+++++]F4$"+++++_F4$"+++ ++aF4$"+++++cF4$"+++++eF4$F/F4$"+++++iF4$"+++++kF4$"+++++mF4$"+++++oF4$"+++++qF 4$"+++++sF4$"+++++uF4$"+++++wF4$"+++++yF4$F1F4$"+++++#)F4$"+++++%)F4$"+++++')F4 $"+++++))F4$"+++++!*F4$"+++++#*F4$"+++++%*F4$"+++++'*F4$"+++++)*F4$F3!")F& } {RTABLE M7R0 I5RTABLE_SAVE/17562160X*%)anythingG6"6"[gl!#%!!!"aq"aq=F&%&falseGE\[lT"""$"#?"" !""#$"+^9\cA!")""$$"+tQ"ya#F1""%$"+!)fGwGF1""&$"+;I>VKF1""'$"+/P/[OF1""($"+)H;v 3%F1"")$"+"RFSb%F1""*$"+$)4"R.&F1"#5$"+u)[`]&F1"#6$"+y#Hl$fF1"#7$"+hS5&G'F1"#8$ "+]ry+lF1"#9$"+FrKLlF1"#:$"+&>qnM'F1"#;$"+2pcOfF1"#<$"+*R$>S`F1"#=$"+h\;JYF1"#> $"+S"fh*QF1F,$"+!4K$3KF1"#@$"+q-">h#F1"#A$"+\P>A@F1"#B$"+'*ohMj?"F1"#E$"+ibML5F1"#F$"+Ka(>.*!"*"#G$"+!)yAf!)Ffp"#H$"+=bZStFfp"#I$"+, @"3#oFfp"#J$"+H$[0Y'Ffp"#K$"+gWYJiFfp"#L$"+'GmR6'Ffp"#M$"+y(4]4'Ffp"#N$"+KnimhF fp"#O$"+:?!\K'Ffp"#P$"+YUGplFfp"#Q$"+>:9-pFfp"#R$"+Kv\GtFfp"#S$"+N3!f&yFfp"#T$" +$["R%\)Ffp"#U$"+6i`c#*Ffp"#V$"+w0v:5F1"#W$"+[m^@6F1"#X$"+Z_,X7F1"#Y$"+\te)Q"F1 "#Z$"+y&R[b"F1"#[$"+n1hYF1"#]$"+(yj(=AF1"#^$"+>630DF1=F&F(E\[lTF *$F5F-F.$"+')3_SRFfpF2$"+3rlCRFfpF5$"+-xedRFfpF8$"+GwnYSFfpF;$"+"zf@?%FfpF>$"+8 W+QWFfpFA$"+j=&Hx%FfpFD$"+UncJ_FfpFG$"+)p5\%eFfpFJ$"+@%Q*\mFfpFM$"+N2'fo(FfpFP$ "+&Hjd)*)FfpFS$"+*pFf0"F1FV$"+2>6P7F1FY$"+!3:@V"F1Ffn$"+m<]C;F1Fin$"+t%RZz"F1F\ o$"+'oKe#>F1F,$"+X@73?F1Fao$"+F1Fjo$"+C&33">F1F ]p$"+`aPA=F1F`p$"+%G9Us"F1Fcp$"+v>a@;F1Fgp$"+Nn5=:F1Fjp$"+x'*\;9F1F]q$"+OTW=8F1 F`q$"+"GK]A"F1Fcq$"+Mn!p8"F1Ffq$"+C_Ra5F1Fiq$"+Db3w(*FfpF\r$"+;76l!*FfpF_r$"+L\ r4%)FfpFbr$"+je)z!yFfpFer$"+FlndsFfpFhr$"+F.RcnFfpF[s$"+P'=&FfpFgs$"+&Gmd!\FfpFjs$"+"FJ$F-FJ$"+:*o%>8FJ$"+o2I#Q"FJ$"+@E8X9FJ$"+uW'z]"FJ$"+Fjz q:FJ$"+!=GOj"FJ$"+L+Y'p"FJ$"+')=HfFJ$"+*H>1,#FJ$"+_6X t?FJ$"+0IGO@FJ$"+e[6*>#FJ$"+6n%>E#FJ$"+k&yZK#FJ$"+>vbMFJ$"+sPe=NFJ$ "+EcT"e$FJ$"+zuCWOFJ$"+K$zqq$FJ$F6FJ$"+QIuKQFJ$"+"*[d&*QFJ$"+WnSeRFJ$"+(fQ7-%FJ $"+]/2%3%FJ$"+.B!p9%FJ$"+cTt4UFJ$"+4gcsUFJ$"+jyRNVFJ$F8FJ$"+p:1hWFJ$"+AM*Q_%FJ$ "+v_s'e%FJ$"+Grb\YFJ$"+")*)Q7ZFJ$"+M3AvZFJ$"+(o_!Q[FJ$"+SX)3!\FJ$"+$R;P'\FJ$F:F J$"+*4!Q*3&FJ$"+`>@_^FJ$"+1Q/:_FJ$"+fc(yF&FJ$"+7vqS`FJ$"+l$RNS&FJ$"+=7PmaFJ$"+r I?HbFJ$"+C\.#f&FJ$F%>1fFJ$"+Vg-pfFJ$"+' *y&=.'FJ$"+\(*o%4'FJ$"+-;_dhFJ$"+bMN?iFJ$F*FJF& } {RTABLE M7R0 I5RTABLE_SAVE/19229792X*%)anythingG6"6"[gl!#%!!!"ew"ew=F&%&falseGE\[l`q"""$""!F ,""#$"+;vPp>!#7""$$"+c%=J#yF0""%$"+r,!*RFG"#9$"+iEUf@FG"#:$"+)QX6J#FG"#;$"+Sj$HU#FG"#<$"+*>jr[#FG"#=$"+t%pm\#FG "#>$"+6;"\W#FG"#?$"+Du3EBFG"#@$"+?PDN@FG"#A$"+\6Vo=FG"#B$"+dXpA:FG"#C$"+gsC'4"F G"#D$"+%))zZ)eF7"#E$!%O^F7"#F$!+;(>Ln'F7"#G$!+`lT59FG"#H$!+_/$\A#FG"#I$!+'*[O0J FG"#J$!+kS3XSFG"#K$!+J6SO]FG"#L$!+W8nqgFG"#M$!+r2SQrFG"#N$!+?"p$H#)FG"#O$!+.[wK $*FG"#P$!+@ItV5!"*"#Q$!+)\^J:"Ffr"#R$!+0NOg7Ffr"#S$!+'[.UO"Ffr"#T$!+[]_j9Ffr"#U $!+-r@d:Ffr"#V$!+'e>Uk"Ffr"#W$!+R&QNs"Ffr"#X$!+y-E%z"Ffr"#Y$!+#[kb&=Ffr"#Z$!+mf t1>Ffr"#[$!+g\Ffr"#\$!+bdSw>Ffr"#]$!+mN3%*>Ffr"#^$!+1&*****>Ffr"#_$!+5O3%*>F fr"#`$!+VeSw>Ffr"#a$!+"4vr%>Ffr"#b$!+Sht1>Ffr"#c$!+*pkb&=Ffr"#d$!+O0E%z"Ffr"#e$ !+R)QNs"Ffr"#f$!+F*>Uk"Ffr"#g$!+#[FG"#!*$"+K&>#oF0"$,"$!+\!pN&Q!#<=F&F(E \[l`qF*F+F-$"+=a4R7!#8F1$"+GS.$))*F\_lF4$"+6y0>LF0F8$"+S]28yF0F;$"+W(RC^"F7F>$" +)=.^e#F7FA$"+s.F_SF7FD$"+]W(*efF7FH$"+b*H8M)F7FK$"+0-dA6FGFN$"+UQ!GY"FGFQ$"+Gv Lb=FGFT$"+MEb*H#FGFW$"+Zgp$z#FGFZ$"+%3()[L$FGFgn$"+YN9>RFGFjn$"+TlUTXFGF]o$"+#[ .d>&FGF`o$"+#HG](eFGFco$"+[wjrlFGFfo$"+n<1xsFGFio$"+]![A)zFGF\p$"+ocpx')FGF_p$" +,!*f`$*FGFbp$"+"p(******FGFep$"+'4$pg5FfrFhp$"+"Qfk6"FfrF[q$"+]#\j;"FfrF^q$"+6 (f%47FfrFaq$"+\)[\C"FfrFdq$"+V(\?F"FfrFgq$"+zL3!H"FfrFjq$"+_.Z)H"FfrF]r$"+q6u'H "FfrF`r$"+&)\a%G"FfrFcr$"+?mlh7FfrFgr$"+l<)zA"FfrFjr$"+r-c$="FfrF]s$"+[ucG6FfrF `s$"+zLJj5FfrFcs$"+(H-C))*FGFfs$"+&)e (4HgFGFbt$"+f(oc*[FGFet$"+rzV9PFGFht$"+%Q%y'\#FGF[u$"+GXda7FGF^u$"&.r$F7Fau$!+> rca7FGFdu$!+7qx'\#FGFgu$!+f1V9PFGFju$!+K:m&*[FGF]v$!+s+4HgFGF`v$!+UX(R5(FGFcv$! +!G_.6)FGFfv$!+O#p"R!*FGFiv$!+]eR#))*FGF\w$!+dFJj5FfrF_w$!+^ocG6FfrFbw$!+-(fN=" FfrFew$!+E7)zA"FfrFhw$!+6hlh7FfrF[x$!+5Xa%G"FfrF^x$!+J2u'H"FfrFax$!+]*p%)H"FfrF dx$!+:I3!H"FfrFgx$!+?%\?F"FfrFjx$!+p&[\C"FfrF]y$!+u%f%47FfrF`y$!+d!\j;"FfrFcy$! +N#fk6"FfrFfy$!+(*Hpg5FfrFiy$!+sr******FGF\z$!+l*)f`$*FGF_z$!+>hpx')FGFbz$!+!** [A)zFGFez$!+(>jqF(FGFhz$!+m&R;d'FGF[[l$!+%pI](eFGF^[l$!+jjq&>&FGFa[l$!+%*)H9a%F GFd[l$!+et9>RFGFg[l$!+V8*[L$FGFj[l$!+P2q$z#FGF]\l$!+Sxb*H#FGF`\l$!+IIMb=FGFc\l$ !+?(4GY"FGFf\l$!+OkdA6FGFi\l$!+LbRT$)F7F\]l$!+HI/ffF7F_]l$!+n;M_SF7Fb]l$!+-p<&e #F7Fe]l$!+pb^7:F7Fh]l$!+l6&Q"yF0F[^l$!+#[[)>LF0F^^l$!+y@0"*)*F\_lFa^l$!+w.>Z7F\ _lFd^l$!+%)HvO")!#;F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F, F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F ,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F, F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F ,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F,F, F,F,F,F,F,F,F,F,F,F,F,F,F& }