for the uncertainties I was using the values obtai...
# mathematics
b
for the uncertainties I was using the values obtained from
np.sqrt(np.diag(pcov))