{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 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 4 "?Svd" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new defin ition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "A:=randmatrix (3,2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7$!#&)!# b7$!#P!#N7$\"#(*\"#]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "eva lf(Svd(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7$$\"+Paan: !\"($\"+e3TX8!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "read(` a:svd.txt`);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$SVDGR6&'%\"AG%'matr ixG'%\"UG%&evalnG'%\"SGF,'%\"VGF,6&%\"mG%\"nG%\"sG%\"iG6\"F6C(>8$-&%'l inalgG6#%'rowdimG6#9$>8%-&F<6#%'coldimGF?>8&-%&evalfG6#-%$SvdG6%F@9%9' >9&-F)6%F9FB7#-%$seqG6$\"\"!/8';\"\"\"*&F9FgnFBFgn?(FenFgnFgn-%$minG6$ F9FB%%trueG>&FR6$FenFen&FH6#Fen-%&printG6#%\"~GF6F6F6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "SVD(A,U,Sigma,V);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%\"~G" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 " evalm(U),evalm(Sigma),evalm(V);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6%-%' matrixG6#7%7%$!+Wk*fX'!#5$!++27d@F*$\"+GBxDtF*7%$!+YsPwJF*$!+AJ8lzF*$! +JijW^F*7%$\"+7`$[%pF*$!+upJ[cF*$\"++U5dWF*-F$6#7%7$$\"+Paan:!\"(\"\"! 7$FD$\"+e3TX8!\")7$FDFD-F$6#7$7$$\"+_)zza)F*$!+\"*)3'*=&F*7$$\"+\"*)3' *=&F*FN" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "u1:=col(U,1); u2 :=col(U,2); u3:=col(U,3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u1G-%' vectorG6#7%$!+Wk*fX'!#5$!+YsPwJF+$\"+7`$[%pF+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u2G-%'vectorG6#7%$!++27d@!#5$!+AJ8lzF+$!+upJ[cF+" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u3G-%'vectorG6#7%$\"+GBxDt!#5$!+Jij W^F+$\"++U5dWF+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "v1:=col( V,1); v2:=col(V,2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#v1G-%'vector G6#7$$\"+_)zza)!#5$\"+\"*)3'*=&F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %#v2G-%'vectorG6#7$$!+\"*)3'*=&!#5$\"+_)zza)F+" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 25 "multiply(transpose(V),V);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7$$\"+**********!#5\"\"!7$F+F(" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "B1:=scalarmul(multiply(u1,tr anspose(v1)),Sigma[1,1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#B1G-%' matrixG6#7%7$$!+QNh]')!\")$!+n$>>D&F,7$$!+K!QhD%F,$!+*GmRe#F,7$$\"+!pD cI*F,$\"+Ise\\cF," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "B2:=sc alarmul(multiply(u2,transpose(v2)),Sigma[2,2]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#B2G-%'matrixG6#7%7$$\"+c`81:!\"*$!+[j!3[#F,7$$\"+A.Q hbF,$!+;rLg\"*F,7$$\"+5JuVRF,$!+,B(e\\'F," }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "matadd(B1,B2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#- %'matrixG6#7%7$$!+-+++&)!\")$!+-+++bF*7$$!+++++PF*$!+,+++NF*7$$\"+,+++ (*F*$\"+++++]F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "B:=randm atrix(2,3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-%'matrixG6#7$7% \"#z\"#c\"#\\7%\"#j\"#d!#f" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "SVD(B,U1,S1,V1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%\"~G" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "evalm(U1),evalm(S1),evalm(V1 );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6%-%'matrixG6#7$7$$!+fG)=U(!#5$!+m U\">q'F*7$F+$\"+fG)=U(F*-F$6#7$7%$\"+n[B'G\"!\"(\"\"!F77%F7$\"+3hf-x! \")F7-F$6#7%7%$!+]y4TyF*$!+>_pK!)!#6$!+5%=S:'F*7%$!+:MJ,iF*$\"+Nk#z>'F D$\"+j$e/#yF*7%$\"+70snCFD$!+ " 0 "" {MPLTEXT 1 0 30 "multiply(U1,S1,transpose(V1));" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7%$\"+.+++z!\")$\"+.+++cF*$\"+**** ****[F*7%$\"+.+++jF*$\"+-+++dF*$!+********eF*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "18 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }