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