We introduce techniques to treat numerically Mellin-Barnes integrals in physical regions, which arise in the need of the computation of Feynman integrals for the electroweak two-loop corrections to the pseudo observables at the Z-boson resonance.