N2-time-dependent SCF scheme

A new procedure for the Fock matrix operator construction is proposed. Its application for RHF calculations on diatomic molecules using Slater orbital basis sets shows that the computation time for the new SCF procedure is proportional to the square of the basis set size.