{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 } {CSTYLE "" -1 256 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 266 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 271 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 273 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Tim es" 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 "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 "Normal" -1 256 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 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 65 "\324\340\351\353: ..\312 \356\354\357\322\345\365\355\356\353\\\312\361\345\355\356\364\356\355 \362\356\342\\\321\356\342\354\321\362\340\362\374\350\\\324\340\351 \353\373\312_\316\357\362\317\356\324\353\356\362\316\367\302\356\344 \373\\" }}{PARA 0 "" 0 "" {TEXT -1 44 "\321\362\316\357\362\324\353 \316\367\302\356\344 \302\345\360\3612.mws. \322\350\362\356\342 \312. \302. 14.02.17\343" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 27 "\307\340\344\340\367\350 \350\347 \305\364\350\354\356 \342 \300.\302., \362.2." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 269 48 "\320\345\370\345\355\350\345 \361\350\361\362\345\354 \344\350\364 \364\345\360\345\355\366\350\340\353\374\355\373\365 \363\360\340\342 \355\345\355\350\351 (\321\304\323)." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plo ts):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords ha s been redefined\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 5 "\304\340\355 \356:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 77 "\307\340\344\340\344\350 \354 \357\356\360\377\344\356\352 \361\350\361\362\345\354\373 \344 \350\364\364\345\360\345\355\366\350\340\353\374\355\373\365 \363\360 \340\342\355\345\355\350\351 (\321\304\323), \356\357\350\361\373\342 \340\376\371\345\351 \357\360\356\366\345\361\361" }}{PARA 0 "" 0 "" {TEXT -1 26 "\364\353\356\362\340\366\350\356\355\355\356\351 \356\367 \350\361\362\352\350 \342\356\344\373 " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "n:=3: " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "\350 \+ \354\340\362\360\350\366\363 \345\270 \352\356\375\364\364\350\366\350 \345\355\362\356\342" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "A:= matrix(n,n,[seq(seq(a[i,j],j=1..n),i=1..n)]):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 16 "\302\345\352\362\356\360 \360\345\370\345\355\350\351 \+ :" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "X:=vector(n,[Ca(t),Cb( t),Cc(t)]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"XGK%'vectorG6#7%-%# CaG6#%\"tG-%#CbGF+-%#CcGF+Q)pprint276\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 24 "\302\345\352\362\356\360 \355\340\367\340\353\374\355\373 \365 \363\361\353\356\342\350\351" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "X0:=vector(n,[100,0,0]): X0[1]:" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 82 "\317\356 \354\340\362\360\350\366\345 \300 \350 \+ \342\345\352\362\356\360\363 \360\345\370\345\355\350\351 \325 \347 \340\357\350\370\345\354 \361\340\354\363 \361\350\361\362\345\354\363 \344\350\364\364\345\360\345\355\366\350\340\353\374\355\373\365 \363 \360\340\342\355\345\355\350\351," }}{PARA 0 "" 0 "" {TEXT -1 34 "\350 \354\345\376\371\363\376 \342\345\352\362\356\360\355\356-\354\340\362 \360\350\367\355\363\376 \347\340\357\350\361\374:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "evalm(Diff(X,t)=A&*X); " }}{PARA 11 "" 1 " " {XPPMATH 20 "6#/K%'vectorG6#7%-%%DiffG6$-%#CaG6#%\"tGF.-F)6$-%#CbGF- F.-F)6$-%#CcGF-F.Q)pprint286\"KF%6#7%,(*&&%\"aG6$\"\"\"FAFAF+FAFA*&&F? 6$FA\"\"#FAF1FAFA*&&F?6$FA\"\"$FAF5FAFA,(*&&F?6$FEFAFAF+FAFA*&&F?6$FEF EFAF1FAFA*&&F?6$FEFIFAF5FAFA,(*&&F?6$FIFAFAF+FAFA*&&F?6$FIFEFAF1FAFA*& &F?6$FIFIFAF5FAFAQ)pprint29F8" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 256 49 "\317\345\360\345\351\344\345\354 \352 \367\350\361\353\345\355\355 \356\354\363 \360\345\370\345\355\350\376 \357\356\361\362\340\342\353 \345\355\355\356\351 \347\340\344\340\367\350" }{TEXT -1 1 "." }}} {EXCHG {PARA 0 "" 0 "" {TEXT 264 10 "\317\360\350\354\345\367\340\355 \350\345" }{TEXT -1 72 ". \304\353\377 \343\360\340\364\350\367\345 \361\352\350\365 \357\356\361\362\360\356\345\355\350\351 \357\360\356 \342\345\344\270\354 \367\350\361\353\345\355\355\356\345 \360\345\370 \345\355\350\345 \362\340\352 \352\340\352 \347\340\357\350\361\340 \355\356" }}{PARA 0 "" 0 "" {TEXT -1 58 "\355\350\346\345, \350\347 \354\345\355\350\342 \350\354\345\355\340 \357\345\360\345\354\345\355 \355\373\365 \350 \364\363\355\352\366\350\351 (\355\345\356\341\377 \347\340\362\345\353\374\355\356). " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 54 "\321\344\345\353\340\345\354 \357\356\344\361\362\340\355\356 \342\352\363 \367\350\361\353\345\355\355\373\365 \347\355\340\367\345 \355\350\351 \352\356\375\364\364\350\366\350\345\355\362\356\342 " } {TEXT 257 1 "k" }{TEXT -1 11 " \342 \354\340\362\360\350\366\363 " } {TEXT 258 1 "A" }{TEXT -1 16 ". \310\354\345\345\362 \347\355\340\367 \345\355\350\345" }}{PARA 0 "" 0 "" {TEXT -1 69 "\352\340\352 \347\340 \357\350\361\340\362\374 \352\356\375\364\364\350\366\350\345\355\362 \373: \342 \344\345\361\377\362\350\367\355\356\354 \342\350\344\345 \+ \350\353\350 - \355\345\362, \355\340\357\360\350\354\345\360, \362 \340\352" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "k1:=1/2*10^(-3) :k2:=1/10*10^(-6):k3:=1/10*10^(-3):k4:=1/10*10^(-5):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "k5:=1/10*10^(-5):k6:=1/2*10^(-6):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "A[1,1]:=-k1-k5:A[1,2]:=k2:A[ 1,3]:=k6:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "A[2,1]:=k1:A[2 ,2]:=-k2-k3:A[2,3]:=k4:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 " A[3,1]:=k5:A[3,2]:=k3:A[3,3]:=-k4-k6:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 20 "\320\340\361\357\345\367\340\362\340\345\354 \354\340\362 \360\350\366\363 " }{TEXT 265 1 "A" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "`evalm(A); eval(A[1,2])`:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 9 "row(A,1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vect orG6#7%#!$,&\"(+++\"#\"\"\"\")+++5#F+\"(+++#Q)pprint306\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 59 "\307\340\357\350\370\345\354 \361\350\361 \362\345\354\363 \344\350\364\364\345\360\345\355\366\350\340\353\374 \355\373\365 \363\360\340\342\355\345\355\350\351 \342 \367\350\361 \353\345\355\355\356\354 \342\350\344\345" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 22 "evalm(Diff(X,t)=A&*X);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/K%'vectorG6#7%-%%DiffG6$-%#CaG6#%\"tGF.-F)6$-%#CbGF-F. -F)6$-%#CcGF-F.Q)pprint316\"KF%6#7%,(*&#\"$,&\"(+++\"\"\"\"F+FA!\"\"*& #FA\")+++5FAF1FAFA*&#FA\"(+++#FAF5FAFA,(*&#FA\"%+?FAF+FAFA*&#\"%,5FEFA F1FAFB*&#FAF@FAF5FAFA,(*&FQFAF+FAFA*&#FA\"&++\"FAF1FAFA*&#\"\"$FHFAF5F AFBQ)pprint32F8" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 80 "\315\340\351 \344\345\354 \361\356\341\361\362\342\345\355\355\373\345 \367\350\361 \353\340 \354\340\362\360\350\366\373 \300 (\356\355\350 \346\345 \352 \356\360\355\350 \365\340\360\340\352\362\345\360\350\361\362\350\367 \345\361\352\356\343\356 \363\360\340\342\355\345\355\350\377)." }} {PARA 0 "" 0 "" {TEXT -1 25 "\304\356\357\363\361\362\350\354\356 eige nvalues(A)." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "ke:=eigenval s(A); `eigenvalues(A,'radical'); eigenvects(A)`:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#keG6%\"\"!,&#\"%8I\")+++5!\"\"*(\"\"$\"\"\"\"(+++&F+ \"'z26#F.\"\"#F.,&#F)F*F+*(F-F.F/F+F0F1F+" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "eigenvects(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6%7 %\"\"!\"\"\"<#K%'vectorG6#7%#\"%.5\"&?]\"F%#\"',5]\"%5vQ)pprint336\"7% ,&#\"%8I\")+++5!\"\"*(\"\"$F%\"(+++&F8\"'z26#F%\"\"#F%F%<#KF(6#7%,&#\" $6\"\"#bF%*&\"$l\"F8F " 0 "" {MPLTEXT 1 0 17 "C:=vector(n,[]): " }{TEXT -1 26 "- \342\345\352\362 \356\360 \355\340\367\340\353\374\355\373\365 \363\361\353\356\342\350 \351" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 297 "sys:=diff(Ca(t),t) =evalm(row(A,1)&*X), diff(Cb(t),t)=evalm(row(A,2)&*X), diff(Cc(t),t)=e valm(row(A,3)&*X): \+ fcns:=\{Ca(t),Cb(t),Cc(t)\}; \+ RA:=dsolve(\{sys,seq(subs(t=0,X[j])=C[j],j=1..n)\},fcn s):" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%fcnsG<%-%#CaG6#%\"tG-%#CbGF( -%#CcGF(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "`C[1]:=100:C[2] :=0: C[3]:=0`:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "plot([rhs (RA[1]),rhs(RA[2]),rhs(RA[3])],t=0..30000,Y=0..100);" }}{PARA 13 "" 1 "" {GLPLOT2D 494 179 179 {PLOTDATA 2 "6'-%'CURVESG6$7]o7$$\"\"!F)$\"3K _7+(>%y\")))!#K7$$\"3-++]7G$R<)!#;$\"3wWN&4flz)R!#<7$$\"3+++]ilyM;!#:$ \"3S)ep9)\\U$y(F37$$\"39++vV)z@X#F7$\"3u!H9j5P%R6F07$$\"3++++DJdpKF7$ \"33F?gzt%G[\"F07$$\"3H++](ofV!\\F7$\"3GHjAt\"H&>@F07$$\"3,+++]i9RlF7$ \"3Sx@s1g#Rp#F07$$\"3'4+voHl:'zF7$\"3iFk\\\\7,ZJF07$$\"3y++vVV)RQ*F7$ \"3e7jOR!=*fNF07$$\"31+D1R.k!3\"!#9$\"3$\\x29JRb$RF07$$\"30++vVA)GA\"F en$\"3eFfRfrcwUF07$$\"3'***\\iSS\"Ga\"Fen$\"3TDc!)e\")))H\\F07$$\"37++ ]Peui=Fen$\"3+[ee1.+YaF07$$\"3/+++]$)z%=#Fen$\"3%)pZnN)H)[eF07$$\"3S++ ]i3&o]#Fen$\"33vm.%4aK:'F07$$\"3>+]i:FPFGFen$\"3%>_$3%y*GujF07$$\"3'** **\\(oX*y9$Fen$\"3@0pn0f!p_'F07$$\"3%)***\\iNf]W$Fen$\"3`$3\"yW'[rh'F0 7$$\"3r***\\P9CAu$Fen$\"3y\\E\"4H!)fm'F07$$\"3!)***\\igsg*QFen$\"3q&fG `A.sn'F07$$\"3!****\\(o5#*\\SFen$\"3u1R$4%*>*zmF07$$\"3+++DJ&pP?%Fen$ \"35dpVk/)[n'F07$$\"35++v$*zhdVFen$\"3iO:#yLxFm'F07$$\"3m++v$fQen%Fen$ \"3fS*\\Fen$\"31&R&y!QT2b'F07$$\"3^***\\(=$f%Gc Fen$\"3g/&>6[6cO'F07$$\"3#3++]#y,\"G'Fen$\"32Om$>W)[EhF07$$\"3O++]7XF07$$\" 3++++!R5'f5Fiu$\"3qE\"zW4%o)G%F07$$\"3++]P/QBE6Fiu$\"33ooLtr8JSF07$$\" 3#******\\\"o?&=\"Fiu$\"3$=5=HRCS\"QF07$$\"3;+]Pa&4*\\7Fiu$\"37=%o(*fq xe$F07$$\"3;+]7j=_68Fiu$\"3W&*H%3<+RQ$F07$$\"3/++vVy!eP\"Fiu$\"3%GlafW %)H=$F07$$\"33+](=WU[V\"Fiu$\"3d/K.j7y3IF07$$\"3%****\\7B>&)\\\"Fiu$\" 3,NtIFy^JGF07$$\"3++]P>:mk:Fiu$\"3S&pc'R%o&eEF07$$\"3)***\\iv&QAi\"Fiu $\"3XBpC.d%o^#F07$$\"3))***\\PPBWo\"Fiu$\"3G%*>rUlasBF07$$\"3w*****\\N m'[F07$$\"37++D6W%)R>Fiu$\"3)*3A(y/:`'=F07$ $\"3w*****\\@80+#Fiu$\"3a_2(4H#yi+m\" F07$$\"3u**\\P4w)R7#Fiu$\"3kM(3y)o`s:F07$$\"39++]x%f\")=#Fiu$\"3WXCwGq l#[\"F07$$\"3m**\\P/-a[AFiu$\"3D!z!>pyD.9F07$$\"3/+](=Yb;J#Fiu$\"37:d$ yR'HD8F07$$\"3)*****\\i@OtBFiu$\"3Ie*e0m\\PD\"F07$$\"3r**\\PfL'zV#Fiu$ \"3G1%Hy3%\\$=\"F07$$\"3A+++!*>=+DFiu$\"3*o5J+UZ+7\"F07$$\"3#****\\i_4 Qc#Fiu$\"3/!=(y%=y\"f5F07$$\"3****\\P%>5pi#Fiu$\"30\\yDsUd-5F07$$\"37+ ++bJ*[o#Fiu$\"3Of*G[8Pm`*F37$$\"35++Dr\"[8v#Fiu$\"3R\"Gi)fyO5!*F37$$\" 34+++Ijy5GFiu$\"3-Zd\\sS$)o&)F37$$\"3-+]P/)fT(GFiu$\"3&4=\\\"fK_E\")F3 7$$\"3&***\\i0j\"[$HFiu$\"3@G363o-HxF37$$\"&++$F)$\"3[+M7#>,$GtF3-%'CO LOURG6&%$RGBG$\"#5!\"\"F(F(-F$6$7Z7$F($\"$+\"F)7$F5$\"3#)*=%Gu:i8#*F07 $F@$\"3Q5#p;]&4*[)F07$FE$\"3i-W!**ed:#yF07$FJ$\"3\\!epN;El?(F07$FT$\"3 e=!G(RvM\\iF07$Fin$\"3HQrGx0N>aF07$F^o$\"3dzZfW)Qph%F07$Fco$\"38^RMQZS LRF07$Fho$\"3S-BlbYeZLF07$F]p$\"3uf\")*)oG4\\GF07$Fbp$\"3g#QqB@unU#F07 $Fgp$\"3S+ZE>H9n?F07$Faq$\"3N*)o#zomb`\"F07$Fer$\"3D)*G1)zX!H6F07$F_s$ \"3pc,d\"*R:=#)F37$Fds$\"3O:%RC?x<*fF37$Fis$\"3n?b)3!yQLVF37$F^t$\"3i7 a\"p[s6E$F37$Fct$\"3#*=]gi`RsBF37$Fht$\"3U*ebX02$Gk9bF[el 7$Fbv$\"3O_-DDlWPTF[el7$Fgv$\"3eH`cugV`KF[el7$F\\w$\"3Mb;i!GUla#F[el7$ Faw$\"3ybx\"z6!*31#F[el7$Ffw$\"3#)\\ad@Xm&p\"F[el7$F[x$\"3w>Vs#pcMX\"F [el7$F`x$\"3'GHq+xF^E\"F[el7$Fex$\"3`g5]vuvF6F[el7$Fjx$\"3+M3WKGGV5F[e l7$F_y$\"3O1gX,aS(y*!#>7$Fdy$\"3!*y***y1vGL*F`gl7$Fiy$\"3\"Q()\\,Bt%Q! *F`gl7$F^z$\"3I/sS$ya_&))F`gl7$Fcz$\"3'Q()o1bF^t)F`gl7$Fhz$\"3'HFn4N$= !o)F`gl7$F][l$\"3YR=\"Gm&3g')F`gl7$Fb[l$\"3sB/@,e#fm)F`gl7$Fg[l$\"37c+ HgU:!p)F`gl7$F\\\\l$\"3ok6=$ojVs)F`gl7$Fa\\l$\"3!fwb8e\\ww)F`gl7$Ff\\l $\"3dij3>iS9))F`gl7$F[]l$\"3iVJI<[ul))F`gl7$F`]l$\"3j_si%fQf\"*)F`gl7$ Fe]l$\"3%>ZO2&4#p'*)F`gl7$Fj]l$\"3y4Mta/R;!*F`gl7$F_^l$\"3oOCc'=m/1*F` gl7$Fd^l$\"3dXV(>EM!4\"*F`gl7$Fi^l$\"3)yG&=UTa]\"*F`gl7$F^_l$\"3[f\\vX Wq#>*F`gl7$Fc_l$\"3pYnO'3/5B*F`gl7$Fh_l$\"3k1m\")e7%*p#*F`gl-F]`l6&F_` lF(F``lF(-F$6$7W7$F($\"3#4]+)yOr_N!#I7$F5$\"3'f&yfov#f.)F`gl7$F@$\"3oJ >zs=r0GF[el7$FE$\"3kAl)ooB8*eF[el7$FJ$\"3qqzuqHya**F[el7$FT$\"3q\"4c!4 UM2>F37$Fin$\"3\"RwnJjA3/$F37$Fco$\"3srAq]&\\f?'F37$F]p$\"3e*4_1PIl(** F37$Fgp$\"3&)z&e]<^fS\"F07$Faq$\"3W-/;@IX)z\"F07$Fer$\"3F-b6ko<3AF07$F _s$\"3,$fd+AVui#F07$Fds$\"3$4dO')z5_.$F07$Fis$\"3V%yuzxs,W$F07$F^t$\"3 ]6OKKys#y$F07$Fct$\"3g9P,]vD^TF07$Fht$\"3K5$z^%RX-XF07$F]u$\"3*)=$y9^< H#[F07$Fbu$\"3aeA10x#))4&F07$Fgu$\"3v\\Rkb]n3aF07$F]v$\"3#eBcd[phl&F07 $Fbv$\"3nG1Th$)[FfF07$Fgv$\"3uXiLY7W`hF07$F\\w$\"3Ql`UxRw'Q'F07$Faw$\" 3MEy(z#4\\&f'F07$Ffw$\"3]#fH)3*e+!oF07$F[x$\"3q^&R+<%owpF07$F`x$\"31i> *\\*3$e:(F07$Fex$\"3h$H)e&)R:ItF07$Fjx$\"3!zmG%o9ssuF07$F_y$\"3N]l)=0m wh(F07$Fdy$\"3[!*=glq8exF07$Fiy$\"31kpV%)y*p)yF07$F^z$\"3?WWI]d5/!)F07 $Fcz$\"3jArdC)\\f7)F07$Fhz$\"3(*z#yc_P&G#)F07$F][l$\"3u[!*oB.KJ$)F07$F b[l$\"3Ka+R'=(z=%)F07$Fg[l$\"3)[Exp\"Gl3&)F07$F\\\\l$\"3YGg7nx,)e)F07$ Fa\\l$\"3!)GHe_f$fm)F07$Ff\\l$\"3^b>ALfVP()F07$F[]l$\"3O\"H`t;Sw!))F07 $F`]l$\"3?mUPTm.r))F07$Fe]l$\"3k$3iU*[&=$*)F07$Fj]l$\"3g " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plo ts):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords ha s been redefined\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 40 "\307\340 \357\350\370\345\354 \321\304\323 \361 \357\345\360\345\354\345\355 \355\373\354\350 \352\356\375\364\364\350\366\350\345\355\362\340\354 \350" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "sys_ode := diff(y(x ),x)=-y(x)/x+x*z(x), diff(z(x),x)=-2*y(x)/x^3+z(x)/x;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(sys_odeG6$/-%%diffG6$-%\"yG6#%\"xGF-,&*&F*\"\" \"F-!\"\"F1*&F-F0-%\"zGF,F0F0/-F(6$F3F-,&*(\"\"#F0F*F0F-!\"$F1*&F3F0F- F1F0" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 29 "\310\371\345\354 \360\345 \370\345\355\350\345 \321\304\323 \342 \356\341\371\345\354 \342\350 \344\345" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "R:=dsolve([sys_ ode]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"RG<$/-%\"zG6#%\"xG,&*&%$ _C1G\"\"\")F*,&*$\"\"##F.F2F.F.!\"\"F.F.*&%$_C2GF.)F*,&F1F4F.F4F.F./-% \"yGF),**&F-F.)F*,&F1F.F.F.F.F.**F2F4F-F.F>F.F2F3F4*&)F*,&F.F.F1F4F.F6 F.F.**F2F4FBF.F6F.F2F3F." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 82 "\335 \362\356 \361\356\342\357\340\344\340\345\362 \361 \356\362\342\345 \362\356\354 \347\340\344\340\367\350 9.427. \320\340\361\357\345\367 \340\362\340\345\354 \357\360\340\342\373\345 \367\340\361\362\350 \+ \360\345\347\363\353\374\362\340\362\356\342 \360\345\370\345\355\350 \377" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "rhs(R[1]);rhs(R[2]) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&*&%$_C1G\"\"\")%\"xG,&*$\"\"##F &F+F&F&!\"\"F&F&*&%$_C2GF&)F(,&F*F-F&F-F&F&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,**&%$_C1G\"\"\")%\"xG,&*$\"\"##F&F+F&F&F&F&F&**F+!\"\" F%F&F'F&F+F,F.*&)F(,&F&F&F*F.F&%$_C2GF&F&**F+F.F0F&F2F&F+F,F&" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 85 "\317\360\356\342\345\344\345\354 \+ \357\360\356\342\345\360\352\363 \360\345\370\345\355\350\377, \357 \356\344\361\362\340\342\350\342 \345\343\356 \342 \350\361\365\356 \344\355\363\376 \321\304\323, \361\362\360\356\343\356 \361\356\356 \362\355\356\361\377 y(x) \350 z(x)" }}{PARA 0 "" 0 "" {TEXT -1 31 " \361 \355\356\354\345\360\356\354 \352\356\354\357\356\355\345\355\362 \373 \342 \347\340\357\350\361\350 R" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "Pr:=subs(y(x)=rhs(R[2]),z(x)=rhs(R[1]),sys_ode); simp lify(Pr);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#PrG/-%%diffG6$,&*&%$_C 1G\"\"\")%\"xG,&*$\"\"##F,F1F,F,!\"\"F,F,*&%$_C2GF,)F.,&F0F3F,F3F,F,F. ,&*(F1F,,**&F+F,)F.,&F0F,F,F,F,F,**F1F3F+F,FF(F2F( F-F.F(F(F*!\"$F/" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 16 "\317\360\356 \342\345\360\352\340 \355\340 \355\356\353\374" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 26 "simplify(lhs(Pr)-rhs(Pr));" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 55 "\307 \340\344\340\344\350\354 \355\340\367\340\353\374\355\373\345 \363\361 \353\356\342\350\377 (\354\356\346\355\356 \352\360\340\345\342\373 \345, \347\344\345\361\374 x <>0):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "ics := y(1)=1, z(1)=0; fcs:=\{y(x),z(x)\};" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$icsG6$/-%\"yG6#\"\"\"F*/-%\"zGF)\"\"!" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%$fcsG<$-%\"yG6#%\"xG-%\"zGF(" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "Rr:=dsolve([sys_ode, ics]); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#RrG<$/-%\"yG6#%\"xG,&*&,&*&\"\" #!\"\"F/#\"\"\"F/F0F1F2F2)F*,&*$F/F1F2F2F2F2F2*&)F*,&F2F2F5F0F2,&*&F/F 0F/F1F2F1F2F2F2/-%\"zGF),**(F-F2)F*,&F5F2F2F0F2F/F1F2*(F/F2F-F2F@F2F2* (F/F2F9F2)F*,&F5F0F2F0F2F2*(F9F2FDF2F/F1F0" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 24 "\317\356\361\362\360\356\350\354 \343\360\340\364\350\352 \350 \360\345\370\345\355\350\377" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "plot([rhs(Rr[1]),rhs(Rr[2])],x=.1..2,Y=0..20);" }} {PARA 13 "" 1 "" {GLPLOT2D 492 143 143 {PLOTDATA 2 "6&-%'CURVESG6$7U7$ $\"3/+++++++5!#=$\"3b_GF'z$>KJ!#<7$$\"3NLLekH227F*$\"3#4*oox8y'*GF-7$$ \"3lmm;Hf999F*$\"3OdW$z.PAr#F-7$$\"3-+vV3!>Vf\"F*$\"3I(o8I*f4!e#F-7$$ \"3QL$3x3#\\uNE#F-7$$ \"3#om;H@swe#F*$\"3P<%Hl45_5#F-7$$\"3UL$3-cmO*HF*$\"3w+++#*)y!y>F-7$$ \"3[m;/\"Hv+P$F*$\"31:gZ4L6z=F-7$$\"3y**\\PHZ#)fPF*$\"3w!p1%fnj!z\"F-7 $$\"3'omTg:/H;%F*$\"3Y;cD;QU5:F- 7$$\"3')****\\pF*$\"3hRq2^%G3 K\"F-7$$\"3OLL$esH\"[tF*$\"3e`lX\"G8IF\"F-7$$\"3FLLL.e'3r(F*$\"3HZBHd* pPB\"F-7$$\"3O**\\P%49G8)F*$\"3c4*G`[^#*=\"F-7$$\"3rKLLG)4j])F*$\"3-tI %y$Hi]6F-7$$\"3I**\\PWQ4;*)F*$\"3rg.v!>`)36F-7$$\"3=+]7L^I1$*F*$\"3-,! p&=p[p5F-7$$\"3SLL3x'\\Mr*F*$\"3^r,8[ImG5F-7$$\"3OL3_YNt35F-$\"3tiwsGB m7**F*7$$\"3gm;zz@1\\5F-$\"3AsjCR%)**3&*F*7$$\"3imTg&H_44\"F-$\"3c;L:) pF\")3*F*7$$\"34+DckxTF6F-$\"3sf4![X?&>()F*7$$\"3om;/q9!o;\"F-$\"3t)G= xc2#=$)F*7$$\"3/++]\"p)[27F-$\"3aPS`+#o#**yF*7$$\"3))**\\()*>$HZ7F-$\" 3#o>h(RnO%[(F*7$$\"35+vV3h!eG\"F-$\"39,85@&3v2(F*7$$\"30+]7F\"o&G8F-$ \"3$GUcnz.)=mF*7$$\"3;LL$Gq\"*pO\"F-$\"3)QDY#)o4)*>'F*7$$\"3(****\\72< !39F-$\"3Y7O'RhOZu&F*7$$\"3mmTg_@>X9F-$\"3KJTvud/D`F*7$$\"37++vNV$e[\" F-$\"3!Gy>`Fux&[F*7$$\"3GL3xia2C:F-$\"3()*3>6Z_'4WF*7$$\"37+v=f%[Sc\"F -$\"39LU?ks1KRF*7$$\"3ILLep$HJg\"F-$\"3%ehjsr3dX$F*7$$\"3()*\\P4YVSk\" F-$\"3O$fHGnMm%HF*7$$\"3wmmm$f[Mo\"F-$\"3]ZoadC'fW#F*7$$\"3wm;HLguBF*7$$\"3QL3x*y4Pw\"F-$\"3;WE%G\"))Q$R\"F*7$$\"3:++]JBV+ =F-$\"3AP/y1'3f'*)!#>7$$\"3PL$e%30_U=F-$\"3EWR4P^m]JFiy7$$\"3cmmmvY;!) =F-$!39u![0yuG;#Fiy7$$\"3/+vV47I?>F-$!3-S3C`oeZzFiy7$$\"3'**\\iN*pre>F -$!3sHS:`q>g8F*7$$\"\"#\"\"!$!3/,rnLd\"4)>F*-%'COLOURG6&%$RGBG$\"#5!\" \"$Fa[lFa[lF[\\l-F$6$7[o7$F($\"3$4LkgnRD$=!#:7$$\"3v;Hd?T)e-\"F*$\"3)* *Gm@a1Fs\"Fb\\l7$$\"3WLe9T#o<0\"F*$\"3ya$=\\MR>i\"Fb\\l7$$\"3+](=i6Fb\\l7$$\"3um\"Hd?T) e7F*$\"3o\"Rz4f1)\\5Fb\\l7$$\"39+](oW41J\"F*$\"3uaKJ@e\\@&*!#;7$$\"3EL 3-)oxBO\"F*$\"3=K4z^)f\"o')Fi^l7$F4$\"3?D5H`\\m=zFi^l7$$\"3L$3-)oCB/:F *$\"3)HEN-mgm\"oFi^l7$F9$\"3pA=zzuh=fFi^l7$$\"3q;H2[bS%o\"F*$\"3k\"H+Y zI#y^Fi^l7$F>$\"3EL->%Rc9c%Fi^l7$$\"3/+D1fb6x>F*$\"3!)[;u]_&R]$Fi^l7$F C$\"3g0?6Q7afFFi^l7$$\"3)omm;i0PQ#F*$\"3VR+\"=^_:&F*7$F]r$\"3PrIF!3hp$QF*7$Fbr$\"33Fs@00z%e#F *7$Fgr$\"3OsO%*\\2qZ:F*7$F\\s$\"3]kb$Q&Q%z)fFiy7$Fas$!3[$=EC;7Ts\"Fiy7 $Ffs$!3!4nt**[-$Q\"*Fiy7$F[t$!3s46#*)f0**f\"F*7$F`t$!3n1$H$3)Gx8#F*7$F et$!3U7-!HF%GlEF*7$Fjt$!3GcgB.n-gJF*7$F_u$!3Ozs'yYS7g$F*7$Fdu$!3oEBR`n 0$*RF*7$Fiu$!3U7;QY'4GR%F*7$F^v$!3X*\\Svd$=CZF*7$Fcv$!37*><)3FU_]F*7$F hv$!3n6y/!4$fH`F*7$F]w$!3Qjw!f;+Ih&F*7$Fbw$!3)[9\"\\*H:H'eF*7$Fgw$!3ye >,%>(f3hF*7$F\\x$!3uOnp()R'\\L'F*7$Fax$!3;I@JM,velF*7$Ffx$!3')H*pC)fui nF*7$F[y$!39S82E%f2'pF*7$F`y$!3)R[QW66v9(F*7$Fey$!3l1VHLOS6tF*7$F[z$!3 o\\%p%p&Q4\\(F*7$F`z$!3n^8fs)>Yk(F*7$Fez$!3fAkluQ*=!yF*7$Fjz$!3oS--B&e l%zF*7$F_[l$!3cW5Ee[7'4)F*-Fe[l6&Fg[lF[\\lFh[lF[\\l-%+AXESLABELSG6$Q\" x6\"Q\"YF\\jl-%%VIEWG6$;$\"\"\"Fj[lF_[l;F[\\l$\"#?Fa[l" 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 90 "---------------------------------- --------------------------------------------------------" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 267 8 "\271 9.427." }{TEXT -1 28 " \342 \342\345 \352\362\356\360\355\356-\354\340\362\360\350\367\355\356\351 \364\356 \360\354\345." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; w ith(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected na mes norm and trace have been redefined and unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "n:=2:A:=matrix(n,n,[-1/x,x,-2/x^3,1 /x]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AGK%'matrixG6#7$7$,$*&\" \"\"F,%\"xG!\"\"F.F-7$,$*&\"\"#F,F-!\"$F.F+Q)pprint366\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "X:=vector(n,[y(x),z(x)]); `map(D,X) ;`:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"XGK%'vectorG6#7$-%\"yG6#%\" xG-%\"zGF+Q)pprint376\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 " sys:=evalm(Diff(X,x)=A&*X); sys1:=evalm(A&*X); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$sysG/K%'vectorG6#7$-%%DiffG6$-%\"yG6#%\"xGF0-F+6$-% \"zGF/F0Q)pprint386\"KF'6#7$,&*&F0!\"\"F-\"\"\"F<*&F0F=F3F=F=,&*(\"\"# F=F0!\"$F-F=F<*&F0F%%sys1GK%'vectorG6#7$,&*&%\"xG!\"\"-%\"yG6#F+\"\"\"F,*&F+F0-%\" zGF/F0F0,&*(\"\"#F0F+!\"$F-F0F,*&F+F,F2F0F0Q)pprint406\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 87 "\317\360\350\354\345\367\340\355\350\345. \307\340\357\350\361\340\355\355\340\377 \342 \362\340\352\356\354 \+ \342\350\344\345 \321\304\323 sys \342 dsolve \355\345 \360\340\341 \356\362\340\345\362. \317\356\375\362\356\354\363 \347\340\357\350 \370\345\354 sys2." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 70 "sys2: =diff(y(x),x)=evalm(row(A,1)&*X),diff(z(x),x)=evalm(row(A,2)&*X):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "fcns:=\{y(x),z(x)\};" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%fcnsG<$-%\"yG6#%\"xG-%\"zGF(" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 30 "\307\340\357\350\370\345\354 \347 \340\344\340\355\350\345 \355\340 \360\345\370\345\355\350\345 \321 \304\323" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "R:=dsolve(\{sys 2\},fcns);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"RG<$/-%\"zG6#%\"xG,* *(\"\"#\"\"\"%$_C1GF.)F*,&*$F-#F.F-F.F.!\"\"F.F.*(F/F.F0F.F-F3F.*(F-F. %$_C2GF.)F*,&F.F4F2F4F.F.*(F7F.F8F.F-F3F4/-%\"yGF),&*&F/F.)F*,&F2F.F.F .F.F.*&F7F.)F*,&F2F4F.F.F.F." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 77 " \317\360\350\354\345\367\340\355\350\345. \310\355\362\345\343\360\350 \360\356\342\340\355\350\345 \321\304\323 \361 \357\345\360\345\354 \345\355\355\373\354\350 \352\356\375\364\364\350\366\350\345\355\362 \340\354\350 \356\341\373\367\355\356 \357\360\356\342\356\344\350\362 \361\377" }}{PARA 0 "" 0 "" {TEXT -1 19 "\354\345\362\356\344\356\354 \+ \350\361\352\353\376\367\345\355\350\377." }}}{EXCHG {PARA 0 "" 0 "" {TEXT 268 19 "\271 \361\362\360.101 \317\360\350\354\345\360 10" } {TEXT -1 28 " \342 \342\345\352\362\356\360\355\356-\354\340\362\360 \350\367\355\356\351 \364\356\360\354\345." }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 64 "\302 \361\353 \363\367\340\345 \361\350\361\362\345\354 \361 \357\356\361\362\356 \377\355\355\373\354\350 \352\356\375\364\364\350\366\350\345\355\362 \340\354\350 \350\354\345\345\354 \354\340\362\360\350\366\363 \367 \350\361\345\353" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "n:=3:A: =matrix(n,n,[4,1,0, 3,2,0, 2,3,4]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%\"AGK%'matrixG6#7%7%\"\"%\"\"\"\"\"!7%\"\"$\"\"#F,7%F/F.F*Q)pprint4 16\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 22 "\307\340\357\350\370\345 \354 \342\345\352\362\356\360 \360\345\370\345\355\350\351" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "X:=vector(n,[x(t),y(t),z(t)]); " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"XGK%'vectorG6#7%-%\"xG6#%\"tG-%\"y GF+-%\"zGF+Q)pprint426\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 37 "\317 \360\356\342\345\344\345\354 \360\345\370\345\355\350\345 \361 \357 \356\354\356\371\374\376 \330\317\316 Maple" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "`sys:=diff(x(t),t)=evalm(row(A,1)&*X),diff(y(t),t )=evalm(row(A,2)&*X),diff(z(t),t)=evalm(row(A,3)&*X);`: " }{TEXT -1 43 "- \347\340\357\350\361\374 \344\353\350\355\355\340\377, \340 \354 \356\346\355\356 \345\270 \361\344\345\353\340\362\374 \352\356\360 \356\367\345" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "sys2:=seq(d iff(X[k],t)=evalm(row(A,k)&*X),k=1..n);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%sys2G6%/-%%diffG6$-%\"xG6#%\"tGF-,&*&\"\"%\"\"\"F*F1F1-%\"yGF ,F1/-F(6$F2F-,&*&\"\"$F1F*F1F1*&\"\"#F1F2F1F1/-F(6$-%\"zGF,F-,(*&F;F1F *F1F1*&F9F1F2F1F1*&F0F1F?F1F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "fcns:=\{x(t),y(t),z(t)\};" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%% fcnsG<%-%\"xG6#%\"tG-%\"yGF(-%\"zGF(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 30 "\307\340\357\350\370\345\354 \347\340\344\340\355\350\345 \355 \340 \360\345\370\345\355\350\345 \321\304\323" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 23 "R:=dsolve(\{sys2\},fcns);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"RG<%/-%\"yG6#%\"tG,&*(\"\"$\"\"\"%$_C2GF.-%$expGF)F .!\"\"*&%$_C3GF.-F16#,$*&\"\"&F.F*F.F.F.F./-%\"xGF),&*&F/F.F0F.F.F3F./ -%\"zGF),(*&#\"\"(F-F.F>F.F.*(F9F.F4F.F5F.F.*&-F16#,$*&\"\"%F.F*F.F.F. %$_C1GF.F." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 41 "\317\360\356\342 \345\344\270\354 \357\360\356\342\345\360\352\363 \357\356\353\363\367 \345\355\355\356\343\356 \360\345\370\345\355\350\377 \321\304\323" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 23 "\315\350\346\345 \347\340\357\350 \361\374 \355\345 \357\360\356\365\356\344\350\362" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 78 "`Prov:=subs(seq(X[j]=rhs(R[j]),j=1..n),sys2 ); simplify(lhs(Prov)-rhs(Prov));`:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 75 "\317\360\350\361\342\356\345\355\350\345 \352\356\354\357\356 \355\345\355\362 \342\345\352\362\356\360\340 \360\345\370\345\355\350 \377 R \352\356\354\357\356\355\345\355\362\340\354 X \357\360\356\342 \345\344\270\354 \342 \340\353\343\356\360\350\362\354\345 " }}{PARA 0 "" 0 "" {TEXT -1 37 "\356\344\355\356\342\360\345\354\345\355\355 \356 \361 \342\373\341\356\360\352\356\351 \361\356\356\362\342\345 \362\361\362\342\350\377." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "if lhs( R[i])=x(t) then Pr:=subs(X[1]=rhs(R[i]),sys2):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "elif lhs(R[i])=y(t) then Pr:=subs(X[2]=rhs(R[i]),sys2 ):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "else Pr:=subs(X[3]=rhs(R[i]), sys2) fi:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "sys2:=Pr:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "simplify(Pr);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/,(*&#\"\"(\"\"$ \"\"\"*&%$_C2GF)-%$expG6#%\"tGF)F)F)*(\"#DF)%$_C3GF)-F-6#,$*&\"\"&F)F/ F)F)F)F)*(\"\"%F)-F-6#,$*&F9F)F/F)F)F)%$_C1GF)F)F$" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 16 "\317\360\356\342\345\360\352\340 \355\340 \355 \356\353\374" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "simplify(lh s(Pr)-rhs(Pr));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 39 "\317\360\356\342\345\360\352\340 \357\360 \356\370\353\340 \363\361\357\345\370\355\356 - \356\362\342\345\362 \+ \342\345\360\355\373\351." }}{PARA 0 "" 0 "" {TEXT 270 71 "\320\345 \370\350\354 \362\363 \346\345 \347\340\344\340\367\363, \350\361\357 \356\353\374\347\363\377 \362\345\356\360\350\376 \350 \354\345\362 \356\344\373 \360\345\370\345\355\350\377 \321\304\323 \361 \357\356 \361\362\356\377\355\355\373\354\350" }}{PARA 0 "" 0 "" {TEXT 272 15 " \352\356\375\364\364\350\366\350\345\355\362\340\354\350." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 79 "\315\340\365\356\344\350\354 \361\356\341 \361\362\342\345\355\355\373\345 \367\350\361\353\340 L (\352\356\360 \355\350) \350 \361\356\356\362\342\345\362\361\362\342\363\376\371 \350\345 \350\354 \361\356\341\361\362\342\345\355\355\373\345 \342 \345\352\362\356\360\373 VL" }}{PARA 0 "" 0 "" {TEXT -1 87 "(\367\340 \361\362\355\373\345 \360\345\370\345\355\350\377) VL*exp(L*t), \343 \344\345 VL - \342\345\352\362\356\360-\361\362\356\353\341\345\366 \+ \367\350\361\345\353 \362\340\352\350\365, \367\362\356 (A - L*E)*VL=0 ." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "V:=eigenvects(A);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"VG6%7%\"\"\"F'<#K%'vectorG6#7%F'! \"$#\"\"(\"\"$Q)pprint436\"7%\"\"%F'<#KF*6#7%\"\"!F9F'Q)pprint44F27%\" \"&F'<#KF*6#7%F'F'F " 0 "" {MPLTEXT 1 0 11 "V[1][3][1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vectorG6#7%\"\"\"!\"$#\"\"(\"\"$Q)pprint466\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 79 "\363\354\355\356\346\345\355\350\345 \352 \356\362\356\360\356\343\356 \355\340 \361\356\341\361\362\342\345\355 \355\363\376 \354\340\362\360\350\366\363 (A - L*E) \344\356\353\346 \355\356 \344\340\362\374 \355\356\353\374. \315\340\357\360\350\354 \345\360:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "M:=charmat(A,V [1][1]): A1:=eval(M); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A1GK%'mat rixG6#7%7%!\"$!\"\"\"\"!F)7%!\"#F*F*Q)pprint476\"" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 22 "evalm(A1&*V[1][3][1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vectorG6#7%\"\"!F'F'Q)pprint486\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 6 "\327.\362.\357." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 66 "\322\340\352\350\354 \356\341\360\340\347\356\354, \360 \345\370\345\355\350\345 \347\344\345\361\374 \341\363\344\345\362 \+ \357\360\345\344\361\362\340\342\353\345\355\356 \361\362\360\363\352 \362\363\360\355\356 \342 \342\350\344\345:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "RR:=sum(C[k]*V[k][3][1]*exp(V[k][1]*t), k=1..n); \+ " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#RRG,(*(&%\"CG6#\"\"\"F*K%'vecto rG6#7%F*!\"$#\"\"(\"\"$Q)pprint496\"F*-%$expG6#%\"tGF*F**(&F(6#\"\"#F* KF,6#7%\"\"!F@F*Q)pprint50F4F*-F66#,$*&\"\"%F*F8F*F*F*F**(&F(6#F2F*KF, 6#7%F*F*\"\"&Q)pprint51F4F*-F66#,$*&FMF*F8F*F*F*F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 78 "\317\360\350\354\345\367\340\355\350\345. \317 \356\361\353\345\344\356\342\340\362\345\353\374\355\356\361\362\374 \+ \361\356\341\361\362\342\345\355\355\373\365 \367\350\361\345\353 \350 \342\345\352\362\356\360\356\342 \356\357\360\345\344\345\353\377\345 \362\361\377 \342 V \350" }}{PARA 0 "" 0 "" {TEXT -1 21 "\354\356\346 \345\362 \341\373\362\374 \361\353\363\367\340\351\355\356\351." }} {PARA 0 "" 0 "" {TEXT -1 51 "\307\340\357\350\370\345\354 \360\345\370 \345\355\350\345 \342 \357\356\352\356\354\357\356\355\345\355\362\355 \356\351 \364\356\360\354\345 \342 \342\345\352\362\356\360 XX " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "XX:=vector(n,[]):" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 18 "\350\361\357\356\353\374\347\363 \377 \340\353\343\356\360\350\362\354" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "XX[i]:=sum(C[k]*V[k][3][1][i]*exp(V[k][1]*t), k=1..n)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 42 "\302\373\357\350\370\345\354 \357\356\353\363\367\345\355\355\373 \351 \362\340\352\350\354 \356\341\360\340\347\356\354 \342\345\352 \362\356\360 XX" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "eval(XX); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vectorG6#7%,&*&&%\"CG6#\"\"\"F, -%$expG6#%\"tGF,F,*&&F*6#\"\"$F,-F.6#,$*&\"\"&F,F0F,F,F,F,,&*(F4F,F)F, F-F,!\"\"F1F,,(*&#\"\"(F4F,F(F,F,*&&F*6#\"\"#F,-F.6#,$*&\"\"%F,F0F,F,F ,F,*(F9F,F2F,F5F,F,Q)pprint526\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 83 "\317\360\350\354\345\367\340\355\350\345. \307\344\345\361\374 \+ \357\356\347\350\366\350\377 N \342 V[k][3][N][i] \360\340\342\355\340 \345\344\350\355\350\366\345, \345\361\353\350 \355\345\362 \352\360 \340\362\355\373\365 \352\356\360\355\345\351." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 25 "\317\360\356\342\345\344\270\354 \357\360\356\342 \345\360\352\363 \360\345\370\345\355\350\377" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "Pro:=subs(X[i]=XX[i],sys2): sys2:=Pro:" }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "simplify(Pro);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/,(*&#\"\"(\"\"$ \"\"\"*&%$_C2GF)-%$expG6#%\"tGF)F)F)*(\"#DF)%$_C3GF)-F-6#,$*&\"\"&F)F/ F)F)F)F)*(\"\"%F)-F-6#,$*&F9F)F/F)F)F)%$_C1GF)F)F$" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 19 "\342 \362\356\354 \367\350\361\353\345 \355\340 \+ \355\356\353\374" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "simplif y(lhs(Pro)-rhs(Pro));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "\317\360\356\342\345\360\352\340 \+ \357\360\356\370\353\340 \363\361\357\345\370\355\356 - \355\356\353 \374." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 74 "\317\360\350\354\345\367 \340\355\350\345. \302\356\356\341\371\345 \357\356\360\377\344\356 \352 \361\353\345\344\356\342\340\355\350\377 \361\356\341\361\362\342 \345\355\355\373\365 \342\345\352\362\356\360\356\342 \342 V \354\356 \346\345\362 \341\373\362\374 " }}{PARA 0 "" 0 "" {TEXT -1 62 "\361 \353\363\367\340\351\355\373\354, \355\356 \355\340 \360\345\370\345 \355\350\345 \375\362\356 \355\345 \342\353\350\377\345\362. \302\340 \346\355\340 \357\360\350\342\377\347\352\340 \342 \357\340\360\345 " }}{PARA 0 "" 0 "" {TEXT -1 82 "\"\361\356\341\361\362\342\345\355\355 \356\345 \367\350\361\353\356 - \361\356\341\361\362\342\345\355\355 \373\351 \342\345\352\362\356\360\", \367\362\356 \344\356\361\362\350 \343\340\345\362\361\377 \342 V. \317\356\375\362\356\354\363 \342\361 \345\343\344\340 1-\356\345" }}{PARA 0 "" 0 "" {TEXT -1 67 "\360\345 \370\345\355\350\345 X[1] \361\356\341\350\360\340\345\362\361\377 \+ \350\347 \357\345\360\342\373\365 \352\356\354\357\356\355\345\355\362 \342 \361\356\341\361\362\342\345\355\355\373\365 \342\345\352\362 \356\360\340\365 " }}{PARA 0 "" 0 "" {TEXT -1 89 "\{V[k][3][1][1], k=1 ..n\}, \362\356 \345\361\362\374 X[i] = XX[i], \350 \362\340\352 \344 \340\353\345\345. \317\356\361\353\345\344\355\345\345 \357\356\347 \342\356\353\377\345\362 \361\356\356\362\355\345\361\362\350" }} {PARA 0 "" 0 "" {TEXT -1 35 "\361\362\360\356\343\356 X[i] \350 \356 \362\342\345\367\340\376\371\345\345 \345\354\363 XX[k]." }}}{EXCHG {PARA 0 "" 0 "" {TEXT 271 36 "\307\340\344\340\355\350\345 \344\353 \377 \361\340\354\356\361\362\356\377\362\345\353\374\355\356\343\356 \+ \360\345\370\345\355\350\377" }{TEXT -1 19 ". \305\364\350\354\356\342 \300.\302., \362.2." }}{PARA 0 "" 0 "" {TEXT -1 25 "\271\271 9.428., \+ 9.429., 9.430." }}{PARA 0 "" 0 "" {TEXT -1 37 "\315\340\351\362\350 \+ \360\345\370\345\355\350\377 \350 \357\356\361\362\360\356\350\362\374 \350\365 \343\360\340\364\350\352\350." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 65 "--------------------------------------------------------- --------" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 61 "\321 \352\360\340\362 \355\373\354\350 \352\356\360\355\377\354\350: \361\354. 2.4., \343 \353.3, \"\312\356\354\357\374\376\362\345\360\355\340\377 \354\340 \362\345\354\340\362\350\352\340\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 273 19 "\271 \361\362\360.102 \317\360\350\354\345\360 12" }{TEXT -1 28 " \342 \342\345\352\362\356\360\355\356-\354\340\362\360\350\367 \355\356\351 \364\356\360\354\345." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 64 "\302 \361\353 \363\367\340\345 \361\350\361\362\345\354 \361 \357\356\361\362\356 \377\355\355\373\354\350 \352\356\375\364\364\350\366\350\345\355\362 \340\354\350 \350\354\345\345\354 \354\340\362\360\350\366\363 \367 \350\361\345\353" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "n:=2:A: =matrix(n,n,[2,-1,4, 6]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AGK%' matrixG6#7$7$\"\"#!\"\"7$\"\"%\"\"'Q)pprint406\"" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 22 "\307\340\357\350\370\345\354 \342\345\352\362\356 \360 \360\345\370\345\355\350\351" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "X:=vector(n,[]); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"XG-%&arrayG6$;\"\"\"\"\"#7\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "V:=eigenvects(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"VG7%\"\"%\"\"#<#K%'vectorG6#7$\"\"\"!\"#Q)pprint416\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "\323\352\340\346\345\354 \361\362\360\363 \352\362\363\360\363 \342\345\352\362\356\360\340 V" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "V[1];V[2];V[3];V[3][1];V[3][1][2];" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<#K%'vectorG6#7$\"\"\"!\"# Q)pprint426\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vectorG6#7$\"\"\"! \"#Q)pprint436\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#!\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 51 "\307\340\357\350\370\345\354 \360\345\370 \345\355\350\345 \342 \357\356\352\356\354\357\356\355\345\355\362\355 \356\351 \364\356\360\354\345 \342 \342\345\352\362\356\360 XX " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "XX:=vector(n,[]):" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 18 "\350\361\357\356\353\374\347\363 \377 \340\353\343\356\360\350\362\354" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 78 "XX[i]:=sum((C[i]*V[3][1][1]+C[i+1]*t*V[3][1][1])*exp(V[1]*t), k= 1..n-(V[2]-1))" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 42 "\302\373\357\350\370\345\354 \357\356\353 \363\367\345\355\355\373\351 \362\340\352\350\354 \356\341\360\340\347 \356\354 \342\345\352\362\356\360 XX" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "eval(XX);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vector G6#7$*&,&&%\"CG6#\"\"\"F,*&&F*6#\"\"#F,%\"tGF,F,F,-%$expG6#,$*&\"\"%F, F1F,F,F,*&,&F.F,*&&F*6#\"\"$F,F1F,F,F,F2F,Q)pprint456\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "M:=charmat(A,V[1]): A1:=eval(M); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A1GK%'matrixG6#7$7$\"\"#\"\"\"7$! \"%!\"#Q)pprint466\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "eva lm(A1&*V[3][1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#K%'vectorG6#7$\"\" !F'Q)pprint476\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "0 3 0" 5 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 1 0 2 33 1 1 }