The computation of fourth virial coefficients for pairwise-additive spherically symmetric interaction potentials