dsolve(diff(y(x),x) = 1/2*(-sin(2*y(x))+cos(2*y(x))*x^3+x^3)/x,y(x))