{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 7 0 0 1 0 0 0 0 0 0 0 } {CSTYLE "" -1 257 "" 0 1 0 0 2 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 "Co urier" 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 "Maple Plot" 0 13 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 256 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 257 1 {CSTYLE "" -1 -1 " " 0 1 56 1 204 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 182 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 0 0 0 1 0 0 0 0 0 0 0 }0 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 256 32 "Global Stability Investig ation 2" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 63 "You are going to draw the phase portrait of a nonlinear system " } }{PARA 260 "" 0 "" {OLE 1 4096 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N :F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::fyyyyya:nYf::G:jy;::::::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::JcvGYMt>^:fBWMtNHm=;:::: :::n:;JZC:bKB::>BmwkoiJwW:A:LZ:>Z:j:vCSmlJ::::::::::OJ;@ jyyyyyy;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN:ry:>:<::::::=J:^ J>:F:AlqfG[maNFO=;::::::::_J;vyyuy:>:<::::::AZ:^E>:nYN:J:[B::D]Hb:R:[DZATRWuREEUv:I J:>:::::FZ:nyyYZDjysyQj;J:>R;B:bK_J:fc[_hb_ds?h_?^?GhoGfnGgioGM:ETV::: ^xsNpk?:;B::::::::vYxy;J<<:=J:vYxY;J:JMQ:sS;Nk:Mb: >Z:^:NZ;F:E:=b:yyyyI:E:M::Gc;YJCv YxY:JZyyRyYjMjrLj?JMJ@fc[;>:@jJ;D_mlVH[KRJ:<:::::::>=?R:>: ?J:@>Z::::::::kJ:vYxI:;Z::::::JywYB: ::::::::::::yay=J:B:::::::::::::::::::jysy:>:<::::::::Fjx]:JBAZ:b :DJ:f<l;F:>ZGC:Uk:^:>x;j>>:_kKAJ`F:I:_kELJfF:FIO?N@Fn=F:FI><X=j>>:_cZ:V@[m>JSVg=nf>F:V@N[:B:Gc;N;yyyxy: J:><:sg:B:=b:?bBaTXaEWEUUtP:C:[q:VZ:B:;:::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::5:" }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 44 "and the level curves of a Lyapu nov function " }}{PARA 261 "" 0 "" {OLE 1 4096 1 "[xm]Br=WfoRrB:::wk;n yyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::fyyyyya:nYf::G:jy;::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::JcvGYMt> ^:fBWMtNHm=;:::::::n:;JZC:bKB::>BmwkoiJwW:A:;F;>Z:j:vCSm lJ::::::::::OJ;@jyyyyyy;jysy;Z:::::::^<>:F::]KRnC=MtFGgml>:;::::::JGN: ry:>:<::::::=J:^J>:F:AlqfG[maNFO=;::::::::_J;vyyuy:>:<::::::AZ:^E>:nYN :J:W:YJ:>\\:B:]:wAyA::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::j:b:=b:KfFN[d:Xb[og]gffk=v:;J:::::j:Z::::::::jysyA:CB:F:;jysy? :;:[u<D:;B:C:?R:=jZ:^:n_;>=E :]c:=Z:f:V[f<>fDn;n^@v;_jysy;B:KJ:>ryapyA^D^\\:; E:EZ:F;NZ:vCS=[LsfFaMR>@>Z::::::::kJ;@:;J;>Z:vYxY:>Z::::::jD_=a=[;;B:::::::JF>:yay=J:B::::::nYyA<::::::::::::jysy:>:<:::::::::: :::::::::vYxI:;Z:::::::::eZ:Vy<>jx]:JBAj:J:DZJ^dcgg_WhZnc_whZNdigg[oGk iZC:US:F[:>Z:N`D>n;>\\:FZ:>Z:>:::::::::::nDOq>JSj_Aj:JvNOM:_;Us:=:mMO; _;;b;=J:>I:_;uH<=:m=e:qAB:>l;>Z:B:DJ:<;::FEC:Uk:^:>x;j>>:_;Ef:=J:F [Z:VY;RyB:>l;Z:b:aTXDpql`W^:f?=JNoG;=J:YC:UK:^:>X=j>>:_Ktgr@F:HM:d:?W<=:<:OB: " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 13 "with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warni ng, new definition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning, n ew definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 " with(DEtools):" }}{PARA 7 "" 1 "" {TEXT -1 35 "Warning, new definition for adjoint" }}}{EXCHG {PARA 259 "" 0 "" {TEXT -1 56 "To plot the pha se plane diagram for the nonlinear system" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 12 "with(plots):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "A:=diff(x(t),t)=-x^3;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "B:=diff(y(t),t)=-y^3;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 156 "p1:=DEplot(\{A,B\},\{x(t),y(t)\},t =0..50,[[x(0)=1,y(0)=1],[x(0)=-.25,y(0)=1],[x(0)=-1,y(0)=-.4],[x(0)=.3 ,y(0)=-1]],x=-1..1,y=-1..1,stepsize=.1,linecolour=blue):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "p1;" }}{PARA 13 "" 1 "" {TEXT -1 0 " " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 13 " Verify that " }}{PARA 256 " " 0 "" {OLE 1 3608 1 "[xm]Br=WfoRrB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]:: yyyyyy:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: :::::::fyyyyya:nYf::wyyyqy;::::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::;lVYVwENYAj;JZI:<:=ja^GE=;:::::::::N;?R:yyyyyyA:ya yA:<::::::JDB:;:=:jR>@Wlj^HMMufF;J:::::::N=?:xI:;Z::::::j:>:C=;:=j[vGU MrvC?MoJ::::::::JCN:yyyxI:;Z::::::j;B:s<;:wA?Z::C:J=j=B:K:OJ:vYxy?J@>: UJ:n;nYvY::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::::::::::::=ZAFXAr=V[KFZ<>klJ?DJ=iLCHrv HjA:Z:>:jysy::;J:<::::::::::=Jyyy;d:yayY;A:;JZC:bK_J:fc[_hb_ds?h_?^?Gh oGfnGgioGM:ETV:::^xsNpkK:<::::::::yayY:^Z:j:>:yayQ:>:F]?jwD:oX:F:MZ=fW N:eZ:>:EJ:F[DEZ:F;NZ:vCS=[LsfFaMR>@>Z ::::::::kJ;@:;J;>Z:vYxY:>Z::::::j@>Z::::::::kJ:vYx I:;Z::::::JywYB:::::::::::::yay=J:B:::::::::::::::::::jysy:>:<:::::::: f<Z::::::::::>IKG:M:_;Ms:=J:FIKoDjw?^y]:JBAJ:x;j>>:_KtG\\;FZ:>Z:B:[m>ZDJqIj:Z:JBB:qQ:[:JBA:DZaTXDpql`M^Z:>:UK :^:>X=j>>:_c<;W:YY:=J:UOZ:n>N;yyyxy:J:><:sg:B:=J;Dlc`qsLq lp`h_:f?=JOM:_k>Ij=Ij:B:k=[Z:VY;><:[V:b:DZJ:Y]N^:f??Jfjysy=:;JBB:qQBv:>:s?<::::5:" }{MPLTEXT 1 0 0 " " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 38 "is a Lyapunov function for th e system." }}{PARA 0 "" 0 "" {TEXT -1 35 "If it is a Lyapunov functio n then " }}{PARA 257 "" 0 "" {OLE 1 3612 1 "[xm]Br=WfoRrB:::wk;nyyI;G: ;:j::>:B>N:F:nyyyyy]::yyyyyy:::::::::::::::::::::::::::::::::::::::::: :::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::NDYmq^H;C:E Lq^H_mvJ::::::::gjR<:T><::;lVYVwENYAj;JZI:<:=ja^GE=;:: :::::::N;?R:yyyyyyA:yayA:<::::::JDB:;:=:jR>@Wlj^HMMufF;J:::::::N=?:xI: ;Z::::::j:>:K=;:=j[vGUMrvC?MoJ::::::::JCN:yyyxI:;Z::::::j;B:s<;:wA?Z:: C:J=j=B:K:M:OJ:V;;J@>:UJ:n;nYvY::::::::::::::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::=ZZ;N Z@ppLVjPS:?:LA?b[Ag:KFFFZ:VZJ?j`Fe_WHt]_bFtZA:[DZGcM>Z::::::kRTq<;:::: =B:wyyAb:EJ:F[Z:^:n_;>=E:]c:=Z:f:V[v=>r>n;;JMTjAN;;JZyyTyY:MjDIj?JMJ@fc[;>:;H:MB:@[C:>Z::::::::kB:? R:>:?J:@>Z::::::::kJ:vYxI:;Z::::::Jy wYB:::::::::::::yay=J:B:::::::::::::::::::jysy:>:<::::::::[B:<:N>C:US: f::Mb:B:[J:IjLG:nP:_KjDJPEj:>:AOfg:;b:;Jd``_WhZ nc_whZNdigg[oGKQ>:_kov\\;F:;Je^VM:_k>gf;F:^HSmDjw;<:[V:;Z :b::::SBBC:Uk:^:>X?B:MJ:N@;mdIj:@;P?B:G;OjyyiyI:>:[Z:VY[j=J:^q:UK;^:>X=j;::::::::::::::::::::::::::5:" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "V :=x^2+y^2;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "Vx:=diff(V,x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "Vy:=d iff(V,y);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 " f:=-x^3;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "g :=-y^3;" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "V t:=Vx*f+Vy*g;" }}{PARA 11 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 5 "Thus " }}{PARA 258 "" 0 "" {OLE 1 3612 1 "[xm]Br=WfoR rB:::wk;nyyI;G:;:j::>:B>N:F:nyyyyy]::yyyyyy::::::::::::::::::::::::::: ::::::::::::::::::::::::::::::::::::::::::::fyyyyya:nYf::wyyyqy;:::::: :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: ::::NDYmq^H;C:ELq^H_mvJ::::::::gjR<:T><::;Ub;xkoiJwW: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:::::::::::::::::::::::: :::::::::::::::::::::::::::::::::::::::::::::::::::=ZAfO@J;ZnK;TrKEJ>===B:ABN;fCmTSa]xRCD]Hr;Jb ;B::::::::::FZ:nyyYZDjysyQj;:;`:Z@O<;j`@Pt\\Pd`QrPPJPnrPqjLqnxPqF; fbk;::JtaMSAA;B::::::::vYxy;J<JGD:=j>r:UR:VZ:f<D:;B:CZ:NZ;F:E:=b:yyyyI:E:M: :Gc;YJCvYxY:>:;xyewy;F;ev:Q:G;Sj`@@J:>r:F[:B:oi:NZ:vCS=[LsfFaMR>`:J :<:::::::>=?R:>:?J:@>Z::::::::kJ:vYx I:;Z::::::JywYB:::::::::::::yay=J:B:::::::::::::::::::jysy:>:<:::::::: [B:<:N>C:US:f::Mb:B:[J:IjLG:nP:_KjDJPEj:>:AOfg: ;b:;Jd``_WhZnc_whZNdigg[oGSBBC:UK;^:>x;F:MJ:N@QmE@j:>:gLlG;N@Mnl@j:Jt^ Ge:qAB:>l;>:X?B:MJ:N@;mdIj:@;P?B:G;OjyyiyI:>:[Z:VY [j=J:^q:UK;^:>X=j;::::::::::::::::::::::::::5 :" }}{PARA 0 "" 0 "" {TEXT -1 69 "for all x and y so V is a Lyapunov f unction for the nonlinear system." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 257 29 "To plot the level curves of V" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "p2:=implicitplot(\{x^2+y^2=1 /4,x^2+y^2=1,x^2+y^2=.5\},x=-1..1,y=-1..1):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "p2;" }}{PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 122 "The two diagrams are obviously not the s ame. to examine them more closely look at the two graphs plotted on th e same axes." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "display(\{p 1,p2\});" }}{PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 109 "Notice that the solution curves for the nonlinear system all cross the level curves of the Lyapunov function." }}}}{MARK "24 0 0" 77 }{VIEWOPTS 1 1 0 1 1 1803 }