rugo.mws

>    deq:=m*diff(x(t),t,t)=-d*x(t)-k*diff(x(t),t)+A*sin(omega*t);

deq := m*diff(x(t),`$`(t,2)) = -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;

m := 1

d := .1

k := .1

A := .1

omega := 1

>    DEplot(deq,x(t),t=0..100,[[x(0)=1,D(x)(0)=0]],stepsize=0.01);

[Maple Plot]

>    dsolve(deq,x(t));

x(t) = exp(-1/2*(k-(k^2-4*d*m)^(1/2))/m*t)*_C2+exp(-1/2*(k+(k^2-4*d*m)^(1/2))/m*t)*_C1+A*((d-omega^2*m)*sin(omega*t)-omega*cos(omega*t)*k)/(omega^4*m^2+(k^2-2*d*m)*omega^2+d^2)

>