{VERSION 6 0 "IBM INTEL NT" "6.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 61 "with(PDEtools): diff usion:=diff(u(t,x),t)=2*diff(u(t,x),x,x);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 52 "pdsolve(diffusion, u(t,x), HINT=X(x)*T(t), 'build') ;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 17 "Neumann boundary " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "IBC:=\{u(0,x)=exp(-0.2*x), D[2](u)(t,0)=0 , D[2](u)(t,1)=0\};" }{TEXT -1 0 "" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 54 "pds:=pdsolve(diffusion, IBC, numeric, spaces tep=1/40);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "pds:-animate( t=0.2,frames=100);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 14 "Robin Bound ary" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "IBC1:=\{u(0,x)=exp(-0.2*x), \+ D[2](u)(t,0)=-2, D[2](u)(t,1)+u(t,1)=0\};" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 56 "pds1:=pdsolve(diffusion, IBC1, numeric, spacestep=1 /40);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "pds1:-animate(t=0. 2,frames=500);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 17 "Periodic Bounda ry" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "IBC2:=\{u(0,x)=exp(-0.2*x), D [2](u)(t,0)=D[2](u)(t,1), u(t,0)=u(t,1)\};" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 57 "pds2:=pdsolve(diffusion, IBC2, numeric, spacestep=1 /160);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "pds2:-animate(t=0 .2,frames=500);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "Smoothering ef fect of diffusion" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 169 "p1:=pds2:-plo t(t=0,numpoints=50, color=red):\np2:=pds2:-plot(t=0.01,numpoints=50,co lor=blue):\np3:=pds2:-plot(t=0.02,numpoints=50,color=green):\nplots[di splay](\{p1,p2,p3\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "IB C3:=\{u(0,x)=sin(20*x), D[2](u)(t,0)=D[2](u)(t,1), u(t,0)=u(t,1)\};" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "pds3:=pdsolve(diffusion, I BC3, numeric, spacestep=1/160);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 169 "p4:=pds3:-plot(t=0,numpoints=50, color=red):\np5:=pds3:-plot( t=0.01,numpoints=50,color=blue):\np6:=pds3:-plot(t=0.02,numpoints=50,c olor=green):\nplots[display](\{p4,p5,p6\});" }}}{EXCHG {PARA 0 "" 0 " " {TEXT -1 36 "Anti-smoothing of backward diffusion" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 171 "p4:=pds3:-plot(t=0,numpoints=50, color=red):\np5: =pds3:-plot(t=-0.01,numpoints=50,color=blue):\np6:=pds3:-plot(t=-0.02, numpoints=50,color=green):\nplots[display](\{p4,p5,p6\});" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "15 1 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }