Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order Fractional Differential Equations

Mohammed S. Mechee; Mohammed Mahmood Salih

doi:10.12688/f1000research.172397.1

Home Browse Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order...

ALL Metrics

-

Views

-

Downloads

Get PDF

Get XML

Export

▬

✚

Research Article

Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order Fractional Differential Equations

[version 1; peer review: awaiting peer review]

Mohammed S. Mechee ¹, Mohammed Mahmood Salih²

PUBLISHED 24 Dec 2025

Author details Author details

¹ TRDC, university of Kufa, Najaf, Najaf-Iraq, 00964, Iraq
² Faculty of Information Technology Ninevah, Ninevah University, Mosel, 00964, Iraq

Mohammed S. Mechee
Roles: Formal Analysis, Methodology, Software, Supervision, Validation

Mohammed Mahmood Salih
Roles: Conceptualization, Data Curation, Formal Analysis, Funding Acquisition, Resources, Software, Writing – Original Draft Preparation, Writing – Review & Editing

OPEN PEER REVIEW

REVIEWER STATUS AWAITING PEER REVIEW

This article is included in the Fallujah Multidisciplinary Science and Innovation gateway.

Abstract

Background

In this paper, we focus on deriving an efficient method for solving ordinary differential equations (ODEs) of sixth order, and then, we modify the proposed method for solving fractional differential equations (FDEs).

Methods

The methodology of this paper used the approach of derivation of implicit numerical method. However, the purpose of this study is to derive a direct implicit block and multi-step strategy for solving sixth-order ODEs.

Results

The main novelty of this work lies in the investigation of constructing a novel method for solving FDEs. This paper presents a general block two-point approach, known as GIBM2P, for solving sixth-order initial-value problems (IVPs) using Hermite-interpolating polynomials. Consequently, the FDE problem has been transformed into ODE using the α-fractional-derivative transform, and then, the sixth-order IVP is solved using the proposed GIBM2P method.

Conclusions

In this paper, we offered a method for solving sixth-order FDEs. The proposed method has been shown to be an accurate and efficient method; also, one can obtain the approximate efficient solutions of FDEs problems. The numerical implementations are used to show the high accuracy and efficacy of the proposed GIBM2P approach. However, the GIBM2P method is effective in terms of accuracy, processing time, and the number of function calls.

Keywords

Implicit Method; Block Method; Multi-Step Method; ODEs; IVPs; Sixth-Order; FDEs.

Corresponding author: Mohammed S. Mechee

Competing interests: No competing interests were disclosed.

Grant information: The author(s) declared that no grants were involved in supporting this work.

Copyright: © 2025 Mechee MS and Mahmood Salih M. 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 work is properly cited.

