  A next-to-leading order Monte Carlo program for the calculation of heavy
quark cross sections in deeply inelastic scattering is described. Concentrating
on charm quark and D^*(2010) production at HERA, several distributions are
presented and their variation with respect to charm quark mass, parton
distribution set, and renormalization-factorization scale is studied.
