Computational Solutions of Matrix Problems Over an Integral Domain