  This article describes the code SusHi (for "Supersymmetric Higgs") which
calculates the cross sections $pp/p\bar{p}\rightarrow\phi+X$ in gluon fusion
and bottom-quark annihilation in the SM and the MSSM, where $\phi$ is any of
the neutral Higgs bosons within these models. Apart from inclusive cross
sections up to NNLO QCD, differential cross sections with respect to the Higgs'
transverse momentum $p_T$ and (pseudo-)rapidity $y(\eta)$ can be calculated
through NLO QCD. In case of gluon fusion, SusHi contains NLO QCD contributions
from the third family of quarks and squarks, NNLO corrections due to
top-quarks, approximate NNLO corrections due to top-squarks, and electro-weak
effects. It supports various renormalization schemes for the sbottom sector and
the bottom Yukawa coupling, as well as resummation effects of higher order
$\tan\beta$-enhanced sbottom contributions. SusHi provides a link to FeynHiggs
for the calculation of the Higgs masses.
