solvePoly
Timing:
Utilities.py
NAM: solvePoly (function)
DES: perform polyomial interpolation: solve linear system
dataY = P_n(dataX)
INP: (int) order : polynomial degree
(flt arrays) dataX/Y : system to solve
OUT: (flt array) coeff : polynomial coefficients