<p>The Jiles-Atherton model has been widely used in describing the hysteretic property of a magnetic material or device. However, the calculation errors are not so easily discovered. With a complex expression, the frequency-dependent Jiles-Atherton model should be solved numerically with appropriate settings. This paper proposes an effective solving method for this model and describes some necessary analysis built on the numerical results. In the numerical method proposed in this manuscript, the anhysteretic magnetization was calculated by the secant method, and the trapezoidal rule was utilized to form the implicit function, which can be calculated by the fixed-point iteration. Compared to the other common methods, the proposed one has a friendly expression and fast computation speed. The Terfenol-D material was taken as an example for the numerical analysis. The feasible region was determined and the commonly used approximation that neglects the term of the magnetic field when calculating the magnetic induction intensity was tested. At last, the required number of sampling points per period was reached to guarantee high precision from analyzing its influence on the computation precision. The proposed numerical method is helpful for high-precision solutions of the frequency-dependent Jiles-Atherton model. The results from the numerical analysis can also help users avoid some incorrect calculations when employing this hysteresis model.</p>