# Computational Procedure Based on the Block Gauss-Seidel Algorithm

## Previous Calculation Procedure

In [23], a calculation procedure to solve the equation (5-36) has been proposed:

- 1. First, the fundamental harmonic
*A*is calculated based on the assumption of setting A_{1}_{2}to A_{n}to zero. If the convergence of*A*_{1}is satisfied, turn to step 2. If not, update*A*_{1}and repeat step 1. - 2. The other harmonics,
*A*_{2}to*A*_{N}, are calculated respectively with the known*A*_{1}. If the convergence of A_{2}to*A*is satisfied, the calculation procedure stops. If not, turn to step 1._{N}

However, the computational procedure above is not dependable nor efficient. The convergent criterion in each harmonic computation and the lack of full use of updated solutions in nonlinear iterations leads to convergence uncertainty and the subsequent inaccuracy of harmonic solutions, especially when the nonlinearities of the magnetic material are strong. Yamada [23] also pointed out that there is convergence uncertainty in the iterative approach. In fact, a compulsory stop is usually done, by setting a maximum number of iterative steps.