Abstract
This paper presents a class of (p + 2)-step backward differentiation formulas of orderp. The two extra degrees of freedom obtained by limiting the order of a (p + 2)-step formula top are used to extend the region of absolute stability. A new formula of orderp has a region of absolute stability very similar to that of a classical backward differentiation formula of orderp - 1 forp being in the range 4–6. The backward differentiation formulas with extended regions of absolute stability are constructed by appending two exponential-trigonometric terms to the polynomial basis of the classical formulas. Besides the absolute stability, the paper discusses relative stability and contractivity. The principles of an experimental implementation of the new formulas are outlined, and a linear problem integrated with this computer program indicates that the extended regions of absolute stability can actually be exploited in practice.
Similar content being viewed by others
References
J. D. Lambert,Computational Methods in Ordinary Differential Equations, John Wiley & Sons, New York, 1973.
G. K. Gupta,Some new high-order multistep formulae for solving stiff equations, Math. Comp., vol. 30 (1976), pp. 417–432.
R. W. Klopfenstein,Numerical differentiation formulas for stiff systems of ordinary differential equations, RCA Review, Vol. 32 (1971), pp. 447–462.
J. M. Varah,Stiffly stable linear multistep methods of extended order, SIAM J. Numer. Anal., vol. 15 (1978), pp. 1234–1246.
Stig Skelboe,Numerical methods for the analysis of mildly nonlinear electrical circuits, Ph. D. Thesis, Report No. 75-08, Oct. 1975, Inst. for Numerical Analysis, Tech. University of Denmark.
T. Bar-David,Exponentially fitted linear multi-value methods, Manuscript, Lockheed Missiles and Space Company, Sunnyvale, California.
K. Madsen,Minimax solution of nonlinear equations without calculating derivatives, Mathematical Programming Study 3 (1975), pp. 110–126.
R. K. Brayton and C. C. Conley,Some results on the stability and instability of the backward differentiation methods with non-uniform time steps, IBM Research, July 28, 1972, RC 3964.
O. Nevanlinna and W. Liniger,Contractive methods for stiff differential equations, Part I, BIT, vol. 18 (1978), pp. 457–474.
G. Dahlquist,Some contractivity questions for one-leg- and linear multistep methods, Report TRITA-NA-7905, Dept. of Numerical Analysis and Computer Science, The Royal Institute of Technology, Stockholm, Sweden.
R. K. Brayton and C. H. Tong,Stability of dynamical systems: A constructive approach, IEEE Trans. Circuits and Systems, vol. CAS-26 (1979), pp. 224–234.
S. Skelboe,Implementation of Chebyshevian linear multistep formulas, BIT, vol. 20 (1980), pp. 356–366.
W. M. G. van Bokhoven,Linear implicit differentiation formulas of variable step and order, IEEE Trans. Circuits and Systems, CAS-22 (1975), pp. 109–115.
B. Christensen,Numerical solution of stiff systems of differential equations (Master's thesis in Danish), Institute for Numerical Analysis, Technical University of Denmark, Lyngby, Denmark, 1980.
S. Skelboe,The control of order and steplength for backward differentiation methods, BIT, vol. 17 (1977), pp. 91–107.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Skelboe, S., Christensen, B. Backward differentiation formulas with extended regions of absolute stability. BIT 21, 221–231 (1981). https://doi.org/10.1007/BF01933167
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF01933167