  We provide nonperturbative fragmentation functions for B mesons, both at
leading and next-to-leading order in the MS-bar factorization scheme with five
massless quark flavors. They are determined by fitting the fractional energy
distribution of B mesons inclusively produced in e^+ e^- annihilation at CERN
LEP1. Theoretical predictions for the inclusive production of B mesons with
high transverse momenta in p p-bar scattering obtained with these fragmentation
functions nicely agree, both in shape and normalization, with data recently
taken at the Fermilab Tevatron.
