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