function dydt=oscill(t,y)
dydt=[y(2); -3*y(1)-4*y(2)+sin(2.2*t)];