In this paper we compute the QCD corrections for the triple Higgs boson production cross section via gluon fusion, within the heavy-top approximation. We present, for the first time, analytical results for the next-to-leading order corrections, and also compute the soft and virtual contributions of the next-to-next-to-leading order cross section. We provide predictions for the total cross section and the triple Higgs invariant mass distribution. We find that the QCD corrections are large at

more »

... perturbative orders, and that the scale uncertainty is substantially reduced when the second order perturbative corrections are included.