Accession Number : AD0787588

Title :   A Stable Variant of the Secant Method for Solving Nonlinear Equations.

Descriptive Note : Interim rept.,

Corporate Author : CARNEGIE-MELLON UNIV PITTSBURGH PA

Personal Author(s) : Gragg,W. B. ; Stewart,G. W.

Report Date : APR 1974

Pagination or Media Count : 59

Abstract : The usual successive secant method for solving systems of nonlinear equations suffers from two kinds of instabilities. First the formulas used to update the current approximation to the inverse Jacobian are numerically unstable. Second, the directions of search for a solution may collapse into a proper affine subspace, resulting at best in slowed convergence and at worst in complete failure of the algorithm. In this report it is shown how the numerical instabilities can be avoided by working with factorizations of matrices appearing in the algorithm. Moreover, these factorizations can be used to detect and remedy degeneracies among the directions. A second part of this report documents and lists a program implementing the algorithm described in the first part. (Author)

Descriptors :   *Nonlinear algebraic equations, *Iterations, Matrices(Mathematics), Computations, Algorithms, Computer programs

Subject Categories : Theoretical Mathematics

Distribution Statement : APPROVED FOR PUBLIC RELEASE