      function fermi(x)
      common/pawpar/par(3)
      fermi=par(1)/(1.+exp((x-par(2))/par(3)))
      return
      end
