{VERSION 5 0 "IBM INTEL LINUX" "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 14 0 
0 0 0 2 2 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R
3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 14 0 0 0 0 2 1 2 0 
0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 2" -1 
257 1 {CSTYLE "" -1 -1 "Courier" 1 14 0 0 0 0 2 2 2 0 0 0 0 0 0 0 }0 
0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots);" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "phi:=x->x*(Pi-x)*(sin(12*x))
^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 17 " b:=array(1..30);" }}}{EXCHG {PARA 0 "> " 0 
"" {MPLTEXT 1 0 23 "for n from 1 by 1 to 30" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 5 "   do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "      b[n]
:=(2/Pi)*int(phi(x)*sin(n*x),x=0..Pi);" }}{PARA 0 "> " 0 "" {MPLTEXT 
1 0 5 "  od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n
');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 49 "u:=(t,x)->sum(b[k]*sin(k*x)*exp(-k^2*t),k=1.
.30);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "plot3d(u(t,x),x=0.
.Pi,t=0..0.02,grid=[100,100]);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 44 "animate(u(t,x),x=0..Pi,t=0..0.5,frames=400);" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 8 "phi:=10;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 " b:=array(
1..30);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "for n from 1 by \+
1 to 30" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "   do" }}{PARA 0 "> " 0 "
" {MPLTEXT 1 0 88 "      b[n]:=(2/Pi)*int(phi(x)*sin(n*x),x=1.5..1.7)+
(2/Pi)*int(phi(x)*sin(n*x),x=0.8..1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 5 "  od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n')
;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 49 "u:=(t,x)->sum(b[k]*sin(k*x)*exp(-k^2*t),k=1..3
0);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "plot3d(u(t,x),x=0..P
i,t=0..2);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "animate(u(t,x
),x=0..Pi,t=0..2,frames=400);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 0 "" }}}}{MARK "19 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }
{PAGENUMBERS 0 1 2 33 1 1 }