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