Numerical simulation is a valuable tool for the study of magnetic susceptibility artifacts from metallic implants. A major difficulty in the simulation lies in the computation of the magnetic field induced by the metallic implant. A new method has been designed and implemented to compute the magnetic field induced by metallic objects of arbitrary shape. The magnetic field is expressed pointwise in terms of a surface integral. In magnetic resonance imaging (MRI) even minor perturbations of the gradients can disturb the imaging process and may render the clinical image inaccurate or useless. Common causes of magnetic field perturbation are changes of magnetic properties in the sample due to metallic implanted objects, such as a dental prosthesis, orthopedic apparatus, etc. Metallic surgical instruments used in interventional MRI are also responsible for magnetic field perturbation. In the case of guided stereotactic surgery, image distortions are especially harmful since geometrical accuracy is required.In the past, several studies involving numerical simulation have been undertaken to help in the understanding of magnetic field perturbation and to establish experimental conditions that determine susceptibility artifacts (1-3). However, these simulation attempts have been restricted to simple test objects (cylinders, spheres, and ellipsoids) for which an analytical expression for the induced magnetic field is known. In general, a precise calculation of the magnetic field involves a boundary value problem with partial differential equations (PDE) derived from Maxwell's equations, and requires the use of PDE approximation schemes. Among the classical numerical methods are the finite element method (FEM), the finite difference method (FDM), and the boundary element method (BEM). For a comprehensive treatment of these numerical methods in electromagnetism see Ref. 4. Our approach, which is based on a surface integral representation formula for the magnetic flux density, is most similar to the BEM. In short, our method consists of expressing the problem in an integral form over the boundary of the implant, and dividing the boundary into elements in which the integral is numerically computed. It is able to deal with objects of any shape, provided that a mesh of the object boundary is available. This requirement is not at all restricting since mesh generation tools are widely used in computer assisted design (CAD). The advantages of our method compared to the FEM (5) or the FDM (6) are numerous. First, the computation depends only on the object boundary, so the discretization of space is reduced from 3D to 2D. Furthermore, in the FDM/FEM the exterior domain must be truncated and an approximation of the behavior of the field at infinity must be introduced on an artificial boundary. With the BEM, the behavior at infinity is always exactly satisfied. Finally, for 3D problems the FEM and the FDM lead to large linear systems to be solved, whereas in the proposed method the solution is obtained pointwise by evaluating a surfac...