{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 56 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 29 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 "Maple Output" 0 11 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 2 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 "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 179 0 4 0 0 
0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 257 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 2 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 
0 0 0 0 -1 0 }{PSTYLE "" 0 258 1 {CSTYLE "" -1 -1 "" 0 1 74 0 29 0 0 
0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 260 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 
0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "Matrix Algebra worked exam
ple 5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 19 "Consider the system" }}
{PARA 256 "" 0 "" {OLE 1 3584 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:
F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::::::::::::::
::::::::::::::::::fyyyyya:nYf::wyyyqy;::::::::::::::::::::::::::::::::
::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_mvJ:::
:::::gj<vyyuyJ<j?V:;J:>R<:T><::;Xr[TM;t_wW:A:;B;B:F:YLpfF>:::::::::J?N
Z;vyyyyyY:vYxY:B:::::::c:;:=:jR>@Wlj^HMMufF;J:::::::N=?:xI:;Z::::::j:>
:C=;:=j[vGUMrvC?MoJ::::::::JCN:yyyxI:;Z::::::j;B:s<;:wA?Z::C:J=j=B:K:M
:OJ:V;;J@>:UJ:n;v;;Jyky;::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::::::::::::::::::j:b:<Z::::::::::JjB:CR:Hj<dJ<@j:@JSHRvTJ;NZ
@ha<>klN\\=@i:FZlVjsoZ<:Zs?:<ZxRCD]Hr;R:[DZGcMB::::::::::FZ:nyyYZDjysy
Qj;J:>R<:TNC>:UTRcETcTX[US;SK]UW=EWMuUWm>j\\@A:::cYH_WV>Z::::::::jysyA
:CB:F:;jysy?:;:gI;vo>JmDZ:F:MZ=V`;B:?jDJ:f:;j>D:;B:CZ:NZ;F:E:=b:yyyyI:
E:M:<J<JQ@JFf:F`<F:<j<j?D:<j<J@DJZLJ:MJA>:Gc;YJCvYxY:<J>JZyyayYZLjALj?
JMJ@fc[;>:;I:EZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:vYxY:>Z::::::j
<B:Wf:B:N:Y<>D_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jysy:>:
<:::::::::::::::::::vYxI:;Z:::::::::]:qi:;fy>Z:JBAZ:b:DJ:<k`LOJPNvLQs@
`I^:f_;j><J:<JSdj\\@JrFZ:B:GC:Cn>JSfj=nH=:G;eZ:VY;><:[V:=:;b:^dcgg_WhZ
nc_whZNdigg[oGuM<jPF:C:[Y:F;;JSFd;^C=Z:>IKG:M:_;Ku:=:k=JSFl=>E=B:FIK?N
@>j=FZ:jvJB<Z:VYZ:JBAJ:b:DZJVdscRYEU`K<jP>:C:[q:F;;JSdjb^e;F:;B:uKfG;N
@]Nn@j:jXjDjw;<:[V:B:D:c\\_;:^?CB:f?=J<Jr?:MJ:N@]L?Ej:>Z:B:;H:M:djBol<
F:B:N;<JMT:OjyyiyI:>:[Z:VY[j=J:^q<Z:FZ<NZ\\VdsWhngfgcUC:UK:^:>X=j;<Z:>
:::::::::3:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 29 "find a) the equili
brium point" }}{PARA 0 "" 0 "" {TEXT -1 96 "      b)  the eigenvalues
\015\n      c) the eigenvectors\015\n     d) the equations of the dire
ctrices\015\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }
}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 
-1 29 "To find the equilibrium point" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 7 "X:=2*x;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"XG,$%\"
xG\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "Y:=2*y;" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"YG
,$%\"yG\"\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 40 "The equilibrium p
oint is obviously (0,0)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "
" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 23 "To find the eigenvalues" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "A:=matrix(2,2,[2,0,0,2]);" }
}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'MATRIXG6#7$7$\"\"#\"\"!7$F+F
*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "eigenvalues:=eigenvals
(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%,eigenvaluesG6$\"\"#F&" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 32 "Thus the eigenvalue is 2 twice. " 
}}{PARA 0 "" 0 "" {TEXT -1 117 "Since the matrix is diagonal there is \+
an infinite number of eigenvectors and hence an infinite number of dir
ectrices." }}{PARA 0 "" 0 "" {TEXT -1 24 "To find the eigenvectors" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "eigensystem:=eigenvectors(A
);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%,eigensystemG7%\"\"#F&<$-%'VEC
TORG6#7$\"\"!\"\"\"-F)6#7$F-F," }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 
131 "Taken at face value this imlis that there are only two independen
t eigenvectors.That this is incorrect may be shown as follows. Let" }}
{PARA 257 "" 0 "" {TEXT -1 1 " " }{OLE 1 3584 1 "[xm]Br=WfoRrB:::wk;ny
yI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::::
::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;::::::::::::::::
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^
H;C:ELq^H_mvJ::::::::gj<vyyuyJ<j?V:;J:>R<:T><::;Xr[TM;t_wW:A:;B;B:F:YL
pfF>:::::::::J?NZ;vyyyyyY:vYxY:B:::::::c:;:=:jR>@Wlj^HMMufF;J:::::::N=
?:xI:;Z::::::j:>:C>;:=j[vGUMrvC?MoJ::::::::JCN:yyyxI:;Z::::::j;B:s<;:w
A?Z:F;^:nYn:v:<J>J?>:wAQ:SJ:f;;JAjA>:wAyA:::::::::::::::::::::::::::::
:::::::::::::::::::::::::::::::::::::::::::::=Z<B:<::::::::::;^:^Z;r:E
b<=r:AZ;fZ;FZ;NZ;j:Dk;Ip;xj:V[;NZBH@]B:dj:Dk;;JSF<<J:dJ;TsgDJEq?La_bFt
Z<Z;>d:Xb[og]gffk=v:;J:::::j:<Jyyy;d:yayY;AZ:>:;P:<:TNC>:UTRcETcTX[US;
SK]UW=EWMuUWm>j\\@A:::cYH_WV;Z::::::::jysyA:CB:F:;jysy?:;:Wt:F]?JTD:=j
>r:oR:N:eZ:>:EJ:F[<J:<J<J;@j:f:FZ<vyyyy:f:F;B:C:WS:kj<jRDj:B:E:Qb:<Z:f
:^[<>n>r:WJ:n^@v;_jysy;B:K:;xyAwy;h:AD;Q:G;Sj`@@J:>XEZ:F;NZ:vCS=[LsfFa
MR>@>Z::::::::kJ;@:;J;>Z:vYxY:>Z::::::j<B:sd:B:N:Y<>D_=a=[;;B:::::::JF
>:yay=J:B::::::nYyA<::::::::::::jysy:>:<:::::::::::::::::::vYxI:;Z::::
:::::eB:<jwEJ:uI<:[V:=:;b:<Kd`ppPps<LaPpu<Lcxpp@Pqf<CJ:f_;j><J:<JSdjHA
JmG:;B:;:::::::::::AE:Ko>JSfu=^T=:A=eJ:VY;syB:>l;>:<Z<:::fICJ:f?=J<JrA
:MJ:N@uY:mU:=J:<jR<JrF;N@uC;eU:=B:B:[MBB:qQ:pI<:[V:B:DZ:nb_og]gfmS@C:U
K:^:>X=j>>:_c<;X:;l:<Z:^C[F:MZ:N@yC;;l:Ja<JBG;<JSF`=>B=:UD:[n>B:_KZG]<
F:;ja>l:F;<JSvAmb:=:]LBG;<JSF@mb:=:[LBO[:B:G;OjyyiyI:;J:><<jw?<IZ:>:sg
:>Z:FZ<NZ\\VdsWhngfgcUC:Uk:^:>X;j;<Z:>:::::::1:" }}{PARA 0 "" 0 "" 
{TEXT -1 19 "be any vector. Then" }}{PARA 258 "" 0 "" {OLE 1 5120 1 "[
xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy:::::::::::::::::
::::::::::::::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::G:I
:K:wAyA:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::::::JcvGYMt>^:fBWMtNHm=;:::::::n<Ejyyiy=C:Qj;>:;JZC:bKB::>r
[AdG?bQpq;V:>bCZ:j:vCSmlJ::::::::::OJ;@jyyyyyy;jysy;Z:::::::^<>:F::]KR
nC=MtFGgml>:;::::::JGN:ry:>:<::::::=J:^b;Z:j:VBYmp>HYLkNG>::::::::N<?j
yyiy=J:B::::::V:<JhJ:JyK;B:M:CJ:nYn:v:<J>j@>:OJ:V;;J@>:wAW:YJ:>\\:B:]:
_J:V<^<f<n<;jEJFjFJG>:q:s:u:w:y:;C:<jJ<:?;A;C;wAyA::::::::::::::::::::
::::::::::::::::::::::::::::::::::::::::=Z<Z::::::::::JZ@:CR:Hj<dj:Hj;
R:ER:=R:?b:FZ<>klNZ@h_<^jfS:?JRF::_;J:NZPXe<nlWs:A:E:@:=b>AvRAr==j?@J;
\\rRF\\:b<=b>A>:_C:]Z:>ZDNZPXEb:GVIAr:=Jyyy;d:yayY;AZ:>:;`:Z@O<;B:UTRc
ETcTX[US;SK]UW=EWMuUWm>j\\@A:::cYH_WVB::::::::vYxy;J<<J:FZ:jysy?:;:=W@
F]?JMU:=j>r:at:NZ:f<<J:f:;j>D:<J<J;@j:f:FZ<vyyyy:f:F;B:C:WS:kj<jRDj:j<
j?DZ:B:E:Sb:;G;;t?W:Gc;YJCvYxY:<J>>:;xyAwyXk[Lj?B:G;Sj`@@J:>XEZ:F;NZ:v
CS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:vYxY:>Z::::::j<B:Wh;<J:<:?ja:[LsjsJR
J:<:::::::>=;jysy:>:<::::::wqy[:::::::::::::vYxI:;Z:::::::::::::::::::
:yay=J:B::::::::jD<Z:Vy<>jx]:JBAZ:B:DZJ^dcgg_WhZnc_whZNdigg[oGSM<<:US:
F[:>Z:N`Dfm;^T=Z:FZ?F:DJ;R:?B\\[k:::B:;b:g>V:FZ?:=J:@JC:ZkJ:N[;N<:bF@;
;B:Ob::?jXB:f:N:FZ<BF@:>Z:N:Z:R:iMKDj:<KRntJ^^<J:<Je<LR>Z:nT=L:L:<::::
:::::::::::JS<JZG;b<Vg@F:;B:[KZ=:N@ux:uY:=:[;JSJ[Pj:>:<:ZDfMcV==:[;JSJ
Ryj:>:[;JSfUcV==:[;JSJRyj:JR:N@u^OFZ:JR:N@u`OF:>@f\\:jw?^y]:JBA:<Z<>:c
<::uM<jPJ<JrA:MJ:N@uY:cc<=Z:F`:>HM:_kXLJRdj:<:<Jr:djx?dLF:;jR:N@uK`Dk:
Z:JBB:qQ:[:JBA:DZ<B>aTXDpql`S^:f?;J<JrG:MJ:N`DfMkX:=J:<jP<JfG;N@uP]Tj:
jPJBB:qQ:pI<:[V:B:DZ:nb_og]gfmS@C:US:=J<Jr?:MJ:N@;X:;l:<Z:^C[F:MZ:N@yC
;;l:Ja>LMZ:N@]s:;l:j`<JBG;<JS>RaD<=:YLB=j>B:_kYVdBF:FD[n>B:_kRVdBF:>D[
n>B:_KZWyBF:^CF;<JSvAqI<=:W<MZ:N@]kw]j:j`j>B:_KZ_tDF:vCF;<JSvAch<=:]<M
Z:N@]Kdej:JbJBB:qQ:[:JBA:DZ<B>jq:f?;J<JrG:MJ:<JSdJZEJ]ij:>Z:fAk@<jw?:B
:>l;Z:b:<::jPF:C:[Q:F;;JSR:;R==Z:^CF;<JSvA;R==:W<MZ:N@]s:;R==:U<MZ:N@;
pq<;;jaj>B:_kYvgJF:FDF;<JSF@YE>=:[<MZ:N@;pE@k:>:S<MZ:N@ykE@k:Jaj>B:_kR
v\\KF:fCF;<JS>R[i>=J:<jaj>B:_kY>xL:]<MZ:N@]KrEk:JbJBB:qQ:[:JBA:DZ<B>a<
Y=jP>:C:[q:F;;B:_c<;h:?y>=J:<jX:<jw?:B:>l;Z:b:<::jPF:C:[Q:F;;JSR:sR?=J
:<J`j>B:_kY^]OF:nCF;<JSF`=^]OF:fCF;<JS>RYe?=:Y<MZ:N@ykqTk:jb>LMZ:N@]kq
Tk:JbJBB:qQ:[:JBAj:>Z<b:<;vG:UK:^:>X=j>>Z:N`DfM[i:=B:;B:SMfG;N@uPiXj:J
pjDB:qA>Z:JBA:;J:DJ:^DP@Jc:SM<jPF:C:[Q:F;;JSfM=W<=Z:Vf:>OM:_khwkCF:VF:
_kHgqKF:;jk:N@upOAk:jkJ?B:G;OjyyiyI:>:[Z:VY[j=B:;JXE:;B:=b:?bBaTXaEWEU
UtP<>:UK:^:>X=j;<J:<J:::::::::::::::::::::::::::::::::::::::::5:" }}
{PARA 0 "" 0 "" {TEXT -1 17 "Hence the vector " }}{PARA 260 "" 0 "" 
{OLE 1 3588 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy:
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
fyyyyya:nYf::wyyyqy;::::::::::::::::::::::::::::::::::::::::::::::::::
::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_mvJ::::::::gj<vyyuyJ<j?V
:;J:>R<:T><::;iZkTM;t_wW:A:;B;B:F:YLpfF>:::::::::J?NZ;vyyyyyY:vYxY:B::
:::::c:;:=:jR>@Wlj^HMMufF;J:::::::N=?:xI:;Z::::::j:>:C>;:=j[vGUMrvC?Mo
J::::::::JCN:yyyxI:;Z::::::j;B:s<;:wA?Z:F;^:nYn:v:<J>J?>:wAQ:S:UJ:n;;j
A>:wAyA:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
:::::::::::=Z<B:<::::::::::;^:^Z;r:Eb<=r:AZ;fZ;FZ;NZ;j:Dk;Ip;xj:V[;NZB
H@]B:dj:Dk;;JSF<<J:dJ;TsgDJEq?La_bFtZ<Z;>d:Xb[og]gffk=v:;J:::::j:<Jyyy
;d:yayY;AZ:>:;P:<:TNC>:UTRcETcTX[US;SK]UW=EWMuUWm>j\\@A:::cYH_WV;Z::::
::::jysyA:CB:F:;jysy?:;:Wt:F]?JTD:=j>r:oR:N:eZ:>:EJ:F[<J:<J<J;@j:f:FZ<
vyyyy:f:F;B:C:WS:kj<jRDj:B:E:Qb:<Z:f:^[<>n>r:WJ:n^@v;_jysy;B:K:;xyAwy;
h:AD;Q:G;Sj`@@J:>XEZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:vYxY:>Z::
::::j<B:sd:B:N:Y<>D_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jy
sy:>:<:::::::::::::::::::vYxI:;Z:::::::::eB:<jwEJ:uI<:[V:=:;b:<Kd`ppPp
s<LaPpu<Lcxpp@Pqf<CJ:f_;B:MB:;B:_c<uV:Gq:>Z:>:::::::::::Vf:>OM:_khIJdG
:VFf<;jw?^y]:JBAJ:B:D:::jx^:;jPF:C:[Y:F;;JSfy;Fi;F:;B:]C:[m>JSfa>fh;FZ
:Z:>H[Z:VY;RyB:>l;Z:b:<J]PPqHplHaS^:f?;J<JrG:MJ:N`D>r;>B=B:<J`>l:F;<JS
va>>B=:WD:[n>B:_kRHJZF:fc:>LMZ:N@;pFDj:>:YLB=j>B:_kYF]<F:FD[n>B:_kRF]<
F:>D[N?<Z:n>N;yyyxy:>:;JBB:qQBv:<J:^q<J:<j:DJ;Dlc`qsLqlp`h_:f?=J<Jr?:A
B:<J::::::::1:" }}{PARA 0 "" 0 "" {TEXT -1 50 "is an eigenvector for a
ny values of the constants." }}{PARA 0 "" 0 "" {TEXT -1 78 "As each ei
genvector is a directrix their is an infinite number of directrices." 
}}}}{MARK "16 7 0" 78 }{VIEWOPTS 1 1 0 1 1 1803 }