High-order P-stable multistep methods

From MaRDI portal





A class of two-step methods for solving the initial value problem \(y=f(t,y)\) is presented. The stability polynomial is directly related to the (m,m)-diagonal Padé approximation to the exponential through the \(\Sigma\)-map, see the reviewer and \textit{O. Nevanlinna} [SIAM J. Numer. Anal. 20, 1210-1218 (1983; Zbl 0532.65065)]. Hence the schemes are P- stable, of order 2m and implicit. If the implicit equations are solved by a Newton like algorithm at each iteration step the linear system consists of a matrix polynomial in the Jacobian \(f_ y\) of degree m. The powers of the Jacobian \(f_ y\) are avoided by factoring this polynomial. For the linear test equation with an oscillatory forcing term it is shown that the schemes are always in phase in the sense of \textit{I. Gladwell}, \textit{R. M. Thomas} [Int. J. Numer. Methods Eng. 19, 495-503 (1983; Zbl 0513.65053)]. A numerical comparison with other schemes is presented for linear and nonlinear equations.




Cited in
(only showing first 100 items - show all)








This page was built for publication: High-order P-stable multistep methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q804254)