{VERSION 6 0 "IBM INTEL NT" "6.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 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 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 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plo
t" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 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 48 "Maple can produce companion matrix by a command:" }}{PARA 0 "" 
0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "p := (
3*x^6 - 5*x^3 + 3*x^2 - 8)/3;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"p
G,**$)%\"xG\"\"'\"\"\"F**&#\"\"&\"\"$F**$)F(F.F*F*!\"\"*$)F(\"\"#F*F*#
\"\")F.F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "C := LinearAlg
ebra:-CompanionMatrix(p,x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"CG-
%'RTABLEG6%\"*[1aZ\"-%'MATRIXG6#7(7(\"\"!F.F.F.F.#\"\")\"\"$7(\"\"\"F.
F.F.F.F.7(F.F3F.F.F.!\"\"7(F.F.F3F.F.#\"\"&F17(F.F.F.F3F.F.7(F.F.F.F.F
3F.%'MatrixG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "LinearAlgeb
ra:-Eigenvalues(1.0*C);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'RTABLEG6
%\"*![@w9-%'MATRIXG6#7(7#^$$\"3%oPRdP:#)G\"!#<$\"\"!F17#^$$\"3(4kjnaj7
#f!#=$\"3%fO:A)[qb')F67#^$F4$!3%fO:A)[qb')F67#^$$!3%)fx,u>w/yF6$\"3d5X
@-S`17F/7#^$F?$!3d5X@-S`17F/7#^$$!3u\"\\&)G!p::\"*F6F0&%'VectorG6#%'co
lumnG" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" 
{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 29 "read \"c:/304/RungeKutta.txt\":" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "f := Vector([1]):" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 21 "f[1] := (t,y) -> y^2:" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 12 "a, b := 0,1:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 18 "y0 := Vector([1]):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 8 "n := 30:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 
"t, z := RungeKutta(f,a,b,y0,n):" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%L
~End~of~Runge-Kutta~method~for~mxm~IVP-ODE~G" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 36 "pts := [seq([t[k],z[1][k]],k=1..n)];" }}{PARA 
12 "" 1 "" {XPPMATH 20 "6#>%$ptsG7@7$$\"\"!F($\"\"\"F(7$$\"+ieF[M!#6$
\"+$G9d.\"!\"*7$$\"+C<b'*oF.$\"+oS2u5F17$$\"+fF[M5!#5$\"+0YQ:6F17$$\"+
X.Jz8F:$\"+%)****f6F17$$\"+Jz8C<F:$\"+5LL37F17$$\"+<b'*o?F:$\"+K&p3E\"
F17$$\"+.Jz8CF:$\"+s<==8F17$$\"+!p?'eFF:$\"+=B&4Q\"F17$$\"+w#[M5$F:$\"
+9****\\9F17$$F-F:$\"+rcJE:F17$$\"+[M5$z$F:$\"+Z466;F17$$\"+M5$z8%F:$
\"+D@)eq\"F17$$\"+@'eF[%F:$\"+z'*\\7=F17$$\"+2ieF[F:$\"+wGLL>F17$$\"+$
z8C<&F:$\"+3zUr?F17$$\"+z8C<bF:$\"+[#o2B#F17$$\"+l*o?'eF:$\"+t^m;CF17$
$\"+_l*o?'F:$\"+%Hhjj#F17$$\"+QTs^lF:$\"+'='****GF17$$F4F:$\"+8d@AKF17
$$\"+5$z8C(F:$\"+=#))\\i$F17$$\"+'*o?'e(F:$\"+dT$G9%F17$$\"+$[M5$zF:$
\"+_QGL[F17$$\"+p?'eF)F:$\"+'>y)*z&F17$$\"+b'*o?')F:$\"+Klj\\sF17$$\"+
Ts^l*)F:$\"+gc?l'*F17$$\"+F[M5$*F:$\"+\\E,\\9!\")7$$\"+9C<b'*F:$\"+c*3
%zGF]t7$$\"+++++5F1$\"+tcuxB!\"(" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 10 "plot(pts);" }}{PARA 13 "" 1 "" {GLPLOT2D 407 407 407 
{PLOTDATA 2 "6%-%'CURVESG6$7@7$$\"\"!F)$\"\"\"F)7$$\"3!)*****>'eF[M!#>
$\"3!******HG9d.\"!#<7$$\"3g*****Rs^l*oF/$\"3/+++oS2u5F27$$\"3(*******
eF[M5!#=$\"3.+++0YQ:6F27$$\"3%******\\M5$z8F;$\"3%******R)****f6F27$$
\"3!******4$z8C<F;$\"3********4LL37F27$$\"3')*****p^l*o?F;$\"33+++K&p3
E\"F27$$\"36+++.Jz8CF;$\"3(******>x\"==8F27$$\"3))*******o?'eFF;$\"31+
++=B&4Q\"F27$$\"3%)*****fF[M5$F;$\"3%******R\"****\\9F27$$F.F;$\"35+++
rcJE:F27$$\"3w*****zW.Jz$F;$\"3,+++Z466;F27$$\"3t*****R.Jz8%F;$\"3.+++
D@)eq\"F27$$\"3y*****4ieF[%F;$\"34+++z'*\\7=F27$$\"3u*****p?'eF[F;$\"3
#******f(GLL>F27$$\"3E+++$z8C<&F;$\"3)******z!zUr?F27$$\"3n******y8C<b
F;$\"3&******zCo2B#F27$$\"3=+++l*o?'eF;$\"3;+++t^m;CF27$$\"3o*****>b'*
o?'F;$\"3')*****RHhjj#F27$$\"3?+++QTs^lF;$\"33+++'='****GF27$$F5F;$\"3
?+++8d@AKF27$$\"37+++5$z8C(F;$\"3$******z@))\\i$F27$$\"3k+++'*o?'e(F;$
\"3C+++dT$G9%F27$$\"39+++$[M5$zF;$\"35+++_QGL[F27$$\"3a******o?'eF)F;$
\"3))*****f>y)*z&F27$$\"31+++b'*o?')F;$\"3C+++Klj\\sF27$$\"3Y*****4C<b
'*)F;$\"3!)******fc?l'*F27$$\"3*******p#[M5$*F;$\"3)*******[E,\\9!#;7$
$\"3[*****RTs^l*F;$\"3)******f&*3%zGF^t7$F*$\"3++++tcuxB!#:-%'COLOURG6
&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$Q!6\"Fbu-%%VIEWG6$%(DEFAULTGFgu" 
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 25 "f[1] := (t,y) -> t^2-y^2:" }
}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "t, z := RungeKutta(f,a,b,y
0,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%\"zG6$-%
'RTABLEG6*\"*So'z:%)anythingG&%'VectorG6#%'columnG%,rectangularG%.Fort
ran_orderG7\"\"\"\";F4\"#I%\"xG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 36 "pts := [seq([t[k],z[1][k]],k=1..n)]:" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 10 "plot(pts);" }}{PARA 13 "" 1 "" {GLPLOT2D 407 
407 407 {PLOTDATA 2 "6%-%'CURVESG6$7@7$$\"\"!F)$\"\"\"F)7$$\"3!)*****>
'eF[M!#>$\"3u******)4,om*!#=7$$\"3g*****Rs^l*oF/$\"3_******>t*eN*F27$$
\"3(*******eF[M5F2$\"3S*****p[?g1*F27$$\"3%******\\M5$z8F2$\"3Y*****Rb
3hz)F27$$\"3!******4$z8C<F2$\"3Q+++'=w_a)F27$$\"3')*****p^l*o?F2$\"3\"
******R4!z7$)F27$$\"36+++.Jz8CF2$\"3w+++)>[!)4)F27$$\"3))*******o?'eFF
2$\"3r+++w=c+zF27$$\"3%)*****fF[M5$F2$\"31+++<2%*>xF27$$F.F2$\"3A+++(G
zeb(F27$$\"3w*****zW.Jz$F2$\"3;+++Zb93uF27$$\"3t*****R.Jz8%F2$\"37+++#
fqlF(F27$$\"3y*****4ieF[%F2$\"3=+++q%R5;(F27$$\"3u*****p?'eF[F2$\"3O**
***R'G[hqF27$$\"3E+++$z8C<&F2$\"3u+++![py(pF27$$\"3n******y8C<bF2$\"3[
+++!**)>5pF27$$\"3=+++l*o?'eF2$\"3s*****fU&\\eoF27$$\"3o*****>b'*o?'F2
$\"3r*****f$4!G#oF27$$\"3?+++QTs^lF2$\"3G+++N)pJ!oF27$$F5F2$\"3M+++%*G
m*z'F27$$\"37+++5$z8C(F2$\"3P+++U=M7oF27$$\"3k+++'*o?'e(F2$\"3$*******
\\SEToF27$$\"39+++$[M5$zF2$\"3g*****4Sxk)oF27$$\"3a******o?'eF)F2$\"3K
+++W`,[pF27$$\"31+++b'*o?')F2$\"3G+++[@*e-(F27$$\"3Y*****4C<b'*)F2$\"3
]+++#R)4?rF27$$\"3*******p#[M5$*F2$\"3X*****\\r'fIsF27$$\"3[*****RTs^l
*F2$\"3?+++wyJdtF27$F*$\"3l*****f<d,](F2-%'COLOURG6&%$RGBG$\"#5!\"\"F(
F(-%+AXESLABELSG6$Q!6\"F_u-%%VIEWG6$%(DEFAULTGFdu" 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 20 "Problem 9, page 357:
" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 15 "f := Vector(3):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 90 "f[1] := (t,u1,u2,u3) -> u2:\nf[2] := (t,u1,u2,u3) -> \+
u3:\nf[3] := (t,u1,u2,u3) -> exp(t)*u1:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 22 "y0 := Vector([1,0,1]):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 11 "a,b := 0,2:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 8 "n := 30:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "t, z := Ru
ngeKutta(f,a,b,y0,n):" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%L~End~of~Run
ge-Kutta~method~for~mxm~IVP-ODE~G" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 36 "pts := [seq([t[k],z[1][k]],k=1..n)]:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "plot(pts);" }}{PARA 13 "" 1 "" 
{GLPLOT2D 405 241 241 {PLOTDATA 2 "6%-%'CURVESG6$7@7$$\"\"!F)$\"\"\"F)
7$$\"3g*****Rs^l*o!#>$\"3,+++]PV-5!#<7$$\"3%******\\M5$z8!#=$\"3,+++\"
yl*45F27$$\"3')*****p^l*o?F6$\"31+++!QiH-\"F27$$\"3))*******o?'eFF6$\"
3)******p-?=/\"F27$$\"3!)*****>'eF[MF6$\"31+++0-(p1\"F27$$\"3t*****R.J
z8%F6$\"3#******\\5&)))4\"F27$$\"3u*****p?'eF[F6$\"3,+++sV3Q6F27$$\"3n
******y8C<bF6$\"3++++)3V^=\"F27$$\"3o*****>b'*o?'F6$\"33+++:7qS7F27$$F
.F6$\"3)*******eZZ08F27$$\"3k+++'*o?'e(F6$\"3#*******p*o-Q\"F27$$\"3a*
*****o?'eF)F6$\"3'*******eR*fY\"F27$$\"3Y*****4C<b'*)F6$\"3#******p<$o
j:F27$$\"3[*****RTs^l*F6$\"3++++,`^u;F27$$\"34+++fF[M5F2$\"3/+++r,%)*z
\"F27$$\"3&******fF[M5\"F2$\"3-+++v&47%>F27$$\"3/+++$z8C<\"F2$\"3()***
**R#QT+@F27$$\"3!*******4$z8C\"F2$\"3=+++\"GD&zAF27$$\"32+++G[M58F2$\"
35+++f-&4[#F27$$F5F2$\"39+++15\\2FF27$$\"3-+++ieF[9F2$\"3A+++3*>C'HF27
$$\"3*)******y8C<:F2$\"3!)*****z^o&\\KF27$$\"31+++(*o?'e\"F2$\"3!*****
*RiMMd$F27$$\"3#******RTs^l\"F2$\"3')*****RP6$RRF27$$\"3,+++Jz8C<F2$\"
3I+++'RVMN%F27$$\"35+++[M5$z\"F2$\"3x*****>8;K#[F27$$\"3'******\\'*o?'
=F2$\"3a+++9hQd`F27$$\"3\"******H[M5$>F2$\"3]+++=>OmfF27$$\"\"#F)$\"3%
*******H![Dm'F2-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%+AXESLABELSG6$Q!6\"Fb
u-%%VIEWG6$%(DEFAULTGFgu" 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 "26 1" 1 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 
0 1 2 33 1 1 }