An algorithm to calculate the numerical values of the Clebsch-Gordan (
C-G) coefficients is described. It uses a well elaborated recursive pr
ocedure yielding significantly high computation efficiency and accurac
y. A 60-line Fortran program that implements the algorithm can be appl
ied to the calculation of various symmetry-adapted C-G coefficients.