real FUNCTION HIGGS_POTENTIAL(X,Y)
************************************************************************
* Makes the plot of the SM Higgs potential for a 115 GeV Higgs mass.   *
************************************************************************      
      implicit none
      real x,y,z
      real mu,lambda
      parameter (mu2=-6612.5)     ! For mh=115 GeV and v=250 GeV
      parameter (lambda=0.1058)

      higgs_potential = mu2*(x**2+y**2)+lambda*(x**4+2*x**2*y**2+y**4)

      END