This is a list of changes to higlu.f since October 1995:

February 18 2010: - version 2.400:
                    Electroweak corrections to SM Higgs boson production
                    implemented in factorized form (thanks to G. Passarino
                    et al.). Please, cite the following papers for these
                    corrections:

%\cite{Actis:2008ts}
\bibitem{Actis:2008ts}
  S.~Actis, G.~Passarino, C.~Sturm and S.~Uccirati,
  %``NNLO Computational Techniques: the Cases H -> gamma gamma and H -> g g,''
  Nucl.\ Phys.\  B {\bf 811} (2009) 182
  [arXiv:0809.3667 [hep-ph]].
  %%CITATION = NUPHA,B811,182;%%

%\cite{Actis:2008ug}
\bibitem{Actis:2008ug}
  S.~Actis, G.~Passarino, C.~Sturm and S.~Uccirati,
  %``NLO Electroweak Corrections to Higgs Boson Production at Hadron
  %Colliders,''
  Phys.\ Lett.\  B {\bf 670} (2008) 12
  [arXiv:0809.1301 [hep-ph]].

%\cite{Actis:2008uh}
\bibitem{Actis:2008uh}
  S.~Actis, G.~Passarino, C.~Sturm and S.~Uccirati,
  %``Two-Loop Threshold Singularities, Unstable Particles and Complex Masses,''
  Phys.\ Lett.\  B {\bf 669} (2008) 62
  [arXiv:0809.1302 [hep-ph]].
  %%CITATION = PHLTA,B669,62;%%

November 20 2009: - Link to LHAPDF implemented. Fortran compiler
                    gfortran added to makefile. Links to PDFLIB, MRST2001
                    and MSTW2008 removed (included in LHAPDF).

October 16 2009: - Option to choose between bottom pole mass and bottom
                   MSbar mass (including independent scale choice) added

                 **********************************************************
                 BOTTOM MASS:  0 = POLE MASS   1 = MSBAR MASS
                 ============  SCALE = MU_1*M_HIGGS + MU_2 (FOR MSBAR MASS)

                 BOTMASS  = 0
                 MU_1     = 0.D0
                 MU_2     = 10.D0
                 **********************************************************

                 Careful: Best choices are pole mass or MSbar bottom
                 mass at a scale of O(m_b) due to large logs if the
                 scale is of O(Higgs mass)! (see Fig.7a of NPB453 (1995) 17)

                 --> version 2.300

October 15 2009: - Link to MSTW2008 PDFs
                 - dsigma/p_T and d2sigma/dp_T/dy distributions added (LO):

                 **********************************************************
                 DISTRIBUTIONS: 0=TOT CXN  1=D2SIGMA/DP_T/DY  2=DSIGMA/DP_T
                 =============  (ONLY FOR GG --> H)

                 DISTRIB  = 0
                 P_T      = 100.D0
                 Y        = 0.D0
                 **********************************************************

                 --> version 2.200

April 3 2006: - Adding delta_mb corrections to bottom Yukawa couplings

March 1 2006: - Removing an instability in the subroutine AMHAMA
                (thanks to P. Trueb)

April 12 2005: - Link to CTEQ6 and MRST2001 PDFs
                 (NGROUP = -1: CTEQ6    -2: MRST2001)

November 7 2003: - Lines 919,920 added due to problems with g77 compiler
                   (thanks to M. Schumacher)

November 6 2003: - Lines 491-495 commented out to allow arbitrary Yukawa
                   couplings (thanks to M. Schumacher)

February 28 2003: - Extension of MSSM subroutines such that they are in
                    agreement with HDECAY (version 3.0) now.