dsolve(diff(y(x),x) = 1/4*(-4*cos(x)*x+4*x^2*sin(x)+4*x+4+4*y(x)^2+8*y(x)*cos(x)*x-8*x*y(x)+2*x^2*cos(2*x)+6*x^2-8*x^2*cos(x)+4*y(x)^3+12*y(x)^2*cos(x)*x-12*x*y(x)^2+6*y(x)*x^2*cos(2*x)+18*x^2*y(x)-24*y(x)*cos(x)*x^2+x^3*cos(3*x)+15*x^3*cos(x)-6*x^3*cos(2*x)-10*x^3)/x,y(x))