htrianr - triangularization of polynomial matrix
triangularization of polynomial matrix A.
A is [m,n] , m <= n.
Ar=A*U
Warning: there is an elimination of "small" terms (see function code).
x=poly(0,'x'); M=[x;x^2;2+x^3]*[1,x-2,x^4]; [Mu,U,rk]=htrianr(M) det(U) M*U(:,1:2)