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