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