The equation of state of nuclear matter suggests that at suitable beam energies the disassembling hot system formed in heavy ion collisions will pass through a liquid-gas coexistence region. Searching for the signatures of the phase transition has been a very important focal point of experimental endeavours in heavy-ion collisions, in the last fifteen years. Simultaneously theoretical models have been developed to provide information about the equation of state and reaction mechanisms consistent with the experimental observables.