How to cite: Mechee MS and Mahmood Salih M. Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order Fractional Differential Equations [version 1; peer review: awaiting peer review]. F1000Research 2025, 14:1444 (https://doi.org/10.12688/f1000research.172397.1) First published: 24 Dec 2025, 14:1444 (https://doi.org/10.12688/f1000research.172397.1) Latest published: 24 Dec 2025, 14:1444 (https://doi.org/10.12688/f1000research.172397.1)

1. Introduction

Differential equations (DEs) are very important in many applications of mathematical models in the applied sciences and engineering, particularly in the physics, electrical engineering, economics, chemistry, and biology.¹ However, the mathematical models in the area of applied science, and engineering of actual issues are constructed using DEs techniques. Unfortunately, determining the exact solutions of some kinds of DEs might be somewhat challenging, in spite of numerous efficient common modern or classical analytical and numerical techniques are currently available to the researchers of engineers and scientists. In general, the computation of the analytical solutions of general DEs is sometimes difficult or complicated. Various researchers derived some numerical linear multistep methods for solving IVPs for different orders of ODEs. For the literature reviews (LRs) of this study, we have introduced some studies of numerical methods for solving some classes of ODEs. Firstly, for a review several researchers developed some methods for solving boundary value problems (BVPs) of sixth order ODEs to increase the effectiveness and the accuracy of these approaches^2–7 while Ref. 8 presented a comparison study and a survey of the researches for the generalized RK integrators for solving IVPs of various orders ODEs. Secondly, in likewise, some numerical algorithms of implicit block multistep methods for solving IVPs of ODEs for various orders have been improved by several authors in Refs. 9–15. Meanwhile, for example, the researchers Refs. 16–18 derived block backward methods for solving first-order stiff ODEs, as well as Refs. 19–21 derived one-step, implicit block, block predictor, and 4-step corrector methods for solving a general second-order ODEs. However, the implicit block method and a new strategy of the implicit hybrid block method for solving third-order ODEs are introduced by Refs. 22, 23. In the same manner,²⁴introduced a new special 15-step block method for solving general fourth-order ODEs. In contrast, the general implicit-block method for solving fifth-order ODEs has been derived by Ref. 25 but Turki et al. are constructed a general implicit-block method for solving seventh-order ODEs.²⁶ Moreover, for variable-step block methods, the authors in Refs. 27, 28 have derived a variable-step variable-order block fully implicit method and direct variable-step block multistep method for solving general third-order ODEs. Finally, an explicit multistep block method for solving neutral delay differential equations (DDEs) has been derived by Ref. 29 while Ref. 30 developed three-point diagonally implicit block methods for solving ordinary and fuzzy DEs.

Nowadays, the FDEs and their applications have a significate role in the fields of science and engineering. As a consequence, the researches in finding new approaches for FDEs are increasing.³¹ As a result, numerous researchers have proposed new definitions for fractional derivatives. For instance, the α-fractional-derivative or integral of real functions was defined in a novel way by Mechee et al., which is used in this article.³² Also, Khalil et al. and Z. Zheng et al.^33,34offered more new definitions of the fractional derivative. This article utilizes the fundamental functions of a Hermite-interpolating polynomial to propose a general implicit block two-point multistep approach, called GIBM2P, for solving the IVPs of sixth-order ODEs. Furthermore, this new technique has been used for solving sixth-order FDEs by converting FDE to ODE using α-fractional derivative.

2. Preliminary and background

Several concepts and definitions have been introduced in this section.

2.1 Quasi-Linear Sixth-Order Ordinary Differential Equations

The general formula for the sixth-order ODE’s initial value problem can be expressed as follows⁷:

(1)

Ψ (ζ, u (ζ), u^{'} (ζ), u^{''} (ζ), u^{'''} (ζ), u^{(4)} (ζ), u^{(5)} (ζ), u^{(6)} (ζ)) = 0, a \leq ζ \leq b

with the initial-conditions (ICs),

(2)

u^{(j)} (a) = ζ^{j}, j = 0, 1, 2, \dots, 5 .

while the initial-values problem of the special quasi-linear ODEs of sixth-order can be given as follows:

(3)

u^{(6)} (ζ) = ϕ (ζ, u (ζ), u^{'} (ζ), u^{''} (ζ), u^{'''} (ζ), u^{(4)} (ζ), u^{(5)} (ζ)), a \leq ζ \leq b

with the ICs in Equation (2).

Definition 2.1:

General Class of Sixth-Order FDEs

Below is the definition of the general class of sixth-order FDEs.

(4)

ϕ (ζ, u (ζ), u^{(j)} (ζ), u^{(jα)} (ζ)) = 0, j = 1, 2, \dots, 6; 0 < α < 1, a \leq ζ \leq b,

with the ICs,

(5)

u^{(j α)} (a) = ζ^{j}, j = 0, 1, 2, \dots, 5 .

Definition 2.2:

Special Cass of Quasi-Linear FDEs of Sixth-Order

The following formula illustrates the class of sixth-order FDEs:

(6)

u^{(6 α)} (ζ) = f (ζ, u (ζ), u^{(j)} (ζ), u^{(jα)} (ζ)), a \leq ζ \leq b; 0 < α < 1; j = 1, 2, \dots, 5 .

with the ICs in Equation (5).

2.1 Fractional derivatives

In this subsection, we introduced the α-fractional derivative of a function

Definition 2.3:

α-Fractional-derivative¹⁵

Mechee et al.³² introduced a novel α-fractional-derivative definition for the real function $ϕ (ζ) : [a, \infty) \to R$ by following:

(7)

T_{α} (ϕ (ζ)) = ϕ^{(α)} (ζ) = lim_{ϵ \to 0} \frac{ϕ (ζ + ϵ ζ^{1 - α}) - ϕ (ζ - ϵ ζ^{1 - α})}{2 ϵ},

for

α \in (0, 1) .

From the α-fractional derivative in Definition 2.3, we get the following property of α-fractional derivative of the function $ϕ (ζ) .$

Property 2.1

¹⁵Consider $ϕ (ζ) : [a, b) \to R$ be a real function.

Then, the α-fractional-derivative of the function $ϕ (ζ)$ defined in the Equation (7) has the following property:

(8)

T_{α} (ϕ (ζ)) = ζ^{1 - α} \frac{dϕ (ζ)}{dζ} .

3. Proposed method

We introduce the analysis of GIBM2P method which approximate the solution of IVP using Hermit polynomials in this section.

3.1 Analysis of the proposed numerical GIBM2P method

The form of a Hermit polynomial is defined as follows:

P_{2} (ξ) = \sum_{i = 0}^{m} \sum_{k = 0}^{m_{i} - 1} ϕ_{i}^{(k)} L_{ik} (ξ) .

Consider $h = \frac{b - a}{m}$ and $ξ_{i} = a + ih$ and $f_{i} = f (ξ_{i}, y_{i}, y_{i}^{'}, y_{i}^{''}),$ for $i = 0, 1, \dots, m$ where m = number of partitions in domain interval which is equal to the number of Lagrange’s polynomials and $m_{i} - 1$ is the degree of Lagrange polynomial. In this case, $L_{ik} (ξ)$ is the generalized Lagrange polynomial for m = 0,1,2, …, $m_{i}$ , with a step-size equal to 2h, and $y_{m + 1}$ and $y_{m + 2}$ are the approximation generates the approximated solutions at two-points $ξ_{m + 1}$ and $ξ_{m + 2}$ resp. In the block interval $[ξ_{m}, ξ_{m + 2}],$ we have, $ξ_{m}$ is the starting-point while $ξ_{m + 2}$ is the end-point. However, the initial value for all next iterations should be the approximated solution by $y_{m + 2}$ at the end-point $ξ_{m + 2}$

(9)

P_{2} (ξ) = f_{0} L_{00} (ξ) + f_{1} L_{10} (ξ) + f_{2} L_{20} (ξ) + g_{0} L_{01} (ξ) + g_{1} L_{11} (ξ) + g_{2} L_{21} (ξ),

3.1.1 Lagrange polynomials

Lagrange polynomials are defined as follows in this article:

(10)

L_{00} (λ) = {(λ - λ_{n}) (1 - χ_{1}) (\frac{(λ_{n + 2} - λ) (λ_{n + 1} - λ)}{(λ_{n} - λ_{n + 2}) (λ_{n} - λ_{n + 1})})}^{2},

(11)

L_{10} (λ) = {(λ - λ_{n + 1}) (1 - χ_{2}) (\frac{(λ - λ_{n}) (λ_{n + 2} - λ)}{(λ_{n + 1} - λ_{n}) (λ_{n + 1} - λ_{n + 2})})}^{2},

(12)

L_{20} (λ) = {(λ - λ_{n + 1}) (1 - χ_{3}) (\frac{(λ - λ_{n}) (λ_{n + 1} - λ)}{(λ_{n + 2} - λ_{n}) (λ_{n + 2} - λ_{n + 1})})}^{2},

(13)

L_{01} (λ) = (λ - λ_{n}) {(\frac{λ_{n + 1} - λ}{λ_{n} - λ_{n + 1}})}^{2} {(\frac{λ_{n + 2} - λ}{λ_{n} - λ_{n + 2}})}^{2},

(14)

L_{11} (λ) = (λ - λ_{n + 1}) {(\frac{(λ - λ_{n}) (λ_{n + 2} - λ)}{(λ_{n + 1} - λ_{n}) (λ_{n + 1} - λ_{n + 2})})}^{2},

(15)

L_{21} (λ) = (λ - λ_{n + 2}) {(\frac{(λ_{n} - λ) (λ_{n + 1} - λ)}{(λ_{n + 2} - λ_{n}) (λ_{n + 2} - λ_{n + 1})})}^{2},

where,

χ_{1} = \frac{2}{λ_{n} - λ_{n + 2}} + \frac{2}{λ_{n} - λ_{n + 1}}, χ_{2} = \frac{2}{λ_{n + 1} - λ_{n + 2}} + \frac{2}{λ_{n + 1} - λ_{n}}, χ_{3} = \frac{2}{λ_{n + 2} - λ_{n + 1}} + \frac{2}{λ_{n + 2} - λ_{n + 1}},

The expressions for Lagrange polynomials in the independent variable s are as follows, assuming that $s = \frac{λ - λ_{n + 2}}{h}$ and,

(16)

L_{00} (λ) = \frac{1}{4} {(β_{1} β_{2})}^{2} + 3 β_{3},

(17)

L_{10} (λ) = {(λ β_{3})}^{2},

(18)

L_{20} (λ) = \frac{1}{4} β_{3} {(β_{2})}^{2} β_{4}

(19)

L_{01} (λ) = \frac{h}{4} β_{3} {(β_{1} β_{2}))}^{2},

(20)

L_{11} (λ) = h {(λ β_{3})}^{2} β_{2},

and

(21)

L_{21} (λ) = \frac{h}{4} {(β_{3} β_{2})}^{2} .

where $β_{1} = λ, β_{2} = λ + 1, β_{3} = λ + 2, β_{4} = 1 - 3 λ$ .

3.1.2 Construction of the numerical GIBM2P approach

The approximated solution $u_{n} (ζ)$ of Equation (3) at the initial point $ζ_{n}$ can be obtained by multiple-integrating Equation (3) up to five times with respect to the independent variable $ξ$ along the range $[ζ_{n}; ζ_{n + 1}],$ then, compute these integrations to obtain the formulas in Equation (22).

(22)

u_{n + 1}^{(6 - j)} = \sum_{k = 6 - j}^{5} \frac{h^{6 - j - k}}{(6 - j - k)!} v_{n}^{(k)} + \int_{t_{n}}^{t_{n + 1}} ψ (λ, u (λ) d λ, j = 1, 2, \dots, 6

By assuming that $t_{n + 2} = t - hs$ , the coordinate t is changed to s.

The following six formulas were created by integrating the Equation (22)

u_{n + 1}^{(j)} = ∆_{6 - j} + \sum_{k = 1}^{3} a_{jk} f_{n + k} + a_{jk + 3} g_{n + k}, for j = 5,4,3,2,1,0

where,

(23)

A = (a_{ij}) = (\begin{array}{l} \begin{matrix} \frac{101}{240} & \frac{8}{15} & \frac{11}{240} \\ \frac{13}{42} & \frac{1}{6} & \frac{1}{42} \\ \frac{817}{6720} & \frac{4}{105} & \frac{47}{6720} \end{matrix} & \begin{matrix} \frac{13 h}{240} & - \frac{h}{6} & - \frac{h}{80} \\ \frac{59}{1680} & \frac{8}{105} & \frac{11}{1680} \\ \frac{83}{6720} & \frac{1}{48} & \frac{13}{6720} \end{matrix} \\ \begin{matrix} \frac{67}{2016} & \frac{1}{144} & \frac{1}{672} \\ \frac{1690}{241920} & \frac{1}{945} & \frac{61}{241920} \\ \frac{577}{475200} & \frac{1}{7200} & \frac{17}{475200} \end{matrix} & \begin{matrix} \frac{37}{12069} & \frac{4}{945} & \frac{5}{12096} \\ \frac{143}{241920} & \frac{1}{1440} & \frac{17}{241920} \\ \frac{179}{1900800} & \frac{1}{10395} & 19008 \end{matrix} \end{array})

and,

∆_{j} = \sum_{k = 6 - j}^{5} \frac{h^{6 - j - k}}{(6 - j - k)!} v_{n}^{(k)}, for j = 1, 2, \dots, 6 .

4. Main results

4.1 Analysis of proposed approach for solving of FDEs sixth-order

In this subsection, the analysis of the proposed method for solving FDEs of sixth-order is introduced and studied. In Equation (7), Mechee et al. presented a new definition of the α-fractional-derivative for the real function $f (τ) : [a, \infty) \to R .$ From the α-fractional- derivative in Definition 2.3, we get the property of α-fractional-derivative of the function $ϕ (ζ)$ in Equation (8). Consequently, by deriving the two sides of Equations (8) five times to get the sixth α-fractional-derivative of the function $f (τ)$ . Thus, the α-fractional-derivative of sixth-order for the real function $f (τ)$ in the domain $I_{α}$ , denoted by $T_{6 α}$ , has the following definition:

T_{6 α} (f (ζ)) = ζ^{1 - 6 α} (ζ^{5} f^{(6)} (ζ) + 15 (1 - α) ζ^{4} f^{(5)} (ζ) + 5 (1 - α) (13 - 7 α) ζ^{3} f^{(4)} (ζ) + (1 - α) (15 (5 - 7 α) + 2 (1 - 2 α) (7 - 11 α)) ζ^{2} f^{(3)} (ζ) + 2 (1 - α) (1 - 2 α) (7 - 11 α) (2 - 5 α) ζ f^{''} (ζ) + (1 - α) (1 - 2 α) (1 - 3 α) (1 - 4 α) (1 - 5 α) f^{'} (ζ)),

which is simplified to the following

(24)

T_{6 α} (f (ζ)) = ζ^{1 - 6 α} (ζ^{5} f^{(6)} (ζ) + 15 (1 - α) ζ^{4} f^{(5)} (ζ) + (1 - α) (5 (13 - 7 α) ζ^{3} f^{(4)} (ζ) + (15 (5 - 7 α) + 2 (1 - 2 α) (7 - 11 α)) ζ^{2} f^{(3)} (ζ) + 2 (1 - 2 α) (7 - 11 α) (2 - 5 α) ζ f^{''} (ζ) + (1 - 2 α) (1 - 3 α) (1 - 4 α) (1 - 5 α) f^{'} (ζ))) .

To prove Equation (24), by using $α$ -fractional-derivative five times for two sides of Equation (8) to obtain the formula in Equation (24)

4.2 Proposed analytical method

In this subsection, using the α-fractional derivative property from Equation (24), the FDE in Equation (6) with ICs in Equation (5) has been transformed into an ODE in Equation (3) with ICs in Equation (2) and then, using GIBM2P approach to obtain the numerical solution of the IVP.

4.2.1 The proposed method’s algorithm

The solution of FDE in Equation (6) with the ICs in Equation (5) has been introduced the proposed method using the following algorithm:

4.2.1.1 The Proposed Algorithm for Solving Sixth-Order FDEs

Step I: Convert the FDE in Equation (6) with ICs in Equation (5) into an ODE in Equation (3) with ICs in Equation (2).

Step II: Utilizing an appropriate analytical technique, or the proposed GIBM2P approach, solve the ODE in Equation (3) with the ICs in Equation (2).

Step III: The analytical or numerical solution for the IVP in step II is consequently exactly the same solution of given FDE problem in Equation (6) with the ICs in Equation (5).

5. Implementation

The performance of the numerical GIBM2P approaches was evaluated in this section since we transformed the FDE in Equation (6) to a sixth-order ODE and then solved it using proposed strategy specified in Equation (23) compared with Taylor method. The comparison of numerical solutions of proposed GIBM2P method against Taylor method is shown in Figure 1. The following is how the notations were used:

• GIBM2P: A two-point general implicit block approach is proposed.

Example 5.1:

Linear FDE³¹

η^{(6 α)} (ζ) = η (ζ) + η^{'} (ζ) - 3 η^{''} (ζ) - 11 η^{'''} (ζ); 0 < α < 1; 0 < ζ \leq 1

.

With the ICs: $η^{(κα)} (0) = {(- 1)}^{κ}$ for $κ = 0, 1, 2, 3, 4, 5 .$

Using the formula in Equation (24), we get the following ODE

T_{6 α} (f (ζ)) = ζ^{\frac{- 39}{9}} (ζ^{5} f^{(6)} (ζ) + \frac{15}{9} ζ^{4} f^{(5)} (ζ) \frac{305}{81} ζ^{3} f^{(4)} (ζ) - \frac{1135}{729} ζ^{2} f^{(3)} (ζ) - \frac{7700}{6561} ζ f^{''} (ζ) + \frac{74865}{59049} f^{'} (ζ)) .

Hence,

T_{6 α} (f (ζ)) = 2 ζ^{\frac{2}{3}} f^{(6)} (ζ) + ζ^{\frac{- 39}{9}} (\frac{15}{9} ζ^{4} f^{(5)} (ζ) \frac{305}{81} ζ^{3} f^{(4)} (ζ) - \frac{1135}{729} ζ^{2} f^{(3)} (ζ) - \frac{7700}{6561} ζ f^{''} (ζ) + \frac{74865}{59049} f^{'} (ζ)),

With the ICs: $η^{(k)} (0) = {(- 1)}^{κ}$ for $κ = 0, 1, 2, 3, 4, 5 .$

Example 5.2:

Non-Linear FDE³¹

η^{(6 α)} (ζ) = η (ζ) η^{'} (ζ) - 3 η^{(5)} (ζ) - 11 η^{'''} (ζ); 0 < α < 1; 0 < ζ \leq 1

.

With the ICs: $η^{(κα)} (0) = {(- 1)}^{κ}$ for $κ = 0, 1, 2, 3, 4, 5 .$

Using the formula in Equation (24), we get the following ODE

ζ^{5} η^{(6)} (ζ) + 15 (1 - α) ζ^{4} η^{(5)} (ζ) + 5 (1 - α) (13 - 7 α) ζ^{3} η^{(4)} (ζ) + (1 - α) ((15 (5 - 7 α) + 2 (1 - 2 α) (7 - 11 α))) ζ^{2} η^{(3)} (ζ) + 2 (1 - α) (1 - 2 α) (7 - 11 α) (2 - 5 α) ζ η^{''} (ζ) + (1 - α) (1 - 2 α) (1 - 3 α) (1 - 4 α) (1 - 5 α) η^{'} (ζ) = ζ^{6 α - 1} (η (ζ) η^{'} (ζ) - 3 η^{(5)} (ζ) - 11 η^{'''} (ζ))

with the ICs:

η^{(k)} (0) = {(- 1)}^{κ}

for

κ = 0, 1, 2, 3, 4, 5 .

Figure 1. A comparison of the numerical solutions for (A) Example 5.1 and (B) Example 5.2 using the proposed GIBM2P method with RKM method.

6. Discussion and Conclusion

The two objectives of this article are to convert the FDE to an ODE and to derive a direct-implicit block numerical method for solving a class of ODE. The current paper developed and analyzed the method of converting a sixth-order quasi-linear FDE to an ODE and then solving this ordinary differential equation using a novel direct implicit block technique with two points, which is named GIBM2P. Accordingly, we can conclude that the proposed method GIBM2P is powerful and efficient. From the discussion of the results of this paper of the two examples of implementation, we may conclude that the proposed approach performs efficiently and accurately.

Data availability

No data associated with this article.

Acknowledgment

Our researcher extends his Sincere thanks to the editor and members of the preparatory committee of the journal.

References

1. Farlow SJ: Partial differential equations for scientists and engineers. Courier Corporation; 1993.
2. Twizell E, Boutayeb A: Numerical methods for the solution of special and general sixth-order boundary-value problems, with applications to bénard layer eigenvalue problems. Proc. R. Soc. Lond. A. 1990; 431(1883): 433–450.
3. Boutayeb A, Twizell E: Numerical methods for the solution of special sixth-order boundary-value problems. Int. J. Comput. Math. 1992; 45(3-4): 207–223. Publisher Full Text
4. Twizell E: Numerical methods for sixth-order boundary value problems. Numerical Mathematics Singapore. Springer; 1988; pp. 495–506.
5. Wazwaz A-M: The numerical solution of sixth-order boundary value problems by the modified decomposition method. Appl. Math. Comput. 2001; 118(2-3): 311–325. Elsevier Publisher. Publisher Full Text
6. Akram G, Siddiqi SS: Solution of sixth order boundary value problems using nonpolynomial spline technique. Appl. Math. Comput. 2006; 181(1): 708–720. Publisher Full Text
7. Mechee MS: Generalized RK integrators for solving class of sixth-order ordinary differential equations. J. Interdiscip. Math. 2019; 22(8): 1457–1461. Publisher Full Text
8. Mechee MS, Rajihy Y: Generalized RK integrators for solving ordinary differential equations: A survey & comparison study. Glob. J. Pure Appl. Math. 2017; 13(7): 2923–2949.
9. Allogmany R, Ismail F: Implicit three-point block numerical algorithm for solving third order initial value problem directly with applications. Mathematics. 2020; 8(10): 1771.
10. Turki MY, Ismail F, Senu N, et al.: Second derivative multistep method for solving first-order ordinary differential equations. AIP Conf. Proc. 2016; 1739(1): 2020. AIP Publishing LLC.
11. Turki MY, Salih MM, Mechee MS: Construction of General Implicit-Block Method with Three-Points for Solving Seventh-Order Ordinary Differential Equations. Symmetry. 2022; 14(8): 1605. MDPI publisher. Publisher Full Text
12. Allogmany R, Ismail F: Direct solution of u” = f (t, u, u’) using three-point block method of order eight with applications. J. King Saud Univ., Sci. 2021; 33(2): 101337. Publisher Full Text
13. Ramos H, Rufai, Ajani M: A third-derivative two-step block Falkner-type method for solving general second-order boundary-value systems. Math. Comput. Simul. 2019; 165: 139–155. Elsevier Publisher.
14. Ramos H, Mehta S, Vigo-Aguiar: A unified approach for the development of k-step block Falkner-type methods for solving general second-order initial-value problems in ODEs. J. Comput. Appl. Math. 2017; 318: 550–564. Elsevier Publisher. Publisher Full Text
15. Mechee MS, Senu N: A new numerical multistep method for solution of second order of ordinary differential equations. Asian Transactions on Science Technology. 2012; 02(02): 18–22.
16. Zainal, SITI ZHAFIRAH: Fifth Order 2-point Implicit Block Method with and Off-stage Function for Solving First Order Stiff Initial Value Problems. UPM; 2014. Diss. Master Thesis.
17. Noor NM: Extended two-point and three-point block backward differentiation formulas for solving first order stiff ordinary differential equations.2019.
18. Musa H, et al.: A new fifth order implicit block method for solving first order stiff ordinary differential equations. Malays. J. Math. Sci. 2014; 8: 45–59.
19. Majid A, Zanariah NZM, Suleiman M: Direct Two‐Point Block One‐Step Method for Solving General Second‐Order Ordinary Differential Equations. Math. Probl. Eng. 2012; 1(184253): 2012.
20. Allogmany R, Ismail F, Ibrahim ZB: Implicit two-point block method with third and fourth derivatives for solving general second order ODEs. Math. Stat. 2019; 7(4): 116–123.
21. Oghonyon JG, Okunuga SA, Bishop SA: A 5-step block predictor and 4-step corrector methods for solving general second order ordinary differential equations. Glob. J. Pure Appl. Math. 2015; 11(5): 3386–3847.
22. Omar Z, Kuboye JO: Computation of an accurate implicit block method for solving third order ordinary differential equations directly. Glob. J. Pure Appl. Math. 11(1): 177–186. 20215.
23. Mansor KH, Omar Z: New Strategy of Developing a Direct Two-Step Implicit Hybrid Block Method with Generalized Two Off-Step Points for Third Order Ordinary Differential Equations. Int. J. Interdiscip. Glob. Stud. 2018; 12(1): 1–11.
24. Atabo VO, Adee SO: A new special 15-step block method for solving general fourth order ordinary differential equations. J. Niger. Soc. Phys. Sci. 2021; 308–333. Publisher Full Text
25. Turki MY, Alku SY, Mechee MS: The general implicit-block method with two-points and extra derivatives for solving fifth-order ordinary differential equations. Int. J. Nonlinear Anal. Appl. 2022; 13(1): 1081–1097.
26. Turki MY, Salih MM, Mechee MS: Construction of General Implicit-Block Method with Three-Points for Solving Seventh-Order Ordinary Differential Equations. Symmetry. 2022; 14(8): 1605. Publisher Full Text
27. Majid ZA, Suleiman MB: Variable step variable order two-point block fully implicit method for solving ordinary differential equations. Editorial Advisory Board e. 2005; 21(3): 521–529.
28. Mehrkanoon S: A direct variable step block multistep method for solving general third-order ODEs. Numer. Algorithms. 2011; 57(1): 53–66. Publisher Full Text
29. Ismail NI, Naqiah ZAM, Senu N: Explicit multistep block method for solving neutral delay differential equation. ASM Science Journal (IQRAC2018). 2019; 12(1): 24–32.
30. Ismail N: Three-point diagonally implicit block methods for solving ordinary and fuzzy differential equations. Science. 2014.
31. Singh H, Kumar D, Baleanu D: Methods of mathematical modelling: fractional differential equations. CRC Press; 2019.
32. Mechee MS, Salman ZI, Alkufi SY: Novel Definitions of α -Fractional Integral and Derivative of the Functions.2023; 64(11): 5866–5877.
33. Khalil R, Al Horani M, Yousef A, et al.: A new definition of fractional derivative. J. Comput. Appl. Math. 2014; 264: 65–70. Publisher Full Text
34. Zheng Z, Zhao W, Dai H: A new definition of fractional derivative. Int. J. Non-Linear Mech. 2019; 108: 1–6. Publisher Full Text

Comments on this article Comments (0)

Version 1

VERSION 1 PUBLISHED 24 Dec 2025

Author details Author details

¹ TRDC, university of Kufa, Najaf, Najaf-Iraq, 00964, Iraq
² Faculty of Information Technology Ninevah, Ninevah University, Mosel, 00964, Iraq

Mohammed S. Mechee
Roles: Formal Analysis, Methodology, Software, Supervision, Validation

Mohammed Mahmood Salih
Roles: Conceptualization, Data Curation, Formal Analysis, Funding Acquisition, Resources, Software, Writing – Original Draft Preparation, Writing – Review & Editing

Competing interests

No competing interests were disclosed.

Grant information

The author(s) declared that no grants were involved in supporting this work.

Article Versions (1)

version 1

Published: 24 Dec 2025, 14:1444

https://doi.org/10.12688/f1000research.172397.1

Copyright

© 2025 Mechee MS and Mahmood Salih M. 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 work is properly cited.

Download

Export To

metrics

	Views	Downloads
F1000Research	-	-
PubMed Central Data from PMC are received and updated monthly.	-	-

Citations

0

SEE MORE DETAILS

CITE

how to cite this article

Mechee MS and Mahmood Salih M. Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order Fractional Differential Equations [version 1; peer review: awaiting peer review]. F1000Research 2025, 14:1444 (https://doi.org/10.12688/f1000research.172397.1)

NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.

track

receive updates on this article

Track an article to receive email alerts on any updates to this article.

Open Peer Review

Comments on this article Comments (0)

Version 1

VERSION 1 PUBLISHED 24 Dec 2025

Open Peer Review

Reviewer Status

AWAITING PEER REVIEW

Comments on this article

All Comments(0)

Add a comment

Sign up for content alerts

Browse by related subjects

[1] 1. Farlow SJ: Partial differential equations for scientists and engineers. Courier Corporation; 1993.

[2] 2. Twizell E, Boutayeb A: Numerical methods for the solution of special and general sixth-order boundary-value problems, with applications to bénard layer eigenvalue problems. Proc. R. Soc. Lond. A. 1990; 431(1883): 433–450.

[3] 3. Boutayeb A, Twizell E: Numerical methods for the solution of special sixth-order boundary-value problems. Int. J. Comput. Math. 1992; 45(3-4): 207–223. Publisher Full Text

[4] 4. Twizell E: Numerical methods for sixth-order boundary value problems. Numerical Mathematics Singapore. Springer; 1988; pp. 495–506.

[5] 5. Wazwaz A-M: The numerical solution of sixth-order boundary value problems by the modified decomposition method. Appl. Math. Comput. 2001; 118(2-3): 311–325. Elsevier Publisher. Publisher Full Text

[6] 6. Akram G, Siddiqi SS: Solution of sixth order boundary value problems using nonpolynomial spline technique. Appl. Math. Comput. 2006; 181(1): 708–720. Publisher Full Text

[7] 7. Mechee MS: Generalized RK integrators for solving class of sixth-order ordinary differential equations. J. Interdiscip. Math. 2019; 22(8): 1457–1461. Publisher Full Text

[8] 8. Mechee MS, Rajihy Y: Generalized RK integrators for solving ordinary differential equations: A survey & comparison study. Glob. J. Pure Appl. Math. 2017; 13(7): 2923–2949.

[9] 9. Allogmany R, Ismail F: Implicit three-point block numerical algorithm for solving third order initial value problem directly with applications. Mathematics. 2020; 8(10): 1771.

[10] 10. Turki MY, Ismail F, Senu N, et al.: Second derivative multistep method for solving first-order ordinary differential equations. AIP Conf. Proc. 2016; 1739(1): 2020. AIP Publishing LLC.

[11] 11. Turki MY, Salih MM, Mechee MS: Construction of General Implicit-Block Method with Three-Points for Solving Seventh-Order Ordinary Differential Equations. Symmetry. 2022; 14(8): 1605. MDPI publisher. Publisher Full Text

[12] 12. Allogmany R, Ismail F: Direct solution of u” = f (t, u, u’) using three-point block method of order eight with applications. J. King Saud Univ., Sci. 2021; 33(2): 101337. Publisher Full Text

[13] 13. Ramos H, Rufai, Ajani M: A third-derivative two-step block Falkner-type method for solving general second-order boundary-value systems. Math. Comput. Simul. 2019; 165: 139–155. Elsevier Publisher.

[14] 14. Ramos H, Mehta S, Vigo-Aguiar: A unified approach for the development of k-step block Falkner-type methods for solving general second-order initial-value problems in ODEs. J. Comput. Appl. Math. 2017; 318: 550–564. Elsevier Publisher. Publisher Full Text

[15] 15. Mechee MS, Senu N: A new numerical multistep method for solution of second order of ordinary differential equations. Asian Transactions on Science Technology. 2012; 02(02): 18–22.

[16] 16. Zainal, SITI ZHAFIRAH: Fifth Order 2-point Implicit Block Method with and Off-stage Function for Solving First Order Stiff Initial Value Problems. UPM; 2014. Diss. Master Thesis.

[17] 17. Noor NM: Extended two-point and three-point block backward differentiation formulas for solving first order stiff ordinary differential equations.2019.

[18] 18. Musa H, et al.: A new fifth order implicit block method for solving first order stiff ordinary differential equations. Malays. J. Math. Sci. 2014; 8: 45–59.

[19] 19. Majid A, Zanariah NZM, Suleiman M: Direct Two‐Point Block One‐Step Method for Solving General Second‐Order Ordinary Differential Equations. Math. Probl. Eng. 2012; 1(184253): 2012.

[20] 20. Allogmany R, Ismail F, Ibrahim ZB: Implicit two-point block method with third and fourth derivatives for solving general second order ODEs. Math. Stat. 2019; 7(4): 116–123.

[21] 21. Oghonyon JG, Okunuga SA, Bishop SA: A 5-step block predictor and 4-step corrector methods for solving general second order ordinary differential equations. Glob. J. Pure Appl. Math. 2015; 11(5): 3386–3847.

[22] 22. Omar Z, Kuboye JO: Computation of an accurate implicit block method for solving third order ordinary differential equations directly. Glob. J. Pure Appl. Math. 11(1): 177–186. 20215.

[23] 23. Mansor KH, Omar Z: New Strategy of Developing a Direct Two-Step Implicit Hybrid Block Method with Generalized Two Off-Step Points for Third Order Ordinary Differential Equations. Int. J. Interdiscip. Glob. Stud. 2018; 12(1): 1–11.

[24] 24. Atabo VO, Adee SO: A new special 15-step block method for solving general fourth order ordinary differential equations. J. Niger. Soc. Phys. Sci. 2021; 308–333. Publisher Full Text

[25] 25. Turki MY, Alku SY, Mechee MS: The general implicit-block method with two-points and extra derivatives for solving fifth-order ordinary differential equations. Int. J. Nonlinear Anal. Appl. 2022; 13(1): 1081–1097.

[26] 26. Turki MY, Salih MM, Mechee MS: Construction of General Implicit-Block Method with Three-Points for Solving Seventh-Order Ordinary Differential Equations. Symmetry. 2022; 14(8): 1605. Publisher Full Text

[27] 27. Majid ZA, Suleiman MB: Variable step variable order two-point block fully implicit method for solving ordinary differential equations. Editorial Advisory Board e. 2005; 21(3): 521–529.

[28] 28. Mehrkanoon S: A direct variable step block multistep method for solving general third-order ODEs. Numer. Algorithms. 2011; 57(1): 53–66. Publisher Full Text

[29] 29. Ismail NI, Naqiah ZAM, Senu N: Explicit multistep block method for solving neutral delay differential equation. ASM Science Journal (IQRAC2018). 2019; 12(1): 24–32.

[30] 30. Ismail N: Three-point diagonally implicit block methods for solving ordinary and fuzzy differential equations. Science. 2014.

[31] 31. Singh H, Kumar D, Baleanu D: Methods of mathematical modelling: fractional differential equations. CRC Press; 2019.

[32] 32. Mechee MS, Salman ZI, Alkufi SY: Novel Definitions of α -Fractional Integral and Derivative of the Functions.2023; 64(11): 5866–5877.

[33] 33. Khalil R, Al Horani M, Yousef A, et al.: A new definition of fractional derivative. J. Comput. Appl. Math. 2014; 264: 65–70. Publisher Full Text

[34] 34. Zheng Z, Zhao W, Dai H: A new definition of fractional derivative. Int. J. Non-Linear Mech. 2019; 108: 1–6. Publisher Full Text

Construction an Implicit Block Multi-Steps Approach for Solving Sixth-Order Fractional Differential Equations

Abstract

Background

Methods

Results

Conclusions

Keywords

1. Introduction

2. Preliminary and background

2.1 Quasi-Linear Sixth-Order Ordinary Differential Equations

(1)

(2)

(3)

Definition 2.1:

(4)

(5)

Definition 2.2:

(6)

2.1 Fractional derivatives

Definition 2.3:

(7)

Property 2.1

(8)

3. Proposed method

3.1 Analysis of the proposed numerical GIBM2P method

(9)

(10)

(11)

(12)

(13)

(14)

(15)

(16)

(17)

(18)

(19)

(20)

(21)

(22)

(23)

4. Main results

4.1 Analysis of proposed approach for solving of FDEs sixth-order

(24)

4.2 Proposed analytical method

5. Implementation

Example 5.1:

Example 5.2:

Figure 1. A comparison of the numerical solutions for (A) Example 5.1 and (B) Example 5.2 using the proposed GIBM2P method with RKM method.

6. Discussion and Conclusion

Data availability

Acknowledgment

References

Comments on this article Comments (0)

Open Peer Review

Comments on this article Comments (0)

Open Peer Review

Reviewer Status

Comments on this article

Browse by related subjects

Competing Interests Policy

Stay Updated