Department of Mathematical Sciences, Federal University of Technology, Nigeria
Received Date: February 27, 2017; Accepted Date: April 21, 2017; Published Date: April 29, 2017
Citation: Kayode SJ, Obarhua FO (2017) Symmetric 2-Step 4-Point Hybrid Method for the Solution of General Third Order Differential Equations. J Appl Computat Math 6: 348. doi: 10.4172/2168-9679.1000348
Copyright: © 2017 Kayode SJ, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Applied & Computational Mathematics
This research considers a symmetric hybrid continuous linear multistep method for the solution of general third order ordinary differential equations. The method is generated by interpolation and collocation approach using a combination of power series and exponential function as basis function. The approximate basis function is interpolated at both grid and off-grid points but the collocation of the differential function is only at the grid points. The derived method was found to be symmetric, consistent, zero stable and of order six with low error constant. Accuracy of the method was confirmed by implementing the method on linear and non-linear test problems. The results show better performance over known existing methods solved with the same third order problems. AMS 2010 Subject Classification: 65D05; 65L05; 65L06
Symmetric; Hybrid method; Power series and exponential function; Continuous predictor-corrector method
In this paper, we considered the solution of initial value problems for general third order ordinary differential equations of the form
Where t, y, f∈R.
The numerical and theoretical studies of eqn. (1) have appeared in literature severally. The direct approach for solving this type of ordinary differential equations have been studied and appeared in different literatures [1-7]. This direct approach has demonstrated advantages over the popular approach (reduction to system of first order approach) in terms of speed and accuracy [8,9]. Many authors have focused on direct solution of general second order ivps of odes of the form
Majid et al.  proposed two point four step direct implicit block method for the solution of second order system of ordinary differential equations (ODEs), using variable step size. The method estimated the solutions of initial value problems at two points simultaneously by using four backward steps but with lower order of accuracies. Akinfenwa  presented ninth order hybrid block integrator for solving second order ordinary differential equations. In the paper, the proposed block integrator discretizes the problem using the main and the additional methods to generate system of equations. The resulting system was solved simultaneously in a block-by-block fashion but the order of accuracies is low compare to the order of the method. The authors came up with direct implementation of predictor-corrector methods [3,4,7]. The authors emphasized the need to develop the same order of accuracy of the main predictors and that of the correctors to ensure good accuracy of the method. The order of accuracies in these works improved significantly compare to the existing methods with lower order of their main predictors.
Attempts have also been made by these scholars [12-14,6,7,15]. Olabode  proposed a 5-step block scheme for the solution of special type of eqn. (1). The order of accuracy in Olabode et al.  improves more than that of Olabode et al . Awoyemi et al. , developed a four-point implicit method for the numerical integration of third order ODEs using power series polynomial function . Kuboye and Omar  proposed numerical solution of third order ordinary differential equations using a seven-step block method to improve on Awoyemi et al.  and Olabode  which are of lower order of accuracy. Furthermore, a symmetric hybrid linear multistep method of order six having two off-step points for the solution of eqn. (1) directly was presented by Obarhua and Kayode .
To improve on the study of Obarhua and Kayode  a symmetric of two-step four-point hybrid method for the solution of third order initial value problems of ordinary differential equations directly is therefore proposed using the combination of power series and exponential function as the approximate basis function .
This research work considers the derivation of 2-step 4-point hybrid method for the solution of general third order initial value problems of ordinary differential equations. The approach is to solve eqn. (1) directly without reducing it to a system of first order differential equation. A combination of power series and exponential function is used as the basis function for eqn. (1). The approximate solution eqn. (1) and the resulting differential systems are respectively given as
Where r and s are the number of interpolation and collocation points respectively.
The third derivative of eqn. (3) as compared with eqn. (1) gives
Collocating eqn. (4) at only the grid points, tn+j , j=0(1)2, and interpolating (3) at both grid and off-grid points, , leads to the following system of equations.
where t=[λ0...λ8]T; b=[yn…fn+2]T
ψ=j (j–1) (j–2) as j=5,6.
Solving eqn. (5) for and substituting back into eqn. (3), with some manipulation yields, a linear multistep method with continuous coefficients in the form:
Taking k=2, the coefficients αj(t) and βj(t) are expressed as function of as follows:
Evaluating eqn. (7) at v=1 gives the discrete method:
The first and second derivatives of eqn. (8) are:
Applying the truncation error formula in Awoyemi et al. , associated with eqn. (6) by the difference operator eqn. (9) to determine the order and error constant of the methods:
Where y(x) is assumed to be continuously differentiable of high order. Therefore, expanding eqn. (11) in Taylor’s series and comparing the coefficient of h to give the expression
The linear operator L and the associated methods are said to be of order P if c0=c1=c2=…=cp=…=cp+2=0, cp+3≠0 cp+3. is equal to the error constant. For the purpose of this work, expanding methods (8), (9) and (10) in Taylor’s series and comparing the coefficient of h gives both methods of order p=6 and error constant, cp+3=3.633772 × 10-6, cp+3=1.1804754 × 10-3 and cp+3=1.029201 × 10-4 respectively.
Equations (8-10) are of order six, symmetric, consistent, low error constants and capable of handling oscillatory problems.
To implement the implicit linear 2-step 4-point discrete scheme eqn. (8) and its first and second derivatives eqns.(9) and (10), respectively, the following symmetric explicit schemes and their derivatives are also developed by the same procedure for the evaluation of.
p=6 cp+3=3.107634 × 10-7.
p=6, cp+3=9.0 × 10-6.
p=6, cp+3=1.35 × 10-4.
The methods eqns. (13), (14) and (15) are of order p=6 and error constant, cp+3=3.107634 × 10-7, cp+3=9.0 × 10-6 and cp+3=1.35 × 10-4 respectively.
Other explicit schemes were also generated to evaluate other starting values and Taylor’s series was used to evaluate the values for
Three third order problems out of which one is linear and two are non-linear with exact solutions are solved with our method to test the effectiveness and its accuracy.
Theoretical solution: y(x)=2+2x2–ex.
Table 1 shows the maximum absolute error of our predictorcorrector method and that of Olabode  block method for Problem 1. It reveals that the new method performed creditably well than that of Olabode  of higher order.
|x||yexact||ycomputed||Error in Olabode, (2009), p=7, k=5||Error in new scheme, p=6, k=2|
Table 1: The numerical solution of our methods of order 6 compared with the method of Olabode, (2009), of order 7.
Theoretical solution:. In Table 2, y-exact, the y-computed, the errors of the new method and the time(s) of iteration for Problem 2 are shown.
|x||yexact||ycomputed||Error in new scheme, p=6, k=2||Time(s)|
Table 2: Numerical solution and errors for problem 2.
Table 3 shows the y-exact, the y-computed, and the errors of the new method and the time of iteration for Problem 3.
|x||yexact||ycomputed||Error in new scheme, p=6, k=2||Time(s)|
Table 3: Numerical solution and errors for problem 3.
This paper has produced 2-step 4-point hybrid method for direct solution of higher order ordinary differential equations. The method developed is symmetric, consistent and convergent which can handle oscillatory type of problems. The numerical tests results obtained were compared with block method of Olabode  which was found to perform favorably than the existing method.
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals