The transmission-line-matrix method is a time-domain numerical method for solving wave problems. The method uses a mesh of transmission lines to represent a propagation space, and the losses in the space are accounted for by making the transmission lines lossy. Lossy boundaries are simulated by imperfect boundary reflections on the transmission lines. A FORTRAN program implementing this technique is presented.