We present a fast and robust iterative method for obtaining self-consistent solutions to the coupled system of Schrödinger's and Poisson's equations in quantum structures. A simple expression describing the dependence of the quantum electron density on the electrostatic potential is used to implement a predictor – corrector type iteration scheme for the solution of the coupled system of differential equations. This approach simplifies the software implementation of the nonlinear problem, and provides excellent convergence speed and stability. We demonstrate the algorithm by presenting an example for the calculation ofthe two-dimensional bound electron states within the cross-section of a GaAs-AlGaAs based quantum wire. For this example, six times fewer iterations are needed when our predictor – corrector approach is applied, compared to a corresponding underrelaxation algorithm.