{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 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 14 "phi:=r->1-r^2;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "psi:=r->0;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "MAX:=5;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "k:=array(1..MAX);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 " k[1]:=fsolve(BesselJ(0,x),x=0.1..5);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 " for m from 2 by \+ 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 72 " k[m]:=f solve(BesselJ(0,x),x=k[m-1]..k[m-1]+4);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 22 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " " }}}{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 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "normfac:=array(1..MAX) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for m from 1 by 1 to M AX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 62 " normfac[m]:=2*Pi*int(r*(BesselJ(0,k[m]*r))^2 \+ ,r=0..1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('m');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "b:=array(1..MAX); B:=array(1..MAX); " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for m from 1 by 1 to MAX" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 59 " b[m]:=2*Pi*int(r*BesselJ(0,k[m]*r)*phi(r) ,r=0..1);" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{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 24 "for m from 1 by \+ 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 59 " B[m]:=2*Pi*int(r*BesselJ(0,k[m]*r)*psi(r) \+ ,r=0..1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{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 73 "u:=( t,r)->sum((b[m]/normfac[m])*BesselJ(0,k[m]*r)*cos(k[m]*t) ,m=1..MAX) " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 101 " +s um((B[m]/(normfac[m]*k[m]))*BesselJ(0,k[m]*r)*sin(k[m]*t) ,m=1..MAX) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "animate(u(t,r),r=0..1, t=0..6,frames=25);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "phi:=r->0;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "psi:=r->1 0-1000*r^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "MAX:=20;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "k:=array(1..MAX);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 " k[1]:=fsolve(BesselJ (0,x),x=0.1..5);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 " \+ for m from 2 by 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 " \+ do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 72 " \+ k[m]:=fsolve(BesselJ(0,x),x=k[m-1]..k[m-1]+4);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 22 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " " }}}{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 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "normfac:=array(1..MAX);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for m from 1 by 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 62 " no rmfac[m]:=2*Pi*int(r*(BesselJ(0,k[m]*r))^2 ,r=0..1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('m');" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "b:=array(1..MAX); B:=ar ray(1..MAX); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for m from 1 by 1 to MAX" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 54 " b[m]:=int(r*BesselJ(0,k[m]*r)*phi(r) ,r=0..1);" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{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 24 "for m from 1 by \+ 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 56 " B[m]:=int(r*BesselJ(0,k[m]*r)*psi(r) ,r=0. .0.1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 14 "unassign('m');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "u:=(t,r)->sum((b[m]/normfac[m])*BesselJ(0,k[m]*r)*cos (k[m]*t) ,m=1..MAX)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 76 " +sum((B[ m]/(normfac[m]*k[m]))*BesselJ(0,k[m]*r)*sin(k[m]*t) ,m=1..MAX);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "animate(u(t,r),r=0..1,t=0..2 ,frames=100);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "48 1 0" 76 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }