In this paper, the use of the standard finite difference time domain (FDTD) algorithm is extended to quasi-static electromagnetic field problems. While straightforward application of the standard FDTD algorithm at very low frequencies leads to excessively long simulation times, we show that for linear structures this problem can be circumvented by using a ramp excitation function. The use of appropriate absorbing boundary conditions such as Berenger's perfectly matched layer is also shown to be necessary. By combining two plane waves in opposite directions, a uniform electric or magnetic field can be created so that the electric and magnetic field solutions are decoupled, as required in quasi-static analysis. Calculations of the induced fields and currents in a human model exposed to power line frequency fields provide a realistic example of an application of this novel FDTD technique.¿ 1997 American Geophysical Union