An iterative method for solving problems in neutron transport theory is described. The method uses the result of a low-order spherical harmonies calculation as the initial function and attempts to improve this with one iteration. The results for the one-velocity Boltzmann equation with plane symmetry show that good accuracy is obtained for many practical cases.