The performance of a nonstationary iterative scheme for solving the large system of linear equations arising in neutron diffusion calculations is reported. The method is compared with successive overrelaxation and conjugate gradient methods. Practical aspects of its implementation on computers are considered.