dsolve(diff(y(x),x) = 1/2*(-sin(y(x)/x)*y(x)*x-y(x)*sin(y(x)/x)+y(x)*sin(3/2*y(x)/x)*cos(1/2*y(x)/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)*x+y(x)*cos(1/2*y(x)/x)*sin(1/2*y(x)/x)+2*sin(y(x)/x)*x^4*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/(1+x),y(x))