A general three-dimensional eddy-current probe model, developed by Sabbagh Associates and reported in [1], [2] and [3], has been adapted for the calculation of probe-flaw interactions. The theoretical model, [4] and [5], uses integral equations with dyadic Green’s function kernels, and is applicable to both probe and flaw calculations at arbitrary skin depths and frequencies. Discrete approximations of the integral equations are solved using a highly efficient algorithm based on recent developments in numerical techniques and their application to the solution of large problems in electromagnetic field-theory.