plt.plot(x, y, '.', x, k * x ** 2, 'r', linewidth=3)
plt.legend(('measurement', 'fit'), loc=2)
plt.title('best fit: y = {:.8f} * x**2'.format(k[0]))
plt.show()
numpy.polyfit с адаптированными параметрами