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