An efficient algorithm for calculation of the Luenberger canonical form

An algorithm is suggested to obtain the Luenberger canonical form for multivariable systems. The method computes the canonical form directly without having to compute the transformation matrix. In addition, there is a large reduction in the number of calculations. The reduced computations along with Gaussian techniques lend to greater inherent accuracy and the ability to refine the solution with additional computations.