function rhs=rhs(t,y,dummy,A)

rhs=[y(2)
     y(3)
     -(y(1)^2)*y(2)-cos(t)*y(1)+A*(sin(t)^2)];