We use the S-matrix bootstrap to carve out the space of unitary, analytic, crossing symmetric and supersymmetric graviton scattering amplitudes in nine, ten and eleven dimensions. We extend and improve the numerical methods of our previous work in ten dimensions. A key new tool employed here is unitarity in the celestial sphere. In all dimensions, we find that the minimal allowed value of the Wilson coefficient α, controlling the leading correction to maximal supergravity, is very close but not equal to the minimal value realized in Superstring theory or M-theory. This small difference may be related to inelastic effects that are not well described by our numerical extremal amplitudes. Although α has a unique value in M-theory, we found no evidence of an upper bound on α in 11D.