In this Letter we present an application of the correlation function approach to calculate the Franck-Condon spectra of polyatomic molecules. A coupled-cluster method was used for the time propagation of the wavefunction using a quadratic Hamiltonian. The dispersed fluorescence spectra from various single vibronic excitations of trans-hydroquinone were calculated and compared with the experimentally observed spectra. The agreement between the calculated and experimental spectra was found to be excellent after some empirical adjustments to the potential displacements.