Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2017, Article ID 8510948, 8 pages
https://doi.org/10.1155/2017/8510948
Research Article

Implicit One-Step Block Hybrid Third-Derivative Method for the Direct Solution of Initial Value Problems of Second-Order Ordinary Differential Equations

Department of Mathematics, School of Quantitative Sciences, College of Arts and Sciences, Universiti Utara Malaysia, Sintok, Malaysia

Correspondence should be addressed to Mohammad Alkasassbeh; moc.oohay@hebssasakfm

Received 25 August 2016; Revised 20 October 2016; Accepted 23 October 2016; Published 18 January 2017

Academic Editor: Mehmet Sezer

Copyright © 2017 Mohammad Alkasassbeh and Zurni Omar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

A new one-step block method with generalized three hybrid points for solving initial value problems of second-order ordinary differential equations directly is proposed. In deriving this method, a power series approximate function is interpolated at while its second and third derivatives are collocated at all points in the given interval. The proposed method is then tested on initial value problems of second-order ordinary differential equations solved by other methods previously. The numerical results confirm the superiority of the new method to the existing methods in terms of accuracy.

1. Introduction

Numerous problems such as chemical kinetics, orbital dynamics, circuit and control theory, and Newton’s second law applications involve second-order ordinary differential equations (ODEs). Normally, those equations have no analytical solutions. To approximate the solution of such problems several numerical methods were developed on the hands of many scholars such as [13].

Block methods for solving ODEs were first proposed by Milne ([4]). Later [5] adopted Milne’s methods to provide starting values for predictor-corrector scheme. However, the block methods have some drawbacks and this led to the introduction of hybrid methods. According to [6], hybrid methods were initially introduced to overcome zero-stability barrier that occurred in block methods in Dahlquists ([7]). Besides the ability to change step size, the other benefit of these methods is utilizing data off-step points which contribute to the accuracy of the methods.

To increase the accuracy of the numerical methods further, researchers such as [8, 9] proposed high method derivative to overcome stiffness in ODEs. The former presented another type of hybrid methods called second-derivative methods, while the later proposed a Simpson’s type second-derivative method for the solution of a stiff system of first-order IVPs. These scholars motivated us to develop a new generalized three-hybrid one-step third-derivative implicit method for solving second-order ODEs directly using the approach of interpolation and collocation for the general use to improve the efficiency of the approximate solution.

This article is organized as follows: in the coming section we demonstrate the derivation of the method, where we consider three off-step points through the approach of interpolation and collocation. The details of the analysis of the method are discussed in Section 3 which include zero stability, order, consistency, and convergence. In Section 4 some numerical problems are solved and the performance of the developed method is compared with other methods mentioned in literature. Finally, the conclusion is discussed in Section 5.

2. Development of the Method

An approximate power series basis function taking the formwhere and are the number of interpolation and collocation points, respectively, is considered to be a solution to the following ODE:On derivation of (1) twice and thrice we obtainInterpolating (1) at , and collocating (3) at all points , , where , a system of equations in matrix form is produced as below:whereUsing matrix manipulation to solve (4) for the unknown coefficients s and then substituting them back into (1) yieldwhere , is the constant step size for the partition of the interval which is given by , , , and are undetermined constants listed in Appendix  I in Supplementary Material available online at https://doi.org/10.1155/2017/8510948, , and whose first partial derivative isEvaluating (6) at the noninterpolating points and (7) at all points , , produces the following general equations in block form:where is an identity matrix andTheir entries are listed in Appendix  II in Supplementary Material, while the vectors are defined as follows:

3. Analysis of the Method

3.1. Zero Stability

Definition 1. The hybrid block method formula (8) is said to be zero stable if no root of the first characteristic equation has modulus greater than one; that is, and if then the multiplicity of must not exceed two.

To illustrate that the root of the first characteristic equation satisfies the prior definition we assume that and hence which imply that As a result, the developed method is zero stable.

3.2. Order of the Method

The linear operator associated with the hybrid block methods formula (8) is defined as Expanding the above equation in Taylor series and combining like terms we wind up withAccording to [6, 10] method (8) is said to be of order ifThe term is called the error constant and the local truncation error is given by Comparing like terms of and in (14) produces the coefficients with vector of error constants wherewhich conclude that the order of the developed method is .

3.3. Consistency

Definition 2. A block method is said to be consistent if its order is greater than one.

Consistency property is achieved for the hybrid block method from the above analysis since the order .

3.4. Convergence

Theorem 3 (see [16]). Consistency and zero stability are sufficient conditions for a linear multistep method to be convergent

The hybrid block method equation (8) is convergent since it fulfills both the consistency and zero-stability conditions.

4. Numerical Examples

In this section, the efficiency and the performance of the general three-hybrid one-step implicit hybrid block method (8) with order is investigated on five test problems. The first example is highly stiff linear IVP problem with step size , the second is nonlinear IVP with , the third is linear with , the fourth is a nonlinear system with , and finally the fifth is a nonlinear undamped Duffing equation with . It is worth mentioning that this method works even for large interval and different values of step size. The values mentioned in this article are chosen just for the sake of comparison with the existing methods only.

