Eddy current nondestructive evaluation is a widely used in-plant NDE technique in which the flaw information is extracted from the impedance change of a coil placed above the metal testpiece. To obtain quantitative information about flaw size and shape, we would like to have a theoretical model which is able to predict the impedance change of a coil for different flaws in the test geometry. Because of its importance, this eddy current forward problem has been studied extensively for many years. For simple flaw shapes and geometry, it is possible to obtain analytical solutions. However, for flaws with irregular shapes and complex geometry, an analytical solution usually is not available so we must find a numerical solution. There have been several numerical models in the literature, e.g., the finite element method[l], the boundary element method[2], the volume integral method[3–5] and methods based on variational formulas[6].