In this paper, we set up the numerical S-matrix bootstrap by using the crossing symmetric dispersion relation (CSDR) to write down Roy equations for the partial waves. As a motivation behind examining the local version of the CSDR, we derive a new crossing symmetric, 3-channels-plus-contact-terms representation of the Virasoro-Shapiro amplitude in string theory that converges everywhere except at the poles. We then focus on gapped theories and give novel analytic and semi-analytic derivations of several bounds on low-energy data. We examine the high-energy behaviour of the experimentally measurable rho-parameter, introduced by Khuri and Kinoshita and defined as the ratio of the real to the imaginary part of the amplitude in the forward limit. Contrary to expectations, we find numerical evidence that there could be multiple changes in the sign of this ratio before it asymptotes at high energies. We compare our approach with other existing numerical methods and find agreement, with improvement in convergence.