We present the Maxwell-Bloch equations that are a model for the semi-classical description of laser-matter interactions. After having discussed both models and their numerical coupling in a general context regarding relaxation terms and the number of energy levels, we give examples of simulations that show the capacities of our approach.