Language:   Search:   Contact
Zentralblatt MATH has released its new interface!
For an improved author identification, see the new author database of ZBMATH.

Query:
Fill in the form and click »Search«...
Format:
Display: entries per page entries
Zbl 0916.65081
Vigo-Aguiar, Jesus; Ferrándiz, José M.
A general procedure for the adaptation of multistep algorithms to the integration of oscillatory problems.
(English)
[J] SIAM J. Numer. Anal. 35, No.4, 1684-1708 (1998); erratum ibid. 35, No.5, 2116 (1998). ISSN 0036-1429; ISSN 1095-7170/e

Consider the scalar initial value problem (IPV) in the closed finite interval [a, b] as follows: $$P(D)y(x)=f(x,y), \quad D={d\over{dx}}, \quad y(a)=y_o,\dots,y^{(m-1)}(a)=y_o^{(m-1)},\tag 1$$ \noindent where $f$ is a continuous function with continuous derivatives till a certain order $m+k$ and satisfying all conditions for the existence and uniqueness of solutions. $P(D)$ is the operator $(D-\lambda_1)\dots(D-\lambda_m)$, $\lambda_1,\dots,\lambda_m$ are real or complex conjugate constants. The paper presents a general procedure allowing to construct multistep methods with coefficients depending on the step for the exact integration of both homogeneous linear problem and nonhomegenous IVP (1) without truncation error. $f$ is supposed to be a certain Fourier polynomial in one or more frequencies resulting in oscillator solution. Such an algorithm is denoted as adapted to (1). It combines the integration of both the linear part and a particular solution of the complete problem. \par The study of methods adapted to linear equations is presented in a unified way especially with respect to consistency and convergence. The results are valid for all presented examples and different existing methods. The adaptation of a multistep code to the exact solution of a forced oscillator in one frequency includes a new procedure for the computation of coefficients. The exact integration of oscillators in two or more frequencies which are confluent or non-confluent include general formulae for the calculation of the coefficients and the local truncation error. Numerical examples illustrate efficiency of the derived results. Some of the ideas are related to the adaptation of the Chebyshev methods without multistep character by {\it J. Panovsky} and {\it D. L. Richardson} [J. Comput. Appl. Math. 23, 35-51 (1988; Zbl 0649.65048)].
[Lubomír Bakule (Praha)]
MSC 2000:
*65L06 Multistep, Runge-Kutta, and extrapolation methods
65L05 Initial value problems for ODE (numerical methods)
34A34 Nonlinear ODE and systems, general
34C10 Qualitative theory of oscillations of ODE: Zeros, etc.

Keywords: consistency; multistep methods; convergence; oscillatory problems; forced oscillator; numerical examples; Chebyshev methods

Citations: Zbl 0649.65048

Highlights
Master Server