The calculation of Higgs boson pair production in gluon fusion at next-to-leading order in the strong coupling constant with full top-quark mass dependence is presented. The usage of the programs GoSam, Reduze and SecDec for the computation and numerical evaluation of the most complicated piece, the virtual two-loop amplitude, is explained. A comparison of the cross section of the exact result with previous approximations is shown as well as the invariant mass and pT distribution are shown.