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