The primary application of expand is to distribute products over sums. This is done for all polynomials. For quotients of polynomials, only sums in the numerator are expanded; products and powers are left alone.

The expand(M) function computes the element-wise expansion of M.  The result, R, is formed as R[i,j] = expand(M[i,j]).