RecipePythonBeginner Fit a curve to data with error bars and draw a confidence band J = np.empty((tt.size, popt.size)) # d model / d parameter