A method for the numerical solution of a system of coupled, nonlinear elliptic partial differential equations is described, and the application of the method to the equations governing steady, laminar natural convection ...