Solution of a two-point boundary value model of immobilized enzyme reactions, using an S-system-based root-finding method

Fumihide Shiraishi, Eberhard O. Voit

Research output: Contribution to journalArticle

2 Citations (Scopus)


The analysis of immobilized enzyme reactions leads to a two-point boundary value problem that requires numerical solutions based on root-finding methods. The standard Newton-Raphson (N-R) algorithm is adequate if the root-finding is initiated close to the true solution; in other cases, large numbers of iterations may be needed before the algorithm converges. An alternative is a root-finding method based on the theory of S-systems. This method turns out to be much more robust and efficient than the N-R algorithm. It converges quadratically, with fewer iterations, and in a wider range of relevant initial values. An investigation of a simpler algebraic equation with a similar mathematical structure demonstrates that the rate of convergence in the S-system method is not very sensitive to the distance between the root and initial value, whereas convergence of the N-R algorithm slows down drastically with increasing distance.

Original languageEnglish
Pages (from-to)289-310
Number of pages22
JournalApplied Mathematics and Computation
Issue number2-3
Publication statusPublished - Apr 15 2002
Externally publishedYes


All Science Journal Classification (ASJC) codes

  • Computational Mathematics
  • Applied Mathematics

Cite this