x = linspace(0, 1, 1000)
y = cos(x) * x ** 2 + x + sin(x - 1.0)
p, cov = polyfit(x, y, 2, cov=True)
plot(x, y, 'b')
plot(x, polyval(p, x), 'r')
print(sqrt(diag(cov)))
Найдите неопределенность из полифита