Abstract
In this work it is demonstrated that a straightforward implementation of the commonly used textbook expressions of the magnetic field of a current loop is numerically unstable in a large region of interest. A novel expression based on a modification of Bulirsch’s algorithm for the complete elliptic integral is proposed. It is numerically stable and faster to compute than the classical one. New and old expressions are compared to series expansions and dipole approximation in terms of accuracy, numerical stability and computation performance.
Original language | English |
---|---|
Article number | 1 |
Pages (from-to) | 11-31 |
Number of pages | 20 |
Journal | Magnetism |
Volume | 3 |
Issue number | 2673-8724 |
Publication status | Published - 31 Dec 2022 |