function xdot=cascade(t,x)
xdot=[-3*x(1)+1+cos(t),3*x(1)-x(2)];