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