The Symbolic Computation of Integrability Structures for Partial Differential Equations