{VERSION 5 0 "SUN SPARC SOLARIS" "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 "" 0 1 0 0 0 0 
0 0 0 0 0 0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "Utt = Uxx,     t>0,  0<x<L
" }}{PARA 0 "" 0 "" {TEXT -1 24 "U=0,     t>0,  x=0 & x=L" }}{PARA 0 "
" 0 "" {TEXT -1 26 "U=f,  Ut=g,     t=0, 0<x<L" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 7 "L:=2.0;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 7 "H:=0.1;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "e:=.1;" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 134 "Pluck:= (r,x) -> H*(x-(r*L
-e))/(e)*(Heaviside(x-(r*L-e)) -Heaviside(x-r*L)) + (H-H*(x-r*L)/e)*(H
eaviside(x-r*L)-Heaviside(x-(r*L+e)));" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 7 "r:=1/2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "
plot(Pluck(r,x),x=0..L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "
p:=1/2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "f:= x -> (Pluck(
p,x) + Pluck(1-p,x))/2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "
plot(f(x),x=0..L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "Impul
se:= (s,x) -> (Heaviside(x-s*L+e/2) - Heaviside(x-s*L-e/2))/e;" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "q:= 0.25;" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 47 "g:= x -> 0*(Impulse(q,x) + Impulse((1-q),x)
)/2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "plot(g(x),x=0..L);
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "lam:=(L,n) -> -(n*Pi/L)
^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "Z:=(L,n,x) -> sin(n*
Pi*x/L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "c:=(L,n) -> int
(f(x)*Z(L,n,x),x=0..L)/int(Z(L,n,x)^2,x=0..L);" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 72 "d:=(L,n) -> (L/(n*Pi))*int(g(x)*Z(L,n,x),x=0..
L)/int(Z(L,n,x)^2,x=0..L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
7 "c(L,5);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "d(L,5);" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "T:= (L,n,t) -> c(L,n)*cos(n*
Pi*t/L) + d(L,n)*sin(n*Pi*t/L);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 47 "U:= (m,L,x,t) -> sum(T(L,n,t)*Z(L,n,x),n=1..m);" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "Z(L,5,x);" }}}{EXCHG {PARA 0 "> " 0 
"" {MPLTEXT 1 0 9 "T(L,5,t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 11 "U(3,L,x,t);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(p
lots):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "animate(U(30,L,x,
t),x=0..L,t=0..8,frames=200,color=black);" }}}}{MARK "26 0 0" 58 }
{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }