The production of heavy quarks at large () in hadronic collisions is considered. The analysis is carried out in the framework of perturbative fragmentation functions, thereby allowing a resummation at the NLO level of final state large mass logarithms of the kind . The case of -quark production is considered in detail. The resulting theoretical uncertainty from factorization/renormalization scales at large is found to be much smaller than that shown by the full perturbative calculation.