{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 }} {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 301 "n:=5;\nA:=randmatrix(n,n);A0:=evalm(A):\nb:=r andvector(n); b0:=evalm(b):\n\n#G.A.:\n\nfor k from 1 to n-1 do\n f or i from k+1 to n do\n l:=A[i,k]/A[k,k];\n for j from k+1 t o n do\n A[i,j]:=A[i,j]-A[k,j]*l;\n od;\n b[i]:=b [i]-b[k]*l;\n A[i,k]:=l;\n od;\n od;\n\nprint(A,b);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"\"&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7'7'!#B!#P\"#J!#M!#U7'\"#))!#w!#l\"#D \"#G7'!#h!#g\"\"*\"#H!#m7'!#K\"#y\"#R\"#%*\"#o7'!#%\"bG-%'vectorG6#7'\"\"&!#))!#V!#t\"# D" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$-%'matrixG6#7'7'!#B!#P\"#J!#M!#U7 '#!#))\"#B#!%/]F0#\"%L7F0#!%#\"&\"pF\"%^77'#\"#KF0#!%*[\"\"%-D#!&Ua\"\"&&[N#\"(00h#\"&\"H @#\"(c#yg\"'bk57'#\"##\"%([)\"%(4(#!(siN\"\"(::$y#\"+_[4:J \")vv:R-%'vectorG6#7'\"\"&#!%%e\"F0#!%)\\*\"$R\"#!(4oM&FO#\"*khi/#\"(D 41#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "U:=LUdecomp(A0,L=L); print(L,evalm(L&*U),evalm(A0-L&*U));" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%\"UG-%'matrixG6#7'7'!#B!#P\"#J!#M!#U7'\"\"!#!%/]\"#B#\"%L7F3#!% " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "gausselim(concat(A0,b0));\nl insolve(A0,b0);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7&!#' )\"#B!#%)\"#y7&\"\"!#!%jQ\"#')#\"%')H\"#V#\"%s9F37&F-F-#\"'NA\")\"%jQ# \"'`4IF9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%#!(^<2\"\"'N A\")#!'[P:F)#\"'`4IF)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}}{MARK "3 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }