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