bessel_transform
Graham Dennis
Solve the wave equation on a sphere of radius 1 utilising spherical symmetry by using the Spherical Bessel function transform.
const real T = 10.0; const real mass = 1e-3; const real length = 1.0; const real mu = mass/length;
t
u uDot
u = exp(-100.0*(x-0.25)*(x-0.25)); uDot = 0.0;
100 100
L
L = -T*kx*kx/mu;
main
du_dt = uDot; duDot_dt = L[u];
amp
main
amp = u.Re();
amp
main
amp = u.Re();