Error calculated as from Eq. (31) as a function of spatial (and temporal) resolution. N represents the number of points used to dis-cretize the domain. The points in the plot correspond to N = 4096 (=N_max), 2048, 1024, 512, 256. Different combinations of flux schemes (TVDLF, HLL) and limiters (minmod, cada3, woodward, mp5) are shown by different curves, some of them overlap and are grouped as indicated in the legend.