(SCF) approach. This simplest of wavefunctions
determines many properties with surprising accuracy and recovers over
99% of the total energy.

 <P>The computational kernel of the SCF algorithm is construction of the
Fock matrix <I>F</I> by contracting the integrals <I>(ij|kl)</I>
against the density matrix <I>D</I> 
<PRE>
        F   = h   + (1/2) SUM [ 2 (ij|kl) - (ik|jl) ] D
         ij    ij         kl                           kl
</PRE>

 <P> Symmetry and sparsity in the matrix of integrals <I>(ij|kl)</