Added Boltzmann distribution function BOLTZ_E