A method of computing wave propagation. The region of interest is broken up into a collection of small cubic elements, each of which contains information about the fields and material parameters at that point. Signals are applied to appropriate points within the region and the method proceeds by calculating the effects of these at successive small time steps (hence the term time domain). The method of determining the field variation between spatial and temporal steps is based on mathematical difference techniques applied to finite steps (hence finite difference).