Problem 4. , , , .

Exact Solution. with

Source (see [11]). See Table 1.

Table 1: Comparison of the proposed method with [11].

Problem 5. , , , .

Exact Solution. with

Source (see [12]). See Table 2.

Table 2: Comparison of the proposed method with [12].

Problem 6. , , , .

Exact Solution. with

Source (see [13]). See Table 3.

Table 3: Comparison of the proposed method with [13].

Problem 7 (a system of nonlinear IVP). , , , .
, , , .

Exact Solutions. , , with , , ; see Tables 4 and 5.

Table 4: Exact and approximate solutions for solving using the developed method.
Table 5: Exact and approximate solutions for solving using the developed method.

Problem 8 (the nonlinear undamped Duffing equation). , , ,

Exact Solutions. , where , , , , , , , and . , , ; see Table 6.

Table 6: Comparison of the proposed method with [15].

5. Conclusion

A general three-hybrid one-step block method of order has been proposed for the direct solution of general second-order ODEs. The developed method is tested on five different problems. Numerical analysis shows that the developed method is consistent and zero stable which conclude its convergence. The computed results are then compared with the results of existing methods in terms of error by considering different values of , , and . The new method is found to have superiority over them as shown in Tables 16.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Universiti Utara Malaysia for the financial support of the publication of this article.

References

  1. M. M. Chawla and S. R. Sharma, “Families of three-stage third order Runge-Kutta-Nystrom methods for y′′=fx,y,y,” The Journal of the Australian Mathematical Society. Series B. Applied Mathematics, vol. 26, no. 3, pp. 375–386, 1985. View at Publisher · View at Google Scholar
  2. E. Hairer and G. Wanner, “A theory for Nyström methods,” Numerische Mathematik, vol. 25, no. 4, pp. 383–400, 1976. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  3. J. Vigo-Aguiar and H. Ramos, “Variable stepsize implementation of multistep methods for y = f(x,y,y′),” Journal of Computational and Applied Mathematics, vol. 192, no. 1, pp. 114–131, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  4. W. E. Milne, Numerical Solution of Differential Equations, John Wiley & Sons, New York, NY, USA, 1953. View at MathSciNet
  5. R. J. Barkley, “A Runge-kutta for all seasons,” SIAM Review, vol. 9, no. 3, pp. 417–452, 1967. View at Publisher · View at Google Scholar · View at MathSciNet
  6. J. D. Lambert, Computational Methods in Ordinary Differential Equations, John Wiley & Sons, New York, NY, USA, 1973. View at MathSciNet
  7. G. Dahlquist, “Convergence and stability in the numerical integration of ordinary differential equations,” Mathematica Scandinavica, vol. 4, pp. 33–53, 1956. View at Publisher · View at Google Scholar · View at MathSciNet
  8. F. F. Ngwane and S. N. Jator, “Block hybrid-second derivative method for stiff systems,” International Journal of Pure and Applied Mathematics, vol. 80, no. 4, pp. 543–559, 2012. View at Google Scholar · View at Scopus
  9. R. K. Sahi, S. N. Jator, and N. A. Khan, “A Simpson's-type second derivative method for stiff systems,” International Journal of Pure and Applied Mathematics, vol. 81, no. 4, pp. 619–633, 2012. View at Google Scholar · View at Scopus
  10. S. O. Fatunla, Numerical Methods for Fatunla Initial Value Problems in Ordinary Differential Equations, Academic Press, New York, NY, USA, 1988.
  11. O. Adeniran and B. Ogundare, “An efficient hybrid numerical scheme for solving general second order initial value problems (IVPs),” International Journal of Applied Mathematical Research, vol. 4, no. 2, pp. 411–419, 2015. View at Publisher · View at Google Scholar
  12. A. A. James, A. O. Adesanya, and S. Joshua, “Continuous block method for the solution of second order initial value problems of ordinary differential equation,” International Journal of Pure and Applied Mathematics, vol. 83, no. 3, pp. 405–416, 2013. View at Publisher · View at Google Scholar · View at Scopus
  13. O. O. Olanegan, D. O. Awoyemi, B. G. Ogunware, and F. O. Obarhua, “Continuous double-hybrid point method for the solution of second order ordinary differential equations,” International Journal of Advanced Scientific and Technical Research, vol. 5, no. 2, pp. 549–562, 2015. View at Google Scholar
  14. R. F. Abdelrahim and Z. Omar, “Direct solution of second-order ordinary differential equation using a single-step hybrid block method of order five,” Mathematical and Computational Applications, vol. 21, no. 2, pp. 1–12, 2016. View at Publisher · View at Google Scholar
  15. S. N. Jator, “Block third derivative method based on trigonometric polynomials for periodic initial-value problems,” Afrika Matematika, vol. 27, no. 3-4, pp. 365–377, 2016. View at Publisher · View at Google Scholar · View at MathSciNet
  16. P. Henrici, Discrete Variable Methods in ODEs, John Wiley & Sons, New York, NY, USA, 1962.