Switch to curve_fit wrapper to have a nice covariance matrix (diagonal returned).