The reluctivity matrix is obtained by using the trigonometric function transformation and the form of coefficient matrix is simplified.

该方法适用于任意阶多项式描述磁化曲线的非线电磁场问 的求解，并计算的精度。