ToolPythonBeginner Numerical integration with scipy.integrate: the area under a measured peak net = counts - np.polyval(np.polyfit(lam[edge], counts[edge], 1), lam)line = counts[win] - np.polyval(p, lam[win] - peak)resid = counts[flank] - np.polyval(p, lam[flank] - peak)baseline = np.polyval(p, lam - peak)