The use of a linear combination of Gaussian-type orbitals (CGTO), instead of an individual Gaussiantype orbital (GTO), as a unit of basis functions for large-scale molecular calculations, is discussed. A systematic construction of the CGTO basis functions is attempted and the results for the atoms from Li through Ar are reported.