DSolve[Derivative[1][y][x] == (1 + x*y[x])^3/x^5,y[x],x]