On some 4-Point Spline Collocation Methods for Solving Ordinary Initial Value Problems

This paper discusses collocation schemes based on seventh C 3 -splines with four collocation points s_{i-1+\alpha}=x_{i-1}+\alpha h , x_{i-1+\beta}=x_{i-1}+\beta h , x_{i-1+\theta}=x_{i-1}+\theta h and x_i=x_{i-1}+h in each subinterval [ x i @ 1 , x i ], i = 1(1) N for solving initial value problems in ordinary differential equations including stiff equations. Here 0 < f < g < è < 1 are arbitrarily given. It is shown that the methods are convergent and the order of convergence is seven if: \eqalign {&\beta(\theta -\theta^2)+\beta^2(\theta^2-\theta)+\alpha^2[\theta^2-\theta +\beta^2(1+4\theta -5\theta^2)+\theta(4\theta^2-4\theta -1)]\cr &\qquad \qquad \qquad +\alpha[\theta -\theta^2+\beta (1+4\theta -4\theta^2)+\beta^2(4\theta^2+4\theta -1)]\le 0 and they are unstable if f , g , è < 0.7279115. Moreover, the absolute stability properties of the methods are considered. It shows that with 0.888035 h f < g < è < 1 the methods are A-stable while they are not if f h 0.5; on other hand, the sizes of regions of absolute stability increase when f , g , è M 1 m.

[1]  F. R. Loscalzo,et al.  Spline Function Approximations for Solutions of Ordinary Differential Equations , 1967 .

[2]  S. Sallam,et al.  Numerical solution of general n th-order differential equations via splines , 1990 .

[3]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[4]  H. M. El-Hawary,et al.  A deficient spline function approximation to systems of first-order differential equations: part 2 , 1983 .

[5]  C. D. Boor,et al.  Collocation at Gaussian Points , 1973 .

[6]  H. Meyer,et al.  Modified backward differentiation methods of the Adams-type based on exponential interpolation , 1991 .

[7]  Syvert P. Norsett Splines and Collocation for Ordinary Initial Value Problems , 1984 .

[8]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[9]  B. Hulme Discrete Galerkin and related one-step methods for ordinary differential equations , 1972 .

[10]  H. De Meyer,et al.  Families of backward differentiation methods based on a new type of mixed interpolation , 1990 .

[11]  M. Kerimov,et al.  Modern numerical methods for ordinary differential equations , 1980 .

[12]  S. Sallam,et al.  A quartic C 3 -spline collocation method for solving second-order initial value problems , 1996 .

[13]  F. V. Niekerk Rational one-step methods for initial value problems , 1988 .

[14]  Christoph W. Ueberhuber,et al.  Iterated defect correction for the efficient solution of stiff systems of ordinary differential equations , 1977 .

[15]  J. L. Walsh,et al.  The theory of splines and their applications , 1969 .

[16]  J. C. Butcher,et al.  Parallelism in the Numerical Integration of Initial Value Problems. , 1993 .

[17]  A. Salama HIGHER ORDER AND OPTIMAL SCHEMES FOR STIFF INITIAL-VALUE PROBLEMS , 2000 .