Symbolic Computation for Integro-Differential-Time-Delay Operators with Matrix Coefficients