standard_errors = np.sqrt(np.diag(covariance))
residuals = y - decay_model(x, *params)
rmse = np.sqrt(np.mean(residuals**2))
condition_number = np.linalg.cond(covariance)