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