{VERSION 3 0 "IBM INTEL NT" "3.0" } {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 "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 }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 }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 "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 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 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with (linalg);" }} {PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }} {PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#7^r%.BlockDiagonalG%,GramSchmidtG%,Jor danBlockG%)LUdecompG%)QRdecompG%*WronskianG%'addcolG%'addrowG%$adjG%(a djointG%&angleG%(augmentG%(backsubG%%bandG%&basisG%'bezoutG%,blockmatr ixG%(charmatG%)charpolyG%)choleskyG%$colG%'coldimG%)colspaceG%(colspan G%*companionG%'concatG%%condG%)copyintoG%*crossprodG%%curlG%)definiteG %(delcolsG%(delrowsG%$detG%%diagG%(divergeG%(dotprodG%*eigenvalsG%,eig envaluesG%-eigenvectorsG%+eigenvectsG%,entermatrixG%&equalG%,exponenti alG%'extendG%,ffgausselimG%*fibonacciG%+forwardsubG%*frobeniusG%*gauss elimG%*gaussjordG%(geneqnsG%*genmatrixG%%gradG%)hadamardG%(hermiteG%(h essianG%(hilbertG%+htransposeG%)ihermiteG%*indexfuncG%*innerprodG%)int basisG%(inverseG%'ismithG%*issimilarG%'iszeroG%)jacobianG%'jordanG%'ke rnelG%*laplacianG%*leastsqrsG%)linsolveG%'mataddG%'matrixG%&minorG%(mi npolyG%'mulcolG%'mulrowG%)multiplyG%%normG%*normalizeG%*nullspaceG%'or thogG%*permanentG%&pivotG%*potentialG%+randmatrixG%+randvectorG%%rankG %(ratformG%$rowG%'rowdimG%)rowspaceG%(rowspanG%%rrefG%*scalarmulG%-sin gularvalsG%&smithG%,stackmatrixG%*submatrixG%*subvectorG%)sumbasisG%(s wapcolG%(swaprowG%*sylvesterG%)toeplitzG%&traceG%*transposeG%,vandermo ndeG%*vecpotentG%(vectdimG%'vectorG%*wronskianG" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 41 "A:=matrix(3,3,[2,-1,1/2,-1,3,1,1/2,1,5]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%\"\"#!\"\"#\"\" \"F*7%F+\"\"$F-7%F,F-\"\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "evalf(eigenvals(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6%,&$\"+WOK; a!\"*\"\"\"%\"IG$F'!#5,&$\"+x,$*Q6F&F'F($\"\"*F*,&$\"+zhuWMF&F'F($F&F* " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "Digits:=5;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'DigitsG\"\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "LUAlg:=proc(a,n)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "B:=evalm(a);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "for k to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 25 " u:=LUdecomp(B,L='l');" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 " B:=evalf(evalm(u&*l));" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 " print(B);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "end;" }} {PARA 7 "" 1 "" {TEXT -1 41 "Warning, `B` is implicitly declared local " }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `k` is implicitly declared l ocal" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `u` is implicitly declar ed local" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%&LUAlgGR6$%\"aG%\"nG6%% \"BG%\"kG%\"uG6\"F-C$>8$-%&evalmG6#9$?(8%\"\"\"F79%%%trueGC%>8&-%)LUde compG6$F0/%\"LG.%\"lG>F0-%&evalfG6#-F26#-%#&*G6$F " 0 "" {MPLTEXT 1 0 12 "LUAlg(A,20);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"++++DE!\"*$!+++++v!#5$\"+ ++++]F-7%$!++++v$*F-$\"++++DJF*$\"++++]7F*7%$\"+++]i5F*$\"++++D@F*$\"+ +++]UF*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+&4Q_4$! \"*$!+******\\K!#5$\"+++++]F-7%$!+.(oF-$\"+18R n\\F*$\"+++++:F*7%$\"+%*Qy>jF-$\"+Ibhr5F*$\"+=co0U$!\"*$\"*\\l/s*!#5$\"+++++] F-7%$\"+pOunDF-$\"+Ue<#H&F*$\"+t,*G]\"F*7%$\"+6BgUCF-$\"+BQOzFF-$\"+T) eeG\"F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+C$\\\\Y $!\"*$\"+%HO%G7!#5$\"+++++]F-7%$\"+egg6]F-$\"+PJ-g`F*$\"+@;Pl9F*7%$\"+ \"[etQ)!#6$\"+I-IDgF:$\"+Rv-v6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-% 'matrixG6#7%7%$\"+J-#[\\$!\"*$\"+hg/#G\"!#5$\"+++++]F-7%$\"+\"p&4k!)F- $\"+F;>d`F*$\"+;I0$R\"F*7%$\"+RR&)yF!#6$\"+Ov'3B\"F:$\"+U\"))z9\"F*" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+y$y$GN!\"*$\"+G6k #H\"!#5$\"+++++]F-7%$\"+3VZR7F*$\"+ylJI`F*$\"+14ox7F*7%$\"+Sg\"\\2*!#7 $\"+!yS%=CF:$\"+W]IT6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6 #7%7%$\"+-K2vN!\"*$\"+]Yh%H\"!#5$\"+++++]F-7%$\"+\"3b$f=F*$\"+BFM&G&F* $\"+&GP?5\"F*7%$\"+;m)4$H!#7$\"+-\"))z\\%!#8$\"+vSeR6F*" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+kX\"Gk$!\"*$\"+`R%\\H\"!#5$ \"+++++]F-7%$\"+n)=Xr#F*$\"+Nn1=_F*$\"+![H*>%)F-7%$\"+/j%*Q$*!#8$\"+u, I-v!#9$\"+,(=\"R6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7 %$\"+$=Q%RP!\"*$\"+&y%)\\H\"!#5$\"+++++]F-7%$\"+vac;QF*$\"+K^d@^F*$\"+ @d2%p%F-7%$\"+1s(*>H!#8$\"+[!o9I*!#:$\"+&o')*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+ionM))!\"*$\"+k=h@^F*$\"+@d2%p%!#57 %$!+HGy*y%F*$!+gQd))[F*$\"**y%y+%F/7%$\"+i%eQr)!#9$\"+z#pC\"*)F:$\"+Ow %*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+KK'z0'! \"*$\"+#G.;7&F*$\"+@d2%p%!#57%$\"+B/'\\9\"F*$!+&)4'=6#F*$\"+#H6d%HF/7% $\"+JFPB6!#9$!+DtA#3#F:$\"+\"R[*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#-%'matrixG6#7%7%$\"+k;&f-(!\"*$\"+.zg@^F*$\"+@d2%p%!#57%$!+\\s$4#eF/ $!+?l%)zIF*$\"+=o_e?F/7%$\"+)\\I?6#!#:$\"+$\\P77\"!#9$\"+%\\X*Q6F*" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+M7j,m!\"*$\"+$>1;7 &F*$\"+@d2%p%!#57%$\"+2*z+?#F/$!+So_bEF*$\"+NnUZCF/7%$\"+8,sBM!#;$!+S4 j[T!#:$\"+Wi%*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7% $\"+![:Bx'!\"*$\"+Epg@^F*$\"+@d2%p%!#57%$!+Fmr=%*!#6$!+/2@EGF*$\"+D0*4 H#F/7%$\"+'=qn!f!#<$\"+(R9*yf\"!#6$!+-aE%y#F*$\"+![M%HBF/7%$ \"+BYR)o\"!#=$\"+r)*)Q'H!#;$\"+Wf%*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+%eb#=n!\"*$\"+.ng@^F*$\"+@d2%p%!#57%$\"+5=$!+hPW77!#;$\"+pf%*Q6F*" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%$\"+!RaKs'!\"*$\"+Cng@ ^F*$\"+@d2%p%!#57%$!+E>c5F!#7$!+6(\\rx#F*$\"+&>cfL#F/7%$\"+)f$yV[!#?$ \"+icK\")\\!#<$\"+ff%*Q6F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matri xG6#7%7%$\"+_&*=@n!\"*$\"+;ng@^F*$\"+@d2%p%!#57%$\"+\"G4)=6!#7$!+x[3vF F*$\"+p'[yL#F/7%$\"+b:c0#)!#@$!+TM\"H/#!#<$\"+jf%*Q6F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "C:=hilbert(6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"CG-%'matrixG6#7(7(\"\"\"#F*\"\"##F*\"\"$#F*\"\"%#F* \"\"&#F*\"\"'7(F+F-F/F1F3#F*\"\"(7(F-F/F1F3F6#F*\"\")7(F/F1F3F6F9#F*\" \"*7(F1F3F6F9F<#F*\"#57(F3F6F9F " 0 " " {MPLTEXT 1 0 11 "LUAlg(C,4);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'m atrixG6#7(7($\"&9\\\"!\"%$\"&uL\"F*$\"&)[8F*$\"&I6\"F*$\"&n;'!\"&$\"&n m\"F37($\"&X6\"F3$\"&-I$F3$\"&T;%F3$\"&ot$F3$\"&[:#F3$\"&C&f!\"'7($\"& M\\(!\"($\"&hx#FC$\"&)4_FC$\"&Q\\&FC$\"&DV$FC$\"&1#**FG7($\"&\\(R!\")$ \"&:g\"FG$\"&<`$FG$\"&9a%FG$\"&W>$FG$FQFU7($\"&$)R\"!\"*$\"&L'eF\\o$\" &5S\"FU$\"&#G?FU$\"&Sk\"FU$\"&*ocF\\o7($\"&fQ#!#6$\"&D-\"!#5$\"&jb#F]p $\"&l(RF]p$\"&*yNF]p$\"&:V\"F]p" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%' matrixG6#7(7($\"&%)f\"!\"%$\"&pY\"F*$\"&-X\"F*$\"&l:\"F*$\"&CB'!\"&$\" &nm\"F37($\"&e)=!\"'$\"&ag#F3$\"&8U$F3$\"&Z-$F3$\"&Er\"F3$\"&pq%F97($ \"&2*))!\"*$\"&:j\"!\"($\"&X&=F9$\"&!*Q#F9$\"&Bf\"F9$\"&(yZFJ7($\"&!*y \"!#6$\"&(RO!#5$\"&z-'FG$\"&Z'p!\")$\"&'RsFhn$\"&vn#Fhn7($\"&bD\"!#9$ \"&Oo#!#8$\"&BB&!#7$\"&#H%*FV$\"&wN\"FG$\"%s#)FG7($\"&Xx\"!#=$\"&!)*Q! #<$\"&*)H)!#;$\"&F%=F`o$\"&JP%Fco$\".FV" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7(7($\"&eh\"!\"%$\"&hZ\"F*$\"&TX\"F*$\"&u: \"F*$\"&IB'!\"&$\"&nm\"F37($\"&x)G!\"($\"&IX#F3$\"&+E$F3$\"&1*GF3$\"&# R;F3$\"&.^%!\"'7($\"&F9*!#6$\"&x/\"!\")$\"&nk\"FD$\"&2?#FD$\"&X[\"FD$ \"&?[%F97($\"&)>p!#:$\"&]$))!#8$\"&E4#!#5$\"&w='FK$\"&Sr'FK$\"&#>DFK7( $\"&@))*!#?$\"&&G8!#<$\"&Ou$FW$\"&1y\"!#7$\"&$e7!\"*$\"& " 0 "" {MPLTEXT 1 0 20 "evalf(eigenvals(C));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+&[*z#3\"!#;" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 22 "hatvit:=proc(a,v0,i,n)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "v1:=evalm(v0);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 " for k to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 " v2:=evalm(a&*v 1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 " print(evalf(v2[i]/v1[i]) );" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 18 " v1:=evalm(v2);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "e nd;" }}{PARA 7 "" 1 "" {TEXT -1 42 "Warning, `v1` is implicitly declar ed local" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `k` is implicitly de clared local" }}{PARA 7 "" 1 "" {TEXT -1 42 "Warning, `v2` is implicit ly declared local" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%'hatvitGR6&%\"a G%#v0G%\"iG%\"nG6%%#v1G%\"kG%#v2G6\"F/C$>8$-%&evalmG6#9%?(8%\"\"\"F99' %%trueGC%>8&-F46#-%#&*G6$9$F2-%&printG6#-%&evalfG6#*&&F>6#9&\"\"\"&F2F M!\"\">F2-F46#F>F/F/F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "v 0:=vector(3,[1,1,1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#v0G-%'vect orG6#7%\"\"\"F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "hatvit (A,v0,1,30);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&+]\"!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&n;#!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&#pK!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&f?%!\"%" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&)4Z!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&%o\\!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&;7& !\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&;A&!\"%" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#$\"&')G&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&L L&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&FO&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&>Q&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\" &VR&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&BS&!\"%" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#$\"&uS&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$ \"&1T&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&FT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&ST&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# $\"&\\T&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&aT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&dT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&fT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&hT&!\"%" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&iT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&iT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&jT&! \"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&jT&!\"%" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#$\"&jT&!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&jT &!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&jT&!\"%" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%#\"?D\"GQ \+ " 0 "" {MPLTEXT 1 0 20 "evalf(eigenvals(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6%,&$\"&jT&!\"%\"\"\"%\"IG$\"\"#!\"&,&$\"&)Q6F&F'F($!\"\" F+,&$\"&[W$F&F'F(F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} }{MARK "13 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }