syms q m L F n pi v t w
q=dsolve('D2q+w^2*q=2/(m*L)*F*sin(n*pi*v*t/L)','Dq(0)=0','q(0)=0','t'); F=subs(F);
L=subs(L); q=subs(q); m=subs(m); v=subs(v);
w=subs(w); pi=subs(pi); n=subs(n); t=subs(t);
F=10e4; L=32; EI=3e10; m=3e4; v=5; x=16; n=2;
w=(n*pi)^2*sqrt(EI/(m*L^4));
a=sin(n*pi*x/L);
for t=0:0.5:6.4
q
end