rugo.mws
>
deq:=m*diff(x(t),t,t)=-d*x(t)-k*diff(x(t),t)+A*sin(omega*t);
>
with(DEtools):
>
m:=1; d:=0.1; k:=0.1; A:=0.1; omega:=1;
>
DEplot(deq,x(t),t=0..100,[[x(0)=1,D(x)(0)=0]],stepsize=0.01);
>
dsolve(deq,x(t));
>