The normal form derived by Mansour [I] for single-input / single-output, time-invariant, linear discrete systems will be extended to the multivariable case. An algorithm will be used to achieve simple computation of this form starting with the Luenberger first canonical form. The construction of this form will be illustrated by examples applied to model reduction of multivariable discrete-time ...