vibstring3
Joe Hope
Vibrating string
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;
mom
moment
integrated_u main
moment = mod2(u);
50 50
sparemoment
integrated_u main
sparemoment = mod2(u);
L
L = -T*kx*kx/mu;
main
du_dt = uDot; duDot_dt = L[u /**/];
integrated_u
filter1
dmom_dt = moment;
filter1 main
u = 1.0 * u;
sparemomentagain
integrated_u main
sparemomentagain = mod2(u);
L
L = -T*kx*kx/mu;
amp ke
main filter1
amp = mod2(u + moment); ke = mod2(L[u]);
momR momentR
integrated_u filter1
momR = mom.Re(); momentR = moment.Re();