A method for simulating a laser machining using multi-pulse irradiation of a material by a laser machining system capable of emitting a machining laser beam, said method being characterized by the following steps: a) providing a central unit with: - information relating to the interaction between the machining laser beam and the material to be machined comprising: ∘ a delta δ relative to a penetration depth of the machining laser beam in the material, ∘ a threshold fluence F th relative to a minimum energy density of the laser beam allowing a machining, ∘ an incubation coefficient S, between 0 and 1, relating to a dependence of the threshold fluence F th in function of a number of pulses of the laser beam at a single point, ∘ a complex refractive index n + ik, - information relating to said laser machining system comprising: ∘ information relating to a polarization of said machining laser beam, ∘ a pulse energy of said machining laser beam E P, ∘ a diameter of said machining laser beam at a focal point w, ∘ an order of a Gaussian p of said machining laser beam between 1 and 20, ∘ a pulse repetition rate PRR of pulses n of the machining laser beam, ∘ a wavelength of the machining laser beam; b) determining with said central unit on the basis of: - said information relating to the interaction between the machining laser beam and the material to be machined, and, - said information relating to said laser machining system, a machining profile in two and/or three dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system the method further comprising taking into account the increase in surface depth after each consecutive laser pulse. a) providing a central unit with: - information relating to the interaction between the machining laser beam and the material to be machined comprising: ∘ a delta δ relative to a penetration depth of the machining laser beam in the material, ∘ a threshold fluence F th relative to a minimum energy density of the laser beam allowing a machining, ∘ an incubation coefficient S, between 0 and 1, relating to a dependence of the threshold fluence F th in function of a number of pulses of the laser beam at a single point, ∘ a complex refractive index n + ik, - information relating to said laser machining system comprising: ∘ information relating to a polarization of said machining laser beam, ∘ a pulse energy of said machining laser beam E P, ∘ a diameter of said machining laser beam at a focal point w, ∘ an order of a Gaussian p of said machining laser beam between 1 and 20, ∘ a pulse repetition rate PRR of pulses n of the machining laser beam, ∘ a wavelength of the machining laser beam; - information relating to the interaction between the machining laser beam and the material to be machined comprising: ∘ a delta δ relative to a penetration depth of the machining laser beam in the material, ∘ a threshold fluence F th relative to a minimum energy density of the laser beam allowing a machining, ∘ an incubation coefficient S, between 0 and 1, relating to a dependence of the threshold fluence F th in function of a number of pulses of the laser beam at a single point, ∘ a complex refractive index n + ik, ∘ a delta δ relative to a penetration depth of the machining laser beam in the material, ∘ a threshold fluence F th relative to a minimum energy density of the laser beam allowing a machining, ∘ an incubation coefficient S, between 0 and 1, relating to a dependence of the threshold fluence F th in function of a number of pulses of the laser beam at a single point, ∘ a complex refractive index n + ik, - information relating to said laser machining system comprising: ∘ information relating to a polarization of said machining laser beam, ∘ a pulse energy of said machining laser beam E P, ∘ a diameter of said machining laser beam at a focal point w, ∘ an order of a Gaussian p of said machining laser beam between 1 and 20, ∘ a pulse repetition rate PRR of pulses n of the machining laser beam, ∘ a wavelength of the machining laser beam; ∘ information relating to a polarization of said machining laser beam, ∘ a pulse energy of said machining laser beam E P, ∘ a diameter of said machining laser beam at a focal point w, ∘ an order of a Gaussian p of said machining laser beam between 1 and 20, ∘ a pulse repetition rate PRR of pulses n of the machining laser beam, ∘ a wavelength of the machining laser beam; b) determining with said central unit on the basis of: - said information relating to the interaction between the machining laser beam and the material to be machined, and, - said information relating to said laser machining system, - said information relating to the interaction between the machining laser beam and the material to be machined, and, - said information relating to said laser machining system, a machining profile in two and/or three dimensions corresponding to the simulation of a machining of said material to be machined with said laser machining system the method further comprising taking into account the increase in surface depth after each consecutive laser pulse.