The numerical S-matrix Bootstrap exploits the analyticity, crossing and unitarity properties of scattering amplitudes to bound physical observables. Low energy constants in EFTs represent a natural target to test those numerical methods. In this talk I will focus on a specific class of massless EFTs which include QCD in the chiral limit. I will review the numerical formulation of the problem, I will show our numerical results and compare them with existing experimental data.