{VERSION 6 0 "IBM INTEL NT" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 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 2 0 2 0 2 2 0 1 }{PSTYLE "Headi ng 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 2 0 2 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 56 " \304\350\364\364\345\360\345\355\366\350 \340\353\374\355\373\345 \363\360\340\342\355\345\355\350\377 \361 \+ \360\340\347\360\373\342\355\373\354\350 \357\360\340\342\373\354\350 \+ \367\340\361\362\377\354\350" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "eq:=diff(y(t),t)=-y(t)*Heaviside(t-1);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 54 "plot('Heaviside'(x-1),x=-2..2,thickness=2,axes=BOXE D);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "ini:=y(0)=3;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "dsolve(\{eq,ini\},\{y(t)\}); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "subs(%,y(t));" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "f:=unapply(%,t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "plot(f,0..4);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "sol1:=dsolve(\{eq,ini\},\{y(t)\},type=numeric );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "plots[odeplot](sol1,[ t,y(t)],0..4);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 78 "#======== ===================================================================== " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "eq:=diff(y(t),t)=-y(t)*Dirac(t-1); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "ini:=y(0)=3;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "dsolve(\{eq,ini\},\{y(t)\}); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "f:=unapply(subs(%,y(t)) ,t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "plot(f,0..4,discont =true,view=[0..4,0..4]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "sol2:=dsolve(\{eq,ini\},\{y(t)\},type=numeric,method=classical[rk2]); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "with(plots,odeplot);" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "odeplot(sol2,[t,y(t)],0..4 ,view=[0..4,0..4]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "#=== ==============================================================" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "f:=x->piecewise(1<=x and x<2 ,1,0);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "f(x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "plot(f,0..3,axes=BOXED);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "eq:=diff(y(t),t)=1-y(t)*f(t) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "ini:=y(0)=3;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "dsolve(\{eq,ini\},\{y(t)\}); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "sol3:=dsolve(\{eq,ini\} ,\{y(t)\},type=numeric,method=classical[rk4]);" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 20 "with(plots,odeplot):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "odeplot(sol3,[t,y(t)],0..4);" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 49 " \304\350\364\364\345\360\345\355\366\350\340 \353\374\355\373\345 \363\360\340\342\355\345\355\350\377 \342 \367 \340\361\362\355\373\365 \357\360\356\350\347\342\356\344\355\373\365 " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "wave:=diff(u(x,t),t,t)-c ^2*diff(u(x,t),x,x)=0;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "s ol:=pdsolve(wave,u(x,t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "f1:=xi->exp(-xi^2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 " f2:=xi->piecewise(-1/2 " 0 " " {MPLTEXT 1 0 28 "subs(_F1=f1,_F2=f2,c=1,sol);" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 15 "subs(%,u(x,t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "f:=unapply(%,x,t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "plot3d(f,-8..8,0..5,grid=[60,60],axes=BOXED,style=pat chcontour);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "#----------- -----------------------------------------------------" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "heat:=diff(u(x,t),t)-k*diff(u(x,t), x,x)=0;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "pdsolve(heat,u(x ,t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "eq:=subs(u(x,t)=X( x)*T(t),heat);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "eq/X(x)/T (t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "expand(%);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "sep:=%+(k*diff(X(x),x,x)/X(x )=k*diff(X(x),x,x)/X(x));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "lhs(sep)=C;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "T_sol:=d solve(%,T(t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "rhs(sep)= C;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "X_sol:=dsolve(%,X(x), explicit=true);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "subs(X_s ol,T_sol,X(x)*T(t));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "sol :=simplify(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "subs(C=-k ,k=1,_C1=1,_C2=1,sol);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "S :=convert(%,trig);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "plot3 d(S,x=-5..5,t=0..5,axes=BOXED);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "subs(u(x,t)=sol,heat);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "simplify(\{%\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "#==================================================================== " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "with(PDEtools):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "pde:=diff(u(x,y),x)+cos(2*x)*diff(u(x,y),y) =-sin(y);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "ini:=[0,s,1+s^ 2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "PDEplot(pde,u(x,y),i ni,s=-2..2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "PDEplot(pde ,u(x,y),ini,s=-2..2,basechar=only);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "PDEplot(pde,u(x,y),ini,s=-2..2,basechar=true);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 141 "PDEplot(pde,u(x,y),ini,s=-2 ..2,\n basechar=true,initcolor=white,\n style=patchconto ur,contours=20,\n orientation=[-43,45]);" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 17 " \+ \320\340\341\356\362\340 \361 \364\340\351\353\340\354\350" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 72 "#================================== =====================================" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "#MyDir:=\"D:\\\\ASN\\\\2005\\\\Day_5\\\\\":" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "MyDir:=\"D:/ASN/2005/Day_5/\":" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "L:=readdata(cat(MyDir,\"data.txt\"),3);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "map(u->[u[1],u[3]],L);" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "plot(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "L[5,2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "L[1..nops(L),2];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(stats):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "de scribe[mean](%%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "A:=con vert(L,array);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "transpose (A)&*A;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "P:=evalm(%);" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "linalg[eigenvalues](P);" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "linalg[eigenvectors](P);" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "map(z->linalg[norm](op(z[3]),2),[% ]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "L:=readdata(cat(MyDi r,\"integers.txt\"),integer,2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "#------------------------------------------------------------- -------------------" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "read cat(MyDir,\"maple_text.txt\");" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "read cat(MyDir, \"maple_text.txt\");" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "L:= [3,3.1415,-65,0];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "interf ace(echo=1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "writedata(` terminal`,L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "A:=array([ [1,2,3],[-1.45,0,3/2]]);" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "writedata(`terminal`,A);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 15 "L:=[Pi,exp(9)];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "Lf:=evalf(L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "writedata(`terminal`,Lf);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "L:=[3.7,3.1,-2.2,-1.9];" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "writedata(`terminal`,L,integer);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 13 "map(trunc,L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "writedata(`terminal`,map(round,L),integer);" }}} {EXCHG }{EXCHG }{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "writ edata(cat(MyDir,\"append.txt\"),L,integer);" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "fclose(cat(MyDir,\"append.txt\"));" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "readdata(cat(MyDir,\"append. txt\"),1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "writedata[APP END](cat(MyDir,\"append.txt\"),map(round,L),integer);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "fclose(cat(MyDir,\"append.txt\")); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "readdata(cat(MyDir,\"ap pend.txt\"),1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "#======= ============================================================" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "qbinomial:=(n,k)->product(1- q^i,i=n-k+1..n)/product(1-q^i,i=1..k);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "expr:=qbinomial(10,4);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "nexpr:=normal(expr);" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "save qbinomial,expr,nexpr,cat(MyDir,\"qbinom.m\" );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "expr;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "MyDir:=\"D:/ASN/2005/Day_5/\":" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 27 "read cat(MyDir,\"qbinom.m\");" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "expr;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 69 " \312\356\355\342\345\360\362\340\366\350 \377 \342\373\360\340\346\345\355\350\351 \342 \342\373\360\340\346 \345\355\350\377\n \355\340 \344\360\363\343\350\365 \377\347\373\352 \340\365 \357\360\356\343\360\340\354\354\350\360\356\342\340\355\350 \377" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "#=================== =====================================" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "latex(a/b);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "formula:=Limit(int(f(x),x=-n..n),n=infinity);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "latex(formula);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(codegen):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "f:=convert(taylor(sin(x),x),polynom);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "of:=optimize(f);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "cost(f);" }}}{EXCHG }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "cost(of);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "C([of]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "makeproc (f,x);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "makeproc([of],x); " }}}{EXCHG }}{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG } {EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }{EXCHG }}{MARK "68" 0 } {VIEWOPTS 1 1 1 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }