Previously confined to simplified models or computationally intensive approximations, accurate simulation of resonance Raman spectra for complex molecules is now achievable. A new technique employing Hagedorn wavepackets offers a practical ab initio solution validated against established methods in two dimensions before being applied successfully to anthracene with its 66-dimensional harmonic potential. This advancement enables computation of excitation profiles, including hot bands, for diverse spectral signals without relying on extensive numerical calculation.