MVR4
I'latool=6"%&falseGE\[lC%$pluG:6#%"AG61%"mG%"nG%"PG%"LG%"UG%'pivcolG%&signPG%$t
olG%#rpG%"rG%"pG%"iG%"jG%"kG%%tempG6#%*CopyrightG6"C->8$-%'rowdimG6#9$>8%-%'col
dimGFC>8&-%%bandG6$7#"""F@>8'FK>8(-%'matrixG6%F@FF""!>8)7">8+#FO"*++++">8*FO>8.
FO?(81FOFO-%$minG6$F@FF%%trueGC)>8,-%(findpivG6&-%*submatrixG6%FD;F^oF@;F\oFFF^
oF\oFfn>8-&Feo6#FO>F\o&Feo6#""#@$0F_pF^oC&>FD-%(swaprowG6%FDF^oF_p?(8/FOFO,&F^o
FO!""FOFboC%>82&FQ6$F_pF^q>Fdq&FQ6$F^oF^q>FgqFcq>FJ-F[q6%FJF^oF_p>Fjn,$FjnF`q@$
2Ffn-%$absG6#&FD6$F^oF\oC$>FY7$&FYF=F\o?(F^q,&F^oFOFOFOFOF@FboC$>&FQ6$F^qF^o*&&
FD6$F^qF\oFOFdrF`q?(80F[sFOFFFbo>&FD6$F^qFds,&FfsFO*&F^sFO&FD6$F^oFdsFOF`q?(Fds
F^oFOFFFbo@%2Ffn-Fbr6#Fjs>&FSF[tFjs>FbtFW@$2F\oFF>F\o,&F\oFOFOFO7'-%%evalG6#FJ-
Fjt6#FQ-Fjt6#FS-Fjt6#FY-Fjt6#FjnF=F=%%elimG:F)6(%"EG%"RG%#REG%$eyeG%$augGF-F;F=
C)>FYFA>FQ-FL6$FNFY>FS-%(augmentG6$FDFQ>FJ-%%rrefGF_u>FF-Fjo6%FJ;FOFYF[w>F@-Fjo
6%FJF[w;,&FYFOFOFO,$FYFep7$-Fjt6#F@-Fjt6#FFF=F=%'particG:6$F*%"bG6*F,F-FhuF6F5F
7%"dG%"xGF;F=C)>F@FA>FFFG>FJ-Fgv6#-Fcv6$FD9%>FQ-%)pivotcolGF[u>FS-%%nopsGF]u@%/
&FQF_u,&F@FOFOFO>FfnFZC%>Ffn-%'vectorG6$FFFW>Fjn-%$colG6$-Fjo6%FJ;FOFS;,&FFFOFO
FOF_zFO?(FYFOFOFSFbo>&Ffn6#&FQFau&FjnFau-Fjt6#FfnF=F=%'eigvalG:F)6&F7%(eiglistG
%(evaluesG%(repeatsGF;F=C'>FF7#-%+eigenvectsG6$FD%(radicalG>FF-%%sortG6$FF:6$F]
x%"yGF=6$%)operatorG%&arrowGF=-%&evalbG6#2&FDFdp&FDFapF=F=>FJ7#-%$seqG6$&&FFFdw
Fap/F@;FO-F\yFfw>FQ7#-Fh\l6$&F[]lFdpF\]l7$FitF\uF=F=%)colbasisG:F)6%FhuF,F6F;F=
C&>F@-FgvFC>FFFA>FJ-FixFdw-Fjo6%FD;FOFFFJF=F=%&plot2GFa^l%(linefitG:6$%"tGFjw6/
F-F*%%xhatGF]x%"cGF\x%#tbG%%tseqGF7%#p1G%#p2G%%tminG%%tmaxGF;F=C.>F@-%(vectdimG
FC>FF-Fcv6$-Fey6$F@FOFD>FJ&-%$lsqG6$FFFfxFap>FS&FJFap>FY&FJFdp>Fjn7#-Fh\l6$7$&F
D6#Feo&FfxFg`l/Feo;FOF@>Ffn-Fh\l6$Ff`lFi`l>F^q,&-F`oFgz$"#6F`q-%$maxGFgz$FOF`q>
Fds,&F`alFealFcalFaal>F_p-%%plotG6%Fjn/FQ;F`alFcal/%&styleG%&pointG>F\o-Fjal6$,
&FSFO*&FYFOFQFOFO/FQ;F^qFds-%(displayG6#<$F_pF\oF=F=%)fourbaseG:6#'F*FU6+%$ROWG
%"NG%$COLG%#LNGFguFhu%#ERGF5F1F;F=C,>Ffn-%%rankGFC>Fjn-FduFC>FS&FjnFap>FY&FjnFd
p>Feo-FixFau>F@-%*transposeG6#-Fjo6%FY;FOFfn;FO-FHFau>FF-%)nulbasisGFau>FJ-Fjo6
%FD7#-%"$G6#;FOFAFeo>FQ-Fddl6#-Fjo6%FS;,&FfnFOFOFO-FBF_u;FO-FHF_u7&FcwFewFitF\u
F=F=%(projectG:Fiw6%F.F6%"eGF;F=C'>F@-%(projmatGFC>FF-%)multiplyG6$F@Ffx-%&prin
tG6%F@FfxFF>FJ-%$addG6&FfxFFFOF`q7$FewFitF=F=%&houseG:F=6#%"XGF;F=>F@-FU6%Fep"#
77:!"'F^hl!"(FW""(""'Fahl!"$FbhlFWFWF^hlF_hlFepFO"")FOFepF_hlF_hl!"#FdhlF_hlF_h
lF=F=%)leftnullG:F)F=F;F=-F]el6#-FddlFCF=F=Fhfl:F)6#F.F;F=C$>FD-Fe]lFC>F@-F[gl6
%FD-%(inverseG6#-F[gl6$FihlFDFihlF=F=%$sluG:F)6)F-F/F0F3F7F8F9F;F=C(>F@FA>FF-FU
6%F@F@FW>FJF]jl>FQFgn?(FjnFOFOF@FboC&@$1-Fbr6#&FD6$FjnFjnFQ-%&ERRORG6$%CSmall~p
ivot~encountered~in~column~GFjn>&FFFhjlFO?(FS,&FjnFOFOFOFOF@FboC$>&FF6$FSFjn*&&
FDFd[mFOFgjlF`q?(FYF`[mFOF@Fbo>&FD6$FSFY,&Fi[mFO*&Fc[mFO&FD6$FjnFYFOF`q?(FYFjnF
OF@Fbo>&FJF^\mF]\mFdglF=F=F]el:F)6-FhuF1F,F-F5%*dimensionGFbcl%(freecolGF7F8%"f
GF;F=C,>F@Fj]l>FFF]^l>FJFA>FQFG>FS-F\y6#F6>FY,&FQFOFSF`q>Fjn-FU6%FQFYFW>Ffn-%&m
inusG6$<#-Fcel6#;FOFQ<#&FFF=?(F_pFOFOFYFboC&>F\o-F`o6#&FfnF=>Ffn-Ff]m6$Ffn<#F\o
>&Fjn6$F\oF_pFO?(FeoFOFOFSFbo>&Fjn6$&FFFg`lF_p,$&F@6$FeoF\oF`qFbuF=F=%$slvG:Fiw
6,%#LUGF/F0F-Fh^lF9F8%"sGFe^lF]xF;F=C+>F@-FgilFC>FF&F@Fap>FJ&F@Fdp>FQFA>FS-Fey6
$FQFW?(FYFOFOFQFboC%>FfnFW?(FjnFOFO,&FYFOF`qFOFbo>Ffn,&FfnFO*&&FF6$FYFjnFO&FSFc
uFOFO>&FSFau,&&FfxFauFOFfnF`q>F_pFa`m?(FYFQF`qFOFboC%>FeoFW?(FjnF`wFOFQFbo>Feo,
&FeoFO*&&FJF\amFO&F_pFcuFOFO>&F_pFau*&,&F_amFOFeoF`qFO&FJ6$FYFYF`q-Fjt6#F_pF=F=
%%spluG:F)6.F,F-F.F/F0F2F3F5F7F8F9F:F;F=C,>F@FA>FFFG@$0F@FF-Fjjl6#%7Matrix~must
~be~square.G>FJ-FL6$FNFF>FQF`cm>FS-FU6%FFFFFW>FjnFgn>FYFO?(F\oFOFOFFFboC%@$1-Fb
r6#&FD6$F\oF\oFjnC'?(FfnF\oFOFFFboC$@$2Fjn-Fbr6#&FD6$FfnF\o%&breakG@$/FfnFFC$>F
YFW-%'RETURNG6&FitF\uF^uF`u>FD-F[q6%FDF\oFfn?(FeoFOFO,&F\oFOF`qFOFboC%>F^q&FQ6$
FfnFeo>Fhem&FQ6$F\oFeo>F[fmF^q>FJ-F[q6%FJF\oFfn>FY,$FYF`q?(FeoFgtFOFFFboC$>&FQF
b_m*&&FDFb_mFOF^dmF`q?(F_pFgtFOFFFbo>&FD6$FeoF_p,&F[gmFO*&FffmFO&FDFj^mFOF`q?(F
_pF\oFOFFFbo@%2Fjn-Fbr6#F_gm>&FSFj^mF_gm>FfgmFW7&FitF\uF^uF`uF=F=%)randpermG:6#
F-F=F;F=-&%)combinatG6#FigmFCF=F=%%splvG:Fiw6-%$PLUGF.F/F0F-Fh^lF9Fg_mF8F]xFe^l
F;F=C->F@-FdbmFC>FFF\`m>FJF^`m>FQ&F@6#""$>FSFA>Ffx-F[glF\`l>FY-Fey6$FSFW?(FjnFO
FOFSFboC%>FfnFW?(FeoFOFO,&FjnFOF`qFOFbo>Ffn,&FfnFO*&&FJ6$FjnFeoFO&FYFg`lFOFO>&F
YFcu,&&FfxFcuFOFfnF`q>F_pFaim?(FjnFSF`qFOFboC%>F\oFW?(FeoF`[mFOFSFbo>F\o,&F\oFO
*&&FQF\jmFO&F_pFg`lFOFO>F[bm*&,&F_jmFOF\oF`qFO&FQFhjlF`qFbbmF=F=%)rowbasisG:F)F
=F;F=-Fe]lFhhlF=F=%%initG:F=F=F=F=C$-%%withG6#%'linalgG-Fg[n6$%&plotsGFiblF=F=%
(symmeigG:F)6'%#SLG%"SG%#QRG%"QG%'LambdaGF;F=C$@%1-%%normG6#-Fbgl6&FDFihlFOF`qF
gnC'>F@-%'eigvecGFC>FFF\`m>FSF^`m>FJ-%&gramsGFfw>FQF^`l-Fjjl6#%=The~matrix~is~n
ot~symmetric.G7$F\uF^uF=F=%'eigen2G:F)6(F\xFe^l%#e1G%#e2G%#x1G%#x2GF;F=C0>F@,&*
&&FD6$FOFOFO&FD6$FepFepFOFO*&&FD6$FOFepFO&FD6$FepFOFOF`q>FF,&Ff^nFOFh^nFO>FJ,&F
F#FOFep-%%sqrtG6#,&*$FFFepFOF@!"%Fc_n>FQ,&FFFc_nFd_n#F`qFep@'0F[_nFWC$>FS-Fey6#
7$F[_n,&FJFOFf^nF`q>FY-Fey6#7$F[_n,&FQFOFf^nF`q0F]_nFWC$>FS-Fey6#7$,&FJFOFh^nF`
qF]_n>FY-Fey6#7$,&FQFOFh^nF`qF]_nC$>FS-Fey6#7$FOFW>FY-Fey6#7$FWFO-F^gl6#%^oFor~
this~matrix,~the~polynomial~whose~roots~are~the~eigenvalues~is:G-%'printfG6%%9~
~~~~e^2~-~%g*e~+~%g~=~0GFFF@-F^gl6#%!G-F^gl6#%JThe~first~eigenvalue~and~eigenve
ctor~are:G-F^glF[u-F^glF_u-F^gl6#%KThe~second~eigenvalue~and~eigenvector~are:G-
F^glF]u-F^glFauF=F=%)samespanG:6$%#A1G%#A2G6%%'rankA1G%'rankA2G%)rankbothGF;F=C
'>F@Fhcl>FF-Ficl6#Ffx>FJ-FiclFcx@%3/F@FF/F@FJ-F^gl6#%FA1~and~A2~have~the~same~c
olumn~space.G-F^gl6#%HA1~and~A2~have~different~column~spaces.G-F_em6#%%NULLGF=F
=%(tridiagG:6&%"aGFjwFh^lF-F=F;F=-FL6$7%FDFfx9&9'F=F=Fe]n:F)6)Fc\nFhuF,F-F8F9%%
multGF;F=C(>F@FD>FJ-FBFdw>FQ-FHFdw?(FSFOFOFQFboC$?(FYFOFO,&FSFOF`qFOFboC$>Fjn-%
(dotprodG6$-Fiy6$F@FS-Fiy6$F@FY>F@-%'addcolG6&F@FYFSF^r>F@-%'mulcolG6%F@FS*$-Fi
\n6$F[gnFepF`q>FF-F[gl6$-FddlFdwFDFbwF=F=%'determG:F)6&F7F0%%signG%)splu_outGF;
F=C&>FQFfhm>FF&FQF[im>FJ&FQ6#""%*&FJFO-%$mulG6$&FF6$F@F@/F@;FO-F`o6$-FHFfw-FBFf
wFOF=F=Fgo:6&F*F9F6F36&F7F8F,F-F;F=C(>F@FO>FFFO>FJFA>FQFG?(F=FOFOF=2-Fbr6#&FDFa
oFhen@'2F@FJ>F@F`y2FFFQC$>F@FO>FFF_z-F_em6#7$FfxFgen7$-Fjt6#,(F@FOFfxFOF`qFO-Fj
t6#,(FFFOFgenFOF`qFOF=F=%'cramerG:Fiw6'F7F-F]x%"BG%#b1GF;F=C*>FFFA@$0FFFG-Fjjl6
#%4A~not~square~matrixG@$0FF-Fb_lF`dn-Fjjl6#%:b~not~same~dimension~as~AG>FJFdy>
FQFdcm>FS-FU6%FFFOFfx?(F@FOFOFFFboC%-%)copyintoG6&FDFQFOFO-F_]o6&FSFQFOF@>&FJFd
w*&-%$detGF]uFO-Fg]oFCF`qFitF=F=Fix:6#Fhu6'F7F8F,F-F6F;F=C(>FJFA>FQFG>FSFZ>FFFO
?(F@FOFOFJFboC&?(F=FOFOF=1FFFQC$@$0FcjnFWFidm>FFF_z@$2FQFFFidm>FS7$&FSF=FF>FFF_
zF^uF=F=Fcil:F)6'F,F-F5%%AinvGFeclF;F=C'>F@FA>FFFG>FJFhcl@%3/FJF@/FJFFC$>FSF[dl
>FQ&FSFap>FQFZF\uF=F=%$carG:F)6*FeclFhuFa\nF1%'pivrowGF5Fh^lFcenF;F=C,>F@F[dl>F
FF^`m>FQ-FixFfw>F@-FduFhhl>FJF^`m>FSFhx>FY-Fjo6%FF;FOFhcl;FOFG>Fjn-Fddl6#-Fjo6%
FJF_aoF`ao>Ffn-Fjo6%FDFSFQ-Fjt6#7%FjnFfnFYF=F=F`]n:F)6+F[[lF\[l%)repeats2G%-geo
m_repeatsG%"DGFa\n%)eigvectsGF7F8F;F=C+>F@7#-Fb[lFC>F@-Fg[l6$F@:Fj[lF=F\\lF=F_\
lF=F=>FF7#-Fh\l6$&&F@FgzFap/Ffn;FO-F\yFdw>FJ7#-Fh\l6$-F\y6#&F_coF[imF`co>FQ7#-F
h\l6$-Fh\l6$&FFFgz/Feo;FO&FJFgz/FfnF]]l>FS-%%diagG6#&FQF=>Fjn-Fh\l6$-%#opGFhcoF
`co>FY-FcvFcu7$F`uF^uF=F=F[`l:Fiw6%Fg^lF6FdflF;F=C'@$/-Fg]oFdilFW-Fjjl6#%DColum
ns~of~A~are~linearly~dependentG>F@-Fgw6$Feil-F[gl6$FihlFfx>FF-F[gl6$FDF@>FJFagl
7%FcwFewFitF=F=%(normal2G:F)6(FguF`\nFa\nF0%#URGFd\nF;F=C%>F@-Fbgl6&Feil-F[gl6$
FDFihlFOF`q@%1-Fi\nFdwFgnC'>FFF_]n>FJ&FFFap>FY&FFFdp>FS-Fe]nF[u>FQF]`oC%>FQFZ>F
YFZ-Fjjl6#%:The~matrix~is~not~normal.G7$F\uF`uF=F=6"