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