  We present a calculation of the next-to-leading order cross section for the
inclusive hadroproduction of D and B mesons as a function of the transverse
momentum and the rapidity in a massive variable flavor number scheme. We
compare our numerical results with recent data from the CDF Collaboration at
the Fermilab Tevatron for the production of D^0, D^{*+}, D^+, and D_s^+ mesons
at center-of-mass energy $\sqrt{S} = 1.96$ TeV and find reasonably good
agreement with the measured cross sections.
