Three-dimensional induced polarization and resistivity modeling for buried electrode configurations can be achieved by adapting surface integral techniques for surface electrode configurations to buried electrodes. Modification of. the surface technique is accomplished by considering the additional mathematical terms required to express-the changes in the electrical potential and geometry caused by placing the source and receiver electrodes below the surface. This report presents a listing of a computer program to calculate the resistivity and induced polarization response from a three-dimensional body for buried electrode configurations. The program is designed to calculate the response for the following electrode configurations: (1) hole-to-surface array with a buried bipole source and a surface bipole receiver, (2) hole-to-surface array with a buried pole source and a surface bipole receiver, (3) hole-to-hole array with a buried, fixed pole source and a moving bipole receiver, (4) surface-to-hole array with a fixed pole source on the surface and a moving bipole receiver in the borehole, (5) hole-to-hole array with a buried, fixed bipole source and a buried, moving bipole receiver, (6) hole-to-hole array with a buried, moving bipole source and a buried, moving bipole receiver, and (7) single-hole, buried bipole-bipole array. Input and output examples are given for each of the arrays.