{VERSION 2 3 "SUN SPARC SOLARIS" "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 "" 0 1 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 }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 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "f:=1/x-y*exp(t)/t^2- t:g:=1/y-exp(-t)-exp(t):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "de1:=diff(x(t),t)=f:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "de 2:=diff(y(t),t)=g:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "ics:= x(1)=1,y(1)=exp(-1):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "IVP :=\{de1,de2,ics\}:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "FCNS: =\{x(t),y(t)\}:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "nSoln:=d solve(IVP,FCNS,type=numeric,method=classical[rk4],stepsize=0.001):" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "odeplot(Soln,[[t,x(t)],[t,y(t)]],1. .5);" }}{PARA 13 "" 1 "" {INLPLOT "6$-%'CURVESG6$7T7$$\"\"\"\"\"!F(7$$ \"+`Ej\"3\"!\"*$\"1g%[$>IGX#*!#;7$$\"+1`Ej6F.$\"1!*p\"*G7\\'f)F17$$\"+ fz*[C\"F.$\"1@;P')oyK!)F17$$\"+71`E8F.$\"1S')zR:YQvF17$$\"+lK;39F.$\"1 o&eoF\\95(F17$$\"+=fz*[\"F.$\"1Z$Q(o(GBr'F17$$\"+r&G9d\"F.$\"1u!R`OOOO 'F17$$\"+C71`;F.$\"1fb\"yr#Q\\gF17$$\"+xQpMF.$\"1J)>OR Y:0&F17$$\"+*[C71#F.$\"1^MQ]^[^[F17$$\"+Ur&G9#F.$\"1\"yI&ommmYF17$$\"+ &z*[CAF.$\"1O<*e%GT&\\%F17$$\"+[C71BF.$\"1&3$o(=$GOVF17$$\"+,^v(Q#F.$ \"1F6')*=M!)=%F17$$\"+axQpCF.$\"1!*)o'yne\\SF17$$\"+2/-^DF.$\"1[-z,++? RF17$$\"+gIlKEF.$\"1\\I<9'\\%)z$F17$$\"+8dG9FF.$\"1>Y1G0@%o$F17$$\"+m$ =fz#F.$\"1!f$\\PBkwNF17$$\"+>5bxGF.$\"1;(pmIx^Z$F17$$\"+sO=fHF.$\"1=1^ Y.JzLF17$$\"+Dj\"3/$F.$\"1?xo01f)G$F17$$\"+y*[C7$F.$\"1\"yB2Q9E?$F17$$ \"+J;3/KF.$\"1&>(*)4\">57$F17$$\"+%G9dG$F.$\"1$*[Yi#yM/$F17$$\"+PpMnLF .$\"105FrpppHF17$$\"+!fz*[MF.$\"1'zwbG3%**GF17$$\"+VAhINF.$\"1(>HP%*pB $GF17$$\"+'*[C7OF.$\"1fOV$eh$oFF17$$\"+\\v(Qp$F.$\"1K#=2FF17$$\"+- -^vPF.$\"1%eC,l['[EF17$$\"+bG9dQF.$\"19\\0%f#f#f#F17$$\"+3bxQRF.$\"1:% 3^5g)QDF17$$\"+h\"3/-%F.$\"1-N\"R\"za'[#F17$$\"+93/-TF.$\"1)*zff7=.CF1 7$$\"+nMn$=%F.$\"1\\xoegr0vF17$$\"+?hIlUF.$\"1)oI=QHnx#!#:7$$\"+t(QpM% F.$\"1@'R/hz;w'F_x7$$\"+E9dGWF.$\"1>!>A-szI\"!#97$$\"+zS?5XF.$\"1@%*[% **fP@#Fjx7$$\"+Kn$=f%F.$\"1kp-/#o+W$Fjx7$$\"+&QpMn%F.$\"1aQHB*[0/&Fjx7 $$\"+Q?5bZF.$\"1#>K5RWr2(Fjx7$$\"+\"pMn$[F.$\"181?*py8i*Fjx7$$\"+WtO= \\F.$\"1tk<0Hev7!#87$$\"+(*******\\F.$\"1Uz(*p:ed;Fiz-%'COLOURG6&%$RGB G$\"#5!\"\"F*F*-F$6$7T7$F($\"+7WzyO!#57$F,$\"1ePlj`T!R$F17$F3$\"12l%*> AkCJF17$F8$\"1wB69HqzGF17$F=$\"1Nx!4HkRl#F17$FB$\"1'H,u@@fW#F17$FG$\"1 !zita'=aAF17$FL$\"1+uLs=[x?F17$FQ$\"1Q*Qx**GY\">F17$FV$\"1[&Gr3UXw\"F1 7$Fen$\"1l.IA/AE;F17$Fjn$\"1[6!4tT()\\\"F17$F_o$\"1i\\4PgD\"Q\"F17$Fdo $\"1c:F%**zHF\"F17$Fio$\"107,i;>t6F17$F^p$\"1\\4[#pD73\"F17$Fcp$\"1s5n T*)ok**!#<7$Fhp$\"1T\"[OVhN=*F]_l7$F]q$\"1P4khhmj%)F]_l7$Fbq$\"1B))3eJ ?+yF]_l7$Fgq$\"1z&[Fs[()=(F]_l7$F\\r$\"19(ps#fADmF]_l7$Far$\"1EM!>Vxe5 'F]_l7$Ffr$\"1v\")4L/CFcF]_l7$F[s$\"1eF%\\dBh=&F]_l7$F`s$\"1'*=^qcezZF ]_l7$Fes$\"1'=P74;\\S%F]_l7$Fjs$\"1S\"[Np;'fSF]_l7$F_t$\"1R)yA;&QTPF]_ l7$Fdt$\"1dUyf'*4[MF]_l7$Fit$\"1rFUyY!y<$F]_l7$F^u$\"1!e.@+)pGHF]_l7$F cu$\"1!G^+p=\"*p#F]_l7$Fhu$\"1L*oH+Ov[#F]_l7$F]v$\"1C32*>RDH#F]_l7$Fbv $\"1@,wS\"GG6#F]_l7$Fgv$\"1[hAE[?Z>F]_l7$F\\w$\"1$>/xl$)y0#F]_l7$Faw$ \"1SQ_L\"\\-'=F]_l7$Ffw$!1)=J56'o=YF_x7$F[x$!1Q$[Zh\"=@5Fjx7$Fax$!1`J \\a]TF;Fjx7$Ffx$!1JX;*>M\\G#Fjx7$F\\y$!1n&G$on@)*HFjx7$Fay$!1?\">uTe?x $Fjx7$Ffy$!1(RM8=R;h%Fjx7$F[z$!1Y#3TtoD_&Fjx7$F`z$!1FGqE!G4^'Fjx7$Fez$ !1Lm-z " 0 "" {MPLTEXT 1 0 89 "rkf:=dsolve(IVP,FCNS,type=nu meric,method=rkf45,maxfun=-1):\nst:=time():\nrkf(5):\ntime()-st;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"%5e!\"$" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 77 "mg:=dsolve(IVP,FCNS,type=numeric,method=mgear):\nst :=time():\nmg(5):\ntime()-st;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$5 &!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "J:=matrix(2,2,[diff(f,x),dif f(f,y),diff(g,x),diff(g,y)]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"J G-%'MATRIXG6#7$7$,$*$%\"xG!\"#!\"\",$*&-%$expG6#%\"tG\"\"\"F4F-F.7$\" \"!,$*$%\"yGF-F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "J1:=map (simplify,subs(x=1/t,y=exp(-t),evalm(J)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#J1G-%'MATRIXG6#7$7$,$*$%\"tG\"\"#!\"\",$*&-%$expG6#F ,\"\"\"F,!\"#F.7$\"\"!,$-F26#,$F,F-F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "evs:=eigenvals(J1);" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%$evsG6$,$*$%\"tG\"\"#!\"\",$-%$expG6#,$F(F)F*" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 18 "sr:=evs[2]/evs[1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#srG*&-%$expG6#,$%\"tG\"\"#\"\"\"F*!\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "plot(sr,t=1..5);" }}{PARA 13 "" 1 " " {INLPLOT "6%-%'CURVESG6$7V7$$\"\"\"\"\"!$\"1]1$*)4c!*Q(!#:7$$\"1nmm; ')=(3\"F-$\"1zX!Re_BW(F-7$$\"1LL$e'40j6F-$\"1**)zSE.'ovF-7$$\"1nmm6hO[ 7F-$\"1Ah(f[:(>9F-$ \"1/dQ[(po[)F-7$$\"1nm;>K'*)\\\"F-$\"1aA&>`T2#*)F-7$$\"1++]Kd,\"e\"F-$ \"1j([+=#))[%*F-7$$\"1nm;fX(em\"F-$\"1\"=0(\\%*\\35!#97$$\"1++]U7Y]F-$\"1tG P5p)\\D\"FV7$$\"1+++I,Q+?F-$\"1c&['eFZl8FV7$$\"1+++]*3q3#F-$\"1_6!HTy< \\\"FV7$$\"1+++q=\\q@F-$\"1u94c\\&)H;FV7$$\"1nm;fBIYAF-$\"1&3^fdI3x\"F V7$$\"1LLLj$[kL#F-$\"1SU)=F--'>FV7$$\"1LLL`Q\"GT#F-$\"1x@j*Q!RT@FV7$$ \"1++]s]k,DF-$\"1u\"RK*QIzBFV7$$\"1LLL`dF!e#F-$\"1N[OC,Q/G Tf2Z%FV7$$\"1mm;f`@'3$F-$\"1)e\"G)e*oK]FV7$$\"1++]nZ)H;$F-$\"1T1:)f&\\ 'e&FV7$$\"1mmmJy*eC$F-$\"1%f'fVEahiFV7$$\"1+++S^bJLF-$\"1#***>XkOaqFV7 $$\"1+++0TN:MF-$\"14Lu=:?PzFV7$$\"1++]7RV'\\$F-$\"1X2MEpi1*)FV7$$\"1++ +:#fke$F-$\"1n$o<\"y]85!#87$$\"1LLL`4NnOF-$\"1vGi%Q.&R6Fdu7$$\"1+++],s `PF-$\"1I-pgYv#H\"Fdu7$$\"1mm;zM)>$QF-$\"1u7SI*y1X\"Fdu7$$\"1+++qfaH[3@Fdu7$$\"1LLL$)G[kTF-$\"15+%>cq$)Q#Fdu7$$\"1++]7yh]UF-$\"1n \\m'e^Ns#Fdu7$$\"1mmm')fdLVF-$\"1A*[z`nJ4$Fdu7$$\"1mmm,FT=WF-$\"1:!yyb `d_$Fdu7$$\"1LL$e#pa-XF-$\"1RhMK#4u,%Fdu7$$\"1+++Sv&)zXF-$\"1O'*[7u?KX Fdu7$$\"1LLLGUYoYF-$\"1flI%*)*[2_Fdu7$$\"1mmm1^rZZF-$\"1\\K\"3K'))**eF du7$$\"1LLe*3k**y%F-$\"1xB\\.#GtI'Fdu7$$\"1++]sI@K[F-$\"1bi\")=q&Ru'Fd u7$$\"1+++S2ls[F-$\"1sPqe)y6>(Fdu7$$\"1++]2%)38\\F-$\"1Z*[qI9\"pwFdu7$ $\"1++v.Uac\\F-$\"1zB**)eG%>#)Fdu7$$\"\"&F*$\"1(oAzJ'e5))Fdu-%'COLOURG 6&%$RGBG$\"#5!\"\"F*F*-%+AXESLABELSG6$%\"tG%!G-%%VIEWG6$;F(Fe[l%(DEFAU LTG" 2 368 368 368 2 0 1 0 2 9 0 4 2 1.000000 45.000000 45.000000 10030 10061 10056 10074 0 0 0 20030 0 12020 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "fsolve(sr=1 000,t=1..6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+y+**y]!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "8 1 0" 0 } {VIEWOPTS 1 1 0 1 1 1803 }