bose-einstein distribution function