Fit Polynomial through a set of 2D Points with Rhino/Python

Hi everyone,
is there a possibility to calculate a Polynominal through a set of 2D Points with Python or Rhino Script?
I am used to Python programming and numpy would easily do the job unfortunately it is not available for Rhino 6 and Rhino 5 64 Bit.
Thanks a lot!

Have a look at Math.NET Numerics, Giulio wrote a great guide on how to implement in GHPython. Looks like it has a bunch of fitting routines.

Thanks a lot Anders! I will have a look at it to replace my numpy routines!