function xdot=mfd(t,x)
xdot=zeros(2,1);
global m D K R
F=R*t;
xdot(1)=x(2);
xdot(2)=(1/m).*(-D*x(2)-K*x(1)+F);