A Block Procedure with Linear Multi-Step Methods Using Legendre Polynomials for Solving ODEs

Abstract

In this article, we derive a block procedure for some K-step linear multi-step methods (for K = 1, 2 and 3), using Legendre polynomials as the basis functions. We give discrete methods used in block and implement it for solving the non-stiff initial value problems, being the continuous interpolant derived and collocated at grid and off-grid points. Numerical examples of ordinary differential equations (ODEs) are solved using the proposed methods to show the validity and the accuracy of the introduced algorithms. A comparison with fourth-order Runge-Kutta method is given. The ob-tained numerical results reveal that the proposed method is efficient.

Share and Cite:

Abualnaja, K. (2015) A Block Procedure with Linear Multi-Step Methods Using Legendre Polynomials for Solving ODEs. Applied Mathematics, 6, 717-723. doi: 10.4236/am.2015.64067.

1. Introduction

Numerous problems in Physics, Chemistry, Biology and Engineering science are modeled mathematically by ordinary differential equations (ODEs), e.g., series circuits, mechanical systems with several springs attached in series lead to a system of differential equations [1] . Since most realistic differential equations do not have exact analytic solutions, therefore, approximate and numerical techniques [2] -[14] are used extensively. Recently introduced power series method [15] has been used for solving a wide range of problems. This method has proven rather successful in dealing with both linear as well as nonlinear problems, as it yields analytical solutions and offers certain advantages over standard numerical methods. It is free from rounding off errors since it does not involve discretization, and is computationally inexpensive. Ercan and Mustafa [16] have applied this method to a system of differential-algebraic equations.

In this paper, we present an efficient numerical method to solve numerically the ODEs. The proposed method is a block procedure for some K-step linear multi-step methods (for K = 1, 2 and 3) [17] using Legendre polynomials as the basis functions. In addition, we give discrete methods used in block and implement it for solving the non-stiff IVPs that were the continuous interpolant derived and collocated at grid and off-grid points. In this article, we consider the general form of the first order initial value problem

. (1)

The plan of the paper is as follows: In Section 2, the derivation of the proposed methods is presented. In Section 3, stability and convergence analysis of the block schemes is given. In Section 4, two numerical examples are considered. The paper ends with summary and conclusions in Section 5.

2. The Derivation of the Proposed Methods

In this section, we derive discrete methods to solve (1) at a sequence of nodal points where is the step-length or grid-size defined by and denotes the true solution to (1) while the approximate solution is denoted by, for some positive number N. The proposed method depends on the perturbed collocation method with respect to the power series method is used as the basis for collocation approximation with the Legendre polynomials as the perturbation term.

In the first we consider the approximate solution of the perturbed form of (1) in the following power series [18] -[21]

. (2)

where

. (3)

Substituting from Equation (2) in Equation (1) we have

, (4)

where is the Legendre polynomial of degree K, valid in and is a perturbed parameter. In particular, we shall be dealing with cases K = 1, 2 and 3 in (2) and (3).

The well-known Legendre polynomials are defined on the interval and can be determined with the aid of the following recurrence formula [22] [23]

,

where the first four polynomials are

. (5)

In order to use these polynomials on the interval we define the so called shifted Legendre polynomials by introducing the change of variable

. (6)

Cases Study

Case 1: If K = 1

In this case, we take the polynomial and use (6), then collocate this equation at and, we obtain

.

In addition, from Equation (3) we can deduce that,. Then, the Equation (4) will reduce to the following form

. (7)

We now collocate Equation (7) at, and interpolate (2) at, we get a system of three Equations with and parameter

(8)

Using a suitable method to solve the above system to obtain

From (2), we have

. (9)

Now, the required numerical scheme of the proposed method will be obtained if we collocate the above Equation (9) at and substitute, , as follows

. (10)

Which is the well-known trapezoidal rule.

Case 2: If K = 2

In this case, we take the polynomial and use (6), then collocating this equation at, and, we obtain

In addition, from Equation (3) we can deduce that, ,. Then, the Equation (4) will reduce to the following form

. (11)

We now collocate Equation (11) at, and interpolate (2) at, we get a system of four equations with and parameter

(12)

Using a suitable method to solve the above system to obtain

From (2), we have

. (13)

Now, the required numerical scheme of the proposed method will be obtained if we collocate the above Equation (13) at and substitute as follows

. (14)

Case 3: If K = 3

For case K = 3, we collocate the continuous scheme

, (15)

at grid and off-grid points and and this gives the required block schemes

. (16)

3. Stability and Convergence Analysis of the Block Schemes

In this section, we present a summary on the order, the error constant and the convergence of the proposed block schemes. This summary in given in the following table.

4. Numerical Examples

In this section, we implement the proposed method with K = 2 and K = 3 to solve two first order initial value problems, and compare the obtained numerical results with those obtained from using the fourth-order Runge- Kutta method (RK4).

Example 1.

Consider the following IVP

.

With the exact solution.

The numerical results of this example are presented in Table 1 and Table 2 with the cases K = 2 and K = 3, respectively. In these tables, we presented a comparison the obtained numerical results using the proposed method with the exact solution and those numerical results obtained from using RK4.

Example 2.

Consider the following IVP

With the exact solution.

The numerical results of this example are presented in Figure 1 and Figure 2 with the cases K = 2 and K = 3, respectively. In these figures, we presented a comparison the obtained numerical results using the proposed method with the exact solution and those numerical results obtained from using RK4.

Table 1. A comparison the proposed method at K = 2 with the exact solution and RK4: Example 1.

Table 2. A comparison the proposed method at K = 3 with the exact solution and RK4: Example 1.

Figure 1. A comparison the proposed method at K = 2 with the exact solution and RK4: Example 2.

