Abstract An easily implementable iteration scheme guaranteeing convergence under very lenient conditions is developed and implemented into the MNDO-program package. Calculations done using large atomic clusters obtained from modelling diamond and silicon surfaces show that this scheme is, with respect to the requirements of computer time and memory, superior to some of the other iteration methods successfully used in quantum chemical calculations.