In this talk we discuss computations of NLO virtual corrections to four-point gluon-fusion processes; in particular the production of HH, ZH and ZZ. Recently these processes have been computed numerically, but they are not known analytically. We will discuss how one can perform an expansion of these amplitudes in the high-energy limit, and improve the resulting series through the use of Padé approximants.