A Class of One-Step Hybrid Third Derivative Block Method for the Direct Solution of Initial Value Problems of Second-Order Ordinary Differential Equations
Author Affiliations
1Department of Mathematical Science, Adamawa State University, Mubi-Nigeria
2Department of Mathematics and Statistics, Federal University, Wukari-Nigeria
Corresponding Author
Raymond D, Department of Mathematics and Statistics, Federal University, Wukari-Nigeria
Citation
© 2018 Skwame Y, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 international License.
Abstract
In this paper, we consider the development of a class of one-step hybrid third derivative
block method with three off-grid points for the direct solution of initial value problems
of second order Ordinary Differential Equations. We adopted method of interpolation and
collocation of power series approximate solution to generate the continuous hybrid linear
multistep method, which was evaluated at grid points to give a continuous block method.
The discrete block method was recovered when the continuous block method was evaluated
at selected grid points. The basic properties of the method was investigated and was found
to be zero-stable, consistent and convergent. The efficiency of the method was tested on
some stiff equations and was found to give better approximation than the existing method,
which we compared our result with.
Keywords
One-Step; Hybrid Block Method; Third Derivative; Stiff Odes; Collocation and
Interpolation Method.
Introduction
This paper solves second order initial value problems in the form
where f is continuous within the interval of integration. Solving higher order derivatives method by reducing them to a system of first-order approach involves more functions to evaluate which then leads to a computational burden as in [1-3]. The setbacks of this approach had been reported by scholars, among them are Bun and Vasil’yer and Awoyemi et al [4,5].
The method of collocation and interpolation of the power series approximation to generate continuous linear multistep method has been adopted by many scholars; among them are Fatunla, Awoyemi, Olabode, Vigor Aquilar and Ramos, Adeniran et al, Abdelrahim et al, Mohammad et al to mention a few [5-11]. Block method generates independent solution at selected grid points without overlapping. It is less expensive in terms of number of function evaluation compared to predictor corrector method, moreover it possess the properties of Runge Kutta method for being self-starting and does not require starting values. Some of the authors that proposed block method are: [12-17].
In this paper, we developed a on e-step hybrid third derivative method with three offgrid, which is implemented in block method. The method is self-starting and does not require starting values or predictors. The implementation of the method is cheaper than the predictor-corrector method. This method harnesses the properties of hybrid and third derivative, this makes it e¢cient for sti¤ problems.
The paper is organized as follows: In section 2, we discuss the methods and the materials
for the development of the method. Section 3 considers analysis of the basis properties of
the method which include convergence and stability region, numerical experiments where
the e¢ciency of the derived method is tested on some numerical examples and discussion of
results. Lastly, we concluded in section 4.
Derivation of the Method
We consider a power series approximate solution of the form
where r=2 and s=5 are the numbers of interpolation and collocation points respectively, is considered t o be a solution to (1).
The second and third derivative of (2) gives
Collocating (4) at all points and Interpolating
Equation (2) at
, gives a system of non linear equation
of the form AX = U (5)
Where
AY L = BR1 +CR2 +DR3+ ER4+ GR5 (8)
Analysis of Basic Properties of the Method Order of the Block
According to [7] the order of the new method in Equation (8) is obtained by using the Taylor series and it is found that the developed method has a uniformly order eleven, with an error constants vector of:
Consistency
The hybrid block method [4] is said to be consistent if it has an order more than or equal to one.
Therefore, our method is consistent.
Zero Stability of Our Method
Definition: A block method is said to be zero-stable if as , the root
of the first characteristic polynomial
that is
Satisfies and for those roots with
, multiplicity
must not exceed two. The block method for k=1, with three off grid
collocation point expressed in the form
Hence, our method is zero-stable.
Regions of Absolute Stability (RAS)
Using Mat Lab package, we were able to plot the stability region of the block methods (see Figure 1 below).
The stability polynomial for K=1 with three offstep point using
Scientific Workplace software package we have in the following one:
Using Mat Lab software, the absolute stability region of the new method is plotted and shown in Figure 1.
Numerical Example
Problem I: We consider a highly stiff problem
Problem II:
Exact Solution:
Figure 1: Absolute Stability Region
Conclusion
It is evident from the above tables that our proposed methods are indeed accurate, and can handle stiff equations. Also in terms of stability analysis, the method is A-stable.
Comparing the new method with the existing method, the result
presented in the Tables 1 and 2 shows that the new method performs
better than the existing method, and even the order of new method
is higher than the order of the existing method [11,18]. In this article,
a one-step block method with three off-step points is derived via the
interpolation and collocation approach. The developed method is
consistent, A-stable, convergent, with a region of absolute stability
and order Ten.
Table 1: Comparison of the proposed method
Table 2: Comparison of the proposed new method
References
- Z Omar, M Suleiman.Parallel r-point implicit block method forsolving higher order ordinarydifferential equations directly. JICT , 2004;3(1):53-66.
- SJ kayoed, O
Adeyeye. A 3-step hybrid method for direct solutionof second order initial
value problems. Aust J Basic Appl Sci.2011;5(12):2121-2126.
- A James, A
Adesanya, S Joshua. Continuous block method for the solution of second order
initial value problems of ordinary differential equation. Int J Pure Appl Math.
2013;83(3):405-416.
- RA Bun, TD Vasil’s. A numerical method for solving
Differential equations of any Order Comp. Math Phys. 1992;32(3):317-330.
- DO Awoyemi, AO Adesanya, SN Ogunyebi. Construction of Self
Starting Numeral Method for the Solution of initial value Problem of General
Second Ordinary Differential Equation. Journ Num math. 2008;4(2):267-278.
- SO Fatunla. Block methods
for second order. Int J Comput Maths.2007;41(1-2):55-63.
- BT Olabode. An
accurate scheme by block method for the thirdorder ordinary differential
equation. Pacific journal of scienceand technology. 2009;10(1).
- J Vigo-Aguilar, H Ramos. Variable step size implementation of
multistep methods for ( ) 1 y ' '= f x, y, y . J Comput Appl Math.
2006;192(1):114–131.
- AO Adeniran, BS
Ogundare. An efficient hybrid numericalscheme for solving general second order
Initial Value Problems(IVPs). International Journal of Applied Mathematical
Research.2015;4(2):411-419.
- R Abdelrahim, O
Zurni. Direct Solution of Second-Order OrdinaryDifferential Equation Using a
Single-Step Hybrid Block Methodof Order Five. Mathematical and Computational
Applications.2016;21(2):12.
- A Mohammad, O Zurni.Implicit One-Step Block Hybrid Third Derivative Method for the Direct Solutionof Initial ValueProblems of Second-Order Ordinary DifferentialEquations.Hindawi Journal of Applied Mathematics. 2017.
- SN Jator. A sixorder linear multi step method for the directsolution of ( ) 1 y ' '= f x, y, y. Intern J Pure and Applied of Maths.2007;40(4):457-472.
- Jator SN, Li J. A
self-starting linear multistep method for directsolution general second order
initial value problems. Intern JComput Math. 2009;86(5):827-836.
- ZA Majid, MB Sulaiman, Z Omar. 3
points implicit block methodfor solving ordinary differential equations, Bull.
MalaysiaMathematical Science Soc. 2006;29(1) 23-31.
- DO Awoyemi, EA
Adebile, AO Adesanya, TA Anake. Modified blockfor the direct solution for the
second order ordinary differentialequations. Intern J Applied Math and Comp.
2011;3(3):181-188.
- Anake TA, Awoyeni
DO, Adesanya AO. A one step method for thesolution of general second order
ordinary differential equations.Intern J Sci and Engin. 2012;2(4):159-163.
- Adesanya AO, Odekule MR, Alkali MA. Three steps block predictor-block Corrector method for solution of general second order of ordinary differential equations. Intern J Engin Research and Application. 2012;2(4):2297-2301.
- M Alkasassbeh, Z
Omar. Generalized two-hybrid One-Step ImplicitThird derivatives block method
for the Direct Solution of SecondOrder Ordinary Differential Equations.
International Journal ofPure and Applied Mathematics.
2017;112(3):497-517.