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