We compare the leading and next-to-leading order QCD predictions for the fl
ux of atmospheric muons and neutrinos from decays of charmed particles. We
find that the full NLO lepton fluxes can be approximated to within similar
to 10% by the Born-level fluxes multiplied by an overall factor of 2.2-2.4,
which depends slightly on the PDF. This supports the approach of Thunman,
Ingelman and Gondolo. We also find that their very low lepton fluxes are du
e to the mild slope they used for the gluon distribution function at small
momentum fractions, and that substantially larger lepton fluxes result when
the slope of the gluon distribution function at small momentum fractions i
s larger.