{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 40 "phi:=(x,y)->(x-1.5)*
(2-x)*(y-1.5)*(2-y);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "plo
t3d(phi(x,y),x=0..Pi,y=0..Pi);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 8 "MAX:=20;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "b:=array
(1..MAX,1..MAX);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "for n f
rom 1 by 1 to 20" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "  do" }}{PARA 0 
"> " 0 "" {MPLTEXT 1 0 28 "     for m from 1 by 1 to 20" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 10 "        do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 
86 "            b[m,n]:=(4/Pi^2)*int(int(phi(x,y)*sin(m*x)*sin(n*y) ,x
=1.5..2) ,y=1.5..2);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "        od;
" }}{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 14 "unassign('m');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "u:=(t,x,y)->sum(
sum((b[m,n]*sin(m*x) *sin(n*y)*exp(-(m^2+n^2)*t)  ),n=1..20),m=1..20);
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots);" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "animate3d(u(t,x,y),x=0..Pi,y
=0..Pi,t=0..2,frames=20);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 
"" }}}}{MARK "9 0 0" 12 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }
{PAGENUMBERS 0 1 2 33 1 1 }