Two-dimensional static field problems can be solved by a method based on Fredholm integral equations (equations of the second kind). This has numerical advantages over the mote commonly used integral equation of the first kind. The method is applicable to both magnetostatic and electrostatic problems formulated in terms of either vector or scalar potentials. It has been extended to the solution of eddy current problems with sinusoidal driving functions. The application of the classical Fredholm equation has been extended to problems containing boundary conditions: 1) potential value, 2) normal derivative value, and 3) an interface condition, all in the same problem. The solutions to the Fredholm equations are single or double (dipole) layers of sources on the problem boundaries and interfaces. This method has been developed into computer codes which use piecewise quadratic approximations to the solutions to the integral equations. Exact integrations are used to replace the integral equations by a matrix equation. The solution to this matrix equation can then be used to directly calculate the field anywhere.