  We describe the hexagon function bootstrap for solving for six-gluon
scattering amplitudes in the large $N_c$ limit of ${\cal N}=4$ super-Yang-Mills
theory. In this method, an ansatz for the finite part of these amplitudes is
constrained at the level of amplitudes, not integrands, using boundary
information. In the near-collinear limit, the dual picture of the amplitudes as
Wilson loops leads to an operator product expansion which has been solved using
integrability by Basso, Sever and Vieira. Factorization of the amplitudes in
the multi-Regge limit provides additional boundary data. This bootstrap has
been applied successfully through four loops for the maximally helicity
violating (MHV) configuration of gluon helicities, and through three loops for
the non-MHV case.
