Inspired by the recent success of the numerical approach to the conformal bootstrap, we revisit the S-matrix bootstrap program. We shall explain how to obtain analytic bounds on the interaction strength in 1+1 QFT. In higher dimensions, we propose a numerical algorithm that seems to converge to optimal bounds.