{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 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" }} {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 19 "JAC:=proc(a,b,x0,n)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "d:=diag(a[1,1],a[2,2],a[3,3]);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "l:=matrix(3,3,[0,0,0,a[2,1],0,0,a[3,1],a[3,2],0]);" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "u:=matrix(3,3,[0,a[1,2],a[1,3],0,0 ,a[2,3],0,0,0]);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "di:=inverse(d); " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "B:=evalm((-1)*di&*(l+u));" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "r:=evalm(di&*b);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "y:=evalm(x0);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "for i to n do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "y:=evalf(evalm (B&*y+r));\nprint(y);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "end;" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `d` is implicitly declared local" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `l` is implicitly declared local" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `u` is implicitly declared local" }}{PARA 7 "" 1 "" {TEXT -1 42 "Warning, `di` is implicitly declared local" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `B` is implicitly declared local" }} {PARA 7 "" 1 "" {TEXT -1 41 "Warning, `r` is implicitly declared local " }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `y` is implicitly declared l ocal" }}{PARA 7 "" 1 "" {TEXT -1 41 "Warning, `i` is implicitly declar ed local" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$JACGR6&%\"aG%\"bG%#x0G% \"nG6*%\"dG%\"lG%\"uG%#diG%\"BG%\"rG%\"yG%\"iG6\"F4C*>8$-%%diagG6%&9$6 $\"\"\"F>&F<6$\"\"#FA&F<6$\"\"$FD>8%-%'matrixG6%FDFD7+\"\"!FKFK&F<6$FA F>FKFK&F<6$FDF>&F<6$FDFAFK>8&-FH6%FDFD7+FK&F<6$F>FA&F<6$F>FDFKFK&F<6$F AFDFKFKFK>8'-%(inverseG6#F7>8(-%&evalmG6#-%#&*G6$,$Fhn!\"\",&FFF>FSF>> 8)-F_o6#-Fbo6$Fhn9%>8*-F_o6#9&?(8+F>F>9'%%trueGC$>F_p-%&evalfG6#-F_o6# ,&-Fbo6$F]oF_pF>FhoF>-%&printG6#F_pF4F4F4" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 38 "A:=matrix(3,3,[1,0,0.5,0,-1,0,0,1,1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%\"\"\"\"\"!$\"\"&!\"\"7%F +F.F+7%F+F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "B:=vector( 3,[5,-1,0]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-%'vectorG6#7%\" \"&!\"\"\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "x0:=vector (3,[3,4,6]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x0G-%'vectorG6#7%\" \"$\"\"%\"\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "JAC(A2,B,x 0,40);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+++++g!#5$ \"+++++G!\"*$!+9dG9!\"*$!+q&G9<)!#5$!+++++7F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+H9dw5!\"*$!+Vr&G!p!#5$\"+;3/-NF," }}{PARA 11 " " 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+3/->z!#5$\"+z*[Wn\"F)$\"+ZtOeH F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+Ej@V(*!#5$\"+L ;G(o\"F)$!+5*4i<(!#6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7% $\"+_!)4[5!\"*$!*G@8L$!#5$!+**p?Js!#6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+r(*z25!\"*$!*-%e'4%!#5$\"+?\"4xU\"!#6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+PT^*))*!#5$\"*]o&ej!#6$\"+ m)ycv\"F," }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+gNgx**! #5$\"*'HD$))*!#6$!+'y+^s#!#7" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vec torG6#7%$\"+r;_-5!\"*$!)geQ6!#5$!+S)zcB%!#7" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+>%>1+\"!\"*$!)q+qB!#5$\"+'Ga'z[!#8" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+bSG%***!#5$\")\"3/# =!#6$\"+'G" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+**********!#5$!\"%F)$\"+9d G9Z!#>" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+)********* !#5$\"#D!#6$\"+9dG9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6 #7%$\"+++++5!\"*$\"#6!#6$!+r&G92\"!#>" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"+++++5!\"*$!\"\"!#5$!+9dG9Z!#?" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%'vectorG6#7%$\"+++++5!\"*\"\"!$\"+'G9dG%!#?" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"\"\"\"\"!$\"\"#!#6F) " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"\"\"\"\"!F)$!+r&G 9d)!#@" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"\"\"\"\"!F) F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%$\"\"\"\"\"!F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "A2:=matrix(3,3,[5,-1,1,-1 ,-10,5,0,3,7]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A2G-%'matrixG6#7 %7%\"\"&!\"\"\"\"\"7%F+!#5F*7%\"\"!\"\"$\"\"(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "5 0 0" 13 }{VIEWOPTS 1 1 0 1 1 1803 }