  The next-to-leading order (NLO) QCD corrections to the production and decay
rate of a Higgs boson are computed within the framework of the Minimal
Supersymmetric Standard Model (MSSM). The calculation is based on an effective
theory for light and intermediate mass Higgs bosons. We provide a Fortran
routine for the numerical evaluation of the coefficient function. For most of
the MSSM parameter space, the relative size of the NLO corrections is typically
of the order of 5% smaller than the Standard Model value. We exemplify the
numerical results for two scenarios: the benchmark point SPS1a, and a parameter
region where the gluon-Higgs coupling at leading order is very small due to a
cancellation of the squark and quark contributions.
