The relaxation spectrum of a viscoelastic material holds the key to describing its relaxation mechanisms at a molecular level. It also plays a fundamental role in accessing the molecular weight distribution, and in modelling the dynamics of complex uids. The relaxation spectrum cannot be measured directly, but it may be locally determined from experimental measurements of viscoelastic response at a macroscopic level. In particular, the relaxation spectrum is a continuous distribution of relaxation times which may be recovered, at least locally, from measurements of the complex modulus of the material. Although mathematical expressions for the continuous spectrum have been known for over a century, these were inaccessible to numerical implementation for decades, since they involve inverse operators which are not continuous, resulting in severe instability. Progress was made when regularization methods for approximating discrete line spectra were introduced some two decades ago. It was not until 2012, however, that Davies and Goulding proposed a method of wavelet regularization for recovering continuous spectra in a mathematically rigorous framework. This work was further rened in 2016 by introducing a mathematical form of high-order derivative spectroscopy involving sequences of derivatives of dynamic moduli, termed Maclaurin sequences. In this article, a rigorous justication for the use of Maclaurin sequences is presented. Furthermore, a new sequence is presented, which is termed a wavelet correction sequence, achieving the same accuracy as Maclaurin sequences, but with a reduced order of dierentiation.