We use the numerical S-matrix bootstrap method to obtain bounds on the two leading Wilson coefficients of the chiral lagrangian controlling the low-energy dynamics of massless pions thus providing a proof of concept that the numerical S-matrix bootstrap can be used to derive non-perturbative bounds on EFTs in more than two spacetime dimensions.