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