I need to fit a curve over a set of data and also need the uncertainty - or errors - of the coefficients, for exemple:
fitting ax^2+bx+c, i need the values: a+-da, b+-db and c+-dc. Where da,db and dc are the uncertaints.
I already tried polyfit and optmize.curve_fit, but none of them give me de uncertainty as I want. Some one knows how to do that?
Here is code that works for both linear and nonlinear fitting and should be useful: