A mathematical model has been developed for the time-dependent heat conduction process in a region of arbitrary geometry. An alternating direction implicit (ADI) scheme has been worked out for the computer solution of the problem in three space dimensions. The time and nodal propagation of the errors are depicted graphically. Also included in the computer software is a subroutine for plotting the isotherms in various geometries.