Finding multiple roots of nonlinear algebraic equations using s-system methodology

Michael A. Savageau

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

The solution of nonlinear algebraic equations is a well-known problem in many fields of science and engineering. Several types of numerical methods exist, each with their advantages and disadvantages. S-system methodology provides a novel approach to this problem. The result, as judged by extensive numerical tests, is an effective method for finding the positive real roots of nonlinear algebraic equations. The motivation for this method begins with the observation that any system of nonlinear equations composed of elementary functions can be recast in a canonical nonlinear form known as an S-system. When the derivatives of these ordinary first-order nonlinear differential equations are zero, the resulting nonlinear algebraic equations can be represented as an underdetermined set of linear algebraic equations in the logarithms of the variables. These linear equations, together with a set of simple nonlinear constraints, determine the multiple steady-state solutions of the original system that are positive real. The numerical solution is obtained by approximating the constraints as S-system equations in steady state to yield a determined set of linear algebraic equations, which then can be solved iteratively. This method has been implemented and the experience to date suggests that the method is robust and efficient. The rate of convergence to the solutions is quadratic. A combinatorial method for selecting initial conditions has led to the identification of several, and in many cases all, of the positive real solutions for the set of problems tested. It does so without resorting to analysis in the complex domain, and without having to make random or problem-specific provisions for initializing the procedure. There is an obvious parallel implementation of the algorithm, and there are additional generalizations and efficiencies yet to be realized.

Original languageEnglish (US)
Pages (from-to)187-199
Number of pages13
JournalApplied Mathematics and Computation
Volume55
Issue number2-3
DOIs
StatePublished - 1993
Externally publishedYes

Fingerprint

Multiple Roots
Nonlinear Algebraic Equations
Nonlinear equations
S-system
Linear equations
Methodology
Linear equation
Algebraic Equation
Real Roots
Nonlinear Constraints
Elementary Functions
System of Nonlinear Equations
Steady-state Solution
Numerical methods
First order differential equation
Parallel Implementation
Differential equations
Logarithm
Nonlinear Differential Equations
Derivatives

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics
  • Numerical Analysis

Cite this

Finding multiple roots of nonlinear algebraic equations using s-system methodology. / Savageau, Michael A.

In: Applied Mathematics and Computation, Vol. 55, No. 2-3, 1993, p. 187-199.

Research output: Contribution to journalArticle

@article{cb2e9a1224cb42b5adf6c6a973b24f53,
title = "Finding multiple roots of nonlinear algebraic equations using s-system methodology",
abstract = "The solution of nonlinear algebraic equations is a well-known problem in many fields of science and engineering. Several types of numerical methods exist, each with their advantages and disadvantages. S-system methodology provides a novel approach to this problem. The result, as judged by extensive numerical tests, is an effective method for finding the positive real roots of nonlinear algebraic equations. The motivation for this method begins with the observation that any system of nonlinear equations composed of elementary functions can be recast in a canonical nonlinear form known as an S-system. When the derivatives of these ordinary first-order nonlinear differential equations are zero, the resulting nonlinear algebraic equations can be represented as an underdetermined set of linear algebraic equations in the logarithms of the variables. These linear equations, together with a set of simple nonlinear constraints, determine the multiple steady-state solutions of the original system that are positive real. The numerical solution is obtained by approximating the constraints as S-system equations in steady state to yield a determined set of linear algebraic equations, which then can be solved iteratively. This method has been implemented and the experience to date suggests that the method is robust and efficient. The rate of convergence to the solutions is quadratic. A combinatorial method for selecting initial conditions has led to the identification of several, and in many cases all, of the positive real solutions for the set of problems tested. It does so without resorting to analysis in the complex domain, and without having to make random or problem-specific provisions for initializing the procedure. There is an obvious parallel implementation of the algorithm, and there are additional generalizations and efficiencies yet to be realized.",
author = "Savageau, {Michael A.}",
year = "1993",
doi = "10.1016/0096-3003(93)90020-F",
language = "English (US)",
volume = "55",
pages = "187--199",
journal = "Applied Mathematics and Computation",
issn = "0096-3003",
publisher = "Elsevier Inc.",
number = "2-3",

}

TY - JOUR

T1 - Finding multiple roots of nonlinear algebraic equations using s-system methodology

AU - Savageau, Michael A.

PY - 1993

Y1 - 1993

N2 - The solution of nonlinear algebraic equations is a well-known problem in many fields of science and engineering. Several types of numerical methods exist, each with their advantages and disadvantages. S-system methodology provides a novel approach to this problem. The result, as judged by extensive numerical tests, is an effective method for finding the positive real roots of nonlinear algebraic equations. The motivation for this method begins with the observation that any system of nonlinear equations composed of elementary functions can be recast in a canonical nonlinear form known as an S-system. When the derivatives of these ordinary first-order nonlinear differential equations are zero, the resulting nonlinear algebraic equations can be represented as an underdetermined set of linear algebraic equations in the logarithms of the variables. These linear equations, together with a set of simple nonlinear constraints, determine the multiple steady-state solutions of the original system that are positive real. The numerical solution is obtained by approximating the constraints as S-system equations in steady state to yield a determined set of linear algebraic equations, which then can be solved iteratively. This method has been implemented and the experience to date suggests that the method is robust and efficient. The rate of convergence to the solutions is quadratic. A combinatorial method for selecting initial conditions has led to the identification of several, and in many cases all, of the positive real solutions for the set of problems tested. It does so without resorting to analysis in the complex domain, and without having to make random or problem-specific provisions for initializing the procedure. There is an obvious parallel implementation of the algorithm, and there are additional generalizations and efficiencies yet to be realized.

AB - The solution of nonlinear algebraic equations is a well-known problem in many fields of science and engineering. Several types of numerical methods exist, each with their advantages and disadvantages. S-system methodology provides a novel approach to this problem. The result, as judged by extensive numerical tests, is an effective method for finding the positive real roots of nonlinear algebraic equations. The motivation for this method begins with the observation that any system of nonlinear equations composed of elementary functions can be recast in a canonical nonlinear form known as an S-system. When the derivatives of these ordinary first-order nonlinear differential equations are zero, the resulting nonlinear algebraic equations can be represented as an underdetermined set of linear algebraic equations in the logarithms of the variables. These linear equations, together with a set of simple nonlinear constraints, determine the multiple steady-state solutions of the original system that are positive real. The numerical solution is obtained by approximating the constraints as S-system equations in steady state to yield a determined set of linear algebraic equations, which then can be solved iteratively. This method has been implemented and the experience to date suggests that the method is robust and efficient. The rate of convergence to the solutions is quadratic. A combinatorial method for selecting initial conditions has led to the identification of several, and in many cases all, of the positive real solutions for the set of problems tested. It does so without resorting to analysis in the complex domain, and without having to make random or problem-specific provisions for initializing the procedure. There is an obvious parallel implementation of the algorithm, and there are additional generalizations and efficiencies yet to be realized.

UR - http://www.scopus.com/inward/record.url?scp=0001591382&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001591382&partnerID=8YFLogxK

U2 - 10.1016/0096-3003(93)90020-F

DO - 10.1016/0096-3003(93)90020-F

M3 - Article

AN - SCOPUS:0001591382

VL - 55

SP - 187

EP - 199

JO - Applied Mathematics and Computation

JF - Applied Mathematics and Computation

SN - 0096-3003

IS - 2-3

ER -