Two-Step Hybrid Block Method with Generalized Two Off-Step Points Within Each Step for Solving Second Order Initial Value Problems

Two-Step Hybrid Block Method with Generalized Two Off-Step Points Within Each Step for Solving Second Order Initial Value Problems

Kamarun H. Mansor Oluwaseun Adeyeye Zurni Omar*

School of Quantitative Sciences, Universiti Utara Malaysia, Sintok 06010, Kedah, Malaysia

Corresponding Author Email:
18 October 2022
10 January 2023
25 January 2023
Available online: 
28 April 2023
| Citation

© 2023 IIETA. This article is published by IIETA and is licensed under the CC BY 4.0 license (



This article develops a new two-step hybrid block method for the numerical solution of second order initial value problems with better accuracy. Two off-step points are introduced in generalized form and the resulting block method is developed using interpolation approach with interpolation points at one on-step point and one off-step point. The hybrid points are given in a generalized form to give room for flexibility of the choice of hybrid points which will give more information on which points produces the best solutions. The resultant order seven block method obtained satisfied all basic properties such as order, zero stability, consistency and convergence and produces better accuracy than existing numerical methods for solving second-order initial value problems. Thus, justifying the adoptability of the new block method for solving second-order initial value problems.


two-step, hybrid block method, generalized off-step points, second order initial value problems

1. Introduction

This article focuses on solving second order initial value problems (IVPs) in the form:

$y^{\prime \prime}=f\left(x, y, y^{\prime}\right),\ y(a)=\eta_0,\ y^{\prime}(a)=\eta_1$     (1)

using hybrid block methods with generalized off-step points.

IVPs are a class of ordinary differential equations (ODEs) where a dependent variable and one independent variable are related, and the order being second order results from the value of the highest derivative in the differential equation [1-3]. Numerical approaches are adopted to solve differential equations generally when the exact solution does not exist or is difficult to obtain [4], and for IVPs in the form of Eq. (1), the conventional implementation procedure to solve using a numerical approach was its reduction from a second order IVP to a system of two first order IVPs. This reduction approach is discussed by Awoyemi [5] and said to be saddled with the disadvantage of high computational time and low accuracy. This led to another implementation procedure using predictor-corrector procedure to directly solve IVPs. Studies [6-10] utilised the predictor-corrector procedure. Shokri [6] developed a two-step block method, Panopoulos et al. [7] developed an eight-step symmetric embedded predictor-corrector method, Ndanusa and Tafida [8] introduced predictor-corrector methods of high order, Biasa et al. [9] obtained a predictor-corrector block iteration method, and a predictor-corrector linear multistep method was developed in Kayode and Adegboro [10]. Although, the predictor-corrector procedure was a direct procedure instead of using the concept of reduction, further studies discovered that block methods gave better solutions that predictor-corrector because the order of the predictors are less than the correctors and this impacts the accuracy negatively [11, 12].

Block methods were first proposed by Milne [13] and Sarafyan [14] to obtain starting values for predictor-corrector methods. Adesanya [15] described the formulation of block methods as evaluation of a multistep method at different grid and/or off-points to generate a system of numerical schemes that can be applied to produce approximate solutions at the grid and/or off-grid points simultaneously.

Block methods developed using only grid points are seen in studies [16-21] but the consideration of off-grid points in addition have been seen to result in better accuracy. The use of both grid and off-points in block methods led to the name hybrid block methods. Some studies that have developed hybrid block methods to obtain solutions for models in the form of Eq. (1) include Olukunle and Felix [22] where a one-step hybrid scheme was developed, Shokri et al. [23] which proposed a symmetric two-step semi-hybrid scheme, Ehigie et al. [24] with a generalised two-Step continuous linear multistep method of hybrid type, Fasasi et al. [25] developing a one-step continuous hybrid block method, and Gebremedhin and Jena [26] where the approximate solution of certain ODEs were obtained using a hybrid block approach.

Specifically, this article considers an extension of the previous study by Mansor et al. [27] where a two-step hybrid block method with one generalized off-step point was proposed. In their work, the authors considered one off-step point within each step to find the direct solution of the second order IVPs. Although, the results obtained had good accuracy, there is still room for improvement as only one off-step point between two-step intervals was considered in deriving the method. Therefore, an improvement is introduced in this article by introducing two off-step points within each step for a two-step block method (See Figure 1). The aim of introducing more off-step points is to increase the order of the method, thus, improving the accuracy of the solution. The off-step points are also introduced in a generalised form, such that various suitable values within each interval can be selected as off-step points.

Therefore, this article derives a new block method using interpolation and collocation approach as discussed in the next section. Section 3 tests the basic properties of a numerical method for the developed method, Section 4 shows the results obtained from selected numerical examples, and Section 5 concludes the article.

2. Development of the Two-Step Hybrid Block Method with Generalized Two Off-Step Points within Each Step

The following power series polynomial:

$y(x)=\sum_{j=0}^{i+c-1} a_j\left(\frac{x-x_n}{h}\right)^j$     (2)

is used as an approximate solution of Eq. (1) where $x \in\left(x_n, x_{n+2}\right]$ for n=0, 2, 4, ..., N-2 with $h=x_{\delta+1}-x_\delta$ and $\delta=0,1,2, \ldots, \mathrm{N}$ in the interval [a, b]. i and c denote the number of interpolation and collocation points respectively.

The interpolation-collocation strategy is illustrated in Figure 1 below where two off-step points within each step are denoted by p, q, r, and s for 0<p<q<1<r<s<2.

Figure 1. Interpolation and collocation points for two-step hybrid block method with generalized two off-step points within each step

Referring to Figure 1, i=2 and c=7. Substituting these values in Eq. (2) gives:

$y(x)=\sum_{j=0}^8 a_j\left(\frac{x-x_n}{h}\right)^j$     (3)

whose second derivative is:

$y^{\prime \prime}(x)=f\left(x, y, y^{\prime}\right)=\sum_{j=2}^8 a_j \frac{j(j-1)}{h^2}\left(\frac{x-x_n}{h}\right)^{j-2}$     (4)

Interpolating Eq. (2) at $x_{n+\sigma}(\sigma=1, r)$ yields:

$y_{n+\sigma}=\sum_{j=0}^8 a_j\left(\frac{x-x_n}{h}\right)^j$     (5)

and collocating Eq. (4) at all points $x_{n+\sigma}(\sigma=0, p, q, 1, r, s, 2)$ produces:

$f_{n+\sigma}=\sum_{j=2}^8 a_j \frac{j(j-1)}{h^2}\left(\frac{x-x_n}{h}\right)^{j-2}$     (6)

which form nine simultaneous equations involving the unknown coefficients $a_j^{\prime} s$. The values of $a_j^{\prime} s$ are obtained using Gaussian elimination method and then substituted back in Eq. (2) to produce a continuous implicit hybrid two-step scheme with generalised two off-step points as below:

$\begin{aligned} y(x)= & \sum_{j=0, p} \alpha_j(x) y_{n+j} \\ & +\sum_{j=0}^3 \beta_j(x) f_{n+j}+\sum_{j=p, q, r} \beta_j(x) f_{n+j} \; .\end{aligned}$     (7)

which is then evaluated at the non-interpolating points, $x_{n+\mu}\ (\mu=0, p, q, s, 2)$.

Differentiating Eq. (7) once, we get:

$\begin{aligned} y^{\prime}(x) & =\sum_{j=0, \; p} \frac{\partial}{\partial x} \alpha_j \; (x) y_{n+j} \\ & +\sum_{j=0}^3 \frac{\partial}{\partial x} \beta_j(x) f_{n+j}+\sum_{j=p, q, r} \frac{\partial}{\partial x} \beta_j(x) f_{n+j} \; .\end{aligned}$    (8)

Evaluating Eq. (8) at xn leads to the formation of a discrete implicit hybrid two-step scheme which can be represented in matrix form as:

$\begin{aligned} A^{2[2]_2} Y^{2[2]_2} & =B_1^{2[2]_2} R_1^{2[2]_2}+B_2^{2[2]_2} R_2^{2[2]_2} \\ & +h^2\left[D^{2[2]_2} R_3^{2[2]_2}+E^{2[2]_2} R_4^{2[2]_2}\right] .\end{aligned}$    (9)

Multiplying Eq. (9) by the inverse of $A^{2[4]_2}$ yields:

$\begin{aligned} I_{6 \times 6} Y^{2[2]_2}= & \bar{B}_1^{2[2]_2} R_1^{2[2]_2}+\bar{B}_2^{2[2]_2} R_2^{2[2]_2} \\ & +h^2\left[\bar{D}^{2[2]_2} R_3^{2[2]_2}+\bar{E}^{2[2]_2} R_4^{2[2]_2}\right] .\end{aligned}$    (10)

which is equivalent to:

$\begin{aligned} y_{n+p}= & y_n+h p y_n^{\prime}+h^2\left[\bar{D}_{11} f_n+\bar{E}_{11} f_{n+p}+\bar{E}_{12} f_{n+q}\right. \\ & \left.+\bar{E}_{13} f_{n+1}+\bar{E}_{14} f_{n+r}+\bar{E}_{15} f_{n+s}+\bar{E}_{16} f_{n+2}\right] \\ y_{n+q}= & y_n+h q y^{\prime}{ }_n+h^2\left[\bar{D}_{21} f_n+\bar{E}_{21} f_{n+p}+\bar{E}_{22} f_{n+q}\right. \\ & \left.+\bar{E}_{23} f_{n+1}+\bar{E}_{24} f_{n+r}+\bar{E}_{25} f_{n+s}+\bar{E}_{26} f_{n+2}\right] \\ y_{n+1}= & y_n+h y_n^{\prime}+h^2\left[\bar{D}_{31} f_n+\bar{E}_{31} f_{n+p}+\bar{E}_{32} f_{n+q}\right. \\ & \left.+\bar{E}_{33} f_{n+1}+\bar{E}_{34} f_{n+r}+\bar{E}_{35} f_{n+s}+\bar{E}_{36} f_{n+2}\right] \\ y_{n+r}= & y_n+h r y_n^{\prime}+h^2\left[\bar{D}_{41} f_n+\bar{E}_{41} f_{n+p}+\bar{E}_{42} f_{n+q}\right. \\ & \left.+\bar{E}_{43} f_{n+1}+\bar{E}_{44} f_{n+r}+\bar{E}_{45} f_{n+s}+\bar{E}_{46} f_{n+2}\right] \\ y_{n+s}= & y_n+h s y_n^{\prime}+h^2\left[\bar{D}_{51} f_n+\bar{E}_{51} f_{n+p}+\bar{E}_{52} f_{n+q}\right. \\ & \left.+\bar{E}_{53} f_{n+1}+\bar{E}_{54} f_{n+r}+\bar{E}_{55} f_{n+s}+\bar{E}_{56} f_{n+2}\right] \\ y_{n+2}= & y_n+2 h y_n^{\prime}+h^2\left[\bar{D}_{61} f_n+\bar{E}_{61} f_{n+p}+\bar{E}_{62} f_{n+q}\right. \\ & \left.+\bar{E}_{63} f_{n+1}+\bar{E}_{64} f_{n+r}+\bar{E}_{65} f_{n+s}+\bar{E}_{66} f_{n+2}\right] .\end{aligned}$

Combining Eq. (8) at $x_{n+\mu}(\mu=p, q, 1, r, s, 2)$ and Eq. (10) produces the following first derivative of the main block:

$\begin{aligned} y_{n+p}^{\prime} & =y_n^{\prime}+h\left[D_{11}^{\prime} f_n+E_{11}^{\prime} f_{n+p}+E_{12}^{\prime} f_{n+q}\right. \\ & \left.+E_{13}^{\prime} f_{n+1}+E_{14}^{\prime} f_{n+r}+E_{15}^{\prime} f_{n+s}+E_{16}^{\prime} f_{n+2}\right] \\ y_{n+q}^{\prime} & =y_n^{\prime}+h\left[D_{21}^{\prime} f_n+E_{21}^{\prime} f_{n+p}+E_{22}^{\prime} f_{n+q}\right. \\ & \left.+E_{23}^{\prime} f_{n+1}+E_{24}^{\prime} f_{n+r}+E_{25}^{\prime} f_{n+s}+E_{26}^{\prime} f_{n+2}\right] \\ y_{n+1}^{\prime}= & y_n^{\prime}+h\left[D_{31}^{\prime} f_n+E_{31}^{\prime} f_{n+p}+E_{32}^{\prime} f_{n+q}\right. \\ & \left.+E_{33}^{\prime} f_{n+1}+E_{34}^{\prime} f_{n+r}+E_{35}^{\prime} f_{n+s}+E_{36}^{\prime} f_{n+2}\right] \\ y_{n+r}^{\prime}= & y_n^{\prime}+h\left[D_{41}^{\prime} f_n+E_{41}^{\prime} f_{n+p}+E_{42}^{\prime} f_{n+q}\right. \\ & \left.+E_{43}^{\prime} f_{n+1}+E_{44}^{\prime} f_{n+r}+E_{45}^{\prime} f_{n+s}+E_{46}^{\prime} f_{n+2}\right] \\ y_{n+s}^{\prime}= & y_n^{\prime}+h\left[D_{51}^{\prime} f_n+E_{51}^{\prime} f_{n+p}+E_{52}^{\prime} f_{n+q}\right. \\ & \left.+E_{53}^{\prime} f_{n+1}+E_{54}^{\prime} f_{n+r}+E_{55}^{\prime} f_{n+s}+E_{56}^{\prime} f_{n+2}\right] \\ y_{n+2}^{\prime}= & y_n^{\prime}+h\left[D_{61}^{\prime} f_n+E_{61}^{\prime} f_{n+p}+E_{62}^{\prime} f_{n+q}\right. \\ & \left.+E_{63}^{\prime} f_{n+1}+E_{64}^{\prime} f_{n+r}+E_{65}^{\prime} f_{n+s}+E_{66}^{\prime} f_{n+2}\right] .\end{aligned}$

The coefficients $\left\{\bar{D}_{11}, \bar{D}_{21}, \ldots, \bar{D}_{61}\right\},\left\{\bar{E}_{11}, \bar{E}_{12}, \ldots, \bar{E}_{66}\right\}$, $\left\{D_{11}^{\prime}, D_{21}^{\prime}, \ldots, D_{61}^{\prime}\right\}$, and $\left\{E_{11}^{\prime}, E_{12}^{\prime}, \ldots, E_{66}^{\prime}\right\}$ are obtained with respect to the value of the off-step selected (See Appendix).

The next section discusses the properties of the block method in its generalized form, while in Section 4, comparison will be made with the methods in previous studies for some numerical examples. This is to show the advantage of the improvement introduced in this article over the existing studies.

3. Properties of the Method

Order of the Method

The linear difference operator L associated with Eq. (10) is defined as:

$\begin{aligned} L[y(x) ; h] & =I_{6 \times 6} Y^{2[2]_2}-\bar{B}_1^{2[2]_2} R_1^{2[2]_2}+\bar{B}_2^{2[2]_2} R_2^{2[2]_2} \\ & -h^2\left[\bar{D}^{2[2]_2} R_3^{2[2]_2}+\bar{E}^{2[2]_2} R_4^{2[2]_2}\right] .\end{aligned}$     (11)

where, y(x) is an arbitrary test function continuously differentiable on a. Expanding the components of $Y^{2[2]_2}$ and $R_4^{2[2]_2}$ in Taylors series about x=xn and collecting the terms in powers of h gives:

$\begin{aligned} L[y(x) ; h]= & \bar{C}_0^{2[2]_2} y(x)+\bar{C}_1^{\overline{2}_2} h y^{\prime}(x) \\ & +\bar{C}_2^{2[2]_2} h^2 y^{\prime \prime}(x)+\ldots\end{aligned}$     (12)

Definition 3.1 Hybrid block method (10) and associated linear operator in (11) are said to be of order $d$ if $\bar{C}_0^{2[2]_2}=\bar{C}_1^{2[2]_2}=\bar{C}_2^{2[2]_2}=\ldots=\bar{C}_{d+1}^{2[2]_2}=0$ and  $\bar{C}_{d+2}^{2[2]_2} \neq 0$ with error vector constants $\bar{C}_{d+1}^{2[2]_2}$ [28].

Definition 3.2 A linear multistep method is consistent if it has order $\mathrm{d} \geq 1$ [28].

Expanding all terms in Eq. (11) using Taylor series about xn, the order of method is found to be $[7,7,7,7,7,7]^T$ which implies that the new two-step hybrid block method (10) is consistent since its order is greater than 1.

Zero Stability

The new two-step hybrid block method in (10) is zero-stable if no root of the first characteristic polynomial $\pi(\omega)=$ $\left|\omega I_{6 \times 6}-\bar{B}_1^{2[2] 2}\right|$ is having a modulus greater than one and every root of modulus one is simple, where $I_{6 \times 6}$ is the identity matrix and $\bar{B}_1^{2[2]_2}$ is the coefficients matrix of $y_n$ function.

By setting determinant $\pi(\underline{\omega})=0$, then $\pi(\omega)=\omega^5(\omega-\underline{\underline{1}}=0$ which implies $\omega=0,0,0,0,0,1$. Hence, the newly developed method is zero-stable.

Consistency and Convergence

Theorem 3.1 Consistency and zero stability are sufficient conditions for a linear multistep method to be convergent as seen in study [28].

Based on the above theorem, the developed method satisfies the properties for consistency and zero-stability as shown above, and thus is convergent.

4. Numerical Examples

Problem 1

$\begin{aligned}& y^{\prime \prime}+\left(\frac{6}{x}\right) y^{\prime}+\left(\frac{4}{x^2}\right) y=0, y(1)=1, y^{\prime}(1)=1 \text { with } h=\frac{1}{320} \text {. }\end{aligned}$

Exact Solution: $y(x)=\frac{5}{3 x}-\frac{4}{x^2}$.

Source: Studies [29] and [30].

Problem 2

$y^{\prime \prime}-y^{\prime}=0, \underline{y}(0)=1, y^{\prime}(0)=-1 \text { with } h=\frac{1}{10} \text {. }$

Exact Solution: $y(x)=1-e^x$.

Source: Studies [31] and [32].

Problem 3

$y^{\prime \prime}-2 y^3=0, y(1)=1, y^{\prime}(1)=-1$ with $h=\frac{1}{10}$.

Exact Solution: $y(x)=\frac{1}{x}$.

Source: Study [33].

Tables 1, 2, and 3 have shown the results obtained for solving Problems 1, 2, and 3 respectively, which include linear and nonlinear examples, and the comparison of these results with existing studies. For Problem 1, the solution using the newly developed two-step hybrid block was compared with the one-step hybrid block method developed in Abdelrahim [29] and the continuous implicit hybrid one-step method from Anake [30]. The results in Table 1 show the new method having distinctly better accuracy than Abdelrahim [29] and Anake [30]. For Problems 2 and 3, the existing approaches [27, 31-33] whose solutions were compared with the new hybrid block method compared closely to the new hybrid block method, however there was still a clear improvement in the solutions obtained by the new hybrid block method.

Table 1. Comparison of absolute errors obtained by new method with the previous studies for solving Problem 1


Absolute Error


Absolute Error


Absolute Error New Method

$\begin{gathered}p=4 / 5, q=9 / 10 \\ r=101 / 100, q=11 / 10 .\end{gathered}$









































Table 2. Comparison of absolute errors obtained by new method with the previous studies for solving Problem 2


Absolute Error


Absolute Error


Absolute Error New Method

$\begin{aligned} & p=1 / 4, q=3 / 4 \;,\\ & r=5 / 4, q=7 / 4 \;.\end{aligned}$









































Table 3. Comparison of absolute errors obtained by new method with the previous studies for solving Problem 3


Absolute Error


Absolute Error


Absolute Error New Method

$\begin{array}{ll}p=1 / 4, & q=3 / 4 \;,\\ r=5 / 4, & q=7 / 4 \;.\end{array}$









































5. Conclusions

This study has successfully developed a new two-step hybrid block method with two off-step points within each step for directly solving second order initial value problems of ODEs. The effect of introducing additional off-step points was considered to justify the property that the introduction of more off-step points increases order and thus improves accuracy. The new method was tested on some existing second initial value problems available in literature. These problems included both linear and nonlinear IVPs which covers both types of second order IVPs. The numerical results indicate that the new method produces better accuracy than the previous methods when solving the same linear and nonlinear second order IVPs as considered by existing studies.


This research was supported by Ministry of Higher Education (MOHE) of Malaysia through Fundamental Research Grant Scheme (FRGS/1/2017/STG06/UUM/01/1).


$\begin{aligned} & \bar{D}_{16}=\frac{-p^2}{1680 q r s}\left(5 p^5-560 q r s-8 p^4(3+q+r+s)+\right. \\ & 14 p^3(2+3 s+r(3+s)+q(3+r+ \\ & s))-28 p^2(2 s+r(2+3 s)+q(2+3 s+r(3+s)))+ \\ & 70 p(2 r s+q(2 s+r(2+3 s)))) \\ & \bar{D}_{26}=\frac{q^2}{1680 p r s}\left(2 p\left(4 q^4+280 r s-7 q^3(3+r+s)+\right.\right. \\ & 14 q^2(2+3 s+r(3+s))-35 q(2 s+ \\ & r(2+3 s)))+q\left(-5 q^4-140 r s+8 q^3(3+r+s)-\right. \\ & 14 q^2(2+3 s+r(3+s))+28 q(2 s+ \\ & r(2+3 s)))), \\ & \bar{D}_{36}=\frac{1}{1680 p q r s}(11-20 r-20 s+42 r s+q(42 r- \\ & 20+42 s-112 r s)+2 p(21 r-10+ \\ & 21 s-56 r s+7 q(3-8 r-8 s+35 r s))), \\ & \bar{D}_{46}=\frac{r^2}{1680 p q s}\left(-2 p\left(r\left(70 s-4 r^3+7 r^2(3+s)-\right.\right.\right. \\ & 14 r(2+3 s))+7 q\left(r^3-40 s-2 r^2(3+\right. \\ & s)+5 r(2+3 s)))+r\left(r\left(-5 r^3+56 s+8 r^2(3+s)-\right.\right. \\ & 14 r(2+3 s))+2 q\left(4 r^3-70 s-7 r^2(3+s)+14 r(2+\right. \\ & 3 s)))),\end{aligned}$

$\begin{aligned} & \bar{D}_{56}=\frac{s^2}{1680 p q r}\left(2 p\left(7 q\left(r\left(40-15 s+2 s^2\right)-s(10-\right.\right.\right. \\ & \left.\left.6 s+s^2\right)\right)+s\left(-7 r\left(10-6 s+s^2\right)+\right. \\ & \left.\left.s\left(28-21 s+4 s^2\right)\right)\right)+s\left(-2 q\left(s\left(21 s-28-4 s^2\right)+\right.\right. \\ & \left.7 r\left(10-6 s+s^2\right)\right)+s(s(24 s-28- \\ & \left.\left.\left.\left.5 s^2\right)+r\left(56-42 s+8 s^2\right)\right)\right)\right) \\ & \bar{D}_{66}=\frac{2}{105 p q r s}(4(r+s-2)+q(4-7 r s)+p(4-7 r s+ \\ & 7 q(5 r s-r-s))) \\ & \bar{E}_{11}=\frac{p^2}{840\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(15 p^5-280 q r s-\right. \\ & 20 p^4(3+q+r+s)+28 p^3(2+ \\ & 3 s+r(3+s)+q(3+r+s))-42 p^2(2 s+r(2+ \\ & 3 s)+q(2+3 s+r(3+s)))+ \\ & 70 p(2 r s+q(2 s+r(2+3 s)))) \\ & \bar{E}_{12}=\frac{p^4}{840(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}\left(5 p^4+140 r s-\right. \\ & 8 p^3(3+r+s)+14 p^2(2+3 s+r(3+ \\ & s))-28 p(2 s+r(2+3 s))) \\ & \bar{E}_{13}=\frac{p^4}{840(p-1)(q-1)(r-1)(s-1)}\left(5 p^4+140 q r s-8 p^3(2+\right. \\ & q+r+s)+14 p^2(2 s+r(2+s)+ \\ & q(2+r+s))-28 p(2 r s+q(2 s+r(2+s))))\end{aligned}$

$\begin{aligned} & \bar{E}_{14}=\frac{-p^4}{840 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}\left(5 p^4+140 q s-\right. \\ & 8 p^3(3+q+s)+14 p^2(2+3 s+q(3+ \\ & s))-28 p(2 s+q(2+3 s))), \\ & \bar{E}_{15}=\frac{p^4}{840(q-s) s(s-p)(s-r)\left(2-3 s+s^2\right)}\left(5 p^4+140 q r-\right. \\ & 8 p^3(3+q+r)+14 p^2(2+3 r+q(3+ \\ & r))-28 p(2 r+q(2+3 r))), \\ & \bar{E}_{16}=-\frac{p^4}{1680(p-2)(q-2)(r-2)(s-2)}\left(5 p^4+70 q r s-8 p^3(1+\right. \\ & q+r+s)+14 p^2(r+s+r s+ \\ & q(1+r+s))-28 p(r s+q(r+s+r s))),\end{aligned}$

$\begin{aligned} & \bar{E}_{21}=\frac{q^4}{840 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(8 q^3(3+r+s)-\right. \\ & 5 q^4-140 r s-14 q^2(2+3 s+r(3+ \\ & s))+28 q(2 s+r(2+3 s))), \\ & \bar{E}_{22}=\frac{q^2}{840(p-q)\left(2-3 q+q^2\right)(q-r)(q-s)}\left(2 p\left(10 q^4+140 r s-\right.\right. \\ & 14 q^3(3+r+s)+21 q^2(2+3 s+ \\ & r(3+s))-35 q(2 s+r(2+3 s)))+q\left(20 q^3(3+r+\right. \\ & s)-15 q^4-140 r s-28 q^2(2+3 s+ \\ & r(3+s))+42 q(2 s+r(2+3 s)))), \\ & \bar{E}_{23}=\frac{q^4}{840(p-1)(q-1)(r-1)(s-1)}\left(q\left(5 q^3-56 r s-8 q^2(2+\right.\right. \\ & r+s)+14 q(2 s+r(2+s)))- \\ & \left.2 p\left(4 q^3-70 r s-7 q^2(2+r+s)+14 q(2 s+r(2+s))\right)\right), \\ & \bar{E}_{24}=\frac{q^4}{840 r(-p+r)(r-q)\left(2-3 r+r^2\right)(r-s)}\left(q\left(56 s-5 q^3+\right.\right. \\ & \left.8 q^2(3+s)-14 q(2+3 s)\right)+2 p\left(4 q^3-\right. \\ & \left.\left.70 s-7 q^2(3+s)+14 q(2+3 s)\right)\right),\end{aligned}$

$\begin{aligned} & \bar{E}_{25}=\frac{q^4}{840(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}\left(q\left(5 q^3-56 r-\right.\right. \\ & \left.8 q^2(3+r)+14 q(2+3 r)\right)-2 p\left(4 q^3-\right. \\ & \left.\left.70 r-7 q^2(3+r)+14 q(2+3 r)\right)\right), \\ & \bar{E}_{26}=\frac{q^4}{1680(p-2)(q-2)(r-2)(s-2)}\left(q\left(28 r s-5 q^3+8 q^2(1+\right.\right. \\ & r+s)-14 q(r+s+r s))+ \\ & \left.2 p\left(4 q^3-35 r s-7 q^2(1+r+s)+14 q(r+s+r s)\right)\right), \\ & \bar{E}_{31}=\frac{1}{840 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}(11-20 s+ \\ & r(42 s-20)+q(42 s-20- \\ & 14 r(8 s-3))), \\ & \bar{E}_{32}=\frac{1}{840(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}(r(20- \\ & 42 s)-11+20 s+2 p(10-21 s+ \\ & 7 r(8 s-3))), \\ & \bar{E}_{33}=\frac{1}{840(p-1)(q-1)(r-1)(s-1)}(25-36 r-36 s+56 r s+ \\ & q(56 r-36+56 s-98 r s)+ \\ & p(-36+56 r+56 s-98 r s+14 q(4-7 r-7 s+15 r s))), \\ & \bar{E}_{34}=\frac{1}{840 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}(11-20 s+q(42 s- \\ & 20)+p(42 s-20-14 q(8 s-3))),\end{aligned}$

$\begin{aligned} & \left.\bar{E}_{35}=\frac{1}{840(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}(q)\right), \\ & 20 r+2 p(10-21 r+7 q(8 r-3)), \\ & \bar{E}_{36}=\frac{1}{1680(p-2)(q-2)(r-2)(s-2)}(8 r-5+8 s-14 r s+ \\ & 2 q(4-7 r-7 s+14 r s)-2 p(7 r-4+ \\ & 7 s-14 r s+7 q(1-2 r-2 s+5 r s))), \\ & \bar{E}_{41}=\frac{r^4}{840 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(r\left(56 s-5 r^3+\right.\right. \\ & \left.8 r^2(3+s)-14 r(2+3 s)\right)+2 q\left(4 r^3-\right. \\ & \left.\left.70 s-7 r^2(3+s)+14 r(2+3 s)\right)\right), \\ & \bar{E}_{42}=\frac{r^4}{840(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}\left(r\left(5 r^3-56 s-\right.\right. \\ & \left.8 r^2(3+s)+14 r(2+3 s)\right)-2 p\left(4 r^3-\right. \\ & \left.\left.70 s-7 r^2(3+s)+14 r(2+3 s)\right)\right), \\ & \bar{E}_{43}=\frac{r^4}{840(p-1)(q-1)(r-1)(s-1)}\left(2 p\left(7 q\left(r^2+10 s-\right.\right.\right. \\ & 2 r(2+s))+r\left(7 r(2+s)-4 r^2-\right. \\ & 28 s))+r\left(r\left(5 r^2+28 s-8 r(2+s)\right)+q\left(-8 r^2-56 s+\right.\right. \\ & 14 r(2+s)))),\end{aligned}$

$\begin{aligned} & \bar{E}_{44}=\frac{r^2}{840(p-r)(q-r)\left(2-3 r+r^2\right)(r-s)}\left(2 p\left(r\left(70 s-10 r^3+\right.\right.\right. \\ & \left.14 r^2(3+s)-21 r(2+3 s)\right)+ \\ & \left.7 q\left(2 r^3-20 s-3 r^2(3+s)+5 r(2+3 s)\right)\right)+ \\ & r\left(-2 q\left(10 r^3-70 s-14 r^2(3+s)+\right.\right. \\ & 21 r(2+3 s))+r\left(15 r^3-84 s-20 r^2(3+s)+28 r(2+\right. \\ & 3 s)))) \\ & \bar{E}_{45}=\frac{r^4}{840(r-s) s(-p+s)(s-q)\left(2-3 s+s^2\right)}\left(2 p\left(r\left(21 r-4 r^2-\right.\right.\right. \\ & \left.28)+7 q\left(10-6 r+r^2\right)\right)+r(q(42 r- \\ & \left.\left.\left.56-8 r^2\right)+r\left(28-24 r+5 r^2\right)\right)\right), \\ & \bar{E}_{46}=\frac{r^4}{1680(p-2)(q-2)(r-2)(s-2)}\left(-2 p\left(7 q\left(r^2+5 s-\right.\right.\right. \\ & 2 r(1+s))+r\left(7 r(1+s)-4 r^2-\right. \\ & 14 s))+r\left(2 q\left(4 r^2+14 s-7 r(1+s)\right)+r\left(-5 r^2-14 s+\right.\right. \\ & 8 r(1+s)))), \\ & \bar{E}_{51}=\frac{s^4}{840 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(7 r\left(10-6 s+s^2\right)-\right. \\ & 2 q\left(s\left(21 s-4 s^2-28\right)\right)+s(s(24 s- \\ & \left.\left.\left.28-5 s^2\right)+r\left(56-42 s+8 s^2\right)\right)\right), \\ & \bar{E}_{52}=\frac{s^4}{840(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}\left(2 p\left(s\left(21 s-4 s^2-\right.\right.\right. \\ & \left.28)+7 r\left(10-6 s+s^2\right)\right)+s(r(- \\ & \left.\left.\left.56+42 s-8 s^2\right)+s\left(28-24 s+5 s^2\right)\right)\right),\end{aligned}$

$\begin{aligned} & \bar{E}_{53}=\frac{s^4}{840(p-1)(q-1)(r-1)(s-1)}(-2 p(7 q(2 r(s-5)- \\ & (s-4) s)+s(-7 r(s-4)+2 s(2 s- \\ & 7)))+s(2 q(7 r(-4+s)+2(7-2 s) s)+s(r(28-8 s)+ \\ & s(-16+5 s)))) \\ & \bar{E}_{54}=\frac{s^4}{840 r(r-p)(-q+r)\left(2-3 r+r^2\right)(r-s)}\left(-2 p\left(s\left(21 s-4 s^2-\right.\right.\right. \\ & \left.28)+7 q\left(10-6 s+s^2\right)\right)+s(s(- \\ & \left.\left.\left.28+24 s-5 s^2\right)+q\left(56-42 s+8 s^2\right)\right)\right) \\ & \bar{E}_{55}=\frac{s^2}{840(r-s)(s-p)(s-q)\left(2-3 s+s^2\right)}\left(s\left(s\left(s\left(60 s-15 s^2-\right.\right.\right.\right. \\ & \left.56)+4 r\left(21-21 s+5 s^2\right)\right)- \\ & \left.2 q\left(7 r\left(10-9 s+2 s^2\right)-2 s\left(21-21 s+5 s^2\right)\right)\right)+ \\ & 2 p\left(7 q\left(s\left(9 s-10-2 s^2\right)+r(20-15 s+\right.\right. \\ & \left.\left.\left.\left.3 s^2\right)\right)+s\left(-7 r\left(10-9 s+2 s^2\right)+2 s\left(21-21 s+5 s^2\right)\right)\right)\right) \\ & \bar{E}_{56}=\frac{s^4}{1680(-2+p)(-2+q)(-2+r)(-2+s)}(s(-2 q(7 r(s-2)+ \\ & (7-4 s) s)+s((8-5 s) s+2 r(4 s- \\ & 7)))+2 p(-7 q(r(5-2 s)+(-2+s) s)+s(-7 r(-2+ \\ & s)+s(-7+4 s)))), \\ & \bar{E}_{61}=-\frac{4}{105 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}(q(-4+7 r s)- \\ & 4(-2+r+s)),\end{aligned}$

$\begin{aligned} & \bar{E}_{62}=\frac{4 h^2}{105(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}(p(-4+7 r s)-4(-2+ \\ & r+s)), \\ & \bar{E}_{63}=\frac{4}{105(p-1)(q-1)(r-1)(s-1)}(-4(8 r-10+q(8+7 r(s- \\ & 1)-7 s)+8 s-7 r s)+p(-4(8+ \\ & 7 r(-1+s)-7 s)+7 q(4-4 s+r(-4+5 s)))), \\ & \bar{E}_{64}=-\frac{4 h^2}{105 r(-p+r)(-q+r)\left(2-3 r+r^2\right)(r-s)}(p(-4+7 q s)- \\ & 4(-2+q+s)), \\ & \bar{E}_{65}=\frac{4 h^2}{105(r-s) s(-p+s)(-q+s)\left(2-3 s+s^2\right)}(p(-4+7 q r)- \\ & 4(-2+q+r)), \\ & \bar{E}_{66}=\frac{-2}{105(p-2)(q-2)(r-2)(s-2)}(24 r-40+q(24+7 r(s- \\ & 2)-14 s)+24 s-14 r s+p(24+ \\ & 7 r(-2+s)-14 s+7 q(-2+r+s))),\end{aligned}$

$\begin{aligned} & \acute{D}_{16}=\frac{p}{840 q r s}\left(420 q r s-10 p^5+14 p^4(3+q+r+\right. \\ & s)-21 p^3(2+3 s+r(3+s)+ \\ & q(3+r+s))+35 p^2(2 s+r(2+3 s)+q(2+3 s+ \\ & r(3+s)))-70 p(2 r s+q(2 s+ \\ & r(2+3 s)))) \\ & \acute{D}_{26}=\frac{q}{840 p r s}\left(7 p\left(2 q^4+60 r s-3 q^3(3+r+s)+\right.\right. \\ & 5 q^2(2+3 s+r(3+s))-10 q(2 s+ \\ & r(2+3 s)))+q\left(14 q^3(3+r+s)-10 q^4-140 r s-\right. \\ & 21 q^2(2+3 s+r(3+s))+35 q(2 s+ \\ & r(2+3 s)))) \\ & \acute{D}_{36}=\frac{1}{840 p q r s}(18-28 r-28 s+49 r s-7 q(4-7 r- \\ & 7 s+15 r s)+7 p(7 r-4+7 s- \\ & 15 r s+q(7-15 r-15 s+50 r s))) \\ & \acute{D}_{46}=\frac{r^2}{840 p q s}\left(r\left(70 s-10 r^3+14 r^2(3+s)-21 r(2+\right.\right. \\ & 3 s))+7 q\left(2 r^3-20 s-9 r^2-3 r^2 s+\right. \\ & 5 r(2+3 s)))-7 p\left(r\left(20 s-2 r^3+3 r^2(3+s)-\right.\right. \\ & 5 r(2+3 s))+q\left(3 r^3-60 s-5 r^2(3+\right. \\ & s)+10 r(2+3 s)))),\end{aligned}$

$\begin{aligned} & \acute{D}_{56}=\frac{s}{840 p q r}\left(7 p\left(q\left(s\left(15 s-20-3 s^2\right)+5 r(12-\right.\right.\right. \\ & \left.\left.6 s+s^2\right)\right)+s\left(r\left(-20+15 s-3 s^2\right)+\right. \\ & \left.\left.s\left(10-9 s+2 s^2\right)\right)\right)+s\left(-7 q\left(s\left(-10+9 s-2 s^2\right)+\right.\right. \\ & \left.r\left(20-15 s+3 s^2\right)\right)+s(7 r(10-9 s+ \\ & \left.\left.\left.\left.2 s^2\right)-2 s\left(21-21 s+5 s^2\right)\right)\right)\right), \\ & \acute{D}_{66}=\frac{1}{105 p q r s}(7 q(r+s-2)-2(24+7 r(s-2)- \\ & 14 s)+7 p(q(5 r s-2)-2(r+s-2))), \\ & \acute{E}_{11}=\frac{p}{420\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(60 p^5-420 q r s-\right. \\ & 70 p^4(3+q+r+s)+84 p^3(2+ \\ & 3 s+r(3+s)+q(3+r+s))-105 p^2(2 s+r(2+ \\ & 3 s)+q(2+3 s+r(3+s)))+ \\ & 140 p(2 r s+q(2 s+r(2+3 s))), \\ & \acute{E}_{12}=\frac{p^3}{420(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}\left(10 p^4+140 r s-\right. \\ & 14 p^3(3+r+s)+21 p^2(2+ \\ & 3 s+r(3+s))-35 p(2 s+r(2+3 s))), \\ & \acute{E}_{13}=\frac{p^3}{420(p-1)(q-1)(r-1)(s-1)}\left(10 p^4+140 q r s-\right. \\ & 14 p^3(2+q+r+s)+21 p^2(2 s+r(2+ \\ & s)+q(2+r+s))-35 p(2 r s+q(2 s+r(2+s)))),\end{aligned}$

$\begin{aligned} & \acute{E}_{14}=\frac{p^3}{420 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}\left(14 p^3(3+q+s)-\right. \\ & 10 p^4-140 q s-21 p^2(2+3 s+q(3+ \\ & s))+35 p(2 s+q(2+3 s))) \\ & \acute{E}_{15}=\frac{p^3}{420(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}\left(10 p^4+140 q r-\right. \\ & 14 p^3(3+q+r)+21 p^2(2+3 r+q(3+ \\ & r))-35 p(2 r+q(2+3 r))) \\ & \acute{E}_{16}=\frac{p^3}{840(p-2)(q-2)(r-2)(s-2)}\left(14 p^3(1+q+r+s)-\right. \\ & 10 p^4-70 q r s-21 p^2(r+s+r s+ \\ & q(1+r+s))+35 p(r s+q(r+s+r s))), \\ & \acute{E}_{21}=\frac{q^3}{420 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(14 q^3(3+r+s)-\right. \\ & 10 q^4-140 r s-21 q^2(2+3 s+r(3+ \\ & s))+35 q(2 s+r(2+3 s))), \\ & \acute{E}_{22}=\frac{q}{420(p-q)\left(2-3 q+q^2\right)(q-r)(q-s)}\left(7 p\left(10 q^4+60 r s-\right.\right. \\ & 12 q^3(3+r+s)+15 q^2(2+3 s+ \\ & r(3+s))-20 q(2 s+r(2+3 s)))+q\left(-60 q^4-\right. \\ & 280 r s+70 q^3(3+r+s)-84 q^2(2+ \\ & 3 s+r(3+s))+105 q(2 s+r(2+3 s)))),\end{aligned}$

$\begin{aligned} & \acute{E}_{23}=\frac{q^3}{420(p-1)(q-1)(r-1)(s-1)}\left(-7 p\left(2 q^3-20 r s-\right.\right. \\ & \left.3 q^2(2+r+s)+5 q(2 s+r(2+s))\right)+ \\ & \left.q\left(10 q^3-70 r s-14 q^2(2+r+s)+21 q(2 s+r(2+s))\right)\right) \\ & \acute{E}_{24}=\frac{q^3}{420 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}\left(q\left(70 s-10 q^3+\right.\right. \\ & \left.14 q^2(3+s)-21 q(2+3 s)\right)+ \\ & \left.7 p\left(2 q^3-20 s-3 q^2(3+s)+5 q(2+3 s)\right)\right) \\ & \acute{E}_{25}=\frac{q^3}{420(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}\left(q\left(10 q^3-70 r-\right.\right. \\ & \left.14 q^2(3+r)+21 q(2+3 r)\right)- \\ & \left.7 p\left(2 q^3-20 r-3 q^2(3+r)+5 q(2+3 r)\right)\right) \\ & \acute{E}_{26}=\frac{q^3}{840(p-2)(q-2)(r-2)(s-2)}\left(q\left(35 r s-10 q^3+\right.\right. \\ & 14 q^2(1+r+s)-21 q(r+s+ \\ & r s))+7 p\left(2 q^3-10 r s-3 q^2(1+r+s)+5 q(r+s+\right. \\ & r s))) \\ & \acute{E}_{31}=\frac{2}{420 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}(9-14 s+ \\ & \frac{49 r}{2}\left(s-\frac{4}{7}\right)-\frac{7 q}{2}(4-7 s+ \\ & r(15 s-7))),\end{aligned}$

$\begin{aligned} & \acute{E}_{32}=\frac{2}{420(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}(14 s-9- \\ & \frac{49 r}{2}\left(s-\frac{4}{7}\right)+\frac{7 p}{2}(4-7 s+ \\ & r(15 s-7))), \\ & \acute{E}_{33}=\frac{1}{420(p-1)(q-1)(r-1)(s-1)}(80-98 r-98 s+ \\ & 126 r s-7 q(14-18 r- \\ & 18 s+25 r s)+7 p(-14+18 r+18 s-25 r s+18 q- \\ & 25 q s-25 q r+40 q r s)), \\ & \acute{E}_{34}=\frac{2}{420 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}(9-14 s+ \\ & \frac{7 q}{2}(7 s-4)-\frac{7 p}{2}(4-7 s+ \\ & q(15 s-7))), \\ & \acute{E}_{35}=\frac{2}{420(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}(14 r-9- \\ & \frac{7 q}{2}(7 r-4)+\frac{7 p}{2}(4-7 r+ \\ & q(15 r-7))), \\ & \acute{E}_{36}=\frac{1}{840(p-2)(q-2)(r-2)(s-2)}(14 r-10+14 s-21 r s+ \\ & 7 q(2-3 r-3 s+5 r s)- \\ & 7 p(-2+3 r+3 s-5 r s+q(3-5 r-5 s+10 r s))),\end{aligned}$

$\begin{aligned} & \acute{E}_{41}=\frac{r^3}{420 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}\left(r\left(70 s-10 r^3+\right.\right. \\ & \left.14 r^2(3+s)-21 r(2+3 s)\right)+ \\ & \left.7 q\left(2 r^3-20 s-3 r^2(3+s)+5 r(2+3 s)\right)\right), \\ & \acute{E}_{42}=\frac{r^3}{420(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}\left(-7 p\left(2 r^3-20 s-\right.\right. \\ & \left.3 r^2(3+s)+5 r(2+3 s)\right)+ \\ & \left.r\left(10 r^3-70 s-14 r^2(3+s)+21 r(2+3 s)\right)\right), \\ & \acute{E}_{43}=\frac{2 r^3}{420(p-1)(q-1)(r-1)(s-1)}\left(r\left(2 r\left(5 r^2+21 s-\right.\right.\right. \\ & 7 r(2+s))-7 q\left(2 r^2+10 s-\right. \\ & 3 r(2+s)))+7 p\left(q\left(3 r^2+20 s-5 r(2+s)\right)+r\left(-2 r^2-\right.\right. \\ & 10 s+3 r(2+s)))), \\ & \acute{E}_{44}=\frac{r}{420(p-r)(q-r)\left(2-3 r+r^2\right)(r-s)}\left(7 p\left(r\left(40 s-10 r^3+\right.\right.\right. \\ & 12 r^2(3+s)-15 r(2+ \\ & 3 s))+q\left(12 r^3-60 s-15 r^2(3+s)+20 r(2+\right. \\ & 3 s)))+r\left(-7 q\left(10 r^3-40 s-\right.\right. \\ & \left.12 r^2(3+s)+15 r(2+3 s)\right)+2 r\left(30 r^3-105 s-\right. \\ & \left.\left.\left.35 r^2(3+s)+42 r(2+3 s)\right)\right)\right),\end{aligned}$

$\begin{aligned} & \acute{E}_{45}=\frac{r^3}{420(r-s) s(s-p)(s-q)\left(2-3 s+s^2\right)}(7 p(r(-10+9 r- \\ & \left.\left.2 r^2\right)+q\left(20-15 r+3 r^2\right)\right)+ \\ & \left.r\left(-7 q\left(10-9 r+2 r^2\right)+2 r\left(21-21 r+5 r^2\right)\right)\right) \\ & \acute{E}_{46}=\frac{r^3}{840(p-2)(q-2)(r-2)(s-2)}\left(-7 p\left(q\left(3 r^2+10 s-\right.\right.\right. \\ & 5 r(1+s))+r(3 r(1+s)- \\ & \left.\left.2 r^2-5 s\right)\right)+r\left(7 q\left(2 r^2+5 s-3 r(1+s)\right)+r\left(-10 r^2-\right.\right. \\ & 21 s+14 r(1+s)))) \\ & \acute{E}_{51}=\frac{s^3}{420 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}(-7 q(s(9 s-10- \\ & \left.2 s^2\right)+r(20-15 s+ \\ & \left.\left.\left.3 s^2\right)\right)+s\left(7 r\left(10-9 s+2 s^2\right)-2 s\left(21-21 s+5 s^2\right)\right)\right) \\ & \acute{E}_{52}=\frac{s^3}{420(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}(7 p(s(9 s-10- \\ & \left.\left.2 s^2\right)+r\left(20-15 s+3 s^2\right)\right)+ \\ & \left.s\left(-7 r\left(10-9 s+2 s^2\right)+2 s\left(21-21 s+5 s^2\right)\right)\right) \\ & \acute{E}_{53}=\frac{s^3}{420(p-1)(q-1)(r-1)(s-1)}(-7 p(q(5 r(s-4)+(10- \\ & 3 s) s)+s(r(10-3 s)+ \\ & 2(s-3) s))+s(7 q(-2(s-3) s+r(-10+3 s))+ \\ & 2 s(-7 r(s-3)+s(-14+5 s))))\end{aligned}$

$\begin{aligned} & \acute{E}_{54}=\frac{s^3}{420 r(r-p)(r-q)\left(2-3 r+r^2\right)(r-s)}\left(-7 p\left(s\left(9 s-2 s^2-\right.\right.\right. \\ & \left.10)+q\left(20-15 s+3 s^2\right)\right)+ \\ & \left.s\left(7 q\left(10-9 s+2 s^2\right)-2 s\left(21-21 s+5 s^2\right)\right)\right), \\ & \acute{E}_{55}=\frac{s}{420(q-s)(s-p)(s-r)\left(2-3 s+s^2\right)}\left(7 p\left(q\left(15 r(s-2)^2+\right.\right.\right. \\ & \left.s\left(-40+45 s-12 s^2\right)\right)+ \\ & \left.s\left(r\left(-40+45 s-12 s^2\right)+2 s\left(15-18 s+5 s^2\right)\right)\right)+ \\ & s\left(2 s\left(7 r\left(15-18 s+5 s^2\right)-\right.\right. \\ & \left.3 s\left(28-35 s+10 s^2\right)\right)-7 q\left(-2 s\left(15-18 s+5 s^2\right)+\right. \\ & \left.\left.\left.r\left(40-45 s+12 s^2\right)\right)\right)\right), \\ & \acute{E}_{56}=\frac{s^3}{840(-2+p)(-2+q)(-2+r)(-2+s)}(s(s(2(7-5 s) s+ \\ & 7 r(2 s-3))+7 q(r(5-3 s)+ \\ & s(-3+2 s))+7 p(q(5 r(-2+s)+(5-3 s) s)+s(r(5- \\ & 3 s)+s(-3+2 s)))), \\ & \acute{E}_{61}=-\frac{4}{105 p\left(2-3 p+p^2\right)(p-q)(p-r)(p-s)}(24+7 r(-2+s)- \\ & 14 s+7 q(-2+r+s)), \\ & \acute{E}_{62}=\frac{4}{105(p-q) q\left(2-3 q+q^2\right)(q-r)(q-s)}(24+7 r(-2+s)- \\ & 14 s+7 p(-2+r+s)),\end{aligned}$

$\begin{aligned} & \acute{E}_{63}=\frac{4}{105(p-1)(q-1)(r-1)(s-1)}(80-56 r-56 s+ \\ & 42 r s+7 p(6 r-8+q(6+ \\ & 5 r(s-1)-5 s)+6 s-5 r s)-7 q(8-6 r-6 s+5 r s)), \\ & \acute{E}_{64}=-\frac{4}{105 r(-p+r)(-q+r)\left(2-3 r+r^2\right)(r-s)}(24+7 q(-2+s)- \\ & 14 s+7 p(-2+q+s)), \\ & \acute{E}_{65}=\frac{4}{105(r-s) s(-p+s)(-q+s)\left(2-3 s+s^2\right)}(24+7 q(-2+r)- \\ & 14 r+7 p(-2+q+r)), \\ & \acute{E}_{66}=\frac{1}{105(p-2)(q-2)(r-2)(s-2)}(7 p(q(18+5 r(s-2)- \\ & 10 s)+2(r(9-5 s)-16+\\ & 9 s))+2(200-112 s+7 r(9 s-16)-7 q(16-9 s+ \\ & r(5 s-9))))\end{aligned}$


[1] Jikantoro, D.Y. (2014). Numerical solution of special second order initial value problems by hybrid type methods. Doctoral dissertation, Universiti Putra Malaysia. 

[2] Shampine, L.F. (2018). Numerical solution of ordinary differential equations. Routledge.

[3] Borzì, A. (2020). Modelling with ordinary differential equations: A comprehensive approach. Chapman and Hall/CRC.

[4] Ming, C.Y. (2017). Solution of differential equations with applications to engineering problems. Dynamical Systems-Analytical and Computational Techniques, 15: 233-264.

[5] Awoyemi, D.O. (2003). A p-stable linear multistep method for solving general third order ordinary differential equations. International Journal of Computer Mathematics, 80(8): 985-991.

[6] Shokri, A. (2015). The symmetric two-step P-stable nonlinear predictor-corrector‎‎ methods for the numerical solution of second order‎‎ initial value problems. Bulletin of the Iranian Mathematical Society, 41(1): 201-215. 

[7] Panopoulos, G.A., Anastassi, Z.A., Simos, T.E. (2013). A new eight-step symmetric embedded predictor-corrector method (EPCM) for orbital problems and related IVPs with oscillatory solutions. The Astronomical Journal, 145(3): 1-9.

[8] Ndanusa, A., Tafida, F.U. (2016). Predictor–corrector methods of high order for numerical integration of initial value problems. International Journal of Scientific and Innovative Mathematical Research (IJSIMR), 4(2): 47-55.

[9] Biasa, P.T., Majid, Z.A., Suleiman, M. (2011). Predictor-corrector block iteration method for solving ordinary differential equations. Sains Malaysiana, 40(6): 659-664. 

[10] Kayode, S.J., Adegboro, J.O. (2018). Predictor-corrector linear multistep method for direct solution of initial value problems of second order ordinary differential equations. Asian Journal of Physical and Chemical Sciences, 6(1): 1-9.

[11] Tam, H.W. (1992). Two-stage parallel methods for the numerical solution of ordinary differential equations. SIAM Journal on Scientific and Statistical Computing, 13(5): 1062-1084.

[12] Olabode, B.T. (2009). An accurate scheme by block method for the third order ordinary differential equation. Pacific Journal of Science and Technology, 10(1): 136-142.

[13] Milne, W.E. (1953). Numerical solution of ordinary differential equations. New York: Wiley.

[14] Sarafyan, D. (1965). Multistep methods for the numerical solution of ordinary differential equations made self-starting. Mathematics Research Center, MRC-TSR-495. 

[15] Adesanya, A.O. (2011). Block methods for direct solutions of general higher order initial value problems of ordinary differential equations. Doctoral dissertation, Federal University of Technology, Akure.

[16] Kamoh, N.M., Abada, A.A., Soomiyol, M.C. (2018). A Block procedure with continuous coefficients for the direct solution of general second order initial value problems of (ODEs) using shifted legendre polynomials as basis function. International Journal of Multidisciplinary Research and Development, 5(4): 236-241. 

[17] Kuboye, J.O., Omar, Z. (2014). Numerical solution of third order ordinary differential equations using a seven-step block method. International Journal of Mathematical Analysis, 9(15): 743-754.

[18] Jena, S.R., Gebremedhin, G.S. (2020). Approximate solution of a fifth order ordinary differential equations with block method. International Journal of Computing Science and Mathematics, 12(4): 413-426.

[19] Adeyeye, O., Omar, Z. (2016). A new algorithm for developing block methods for solving fourth order ordinary differential equations. Global Journal of Pure and Applied Mathematics, 12(2): 1465-1471.

[20] Ogunware, B.G., Omole, E.O. (2020). A class of irrational linear multistep block method for the direct numerical solution of third order ordinary differential equations. Turkish Journal of Analysis and Number Theory, 8(2): 21-27.

[21] Sagir, A.M. (2014). Numerical treatment of block method for the solution of ordinary differential equations. International Journal of Bioengineering and Life Sciences, 8(2): 259-263. 

[22] Olukunle, O.E., Felix, M.O. (2020). Bernstein induced one step hybrid scheme for general solution of second order initial value problems. Malaya Journal of Matematik, 8(2): 350-355.

[23] Shokri, A., Mehdizadeh K.M., Mohammad-Sedighi, H., Atashyar, A. (2021). Numerical simulation of second-order initial-value problems using a new class of variable coefficients and two-step semi-hybrid methods. Simulation, 97(5): 347-364.

[24] Ehigie, J.O., Okunuga, S.A., Sofoluwe, A.B., Akanbi, M.A. (2010). On generalized 2-step continuous linear multistep method of hybrid type for the integration of second order ordinary differential equation. Archives of Applied Science Research, 2(6): 362-372. 

[25] Fasasi, K.M., Adesanya, A.O., Adee, S.O. (2014). One step continuous hybrid block method for the solution of. Journal of Natural science Research, 4(10): 55-62.

[26] Gebremedhin, G.S., Jena, S.R. (2019). Approximate solution of ordinary differential equation via hybrid block approach. International Journal on Emerging Technologies, 10(4): 201-211. 

[27] Mansor, K.H, Omar, Z, Rohni, M.A. (2018). Generalized two-step hybrid block method with one off-step point for solving second order ordinary differential equations directly. Journal of Advanced Research in Dynamical and Control Systems, 13: 1-17.

[28] Henrici, P. (1962). Discrete variable methods in ordinary differential equations. New York: Wiley.

[29] Abdelrahim, R. (2016). One step hybrid block methods with generalized off-step points for solving higher order ordinary differential equations directly. PhD Thesis, Universiti Utara Malaysia, Malaysia.

[30] Anake, T.A. (2011). Continuous implicit hybrid one-step methods for the solution of initial value problems of general second-order ordinary differential equations. PhD Thesis, Covenant University, Nigeria.

[31] Kuboye, J. (2015). Block methods for direct solution of higher order ordinary differential equations using interpolation and collocation approach. PhD Thesis, Universiti Utara Malaysia, Malaysia.

[32] Mohammed, U., Jiya, M., Mohammed, A.A. (2010). A class of six step block method for solution of general second order ordinary differential equations. Pacific Journal of Science and Technology, 11(2): 273-277. 

[33] Yahaya, Y.A., Sagir, A.M. (2013). An order five implicit 3-step block method for solving ordinary differential equations. Pacific Journal of Science and Technology, 14(1): 176-181.