This paper describes an application of modern methods of computational electromagnetics to the problem of eddy-current nondestructive evaluation (NDE). Specifically, a volume integral equation is developed that can be used to model eddy-current probes with ferrite cores, or it can be used to compute flaw-field interactions. Both of these problems are of considerable importance in applying electromagnetic techniques to NDE, for they are intimately involved in detecting flaws and inverting eddy-current data for reconstructing flaws.