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