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