The symbolic computation of integrability structures for partial differential equations