Function bose_function¶
-
double o2scl::bose_function(double E, double mu, double T, double limit = 40.0)¶
Calculate a Bose-Einstein distribution function safely.
\( \left[\exp\left(E/T-\mu/T\right)-1\right]^{-1} \)
This function computes a Bose-Einstein distribution function using an expansion to ensure accurate results when \( (E-\mu)/T \) is near zero.
This function is used in the
o2scl::boson_rel
class in .Note
This function may return Inf or NAN if
limit
is too large, depending on the machine precision.