We measure the second-harmonic (SH) response of MoS2 in a broad fundamental photon energy range between 0.85 and 1.7 eV, and present a continuous SH spectrum capturing signatures of both the fundamental A and B excitons in addition to the intense C resonance. Moreover, we interpret our results in terms of the first exciton simulation of the SH properties of multi-layered MoS2 samples, represented by a study of trilayer MoS2. The good agreement between theory and experiments allows us to establish a connection between the measured spectrum and the underlying electronic structure, in the process elucidating the non-linear excitation of electron-hole pairs.In recent years, two-dimensional transition metal dichalcogenides (TMDs), being semiconductors 1-4 with band gaps in the 2-3 eV range 4-7 , have been applied as the basis for a host of novel two-dimensional optoelectronic devices [8][9][10][11][12][13] . Accordingly, a great deal of attention has been devoted to the linear optical properties of these compounds 1-3 . Also, non-linear optical techniques have been demonstrated as particularly powerful probes for the microscopic structure of few-layered TMD crystals, revealing e.g., their crystallographic orientation 7,14-17 , and stacking angle 18 . Moreover, second-harmonic (SH) spectroscopy allows for insights into the electronic structure of TMD flakes 14 , exposing edge-localized effects 19 , or valley-coherent excitations 7 .While experimentally determined linear response functions of various TMDs, dominated by electron-hole pairs bound by several hundred meVs 5,6,20 , are reproduced reasonably well at the Bethe-Salpeter equation (BSE) level of complexity 6,21-23 , similar comparisons for non-linear cases are lacking. Although several proposed theoretical models, based on both independent-particle 19,24 and excitonic 22,23 approaches, have been published, the SH spectrum of MoS 2 has been investigated experimentally 14 only in a narrow fundamental photon energy range between 1.2 and 1.7 eV, probing the so-called C resonance also known from linear optics. Theoretical models have so far been benchmarked by their ability to reproduce this single feature [22][23][24] . However, the parabolic dispersion, and the split valence bands near the K-points of the Brillouin zone, translates into bound electron-hole pairs in the BSE picture. These give rise to sharp peaks in the absorption spectrum due to the so-called A and B excitons 6,25 at photon energies near 1.9 eV, and the question remains how they might affect the SH response. In particular, the relative intensities and lineshapes of such features are of interest.In the present letter, we report an experimental optical SH spectrum generated from many-layered MoS 2 , with fundamental photon energies varied in a broad range between 0.85 and 1.7 eV, capturing both the fundamental A and B excitons in addition to the aforementioned C resonance. We perform our optical experiments on (i) many-layered MoS 2 flakes exfoliated onto fused silica by the well-known scotch tape me...