Numerical treatment of convolution integrals involving distributions with densities having singularities at the origin