{VERSION 5 0 "IBM INTEL NT" "5.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "with(VectorCalculus);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "Digits:=5;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "assume(x,real);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "A:=matr ix([[0,0,1,0],[0,0,0,1],[-3,2,0,0],[2,-3,0,0]]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "dog:=jordan(A,'P');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "print(P);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "z1:=;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "z2:=<0,exp(dog[2,2]*x),0,0>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "z3:=<0,0,exp(dog[3,3]*x),0>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "z4:=<0,0,0,exp(dog[4,4]*x)>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "y1:=;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "y2:=;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "y3:=;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "y4:=;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "ny1:=simplify(convert(y1,trig));" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sol1:=<>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sol2:=<>;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "ny3:=simplify(convert(y3,tri g));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sol3:=<>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sol4:=<> ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "solmat:=augment(sol1,s ol2,sol3,sol4);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "Inmat:=a rray(1..4,1..4);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for i f rom 1 by 1 to 4 do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 29 " for j from 1 by 1 to 4 do " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 " Inmat[i,j] :=eval(solmat[i,j],x=0);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 6 " od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "print(Inmat);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "ics:=<2,4,0,0>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "< ics>;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "c:=linsolve(Inmat, ics);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "c[1];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "yeah:=multiply(solmat,c);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "soly1:=yeah[1];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "soly2:=yeah[2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "dery1:=yeah[3];" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 15 "dery2:=yeah[4];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot(soly1,x=0..20*Pi);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot(soly2,x=0..20*Pi);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "plot(\{soly1-10,soly2-20\},x=0..20*Pi);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "ode1:=diff(w1(x),x$2)=-3*w1( x)+2*w2(x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "ode2:=diff(w 2(x),x$2)=2*w1(x)-3*w2(x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 80 "result:=dsolve(\{ode 1,ode2,w1(0)=2,w2(0)=4,D(w1)(0)=0,D(w2)(0)=0\},\{w1(x),w2(x)\});" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "52" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }