xy'=x^3 sin^2(y) - 2 sin(y) cos(y)