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