In this talk I will 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. I 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.
