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