The application of sparse matrix methods to the numerical solution of nonlinear elliptic partial differential equations