  A method based on the analyzis of the Majorana polynomial roots is introduced
to compute the spectrum of the Lipkin-Meshkov-Glick model in the
thermodynamical limit. A rich structure made of four qualitatively different
regions is revealed in the parameter space whereas the ground state study only
distinguishes between two phases.
