{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 "" -1 256 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 257 "" 0 1 178 0 0 0 1 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 "" 0 256 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 257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 52 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 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 }{PSTYLE "" 0 259 1 {CSTYLE "" -1 -1 "" 0 1 0 0 7 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 }{PSTYLE "" 0 261 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 33 "Global Stability worked ex ample 2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 46 "Find first integrals for the dynamical sy stems" }}{PARA 256 "" 0 "" {OLE 1 3597 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;: j::>:B>N:F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq ^H_mvJ::::::::gjR<:T><::;GXEXDIDNwW:A:;x:B:F:YLpfF>::: ::::::J?NZ;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;nYvY::::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::=Zd:p^H[::::::::::j::;`:Z@O<;j`@Pt\\Pd`QrPPJPnrPqjLqnxPqF;fbk;::JtaMSAA;B:: ::::::vYxy;J<<:=J:vYxY;J:JuM:YG;N]D:;B:CZ:NZ;F:E:= b:yyyyI:E:M::Gc;YJCvYxY:JZyyay YZLjALj?JMJ@fc[;>:;I:EZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:vYxY:> Z::::::jD_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::: :jysy:>:<:::::::::::::::::::vYxI:;Z:::::::::]:qi:;fy>Z:JBAJ:UD O;SKkEW_URXM<>:US:F[:>Z:N`Dfb;>H=B:JSjDIj::C:[q:F;;JSdjb^e; F:;B:uC:kp>JSJCEj:>:UC:k@N@]Nn@j:jXJ?B:Gc;N;yyyxy:J:><:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;:::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::;tyGXDIDNwW:A :;v:B:F:YLpfF>:::::::::J?NZ;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:KJ:F;N;;j?>:S:UJ:nYvY:::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::j:b:<::::::::::=B:wyyAbR<:TNC>:UTRcETcTX[US;SK]UW=EWMuUWm>j\\ @A:::cYH_WV>Z::::::::jysyA:CB:;j:>:yayQ:>:N^>jQM:Kb:F:MZ=^Z;B:?jDJ:f:; j>D:;B:CZ:NZ;F:E:=b:yyyyI:E:M::Gc ;YJCvYxY:<:;xyYxy;v:YB;Q:G;Sj`@@J:>v:f:J;D_mlVH[KRJ:<::::: ::>=?R:>Z:N:;B:yayA:;B::::::f::yay=J:B ::::::nYyA<::::::::::::jysy:>:<:::::::::::::::::::vYxI:;Z:::::::::]:qi :;fy>Z:JBAZ:b:DJ:f<l;F:>Zx;j>>:_kb@J`F:< Jv>o:F;N@Nn C:UK:^:>X=j>>:_c<]Lh@j:>Z:fAkp>JSFLKU:=:u;OZ:n^@J?vyyuy=:;JBB:qQBv:>:s g:B:=b:?bBaTXaEWEUUtP " 0 "" {MPLTEXT 1 0 8 "restart:" }}} {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 6 "f:=-y;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "g:=x;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 91 " To find a first in tegral first divide g by f to obtain a differential equation in x and y" }}{PARA 258 "" 0 "" {OLE 1 3597 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j:: >:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;:::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_ mvJ::::::::gjR<:T><::;tyGXDIDNwW:A:;v:B:F:YLpfF>:::::: :::J?NZ;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:KJ:F;N;;j?>:S:UJ:nYvY:::::::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::j:b:Z;NZ@ppL VjsS:?:\\A?b[us:W:_rZIm:D[`[:RCD]Hr;R:[DZGcM>Z::::::::::j::=J:vYxY;J:j;Q:SF;nZr:AR::Mb:>Z:^::Gc;YJCvYxY:B:;xyewy;H;ex :Q:G;Sj`@@J:>v:F[:B:oi:NZ:vCS=[LsfFaMR>`:J:<:::::::>=?R:>:?J:@>Z::::::::kJ:vYxI:;Z::::::JywYB::::::::::::: yay=J:B:::::::::::::::::::jysy:>:<::::::::[B:<:N^:^:;jP@:EZ:^\\BE Z:V\\f:;j>D:Z:Vy<>jxM:<:[V:;Z>:_cZ<> :cTTPpsJSFKiR:=:kMN=:N @Qm\\Ij:<:IE:;@N@MnhIj:Jm^GOZ:n>N;yyyxy:<:;JBB:qQBv:>:sg:B:=J;Dlc`qsLq lp`h_:f?=J " 0 "" {MPLTEXT 1 0 26 "deq:=diff(y(x),x)=-x/y(x);" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 58 "Use the dsolve command to solve the differential equation." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "dsolve(deq,y(x),implicit);" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 24 "Thus a first integral is" }}{PARA 259 "" 0 "" {OLE 1 3597 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B> N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::fyyyyya:nYf::wyyyqy;:::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_mvJ: :::::::gjR<:T><::;Y[KXDIDNwW:A:;v:B:F:YLpfF>:::::::::J ?NZ;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 J:F;N;;j?>:S:UJ:nYvY:::::::::::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::j:b:=b:Kf FN[AvJbFtZAZ;>d:p^H[::::::::::j:D:;B:C:?R:=jZ:^:n_; >=E:]c:=Z:f:V[b<>fAn;n^@v;_jysy;B:KJ: :[l@[;;B:::::::JFNZ;J:N:;B:yayA:;B::::::f:D_=a=[;;B:::::::JF>:yay=J:B::::::nyyyYB:::::::::::::yay=J:B::::::: ::::::::::::jysy:>:<::::::::f<Z:N`D>n;>E=J:OM:_;Ms:=J:>IKG:JSjBTj:>:ALxgl; >:x;j>>:_Ktg];F:;B:<l;ZC:UK:^:>X=j>>:_c<;owAj:>Z:V?kp>JSjkMj:>:uKfO;:[Z: VY[j=B:;JXE:X;j;:::::::::::::::::::::::: ::::6:" }}{PARA 0 "" 0 "" {TEXT -1 44 "As this expression exists for a ll values of " }{TEXT 256 1 "x" }{TEXT -1 4 " and" }{TEXT 257 3 " y " }{TEXT -1 47 "it is also a conserved quantity for the system." }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 61 "Now enter the second set of equations and repeat the process." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "f1:=x;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "g1:=y(t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "deq1:=diff(y(x),x)=y(x)/x;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "dsolve(deq1,y(x),implicit); " }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 " Thus a first integral is " }}{PARA 260 "" 0 "" {OLE 1 3597 1 "[xm ]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyq y;:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::;Y[KX DIDNwW:A:;t:B:F:YLpfF>:::::::::J?NZ;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@Jyky;:::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::j:b:R<:TNC>:UTRcETcTX[US;SK]UW=EWMuUWm>j\\@A:::cYH_WV>Z::::::::jysy A:CB:F:;jysy?:;:?C;fb>JXA:=j>r:wP:D:Z:^:n_;>=E:]c:=Z:f:V[v=>r=n;n^@v;_jysy;B:KZ:>ryipyAj=fl=V; n>^;UTR:;JZ=j>D_mlVH[KR<:;B:::::::JFNZ;J:N:;B:yayA:;B::: :::f:D_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jy sy:>:<:::::::::::::::::::vYxI:;Z::::::::JBjx]:JBAj:B:;b:JSFk=NH=:kMN=:N@;e:gh:=J:VFKOBB:qAB:>l;>:DJ:DZJVdscRYEUDK< jPN:C:[Y:=:;JSd:KW:=J:fAkP?B:G;OjyyiyI:>:[Z:VY[j=J:^q:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;:::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::;Y[KXDIDNwW:A:; t:B:F:YLpfF>:::::::::J?NZ;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@Jyky;:::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::j:b:AFX?b[uc:KFF>gNGijgFyyYkjLN@>:::::::::::::::j:nyyYZDjysyQj;J:>R;B :bK_J:fc[_hb_dsSS;SK]UW=EWMuUWm>j\\@A:::cYH_WV>Z::::::::jysyA:CB:F:;jy sy?:;:qx:vn;JUAZ:F:MZ=VUN:e:;j<>:Mb:B:C:?R:=jZ:^:n_;>=E: ]c:=Z:f:V[j;>j=n;;JMTjAN;;JZyYjyYZHJj@j?JMJ@fc[;>Z:>^ :f:J;D_mlVH[KR<:;B:::::::JFNZ;J:N:;B:yayA:;B::::::f:=;jysy:>:<::::::wqy[:::::::::::::vYxI:;Z:::::::::: ::::::::::yay=J:B::::::::jDB:qi:;fyB:>L=:;b:F C:US:F[:>Z:N`D>b;^C=J:o:><C:Uk:^:>x;j>>: _c<^e;FZ:>Z:fAkpDjw;<:[V:B:D:c\\_;::f?;J:[Z:VY[j=J:^q