{VERSION 2 3 "IBM INTEL NT" "2.3" } {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 } {CSTYLE "" -1 256 "" 1 14 0 0 10 0 0 1 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 0 }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 "" {TEXT 256 30 "Matrices and matrix opera tions" }}}{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 23 "A:=array(1..2,1..3,[]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%&arrayG6%;\"\"\"\"\"#;F)\"\"$7 \"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "x:=vector(3,[]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%&arrayG6$;\"\"\"\"\"$7\"" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalm(A);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%'MATRIXG6#7$7%&%\"AG6$\"\"\"F+&F)6$F+\"\"#&F)6$F+\" \"$7%&F)6$F.F+&F)6$F.F.&F)6$F.F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalm(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'VECTOR G6#7%&%\"xG6#\"\"\"&F(6#\"\"#&F(6#\"\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "evalm(A &* x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%' VECTORG6#7$,(*&&%\"AG6$\"\"\"F,F,&%\"xG6#F,F,F,*&&F*6$F,\"\"#F,&F.6#F3 F,F,*&&F*6$F,\"\"$F,&F.6#F9F,F,,(*&&F*6$F3F,F,F-F,F,*&&F*6$F3F3F,F4F,F ,*&&F*6$F3F9F,F:F,F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "typ e(\",vector);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%%trueG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "B:=array(1..3,1..2,[]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-%&arrayG6%;\"\"\"\"\"$;F)\"\"#7\"" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalm(B);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%'MATRIXG6#7%7$&%\"BG6$\"\"\"F+&F)6$F+\"\"#7$&F)6$F. F+&F)6$F.F.7$&F)6$\"\"$F+&F)6$F7F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "evalm(A &* B);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%' MATRIXG6#7$7$,(*&&%\"AG6$\"\"\"F-F-&%\"BGF,F-F-*&&F+6$F-\"\"#F-&F/6$F3 F-F-F-*&&F+6$F-\"\"$F-&F/6$F9F-F-F-,(*&F*F-&F/F2F-F-*&F1F-&F/6$F3F3F-F -*&F7F-&F/6$F9F3F-F-7$,(*&&F+F5F-F.F-F-*&&F+FAF-F4F-F-*&&F+6$F3F9F-F:F -F-,(*&FHF-F>F-F-*&FJF-F@F-F-*&FLF-FCF-F-" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "evalm(B &* A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#- %'MATRIXG6#7%7%,&*&&%\"AG6$\"\"\"F-F-&%\"BGF,F-F-*&&F+6$\"\"#F-F-&F/6$ F-F3F-F-,&*&F.F-&F+F5F-F-*&F4F-&F+6$F3F3F-F-,&*&F.F-&F+6$F-\"\"$F-F-*& F4F-&F+6$F3F@F-F-7%,&*&&F/F2F-F*F-F-*&&F/F;F-F1F-F-,&*&F8F-FGF-F-*&F:F -FIF-F-,&*&FGF-F>F-F-*&FIF-FBF-F-7%,&*&&F/6$F@F-F-F*F-F-*&&F/6$F@F3F-F 1F-F-,&*&FSF-F8F-F-*&FVF-F:F-F-,&*&F>F-FSF-F-*&FBF-FVF-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "A:=randmatrix(4,6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'MATRIXG6#7&7(!#&)!#b!#P!#N\"#(*\"#]7(\"#z \"#c\"#\\\"#j\"#d!#f7(\"#X!\")!#$*\"##*\"#V!#i7(\"#x\"#m\"#a!\"&\"#**! #h" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "[rowdim(A),coldim(A)] ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$\"\"%\"\"'" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 58 "for j from 1 to 6 do A[3,j]:=A[1,j]+2*A[2,j] od: evalm(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7&7(!#&) !#b!#P!#N\"#(*\"#]7(\"#z\"#c\"#\\\"#j\"#d!#f7(\"#tF3\"#h\"#\"*\"$6#!#o 7(\"#x\"#m\"#a!\"&\"#**!#h" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "for j from 1 to 6 do A[4,j]:=3*A[1,j]-5*A[2,j] od: evalm(A);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7&7(!#&)!#b!#P!#N\"#(*\"#] 7(\"#z\"#c\"#\\\"#j\"#d!#f7(\"#tF3\"#h\"#\"*\"$6#!#o7(!$]'!$X%!$c$!$?% \"\"'\"$X%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "rank(A);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "N:=nullspace(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %\"NG<&-%'VECTORG6#7(\"\"!F*#\"$b\"\"#))#!$n#\"$;'F*\"\"\"-F'6#7(F1F*# !$&=\"#W#\"$@'\"$3$F*F*-F'6#7(F*F1#!$:#F-#\"#*)F-F*F*-F'6#7(F*F*#\"$z& F7#!%JMF:F1F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "evalf(\",6 );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#<&-%'VECTORG6#7(\"\"!$\"\"\"F($! '=VC!\"&$\"'O65F-F(F(-F%6#7(F(F($\"'Oh " 0 "" {MPLTEXT 1 0 50 "evalf(linsolve(A,vector(4,[0,0,0,0] ),'r',lambda));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'VECTORG6#7(&%'la mbdaG6#\"\"\"&F(6#\"\"#,*F'$\"+XSw))R!\"*F+$\"+nmmm;F1&F(6#\"\"$$!+7/n jSF1&F(6#\"\"%$!+B9'3@$!\")F4F9,*F'$\"+$R&o^YF1F+$\"+LLLLBF1F4$!+0h62B F1F9$!+`u.qDF>" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 93 "nullity:= proc(A) if type(A,matrix)=false then 'error' else RETURN(coldim(A)-ran k(A)) fi; end;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(nullityG:6#%\"AG6 \"F(F(@%/-%%typeG6$9$%'matrixG%&falseG.%&errorG-%'RETURNG6#,&-%'coldim G6#F.\"\"\"-%%rankGF9!\"\"F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "nullity(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"%" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "R:=rowspace(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"RG<$-%'VECTORG6#7(\"\"\"\"\"!#!$B'\"$:%#!$, $\"#$)#!%n&)F.#\"#*)F1-F'6#7(F+F*#\"%U7F.#\"$=&F1#\"&3D\"F.#!$8#F1" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "T:=array(1..2,1..4,[]); " } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"TG-%&arrayG6%;\"\"\"\"\"#;F)\"\"% 7\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "for i from 1 to 2 do for j from 1 to 4 do T[i,j]:=evalm(R[i] &* N[j]) od od: print(T);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7&\"\"!F(F(F(F'" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "C:=colspace(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"CG<$-%'VECTORG6#7&\"\"\"\"\"!F*\"\"$-F'6#7& F+F*\"\"#!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "transpose( A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7(7&!#&)\"#z\"#t!$] '7&!#b\"#c\"#d!$X%7&!#P\"#\\\"#h!$c$7&!#N\"#j\"#\"*!$?%7&\"#(*F/\"$6# \"\"'7&\"#]!#f!#o\"$X%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "N T:=nullspace(\");" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#NTG<$-%'VECTOR G6#7&!\"\"!\"#\"\"\"\"\"!-F'6#7&!\"$\"\"&F-F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "T:=array(1..2,1..2,[]): " }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 87 "for i from 1 to 2 do for j from 1 to 2 do T[i, j]:=evalm(C[i] &* NT[j]) od od: print(T);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$\"\"!F(F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "A:=randmatrix(4,4);" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%\"AG-%'MATRIXG6#7&7&!#`!#h!#B!#P7&\"#J!#M!#U\"#))7&!#w!#l\"#D\"#G7 &F+!#g\"\"*\"#H" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "Eig:=eva lf(Eigenvals(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$EigG-%'VECTORG 6#7&$\"+<[6ls!\"),&$!++1DtaF+\"\"\"%\"IG$\"+\"z)>xbF+,&F-F/F0$!+\"z)>x bF+$\"+uP'Q\"Q!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "det(A -lambda);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,,\"(0>p\"\"\"\"%'lambdaG !'xlV*$F&\"\"#!%()>*$F&\"\"$\"#L*$F&\"\"%F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "charpoly(A,lambda);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,,\"(0>p\"\"\"\"%'lambdaG!'xlV*$F&\"\"#!%()>*$F&\"\"$\"#L*$F&\" \"%F%" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 67 "The eigenvalues are zero s of the characteristic polynomial, however" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "for i from 1 to 4 do subs(lambda=Eig[i],charpoly(A ,lambda)) od;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"#5!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$!#:!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$\"#:!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\" 'r<9!\"(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "Eig:=evalf(Eige nvals(A),20):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "for i from 1 to 4 do subs(lambda=Eig[i],charpoly(A,lambda)) od;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$! \"\"!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$\"\"\"!\"#" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$!%)\\#!\"(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "Eig:=evalf(Eigenvals(A),100):" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 63 "for i from 1 to 4 do subs(lambda=Eig[i],char poly(A,lambda)) od;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$!\"\"!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%\"IG$\"\"\"!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$ !%)\\#!\"(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "0 0 0" 20 }{VIEWOPTS 1 1 0 1 1 1803 }