For teaching it is certainly nice to have numpy.polynomial.polynomial.polyfit
providing modern (vs. traditional) parameter order, but
- it is rather buried
- np.polyfit uses traditional order and has the same name
I recall there was some controversy (?) over all of this,
but might it not be appropriate to have a keyword argument to
both specifying whether the parameter order is to be modern
or traditional (in both polyfits and polyvals)?