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