function dydt=forced(t,y)
dydt=[y(2);-9.0025*y(1)-1/10*y(2)+cos(3*t)];