{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 9 "restart; " }}{PARA 0 
"> " 0 "" {MPLTEXT 1 0 52 "n:='n':\nf:=x->100*Heaviside(x-1)-100*Heavi
side(x-3);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
19 "plot(f(x), x=0..4);" }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "c:=n->0.5*int(f(x)*sin(
n*Pi*x/4),x=0..4):\nc(n);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 64 "uapprox1:=(x,t)->sum(c(n)*sin(n*Pi*x/4)*exp(-4*n^2*t)
, n=1..30);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
30 "tvals:=seq(0.07/20*i,i=0..20):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 27 "plot(uapprox1(x,0),x=0..4);" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 39 "plot(uapprox1(x,1.00),x=0..4,-12..108);" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "toplot:=[seq(uapprox1(x,t),t
=tvals)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "plot(toplot,x=
0..4);" }{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 70 "wi
th(plots):densityplot(uapprox1(x,t),x=0..4,t=0..0.2,colorstyle=HUE);" 
}{TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "animate(ua
pprox1(x,t), x=0..4, t=0..0.2,frames=100);" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}{EXCHG }}{MARK "13" 0 }{VIEWOPTS 1 1 0 1 1 
1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }