The MatrixGCRD(A, B, x) command computes a right matrix GCD of two matrices of polynomials. Both input matrices of polynomials can be square or rectangular, but must have the same number of columns. The entries are either univariate polynomials in x over the field of rational numbers Q, or rational expressions over Q (that is, univariate polynomials in x with coefficients in Q(a1,...,an)).

It is required that the matrix of polynomials  have full column rank.

The fourth argument is optional. It is used to return a list  of matrices. These matrices satisfy the matrix linear diophantine equation U1 . A + U2 . B = MatrixGCRD(A,B).

The fifth argument is optional. It is used to return a list  of matrices. These matrices satisfy the matrix linear diophantine equation . This is useful. If B is square and nonsingular this gives  a left-reduced matrix rational function.

The MatrixGCLD(A, B, x) command returns a left matrix GCD with optional arguments for returning solutions for matrix linear diophantine equations and right-reduced matrix rational functions.  It is required that the matrix of polynomials  have full row rank.

The method used is a fraction-free algorithm by Beckermann and Labahn that computes a matrix GCD using Mahler systems.

Beckermann, B., and Labahn, G. "Fraction-free Computation of Matrix Rational Interpolants and Matrix GCDs." SIAM Journal on Matrix Analysis and Applications. Vol. 22 No. 1, (2000): 114-144.