In this talk I will discuss the computation of the top quark mass dependence of NNLO double Higgs boson production, including both the real and virtual contributions. The dependence is computed in the large top quark mass limit, which gives a good approximation of the exact cross section below the top quark production threshold.