{VERSION 4 0 "IBM INTEL NT" "4.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 1 12 255 0 0 1 0 1 0 
0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 1 12 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 "Geneva" 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 1 }
{CSTYLE "" -1 257 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Normal
" -1 0 1 {CSTYLE "" -1 -1 "Geneva" 1 10 0 0 0 1 2 2 2 0 0 0 0 0 0 1 }
0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE 
"" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 0 1 }1 0 0 -1 -1 -1 
0 0 0 0 0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 
255 1 0 0 0 0 0 0 1 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "M
aple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 
1 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Title" 0 18 1 {CSTYLE "" 
-1 -1 "" 1 18 0 0 0 0 0 1 1 0 0 0 0 0 0 1 }3 0 0 -1 12 12 0 0 0 0 0 0 
19 0 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Monaco" 1 10 0 0 
255 1 2 2 2 0 0 0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R
3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Monaco" 1 9 255 0 0 1 2 2 2 0 0 
0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 258 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 
0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 18 "" 0 "" {TEXT 256 35 "Powers of Matrices and E
igenvectors" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg)
:" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }}
{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "A:= matrix(2,2,[[-1,3],[-2,4
]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'MATRIXG6#7$7$!\"\"\"
\"$7$!\"#\"\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "id := dia
g(1,1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#idG-%'MATRIXG6#7$7$\"\"
\"\"\"!7$F+F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "vects := '
vects';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&vectsGF$" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 27 " evalf(Eigenvals(A,vects));" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'VECTORG6#7$$\"+++++5!\"*$\"+++++?F)
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "print(vects);" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$$!+VH]?$)!#5$!+v7b0O!\"*7$$!
+i>+ZbF*$!+u7b0OF-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "n := 1
0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"#5" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 19 "v := vector([1,0]);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%\"vG-%'VECTORG6#7$\"\"\"\"\"!" }}}{EXCHG {PARA 258 "
" 0 "" {TEXT -1 58 "The power method for finding eigenvalues and eigen
vectors." }{TEXT 257 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 144 "for i \+
from 1 to 10 do\n       v := evalm(A &* v);\n\011\011v_len := evalf(sq
rt(dotprod(v,v))):\n\011\011v := map(evalf,evalm(1/v_len * v)):\n\011
\011print(v):\n\011\n  od:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'VECTO
RG6#7$$!+af8sW!#5$!+3>FW*)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'VEC
TORG6#7$$!+)*R%=S'!#5$!+(z7Ao(F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%
'VECTORG6#7$$!+&*4^/o!#5$!+=\\$zK(F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "
6#-%'VECTORG6#7$$!+q4A]p!#5$!+hP))*=(F)" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#-%'VECTORG6#7$$!+wwM8q!#5$!+)e?$GrF)" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#-%'VECTORG6#7$$!+H.%G/(!#5$!+bI=*4(F)" }}{PARA 11 "" 1 
"" {XPPMATH 20 "6#-%'VECTORG6#7$$!+1s5dq!#5$!+44+&3(F)" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#-%'VECTORG6#7$$!+'=DT1(!#5$!+GO+yqF)" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#-%'VECTORG6#7$$!+Tegnq!#5$!+'3GX2(F)" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#-%'VECTORG6#7$$!+c\"R$pq!#5$!+YgzsqF)" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "8 0 0" 19 }
{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }