An iterative procedure using finite element method without the Lagrange mul
tiplier is proposed for three-dimensional eddy current problems, which is b
ased on an iterative procedure derived from a perturbation problem of the m
agnetostatic problem. To consider the continuity of an excitation current d
ensity, a correction method is also proposed. Numerical results show that t
he BiConjugate Gradient (BiCG) method is applicable to the complex symmetri
c linear systems arising in the iterative procedure, and that approximate p
hysical quantities are suitable.