EFFICIENT FOURTH ORDER P-STABLE FORMULAE.

R. M. Thomas

    Research output: Contribution to journalArticlepeer-review

    Abstract

    A family of fourth order of P-stable methods for solving second order initial value problems is considered. When applied to a nonlinear differential system, all the methods in the family give rise to a nonlinear system which may be solved using a modified Newton method. The classical methods of this type involve at least three (new) function evaluations per iteration (that is, they are 3-stage methods) and most involve using complex arithmetic in factorizing their iteration matrix. We derive methods which require only two (new) function evaluations per iteration and for which the iteration matrix is a true real perfect square. This implies that real arithmetic will be used and that at most one real matrix must be factorized at each step. Also we consider various computational aspects such as local error estimation and a strategy for changing the step size.
    Original languageEnglish
    Pages (from-to)599-614
    Number of pages15
    JournalBIT
    Volume27
    Issue number4
    Publication statusPublished - 1987

    Fingerprint

    Dive into the research topics of 'EFFICIENT FOURTH ORDER P-STABLE FORMULAE.'. Together they form a unique fingerprint.

    Cite this