Figure 2. A comparison the proposed method at K = 3 with the exact solution and RK4: Example 2.

5. Summary and Conclusions

In this paper, we presented three new block-schemes (K = 1, K = 2 and K = 3) that are convergent and absolutely stable. We used the proposed method to solve numerically a wide-range of linear initial value problems. The results of the presented examples show that our method was capable for solving such problems of IVPs and generates the convergence analysis, and closed to their exact solutions. This method is very simple and effective for a wide-range of ODEs. All computations are made by Matlab.

Acknowledgements

We thank the Editor and the referee for their comments.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Hairer, E. and Wanner, G. (1991) Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Springer, Berlin.
http://dx.doi.org/10.1007/978-3-662-09947-6
[2] Atkinson, K.E. (1889) An Introduction to Numerical Analysis. 2nd Edition, John Wiley and Sons, New York.
[3] Khader, M.M. (2013) Numerical Treatment for Solving the Perturbed Fractional PDEs Using Hybrid Techniques. Journal of Computational Physics, 250, 565-573.
http://dx.doi.org/10.1016/j.jcp.2013.05.032
[4] Khader, M.M. (2011) On the Numerical Solutions for the Fractional Diffusion Equation. Communications in Nonlinear Science and Numerical Simulations, 16, 2535-2542.
http://dx.doi.org/10.1016/j.cnsns.2010.09.007
[5] Khader, M.M. (2013) The Use of Generalized Laguerre Polynomials in Spectral Methods for Fractional-Order Delay Differential Equations. Journal of Computational and Nonlinear Dynamics, 8, Article ID: 041018.
[6] Khader, M.M. (2013) An Efficient Approximate Method for Solving Linear Fractional Klein-Gordon Equation Based on the Generalized Laguerre Polynomials. International Journal of Computer Mathematics, 90, 1853-1864.
http://dx.doi.org/10.1080/00207160.2013.764994
[7] Khader, M.M. and Sweilam, N.H. (2013) On the Approximate Solutions for System of Fractional Integro-Differential Equations Using Chebyshev Pseudo-Spectral Method. Applied Mathematical Modelling, 37, 9819-9828.
http://dx.doi.org/10.1016/j.apm.2013.06.010
[8] Khader, M.M. and Hendy, A.S. (2013) A Numerical Technique for Solving Fractional Variational Problems. Mathematical Methods in Applied Sciences, 36, 1281-1289.
http://dx.doi.org/10.1002/mma.2681
[9] Khader, M.M. and Babatin, M.M. (2013) On Approximate Solutions for Fractional Logistic Differential Equation. Mathematical Problems in Engineering, 2013, Article ID: 391901.
[10] Khader, M.M., Sweilam, N.H. and Mahdy, A.M.S. (2013) Numerical Study for the Fractional Differential Equations Generated by Optimization Problem Using Chebyshev Collocation Method and FDM. Applied Mathematics and Information Science, 7, 2011-2018.
http://dx.doi.org/10.12785/amis/070541
[11] Lambert, J.D. (1991) Numerical Methods for ODE. John Wiley and Sons, New York.
[12] Sweilam, N.H. and Khader, M.M. (2010) A Chebyshev Pseudo-Spectral Method for Solving Fractional Integro-Differential Equations. ANZIAM Journal, 51, 464-475.
http://dx.doi.org/10.1017/S1446181110000830
[13] Sweilam, N.H., Khader, M.M. and Nagy, A.M. (2011) Numerical Solution of Two-Sided Space-Fractional Wave Equation Using Finite Difference Method. Journal of Computational and Applied Mathematics, 235, 2832-2841.
http://dx.doi.org/10.1016/j.cam.2010.12.002
[14] Sweilam, N.H., Khader, M.M. and Kota, W.Y. (2013) Numerical and Analytical Study for Fourth-Order Integro-Differential Equations Using a Pseudo-Spectral Method. Mathematical Problems in Engineering, 2013, Article ID: 434753, 7 pages.
[15] Hirayama, H. (2000) Arbitrary Order and A-Stable Numerical Method for Solving Algebraic Ordinary Differential Equation by Power Series. 2nd International Conference on Mathematics and Computers in Physics, Vouliagmeni, Athens, 9-16 July 2000, 1-6.
[16] Çelik, E. and Bayram, M. (2003) On the Numerical Solution of Differential-Algebraic Equations by Padé Series. Applied Mathematics and Computation, 137, 151-160.
http://dx.doi.org/10.1016/S0096-3003(02)00093-0
[17] Onumanyi, P., Awoyemi, D.O., Jator, S.N. and Sirisena, U.W. (1994) New Linear Multistep Methods with Continuous Coefficients for First Order Initial Value Problems. Journal of the Nigerian Mathematical Society, 13, 7-51.
[18] Fatokun, J.O., Onumanyi, P. and Sirisena, U. (1999) A Multistep Collocation Based on Exponential Basis for Stiff Initial Value Problems. Nigerian Journal of Mathematics and Applications, 12, 207-223.
[19] Fatokun, J.O., Aimufua, G.I.O. and Ajibola, I.K.O. (2010) An Efficient Direct Collocation Method for the Integration of General Second Order Initial Value Problem. Journal of Institute of Mathematics & Computer Sciences, 21, 327-337.
[20] Fatunla, S.O. (1998) Numerical Method for Initial Value Problems in ODEs. Academic Press Inc., New York.
[21] Funaro, D. (1992) Polynomial Approximation of Differential Equations. Springer Verlag, New York.
[22] Bell, W.W. (1968) Special Functions for Scientists and Engineers. Butler and Tanner Ltd, Frome and London.
[23] Canuto, C., Hussaini, M.Y., Quarteroni, A. and Zang, T.A. (2006) Spectral Methods. Springer-Verlag, New York.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.