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