We report on a numerical implementation of the QED one-loop 5-point functions. These functions contribute to the NLO corrections to the hard-bremsstrahlung process e⁺ e^- \to \mu ⁺ \mu ^- \gamma , which is an important background process for accurate predictions in experiments at high luminosity meson factories such as DAFNE and PEP-II. The calculation is implemented using publicly available tools and incorporates several numerical and analytical cross-checks. Numerical precision and stability is demonstrated by preliminary test runs with KLOE and BarBar kinematical cuts.

PACS numbers: 12.20.Ds, 13.66.Jn