An efficient Monte Carlo method for order-alpha photon bremsstrahlung
calculations in allowed nuclear, neutron and pion beta decays is prese
nted. This method is useful for radiative correction calculations and
for weighted or unweighted event generations in experimental off-line
data analyses. The most important intermediate steps of the calculatio
ns are also described. Several numerical results of the Monte Carlo pr
ogram are compared with published analytical and numerical quadrature
calculation results.