A Fortran program to calculate the matrix elements of the Coulomb interaction involving hydrogenic wave functions