Numpy read feather. Fit a polynomial p(x) = p[0] * x**deg + .