Approximation of the inverse Langevin function revisited

The main purpose of this paper is to provide an easy-to-use approximation formula for the inverse Langevin function. The mathematical complexity of this function makes it unfeasible for an analytical manipulation and inconvenient for computer simulation. This situation has motivated a series of papers directed on its approximation. The best known solution is given by Cohen. It is used in a lot of statistically based models of rubber-like materials. The formula is derived from rounded Pad, approximation [3/2]. The main idea of the presented approach in this paper relies on improvement of the precision of approximation formula for the inverse Langevin function by using multipoint Pad, approximation method. We focused our study strongly on obtaining a simple and accurate approximation. It is assumed that the proposed approximation formula may be considered a useful tool for verification of the results obtained in other ways. Our results are supported by investigating several numerical examples. The paper also presents a few applications of computer software named Mathematica which can be used to calculate symbolically one point Pad, approximants and numerically multipoint Pad, approximants. Using this software, we showed also how to compute higher order derivatives of the inverse function in a simple and elegant way. This issue was discussed by Itskov et al.