fix Brillouin function and its derivative (in 5th and 7th orders)
fix Langevin function and add macro to define its applicability range
handle special case of infinite spin to msat0 estimation
handle special case of infinite spin in longitudinal susceptibility estimation
handle special case of infinite spin in longitudinal field estimation