DSolve[-2*(x - y[x]^2)*Derivative[1][y][x]^3 + Derivative[1][y][x]*(1 + 4*y[x]*Derivative[1][y][x]) + (x + y[x]^2)*Derivative[2][y][x] == 0,y[x],x]