MNsmithForm(A)

Computes a Smith Form of any polynomial matrix of the size m*n

Function computes a Smith normal form of a polynomia matrix of an arbitrary size n by m. Standard Matlab function 'smithForm' only square matrices, which is rather useless for many linear system related applications.
Besides the Smith Form, the function also provides a vector of Invariant factors and the vector of greatest common dividers of the i-th order minors of the matrix A.