  We solve a special type of linear systems with coefficients in multivariate
polynomial rings. These systems arise in the computation of parametric
Bernstein-Sato polynomials associated with certain hypergeometric ideals in the
Weyl algebra.
