Employment of exact numerical quadratures in the evaluation of matrix elements involving highly accurate wavefunctions of helium (and its isoelectronic congeners) generated with the help of the regularized Krylov sequences of Nakatsuji results in an efficient algorithm for the calculation of natural orbitals and the corresponding natural amplitudes {λnl}. The results of such calculations are presented for the total of 600 natural orbitals pertaining to the ground state of the helium atom. The benchmark-quality values of {λnl} computed for 1 ≤ n ≤ 100 and 0 ≤ l ≤ 5 reveal gross inaccuracies in the previously published data. In particular, the dependence of λnl on n is found to follow very closely a simple power-scaling law λnl≈-Al (n+Bl)-4 with Al that, contrary to previous claims, varies only weakly with l. Even more importantly, the numerical trends observed in the present calculations strongly suggest that in the case of the ground state of the helium atom, the only positive-valued natural amplitude is that pertaining to the strongly occupied orbital, i.e., λ10. The relevance of this finding to the existence of unoccupied natural orbitals pertaining to the ground state wavefunction of the H2 molecule is discussed.