{VERSION 5 0 "SUN SPARC SOLARIS" "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 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 20 "restart;with(plots):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 853 "recreate3d:=proc(kap,ta,a,b,c,d,e,f,g,h)\nloc al sys,p,alph1,alph2,alph3,T1,T2,T3,N1,N2,N3,B1,B2,B3,pl;\nsys:=diff(a lph1(s),s)=T1(s),diff(alph2(s),s)=T2(s),diff(alph3(s),s)=T3(s),\ndiff( T1(s),s)=kap(s)*N1(s),diff(T2(s),s)=kap(s)*N2(s),diff(T3(s),s)=kap(s)* N3(s),\ndiff(N1(s),s)=-kap(s)*T1(s)+ta(s)*B1(s),\ndiff(N2(s),s)=-kap(s )*T2(s)+ta(s)*B2(s),\ndiff(N3(s),s)=-kap(s)*T3(s)+ta(s)*B3(s),\ndiff(B 1(s),s)=-ta(s)*N1(s),diff(B2(s),s)=-ta(s)*N2(s),diff(B3(s),s)=-ta(s)*N 3(s);\nprint(sys);\np:=dsolve(\{sys,alph1(0)=0,alph2(0)=0,alph3(0)=0,T 1(0)=1,T2(0)=0,T3(0)=0,N1(0)=0,N2(0)=1,N3(0)=0,B1(0)=0,B2(0)=0,B3(0)=1 \},\{alph1(s),alph2(s),alph3(s),T1(s),T2(s),T3(s),N1(s),N2(s),N3(s),B1 (s),B2(s),B3(s)\},type=numeric):\npl:=odeplot(p,[alph1(s),alph2(s),alp h3(s)],a..b,numpoints=200,thickness=1,axes=normal,color=red):\ndisplay (pl,scaling=constrained,view=[c..d,e..f,g..h]);\nend:\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "recreate3d(1/5,1/10,-100,100,-10,10 ,-10,10,-10,10);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "kap3d1:=t->t;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "tau3d1:=t->t/10;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "recreate3d(kap3d1,tau3d1,0,10,0,2,0,2,0,1 );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "recreate3d(kap3d1,tau 3d1,-10,10,-2,2,-2,2,-2,2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "kap3d2:=t->t^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "tau3 d2:=t->t^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "recreate3d(k ap3d2,tau3d2,-10,10,-2,2,-2,2,-2,2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "kap3d3: =t->t*sin(t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "tau3d3:=t- >t;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "recreate3d(kap3d3,2, -10,10,-5,5,-2,2,-2,2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 " recreate3d(2,kap3d3,-10,10,-2,2,0,5,-2,2);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 18 "kap3d4:=t->sin(t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "recreate3d(kap3d4,2,-10,10,-2,2,-2,2,-2,2);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "\nrecreate3d(2,kap3d4,-10,10 ,-2,2,-2,2,-2,2);\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} }{MARK "21" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }