A theoretically attractive and computationally fast algorithm is presented for the determination of the coefficients of the determinantal polynomial and the coefficients of the adjoint polynomial matrix of a given three-dimensional (3-D) state space model of Fornasini-Marchesini type. The algorithm uses the discrete Fourier transform (DFT) and can be easily implemented on a digital computer.