Error estimation methods for the measured SED using the VERITAS dataset, as the example case. The solid black lines indicate the upper and lower limits of the error band estimated with the multivariate sampling. These lines represent the 68% containment of 500 spectral realizations (100 shown as grey lines) whose parameters are sampled from a multivariate distribution defined by the fit results.