{VERSION 6 0 "IBM INTEL NT" "6.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 163 "# drumhead problem \n# u_\{tt\} = c^2(laplacian u), u=0 on clamped edge\n# position and v elocity given\n# Separation of variables leads to the following DE to \+ be solved:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "lambda:='lamb da':de:=diff(R(r),r,r)*r+diff(R(r),r)+lambda^2*r*R(r)=0;\ndsolve(de,R( r));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(DEtools):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 75 "lambda:=2:odeadvisor(de); # \+ Advice on what kind of equation is being solved" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 39 "# ?odeadvisor,types\n# ?dsolve\n# ?Bessel" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "# Plots of Bessel functions \+ of first and second kind" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "plot(BesselJ(0,2*r),r=0..40);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "BesselK(0,0);plot(BesselY(0,2*r),r=0..40);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "1 0 0" 17 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }