During inversion we will try to round elememnts to zero.
inv(x, eps = ZERO_EPS)
x
a polynomial matrix
eps
zero threshold
Right now only matrices with numerical determinant are supported.