Interesting, that question would probably have gotten a different response on scicomp, it is a pity we are not attracting more questions there!

I know there are two polyfit modules in numpy, one in numpy.polyfit, the other in numpy.polynomial, the functionality you are suggesting seems to fit in either.

What parameters/functionality are you considering adding?

A