{VERSION 2 3 "IBM INTEL NT" "2.3" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 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 20 "restart;with(plots): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "Lx:=2; Ly:=1;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "g:=sin(Pi*x/Lx)*sin(Pi*y/Ly) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "plot3d(g,x=0..Lx,y=0.. Ly);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "c:=100; # wave spee d on membrane" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "omega:=c*s qrt((Pi/Lx)^2+(Pi/Ly)^2); # lowest mode angular freq" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "Nframes:=10; tfac:=2*Pi/(omega*Nframes); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "for i from 0 to Nframes -1 do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "plt.(i):=plot3d(sin(omega* i*tfac)*g,x=0..Lx,y=0..Ly);od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "display([plt.(0..Nframes-1)],insequence=true);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "9 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }