This article describes the code SusHi (for “Supersymmetric Higgs”) [108] 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 pT 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. Program SummaryProgram title: SusHiCatalogue identifier: AEOY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEOY_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 47725No. of bytes in distributed program, including test data, etc.: 338380Distribution format: tar.gzProgramming language: Fortran 77.Computer: Personal computer.Operating system: Unix/Linux, Mac OS.RAM: A few 100 MBClassification: 11.1.External routines: LHAPDF (http://lhapdf.hepforge.org), FeynHiggs (http://www.feynhiggs.de)Nature of problem:Calculation of inclusive and exclusive Higgs production cross sections in gluon fusion and bottom-quark annihilation in the Standard Model and the MSSM through next-to-leading order QCD, includes next-to-next-to-leading order top-(s)quark contributions and electro-weak effectsSolution method:Numerical Monte Carlo integrationAdditional comments:Code includes ggh@nnlo (http://www.robert-harlander.de/software/ggh@nnlo), bbh@nnlo (http://www.robert-harlander.de/software/bbh@nnlo) and evalcsusy.f (http://www.robert-harlander.de/software)Running time:All the examples provided take a minute or less to run.