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