In [3]:
import numpy as np
x = [1, 3, 4, 5, 6, 7, 8]

y = [i*i for i in x]

np.polyfit(x, y, 2)


Out[3]:
array([ 1.00000000e+00,  5.70919174e-15, -8.14710295e-15])

In [ ]: