{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 "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 Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple P lot" -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 "" {MPLTEXT 1 0 22 "restart: with(linalg ):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm a nd trace have been redefined and unprotected\n" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 39 "A:=matrix(3,3,[-5,2,3,4,-7,3,5,6,-11]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%!\"&\"\"#\"\"$7% \"\"%!\"(F,7%\"\"&\"\"'!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "B:=simplify(exponential(A,t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"BG-%'matrixG6#7%7%,(*&#\"\"$\"#q\"\"\"-%$expG6#,$*&\"#9F/%\"tGF/! \"\"F/F/*&#\"#A\"#XF/-F16#,$*&\"\"*F/F6F/F7F/F/#\"#f\"$E\"F/,(#\"#?\"# jF/*&#F:F;F/F " 0 "" {MPLTEXT 1 0 34 "s:=linsolve(transpose(A),[0,0,0]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%\"sG-%'vectorG6#7%,$*&#\"#f\"#F\"\"\"&%#_tG6#F.F.F.,$*&#\"#SF-F.F/F .F.F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "evalm(s/add(s[k],k =1..3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%#\"#f\"$E\"# \"#?\"#j#\"\"$\"#9" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "plot( \{B[1,2],B[2,2],B[3,2]\},t=0..0.7);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6'-%'CURVESG6$7Z7$$\"\"!F)F(7$$\"3Cmm;z>]9Q!#?$\" 327MOJ0x=A!#>7$$\"3[KLLeR+HwF-$\"3\\f?,8X*HI%F07$$\"3#****\\Pf]V9\"F0$ \"3Q,PfcaJgiF07$$\"3]mmm\"z+e_\"F0$\"3sXrbmv&z4)F07$$\"3#)**\\(oM'f*=# F0$\"3Rl%\\@*[0.6!#=7$$\"3ILL3->R`GF0$\"3esRwk!HFE7$$\"3]m;ay/Gl 6FE$\"35%>L/59[.$FE7$$\"3%***\\PurI88FE$\"3AWkL'QW!HJFE7$$\"3ILL$e#3dl 9FE$\"3%G6X+eu()>$FE7$$\"3imm\"Ht%o*f\"FE$\"3$4E?TxdCC$FE7$$\"3#)**** \\F_m]FE$\"3cMA*QSoyH$FE7$$\"3)** ***\\s2O[?FE$\"3\"RLSlE^(4LFE7$$\"3Om;aG\"H5=#FE$\"3l+PuUn)[J$FE7$$\"3 2LL$ej%yQBFE$\"3BpfMG8%eJ$FE7$$\"3ALLLVUUsCFE$\"3d#fvl[oMJ$FE7$$\"3V** \\(o()yyi#FE$\"3(4Bh]3l\"3LFE7$$\"3%GLL$oD[lFFE$\"3-%Q;w2I>I$FE7$$\"3r **\\(oibk\"HFE$\"3_Av\"3)41%H$FE7$$\"3w**\\i!o<-1$FE$\"3UK?=G)>gG$FE7$ $\"3.LL3-$=-@$FE$\"3ku'*p'>>uF$FE7$$\"3IL$3xplzM$FE$\"3s,)**REg&pKFE7$ $\"3[mm\"H([a'\\$FE$\"3))RBGe?IhKFE7$$\"3?m;ayo(3l$FE$\"3;Z6*z$RSjC$FE7$$\"3TmmT07KIRFE$\"3Sy%\\;Y'[RKFE 7$$\"3K*****\\\\@-3%FE$\"3uV*exq9HB$FE7$$\"3;***\\PopoA%FE$\"3?QF8q>*p A$FE7$$\"3m**\\(oMf(oVFE$\"3l775E7t@KFE7$$\"3\\***\\ii.j_%FE$\"3'zQ/7' GT;KFE7$$\"3=KLLoT'ym%FE$\"3X\"pvHw!37KFE7$$\"3v****\\i-,>[FE$\"32dm?, R*y?$FE7$$\"3!emT&)3rf&\\FE$\"3Y/Woi_Y/KFE7$$\"3m****\\Zaq0^FE$\"3*4og qp%3,KFE7$$\"3/K$3-\"QfY_FE$\"3/evhCsA)>$FE7$$\"3y**\\PWF'QR&FE$\"3;p( 44-Zb>$FE7$$\"3\"GLLe/Xy`&FE$\"3[VEfCI?$>$FE7$$\"3-**\\(=<\"e)o&FE$\"3 5wPy4X,\">$FE7$$\"3'pmmm(zvLeFE$\"3!)pn\"**pQ\"*=$FE7$$\"3#fm;zAAA)fFE $\"3!*>j&3KKu=$FE7$$\"3bK$3-7d%HhFE$\"3-e0]9-$f=$FE7$$\"3#)*****\\p]ZE 'FE$\"3#Gjo_$**p%=$FE7$$\"3RLLe*R7)>kFE$\"3G.>Z&3[M=$FE7$$\"33mmmO9]el FE$\"3#[,`>icC=$FE7$$\"3s**\\(o(GP1nFE$\"3S3Wbg$=:=$FE7$$\"3E**\\78Z!z %oFE$\"3)G#=mIHs!=$FE7$$\"3a**************pFE$\"3%*ew))pv'*zJFE-%'COLO URG6&%$RGBG$\"#5!\"\"F(F(-F$6$7W7$F($\"\"\"F)7$F2$\"3?^_J!G&=([*FE7$F< $\"33;)[f[nW,*FE7$FA$\"3l,6ftOAL')FE7$FG$\"3;YPO7-wM rS$FE7$Fcu$\"3Mz)*Hy+:wLFE7$Fhu$\"37>l0n8g_LFE7$F]v$\"3[P=e8[GILFE7$Fb v$\"3'f#*R!Qf>5LFE7$Fgv$\"3o95l/'zIH$FE7$F\\w$\"3!oYs29)fyKFE7$Faw$\"3 m_:;5v0KFE7$F^z$\"3R[>C!3+>?$FE7$Fcz$\"3'**3 ;(\\aX)>$FE7$Fhz$\"3o!>Glkqa>$FE7$F][l$\"356'>crfI>$FE7$Fb[l$\"3ivkH)Q Q1>$FE7$Fg[l$\"3')z9'>3W()=$FE7$F\\\\l$\"3[fdC)*4(p=$FE7$Fa\\l$\"3MqIk %=$[&=$FE7$Ff\\l$\"3_urB#z$3%=$FE-F[]l6&F]]lF(F^]lF(-F$6$7SF'7$F<$\"3y FKOy8tvHF07$FG$\"3(ptGLi#zEaF07$FQ$\"3)=!GSHW*H,)F07$Fen$\"3Akv(3mlD/ \"FE7$F_o$\"3')H=-(y\"\\j7FE7$Fdo$\"3OLnU2F[^9FE7$Fio$\"3k,MuB6kH;FE7$ F^p$\"3M(zZdasqz\"FE7$Fcp$\"3eNAW\\%\\z%>FE7$Fhp$\"3_yXi$R6w3#FE7$F]q$ \"3i%3[E%=a)>#FE7$Fbq$\"3$o&yPNh'4J#FE7$Fgq$\"3!*>l[*4. " 0 "" {MPLTEXT 1 0 38 "A:=m atrix(3,3,[-3,0,3,4,-6,2,1,3,-4]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"AG-%'matrixG6#7%7%!\"$\"\"!\"\"$7%\"\"%!\"'\"\"#7%\"\"\"F,!\"%" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "B:=simplify(exponential(A,t ));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"BG-%'matrixG6#7%7%,(#\"\"# \"\"&\"\"\"*&#\"\"$F-F.*&-%$expG6#,$*(\"#8F.F,!\"\"%\"tGF.F9F.-%$cosG6 #,$*(F,F9F:F.\"#6#F.F,F.F.F.F.*&#\"\"*\"#bF.*(F@FAF3F.-%$sinGF=F.F.F., (*&#F8FEF.FFF.F9*&#F.F-F.F2F.F9#F.F-F.,(*&#\"\"%FEF.FFF.F.*&#F,F-F.F2F .F9F+F.7%,(*&#F,F-F.F2F.F9*&#\"#9FEF.FFF.F.F+F.,(FNF.*&#FRF-F.F2F.F.*& #\"\")FEF.FFF.F9,(*&#F,F-F.F2F.F9*&#\"\"'FEF.FFF.F9F+F.7%,(*&#F,F-F.F2 F.F9*&#\"#;FEF.FFF.F9F+F.,(*&#F.F-F.F2F.F9*&#\"# " 0 "" {MPLTEXT 1 0 38 "plot(\{B[1,1],B[2,1],B[3,1]\},t=0..0.8);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6'-%'CURVESG6$7X7$$\"\"!F)F(7$$\"3% fmmmmh)=()!#?$\"3E+bc#pi6O$!#>7$$\"3>LLLLBxV5E$F0$\"3\"yorp#>#z8\"!#=7 $$\"3h+++v2<9TF0$\"3!pEN4@2iQ\"F@7$$\"3iLLLLAKn\\F0$\"3GoB4H(Glh\"F@7$ $\"3SLLLL*Gh#eF0$\"3M')*o0M?9$=F@7$$\"3?LLLLc$\\o'F0$\"31#p*RD&>/.#F@7 $$\"3\")*****\\dg'RvF0$\"3$fO\")[fQP@#F@7$$\"3Wmmm;bQ%R)F0$\"3%)HcP.bL $Q#F@7$$\"3KLLL$Qk#z**F0$\"3'Q%*3orqXm#F@7$$\"3%******\\YJ?;\"F@$\"3], V?_V,:HF@7$$\"3ELLL=\"\\=YRF@7$$\"3WmmmEn*Gn#F@$\"3as&G\\>z7+%F@7$$\"3 &ommmqFc#GF@$\"3Q&HCcWd)QSF@7$$\"3\"******\\9!H.IF@$\"3#3R5?^gO2%F@7$$ \"3=mmm1:bgJF@$\"3g$zt*3dl(4%F@7$$\"31+++X@4LLF@$\"3(=h,;jzx6%F@7$$\"3 /+++N;R(\\$F@$\"3u=,;gv&=8%F@7$$\"3)ommm\"4#)oOF@$\"3CP(zBBN@9%F@7$$\" 3jmmm6lCEQF@$\"31Bne)ol#[TF@7$$\"3OLLL$G^g*RF@$\"3!p(4W&HX>:%F@7$$\"37 LLL=2VsTF@$\"3UGlOO\"R_TF@7$$\"3 MLLLjcz\"\\%F@$\"3[Ntd!)f5]TF@7$$\"3g******z-6jYF@$\"3WA17x2YYTF@7$$\" 3%)******4#32$[F@$\"3R^&yd(3#>9%F@7$$\"3q*****\\#y'G*\\F@$\"3DVIXQ%Ro8 %F@7$$\"3u******H%=H<&F@$\"38Y!)f2qiITF@7$$\"3mlmm1>qM`F@$\"3uff/*4)pC TF@7$$\"3^++++.W2bF@$\"3\"[-f$f5:=TF@7$$\"3'HLL$ep'Rm&F@$\"3ahJ1Sq77TF @7$$\"34+++S>4NeF@$\"39^B)='Q`0TF@7$$\"3#emm;@2h*fF@$\"3y&[0:\\#R*4%F@ 7$$\"3S+++lXTkhF@$\"3!)\\(>#\\i4$4%F@7$$\"3bmmmmd'*GjF@$\"3;%4kx32r3%F @7$$\"3>*****\\iN7]'F@$\"396'[3^]53%F@7$$\"3)RLLL(>:nmF@$\"3Ge.gC,XvSF @7$$\"36LLL.a#o$oF@$\"3WVRWo\"y*pSF@7$$\"3?nmm^Q40qF@$\"3uA[Q5'=[1%F@7 $$\"3y******z]rfrF@$\"3K'o>w(yJgSF@7$$\"3/nmmc%GpL(F@$\"3#y_kN5Xa0%F@7 $$\"3!GLLL@Ia\\(F@$\"37kY5)=W80%F@7$$\"3k+++XhUkwF@$\"3[-!fw^Ns/%F@7$$ \"3&******\\\"o_N/%F@7$$\"3U+++++++!)F@$\"3hdsy(Rd)RSF@ -%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-F$6$7SF'7$F2$\"3g\"[Pr/?0\"=F07$F<$\" 3AY2(ojg'oMF07$FG$\"3eirp]G^(Q&F07$FQ$\"3v$)4w1&*yVtF07$Fen$\"330cuTRY (G*F07$Fjn$\"3d_7gxjq16F@7$F_o$\"3i#\\%Q9l:(G\"F@7$Fdo$\"3q1hOM4To9F@7 $Fio$\"3T#oMc;WEk\"F@7$F^p$\"3AU:p`ZP9=F@7$Fcp$\"3AT?K+C!*e>F@7$Fhp$\" 3z$R.b*\\u8@F@7$F]q$\"3OW9Q!Gq1E#F@7$Fbq$\"3ssn]oh5%R#F@7$Fgq$\"3N))>[ lCR3DF@7$F\\r$\"3_slW(*H!fj#F@7$Far$\"3!)oh^=['pt#F@7$Ffr$\"3U^HpzgvYG F@7$F[s$\"3_JgAGgAPHF@7$F`s$\"3ZxA;[kaHIF@7$Fes$\"3[Ir'G'G-6JF@7$Fjs$ \"3cP\"Qum!p*=$F@7$F_t$\"3aM_&fyKlD$F@7$Fdt$\"3z.1$p8pJK$F@7$Fit$\"350 BL'G\\nQ$F@7$F^u$\"3a&\\uX2:xV$F@7$Fcu$\"3u31P;nZ)[$F@7$Fhu$\"3\\s<#pa dl`$F@7$F]v$\"3-o13zwgzNF@7$Fbv$\"3S&*GN(yvxh$F@7$Fgv$\"3qvd\"=mQkl$F@ 7$F\\w$\"3QV!4\"z\\4)o$F@7$Faw$\"3J2]$egD*=PF@7$Ffw$\"3>&G^k2O)pPF@7$F`x$\"33THnX3g\"z$F@7$Fex$\"3vb(*>+CI7QF@7$Fjx$\"3 3c%G+4$oIQF@7$F_y$\"3S*zhsjB\"[QF@7$Fdy$\"3kzt0.qKjQF@7$Fiy$\"3![([#GS )RxQF@7$F^z$\"3WbQ8\"G3+*QF@7$Fcz$\"3XrKA\"3<0!RF@7$Fhz$\"3OG0\\.TS6RF @7$F][l$\"3NL^h5n=?RF@7$Fb[l$\"3k')=%*4&['GRF@7$Fg[l$\"3sf&*=:N&f$RF@7 $F\\\\l$\"3*\\#4W7'=I%RF@-Fa\\l6&Fc\\lF(Fd\\lF(-F$6$7S7$F($\"\"\"F)7$F 2$\"3u'QC(eu\"\\\\*F@7$F<$\"3[\"3FG(H<%3*F@7$FG$\"3'35%Q]V)Hl)F@7$FQ$ \"3KbSi'G6/D)F@7$Fen$\"3-Cw8aeezyF@7$Fjn$\"3iwU*RE35c(F@7$F_o$\"31cL^4 jLbsF@7$Fdo$\"3]*\\(*GvKO'pF@7$Fio$\"3QbB\"p!R2'p'F@7$F^p$\"3&p!Gwn\"o MW'F@7$Fcp$\"38B7Vo)H(QiF@7$Fhp$\"3/1K=8qzEgF@7$F]q$\"3dCR%=/#RKeF@7$F bq$\"3G,lN-:@hcF@7$Fgq$\"3mA)G*Ql^=bF@7$F\\r$\"3[T\"*30uXj`F@7$Far$\"3 ^$or'ekgV_F@7$Ffr$\"3=`=IkOT;^F@7$F[s$\"3^sqGP;A%F@7$Ffw$\"3gMr^.%>&*>%F@7$F[x$\"3_.8&[(oRxTF@7$F`x$\"3e;Vd3HqeTF@7$ Fex$\"3\"*p.>v%\\69%F@7$Fjx$\"37(\\*3mLwDTF@7$F_y$\"3'e*QJ26N6TF@7$Fdy $\"3yTCkMz%*)4%F@7$Fiy$\"3w`J&H^7w3%F@7$F^z$\"3wLPn8CexSF@7$Fcz$\"3m&) o'*zRKpSF@7$Fhz$\"3%)3ssWL(31%F@7$F][l$\"3gMDL&>TT0%F@7$Fb[l$\"397'G7t Lx/%F@7$Fg[l$\"3xwpE'QqA/%F@7$F\\\\l$\"3vYam)o_q.%F@-Fa\\l6&Fc\\lFd\\l Fd\\lF(-%+AXESLABELSG6$Q\"t6\"Q!Fh_m-%%VIEWG6$;F($\"\")Ff\\l%(DEFAULTG " 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "C urve 2" "Curve 3" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "B:=eva lf(exponential(A,17./60));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-% 'matrixG6#7%7%^$$\"+!GCyB&!#5$\"\"!F/^$$\"+Nh]a6F-$\"\"\"!#6^$$\"+'epw g$F-$!\"#F57%^$$\"+e`bSSF-$\"\"#F5^$$\"+H3m%y#F-F.^$$\"+9QyuJF-$!\"\"F 57%^$$\"+U!)*=u#F-F.^$$\"+^M;`CF-F3^$$\"+3&Q\\![F-F." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 20 "Constituent matrices" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "restart:with(linalg):" }} {PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and tra ce have been redefined and unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "A:=matrix(3,3,[-5,2,3,1,-3,2,2,0,-2]); n:=3:" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%!\"&\"\"#\"\"$7% \"\"\"!\"$F+7%F+\"\"!!\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "omega:=[eigenvalues(A)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&ome gaG7%\"\"!,&\"\"&!\"\"*$\"\"##\"\"\"F+F-,&F(F)F*F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "omega:=evalf(omega):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 74 "solve(\{seq(add(c[i]*omega[i]^k,i=1..n)=A^k,k= 0..n-1)\},\{seq(c[i],i=1..n)\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<% /&%\"cG6#\"\"$,&*&$\"+K-.7b!#6\"\"\")%\"AG\"\"#F.F.*&$\"+Dj\\w>!#5F.F0 F.F./&F&6#F1,&*&$\"+>j&)f)*F-F.F/F.!\"\"*&$\"+7CKCjF5F.F0F.F=/&F&6#F., ($F.\"\"!F.*&$\"+(3EyM%F-F.F/F.F.*&$FIF5F.F0F.F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "seq([omega[i],evalm(c[i])],i=1..n);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6%7$$\"\"!F%-%'matrixG6#7%7%$\"*l&p3E!\"*$\"+N/8RC'4W$F07%$!*odJG#F0$\"*!GWP8F0$\"))[rX*F07%$!+v'\\Q w$F0$\"+$47[?#F0$\"+#eP!f:F0" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "evalf(exponential(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'mat rixG6#7%7%$\"+V$>pf#!#5$\"+Lo#y\"=F*$\"+?QD&e&F*7%$\"+X(\\Ba#F*$\"+'o \"fk>F*$\"+q&eI\\&F*7%$\"+(*\\aMEF*$\"+FjVL;F*$\"+\"o=?t&F*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "evalm(add(c[i]*exp(omega[i]) ,i=1..n));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"*M>pf #!\"*$\"+No#y\"=!#5$\"+HQD&e&F-7%$\"+X(\\Ba#F-$\"+!o\"fk>F-$\"+r&eI\\& F-7%$\"+**\\aMEF-$\"+FjVL;F-$\"+u'=?t&F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "evalm(add(c[i]*log(1-omega[i]/10.),i=1..n));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+QW2OQ!#5$!+93E`9F* $!+DO\"GQ#F*7%$!+k1'eU)!#6$\"+9rucDF*$!+[5;9F;F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "restart:wi th(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected nam es norm and trace have been redefined and unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "A:=matrix(3,3,[-3,0,3,4,-6,2,1,3,-4 ]); n:=3:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%! \"$\"\"!\"\"$7%\"\"%!\"'\"\"#7%\"\"\"F,!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "omega:=eigenvalues(A); omega:=evalf(omega);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%&omegaG6%\"\"!,&#\"#8\"\"#!\"\"*&^## \"\"\"F*F/\"#6F.F/,&#F)F*F+*&^##F+F*F/F0F.F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&omegaG6%$\"\"!F'^$$!+++++l!\"*$\"+&R7$e;F+^$F)$!+&R7 $e;F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "solve(\{seq(add(c[ i]*omega[i]^k,i=1..n)=A^k,k=0..n-1)\},\{seq(c[i],i=1..n)\});" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#<%/&%\"cG6#\"\"#,**&$\"+66666!#6\"\"\")%\"AG F(F.!\"\"*&$\"+WWWW9!#5F.F0F.F1*&^#$\"+z&*fYEF5F.F0F.F.*&^#$\"+lQ;bVF- F.F/F.F./&F&6#\"\"$*(^$$!+/&RG<'!#I$!+![N&>C!#HF.F0F.,&*&$\"+++++=\"\" *F.F0F.F.^$$\"++++q6\"#5$!+7B'\\)HFNF.F./&F&6#F.,($F.\"\"!F.*&$\"+@AAA AF-F.F/F.F.*&$\"+))))))))GF5F.F0F.F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "seq([omega[i],evalm(c[i])],i=1..n);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6%7$$\"\"!F%-%'matrixG6#7%7%$\"+,+++S!#5$\"+********>F-$\"+++++SF- 7%$\"+********RF-$\"+)*******>F-F37%$\"+)*******RF-$\"+,+++?F-$\"+-+++ SF-7$^$$!+++++l!\"*$\"+&R7$e;FB-F'6#7%7%^$$\"+********HF-$!+*4-Or#F-^$ $!+**********!#6$\"+yuk>RF-^$$!+********>F-$!+z`/17F-7%^$FU$!*#)e6A%FB ^$F8$\"*w!47CFB^$FU$\"+m!o!4=F-7%^$$!+++++?F-$\"+6:=C[F-^$$!********** F-$!+jGpD^F-^$$\"+)*******HF-$\")M6:IFB7$^$F@$!+&R7$e;FB-F'6#7%7%^$FJ$ \"++@g8FF-^$$!+'*********FQ$!+yuk>RF-^$FU$\"+y`/17F-7%^$F_o$\"+A)e6A%F -^$F3$!+c247CF-^$FU$!+m!o!4=F-7%^$FU$!+6:=C[F-^$$!+)*********FQ$\"+cGp D^F-^$$\"+++++IF-$!+[M6:IFQ" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "evalf(exponential(A));" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "evalm (add(c[i]*exp(omega[i]),i=1..n));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#- %'matrixG6#7%7%^$$\"+,)Rt+%!#5$\"\"$!#6^$$\"+F+$\"\"!F3^$$\"+#=Q T+%F+$!\"\"F.7%^$$\"+l&pJ,%F+$!\"$F.^$$\"+OOs\"*>F+$!\"#F.^$$\"+*z1^*R F+FB7%^$$\"+;a2')RF+$\"\"\"F.^$$\"+lhh:?F+$\"\"#F.^$$\"+>%3$)*RF+F2" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%^$$\"+.)Rt+%!#5$\"\"! F-^$$\"+:?_))>F+$!\"\"!#8^$$\"+#=QT+%F+$\"\"\"F37%^$$\"+j&pJ,%F+F7^$$ \"+MOs\"*>F+F1^$$\"+*z1^*RF+F77%^$$\"+9a2')RF+F1^$$\"+nhh:?F+$\"#B!#9^ $$\"+?%3$)*RF+F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "restart :with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected \+ names norm and trace have been redefined and unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "A:=matrix(3,3,[-2,0,2,1,-2,1,1,1,-2 ]); n:=3:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%!\" #\"\"!\"\"#7%\"\"\"F*F.7%F.F.F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "omega:=eigenvalues(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&o megaG6%\"\"!!\"$F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "solv e(\{seq(c[1]*omega[1]^k+c[2]*omega[2]^k+c[3]*k*omega[3]^(k-1)=A^k,k=0. .n-1)\},\n\{seq(c[i],i=1..n)\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<% /&%\"cG6#\"\"\",(F(F(*&#\"\"#\"\"$F(%\"AGF(F(*&#F(\"\"*F(*$)F.F,F(F(F( /&F&6#F,,&*&#F,F-F(F.F(!\"\"*&#F(F1F(F2F(F:/&F&6#F-,&F.F:*&#F(F-F(F2F( F:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "evalf(exponential(A));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+ArClO!#5$\"+.0nz " 0 "" {MPLTEXT 1 0 71 "evalf(evalm(c[1]*exp(omega[ 1])+c[2]*exp(omega[2])+c[3]*exp(omega[3])));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+ArClO!#5$\"+.0nz " 0 "" {MPLTEXT 1 0 81 "evalf(evalm(c[1]*log(1-omega[1]/10)+c[2]*log (1-omega[2]/10)+c[3]/(omega[3]-10)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+(4&4\\WR)F-$\"+.&*)Rr\"F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "36 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }