FeynHiggs: a program for the calculation of the masses of the neutral   CP-even Higgs bosons in the MSSM
  FeynHiggs is a Fortran code for the calculation of the masses of the neutral CP-even Higgs bosons in the MSSM up to two-loop order. It is based on the complete diagrammatic on-shell results at the one-loop level, the leading diagrammatic two-loop QCD contributions and further improvements taking into account leading electroweak two-loop and leading higher-order QCD corrections. The Higgs-boson masses are calculated as functions of the MSSM parameters for general mixing in the scalar top sector and arbitrary choices of the parameters in the Higgs sector of the model. 
S. Heinemeyer, W. Hollik, G. Weiglein