htrianr

htrianr — triangularization of polynomial matrix

Calling sequence

[Ar,U,rk]=htrianr(A)

`A`	: polynomial matrix
`Ar`	: polynomial matrix
`U`	: unimodular polynomial matrix
`rk`	: integer, normal rank of `A`

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)