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