Top quark mass corrections to single and double Higgs boson production in gluon fusion

In this talk we discuss recent computations of the top quark mass dependence of QCD amplitudes describing Higgs boson production in gluon fusion. We compute terms in the expansion for a large top quark mass, which reduces the Feynman diagrams to products of massless integrals and massive tadpole integrals which contain the top mass dependence. In particular we discuss the real and virtual corrections to double Higgs production at NNLO, and the virtual corrections to single Higgs production at N3LO.