By expanding functions of parton fragmentation into a heavy hadron in
the inverse of the heavy quark mass m(Q) we attempt to factorize them
into perturbative and non-perturbative parts. In our approach the non-
perturbative parts can be defined as matrix elements in heavy quark ef
fective theory, the shape of the functions being predicted by perturba
tive QCD. In this work we neglect the effect at order of m(Q)(-2) and
calculate the perturbative parts at one-loop level for heavy quark and
gluon fragmentation. We compare our results from the leading log appr
oximation with experimental results from e(+)e(-) colliders and find a
deviation below or at the 10% level. Adding the effect of higher orde
rs in alpha(s) it can be expected to reduce this deviation The size of
the matrix elements appearing at the order we consider for several ty
pes of heavy hadrons is determined. (C) 1997 Elsevier Science B.V.