function F=fun_ex(t,w)
global m c k A omega
F=[w(2); 1/m*(A*sin(omega*t)-c*w(2)-k*w(1))];