vibstring_dct
Graham Dennis
Vibrating string using a Discrete Cosine transform (i.e. even boundary conditions on both sides)
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.5)*(x-0.5)); uDot = 0.0;
50 50
moment
main
moment = mod2(u);
L
L = -T*kx*kx/mu;
main
moment_creator
du_dt = uDot; duDot_dt = L[u]; moment;
amp
main
amp = u;
amp
main
amp = u;