Least change secant update strategies are derived that are compatible with algebraic systems of nonlinear equations arising when Galerkin procedures are applied to nonlinear elliptic boundary value problems. The resulting quasi-Newton procedures for solving the algebraic system of equations are observed to result in significant computational savings in experiments involving boundary value problems of one or two spatial variables.