Numerical integrators for stiff and highly oscillatory differential equations
Some L-stable fourth-order explicit one-step numerical integration formulas which requiare no matrix inversion are proposed to cope effictively with systems of ordinary differential equations with large Lipschitz constants (including those having highly oscillatory solutions). The implicit integration procedure proposed in Fatunla (11) is further developed to handle a larger class of stiff systems as well as those with highly oscillatory solutions. The same pair of nonlinear equations as in (11) is solved for the stiffness/oscillatory parameters. However, the nonlinear systems are transformed into linear forms and an efficient computational procedure is developed to obtain these parameters. The new schemes compare favorably with the backward differentiation formula (DIFSUB) of Gear (13), (14), and the blended linear multistep methods of Skeel and Kong (24), and the symmetric multistep methods of Lambert and Watson (17).
- DOE Contract Number:
- EY-76-S-02-2383
- OSTI ID:
- 5346806
- Journal Information:
- Math. Comput.; (United States), Vol. 34:154
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multi-derivative numerical methods for the solution of stiff ordinary differential equations
Blended linear multistep methods. [Blended DIFSUB, in FORTRAN IV]