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