SusHi: a program for the calculation of Higgs production in gluon fusion and bottom-quark annihilation in the standard model and the MSSM. This article describes the code SusHi (for “Supersymmetric Higgs”) which calculates the cross sections pp/pp ¯→ϕ+X in gluon fusion and bottom-quark annihilation in the SM and the MSSM, where ϕ 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(η) can be calculated through NLO QCD. In the 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 β-enhanced sbottom contributions. SusHi provides a link to FeynHiggs for the calculation of the Higgs masses.