am Applied Mathematics 2152-7393 2152-7385 Scientific Research Publishing 10.4236/am.2026.176020 am-151898 Article Physics Mathematics A Decomposition Method via a Fractional Inverse Operator for Fractional Boundary Value Problems Alsulami Mona 1 Al-Mazmumy Mariam 1 Banweer Bashair 1 Department of Mathematics and Statistics, College of Science, University of Jeddah, Jeddah, Saudi Arabia

The authors declare no conflicts of interest regarding the publication of this paper.

 16 06 2026 06 2026 17 06 318 336 29 04 2026 13 06 2026 16 06 2026 © 2026 by the authors and Scientific Research Publishing Inc. 2026 This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( https://creativecommons.org/licenses/by/4.0/ ). https://doi.org/10.4236/am.2026.176020

The current paper numerically tackles the class of fractional ordinary boundary value problems using the modification method of the classical Adomian decomposition method and through a fractional inverse operator. Among the advantages of the adopted fractional inverse operator is that the operator directly incorporates all boundary conditions without introducing unknown constants, thereby making the approximate solution valid and preferred, with the satisfaction of the involved boundary conditions. Moreover, the results of the study show that the proposed method provides exact solutions for problems that have available exact solutions, and on the other hand, gives highly accurate approximate solutions in comparison with other computational techniques, confirming its efficiency and reliability. Indeed, several comparison tables and figures are provided, assessing the effectiveness of the proposed method over existing approaches in the literature.

 Adomian Decomposition Method Fractional Boundary Value Problems Fractional Inverse Operator Dirichlet Boundary Conditions 1. Introduction

Non-integer order calculus is a generalization of classical calculus by extending the notion of integration and differentiation to fractional orders. This field, in recent decades has been of increasing concern owing to its efficiency in perfectly describing a wide range of applications, where classical calculus fails or gives half-backed information. Thus, applications of the new concept have emerged in a wide range of fields in engineering and science, including electromagnetism [1], control engineering [2][3], bio-sciences [4], fluid dynamics and water waves [5], electrochemistry [6], thermodynamical processes [7], continuum and statistical mechanics [8], as well as in viscoelasticity [9]. Building on this foundation, fractional differential equations (FDEs) have emerged as a powerful tool in mathematical modeling. Their significance lies in their ability to accurately represent a variety of physical, scientific, and engineering processes, which makes them highly valuable in many practical and industrial applications, including physics, dynamics, mechanics, and other applied sciences [10]-[13]. In this study, we consider fractional ordinary differential equations with Dirichlet boundary conditions, commonly known as fractional boundary value problems (FBVPs) have been significantly explored in the recent times due to their relevance in modelling different phenomena in many fields such as the modeling of viscoelastic materials, wave propagation in porous media, image processing, various problems in fluid mechanics, control systems, and electrical circuits. Also applied in chemical and biological systems, FBVPs have been used in chemical graph models to represent molecular structures such as cyclopentasilane [14] and enable modeling of nonlocal and multi-dimensional interactions on molecular graphs [15][16], demonstrating their practical significance and the need for effective analytical and numerical solution methods. To solve such problems, scientists and mathematicians developed several methods for obtaining both the approximate closed-form and numerical solutions, including the homotopy analysis method (HAM) [17], homotopy perturbation method (HPM) [18], variational iteration method (VIM) [19], the generalized differential transform method [20], shifted Legendre approximation method [21], shooting method [22], the Chebyshev spectral method [23], the collocation-shooting method [24], the monotone iterative sequences method [25], and the piecewise polynomial collocation approach [26]. Moreover, among these approaches, the Adomian decomposition method (ADM) and its variants [27]-[30], have attracted significant interest. Recently, considerable attention has been devoted to the ADM for solving a wide range of deterministic and stochastic differential equations [31]-[36]. The strength of the ADM lies in its ability to provide analytical approximations for broad classes of differential equations, without resorting to either of discretization, perturbation, linearization, or closure approximations, which often lead to extensive numerical computations.

However, the purpose of the present study is to focus on the ADM and its variants for treating the class of FBVPs for ordinary differential equations, by adapting a fractional inverse operator that directly incorporates all boundary conditions, eliminating the need to introduce unknown constants, in contrast to earlier modified ADM approaches. On the other hand, unlike in [37], where the inverse operator is employed to construct two separate formulations for the solution using the initial and boundary conditions independently, which are then combined through an averaging process to derive the recursive relations, the present work focuses only on the boundary conditions in constructing the recursive scheme. This yields a simpler decomposition procedure with fewer computational steps. The convergence of the Adomian decomposition series and the error analysis for fractional differential equations have been investigated in [38]-[41]. The significance of this work lies in its ability to provide exact solutions for certain problems while producing high-quality approximate solutions for such problems that have no known exact solutions. Additionally, the range of studies applying the classical Adomian method directly to such FBVPs remains limited, let alone the incorporation of the unexplored fractional inversion operator, highlighting the need for further research in this area. The current study is restricted to Caputo fractional ordinary boundary value problems with Dirichlet boundary conditions, and the illustrative examples mainly involve orders 0 < α ≤ 1 and 1 < α ≤ 2 . Thus, the present work is arranged in the following way: Section 2 presents some preliminaries on non-integer order calculus. Section 3 describes the ADM and its modified form suitable for solving FBVPs. Section 4 gives several illustrative examples, highlighting the efficiency and accuracy of the proposed method. Finally, Section 5 summarizes the main results and conclusions.

 2. Fractional Derivatives and Fractional Integrals

This section presents some preliminaries on the non-integer order calculus, which includes the definitions of the Riemann-Liouville (RL) and Caputo fractional integrals and derivatives, and other definitions and properties; see the book by Kilbas et al. [42] for more on the theory of non-integer calculus.

Definition2.1. Assume [ a , b ] to be a finite interval in ℝ . Then, the RL fractional integral for the function f ( x ) of order α ( > 0 ) ∈ ℝ for x > a , and x < b are respectively defined as

( I a + α f ) ( x ) = 1 Γ ( α ) ∫ a x ( x − t ) α − 1 f ( t ) d t ,

and

( I b − α f ) ( x ) = 1 Γ ( α ) ∫ x b ( t − x ) α − 1 f ( t ) d t .

Definition2.2. The RL fractional derivative for the function f ( x ) of order α ( ≥ 0 ) ∈ ℝ for x > a , and x < b when n = [ α ] + 1 are respectively defined as

( D a + α f ) ( x ) = ( d d x ) n ( I a + n − α f ) ( x ) = 1 Γ ( n − α ) ( d d x ) n ∫ a x ( x − t ) n − α − 1 f ( t ) d t ,

and

( D b − α f ) ( x ) = ( − d d x ) n ( I b − n − α f ) ( x ) = 1 Γ ( n − α ) ( − d d x ) n ∫ x b ( t − x ) n − α − 1 f ( t ) d t ,

where [ α ] denotes the greatest integer value for the fractional-order α .

Definition2.3. Assume [ a , b ] is a finite interval in ℝ . The fractional derivatives ( c D a + α f ) ( x ) and ( c D b − α f ) ( x ) of order α ∈ ℝ + in [ a , b ] via (3) and (4) are the left-sided, and right-sided Caputo derivatives of order α , are re-expressed respectively defined as

( c D a + α f ) ( x ) = ( D a + α [ f ( t ) − ∑ j = 0 n − 1 f ( j ) ( a ) j ! ( t − a ) j ] ) ( x ) ,

and

( c D b − α f ) ( x ) = ( D b − α [ f ( t ) − ∑ j = 0 n − 1 f ( j ) ( b ) j ! ( b − t ) j ] ) ( x ) , n = [ α ] + 1

Lemma2.1. Assume α ∈ ℝ + , and n = [ α ] + 1 for α ∉ ℕ 0 ; n = α for α ∈ ℕ 0 . Then, if f ( x ) ∈ C n [ a , b ] or f ( x ) ∈ A C n [ a , b ] , one gets

( I a + α c D a + α f ) ( x ) = f ( x ) − ∑ j = 0 n − 1 f ( j ) ( a ) j ! ( x − a ) j , ( I b − α c D b − α f ) ( x ) = f ( x ) − ∑ j = 0 n − 1 ( − 1 ) j f ( j ) ( b ) j ! ( b − x ) j .

Theorem2.1. [43] For α ∈ ℝ + ( n − 1 < α ≤ n ) , and f ( x ) ∈ C − 1 n , n ∈ N , then

( I α D α f ) ( x ) = f ( x ) − ∑ k = 0 n − 1   f ( j ) ( 0 + ) x j j ! .

Property2.1. For the power function f ( x ) = x m , m > − 1 , x > 0 and α > 0 , one gets

I α x m = Γ ( 1 + m ) Γ ( 1 + m + α ) x m + α .

Property2.2. The Caputo fractional derivative of the function f ( x ) = x m for m > − 1 and n − 1 < α ≤ n satisfies

c D a + α x m = { 0 , for   m < ⌈ α ⌉ , Γ ( 1 + m ) Γ ( 1 + m − α ) x m − α , for   m ≥ ⌈ α ⌉ . 3. Description of the Proposed Method

In proposing the new method for solving the class of FBVPs, the following nonlinear FDE is considered

D α u ( x ) + R u ( x ) + N u ( x ) = g ( x ) ,     a < x < b ,

subject to the following Dirichlet boundary conditions

u ( a ) = γ 1 ,     u ( b ) = γ 2 ,

where D α is the fractional derivative operator in the sense of Caputo, R is a linear differential operator, which may equally feature other fractional operators of order less than α . At the same time, N is an operator that is nonlinear, and g ( x ) is the nonhomogeneous term, also referred to as the source function.

What is more, the governing fractional integral operator I α of much concern in this study has the form [37]

I α ( . ) = ∫ a x ( x − s ) α − 1 Γ ( α ) ( . ) d s − ( x − a b − a ) ∫ a b ( b − s ) α − 1 Γ ( α ) ( . ) d s , = I α ( . ) − ( x − a b − a ) I α ( . ) x = b .

Accordingly, upon applying the inverse operator (14) to both sides of (12), one gets

I α [ D α u ( x ) ] + I α [ R u ( x ) ] + I α [ N u ( x ) ] = I α [ g ( x ) ] ,

such that upon using Lemma 2.1 and the boundary conditions (13) the latter equation becomes

u ( x ) − ∑ k = 0 n − 1   u ( k ) ( a ) ( x − a ) k k ! − ( x − a b − a ) [ u ( b ) − ∑ k = 0 n − 1   u ( k ) ( a ) ( b − a ) k k ! ] +   I α [ R u ( x ) ] + I α [ N u ( x ) ] = I α [ g ( x ) ] ,

or equally,

u ( x ) = u ( a ) + ( x − a b − a ) [ u ( b ) − u ( a ) ] − I α [ R u ( x ) ] − I α [ N u ( x ) ] + I α [ g ( x ) ] .

Next, expressing the solution u ( x ) using an infinite series of the form:

u ( x ) = ∑ n = 0 ∞   u n ( x ) ,

and expanding the nonlinear term N u ( x ) via the use of Adomian polynomials [44]:

N ( u ) = ∑ n = 0 ∞   A n ( u 0 , u 1 , ⋯ , u n ) ,

where the polynomials are plainly determined using the relation

A n = 1 n ! d n d λ n [ N ( ∑ i = 0 n   λ i u i ( x ) ) ] λ = 0 ,     n = 0 , 1 , 2 , ⋯

one thus substitutes the latter series expressions in (16) and (17) into (15) to get

∑ n = 0 ∞   u n ( x ) = u ( a ) + ( x − a b − a ) [ u ( b ) − u ( a ) ] + I α [ g ( x ) ]     − I α R [ ∑ n = 0 ∞   u n ( x ) ] − I α [ ∑ n = 0 ∞   A n ] ,

which thus leads to the resulting recursive scheme as follows

{ u 0 ( x ) = h ( x ) , u n + 1 ( x ) = − I α [ R u n ( x ) ] − I α [ A n ] ,     n ≥ 0.

where

h ( x ) = u ( a ) + ( x − a b − a ) [ u ( b ) − u ( a ) ] + I α [ g ( x ) ] .

Notably, on spitting the latter expression for h , that is, letting h = h 1 + h 2 , the reliable modification of the ADM [29] thus recasts the recursive scheme in (19) to the following

{ u 0 ( x ) = h 1 , u 1 ( x ) = h 2 − I α [ R u 0 ( x ) ] − I α [ A 0 ] , u n + 1 ( x ) = − I α [ R u n ( x ) ] − I α [ A n ] ,     n ≥ 1.

Hence, the overall solution is obtained via the infinite series in (16) by substituting the determined components u 0 , u 1 , u 2 , ⋯ . In fact, truncation yields an approximate solution, while summation - when possible - leads to a closed form.

 4. Applications

The present section shows the efficiency of the proposed new decomposition method through the application of an inverse fractional operator. Several illustrative examples will be extensively examined, portraying Dirichlet boundary conditions and endowed with the Caputo fractional derivative. For each example, the obtained numerical solution will be benchmarked against the exact solution, if available, by computing the consequent absolute error. Moreover, a comparative analysis will be established with the existing numerical methods available in the literature.

Example4.1. Consider the FBVP through the Bagley-Torvik equation

D 3 / 2 u ( x ) + u ″ ( x ) + u ( x ) = 1 + x ,     x ∈ [ 0 , 1 ] , u ( 0 ) = 1 ,     u ( 1 ) = 2 ,

that satisfies the exact solution u ( x ) = 1 + x .

First, with the help of an operator notation, we write (21) as follows

D 3 / 2 u ( x ) = 1 + x − u ″ ( x ) − u ( x ) .

Then, with the application of the inversion operator (14) on the both side of (23) using Theorem 2.1 and boundary conditions (22), one gets

u ( x ) = 1 + x + I 3 / 2 [ 1 + x ] − I 3 / 2 [ u ″ ( x ) + u ( x ) ] ,

where, from the latter equation, the proposed method yields the recursive scheme as follows

{ u 0 ( x ) = 1 + x , u 1 ( x ) = I 3 / 2 [ 1 + x ] − I 3 / 2 [ u ″ 0 ( x ) + u 0 ( x ) ] , u n + 1 ( x ) = I 3 / 2 [ u ″ n ( x ) + u n ( x ) ] ,     n ≥ 1.

In particular, from the above scheme, computing the value of u 1 ( x ) by using I 3 / 2 from (14), one gets

u 1 ( x ) = I 3 / 2 [ 1 + x ] − I 3 / 2 [ u ″ 0 ( x ) + u 0 ( x ) ] = I 3 / 2 [ 1 + x ] − I 3 / 2 [ 1 + x ] = 0.

Hence, u n ( x ) = 0 ∀ n ≥ 1 . Accordingly, the closed-form solution is obtained as u ( x ) = u 0 ( x ) = 1 + x , which coincides with the reported exact solution. Moreover, Table 1 highlights the superior accuracy of the present method in comparison with the Bessel collocation method (BCM) presented in [45]. As shown, the present method reaches the exact solution with two components, demonstrating its high precision and rapid convergence.

Table 1. Comparison of absolute error for Example 4.1.

x BCM [ 45 ]( N = 6 ) Present method( N = 2 )
0.0 0 0
0.1 6.1919 × 10 − 16 0
0.2 1.0292 × 10 − 15 0
0.3 1.4779 × 10 − 15 0
0.4 1.9697 × 10 − 15 0
0.5 2.4941 × 10 − 15 0
0.6 3.0365 × 10 − 15 0
0.7 3.5805 × 10 − 15 0
0.8 4.1090 × 10 − 15 0
0.9 4.6047 × 10 − 8 0
1.0 5.0500 × 10 − 15 0

Example4.2.Consider the FBVP

D 3 / 2 u ( x ) + u ( x ) = 2 x 1 / 2 Γ ( 3 / 2 ) + x 2 − x ,     x ∈ [ 0 , 1 ] , u ( 0 ) = 0 ,     u ( 1 ) = 0 ,

which satisfies the exact solution u ( x ) = x 2 − x .

Now, one begins by re-writing (24) as follows

D 3 / 2 u ( x ) = 2 x 1 / 2 Γ ( 3 / 2 ) + x 2 − x − u ( x ) .

Next, with the application of the inversion operator I 3 / 2 from (14) on (26), one gets after using the property in (7) and the conditions (25) the following

u ( x ) = I 3 / 2 [ 2 x 1 / 2 Γ ( 3 / 2 ) + x 2 − x ] − I 3 / 2 [ u ( x ) ] ,

upon which the application of (14) and the property (10) subsequently yields

u ( x ) = x 2 + 2 Γ ( 9 / 2 ) x 7 / 2 − 1 Γ ( 7 / 2 ) x 5 / 2 − x − [ 2 Γ ( 9 / 2 ) − 1 Γ ( 7 / 2 ) ] x − I 3 / 2 [ u ( x ) ] .

Consequently, employing the proposed decomposition method reveals the overall recursive scheme as follows

{ u 0 ( x ) = x 2 − x , u 1 ( x ) = 2 Γ ( 9 / 2 ) x 7 / 2 − 1 Γ ( 7 / 2 ) x 5 / 2 − [ 2 Γ ( 9 / 2 ) − 1 Γ ( 7 / 2 ) ] x − I 3 / 2 [ u 0 ( x ) ] , u n + 1 ( x ) = − I 3 / 2 [ u n ( x ) ] ,     n ≥ 1.

Computing the value of u 1 ( x ) by using I 3 / 2 from (14) and the property in (10), one gets

u 1 ( x ) = 2 Γ ( 9 / 2 ) x 7 / 2 − 1 Γ ( 7 / 2 ) x 5 / 2 − [ 2 Γ ( 9 / 2 ) − 1 Γ ( 7 / 2 ) ] x − I 3 / 2 [ x 2 − x ] , = 2 Γ ( 9 / 2 ) x 7 / 2 − 1 Γ ( 7 / 2 ) x 5 / 2 − [ 2 Γ ( 9 / 2 ) − 1 Γ ( 7 / 2 ) ] x − 2 Γ ( 9 / 2 ) x 7 / 2     + 1 Γ ( 7 / 2 ) x 5 / 2 + [ 2 Γ ( 9 / 2 ) − 1 Γ ( 7 / 2 ) ] x , = 0.

Hence, u n ( x ) = 0 ∀ n ≥ 1 . Moreover, one gets the solution by summing the components as u ( x ) = u 0 ( x ) = x 2 − x , coinciding with the already reported exact solution. What is more, Table 2 shows that the present technique with two components produces zero absolute error for all values of x , demonstrating that the method attains the exact solution. Compared with Chelyshkov operational matrix method (COMM) [46], Pseudo spectral method (PSM) [47], and Generalized Adams-Bashforth-Moulton method (GABMM) [48], which show small numerical errors. Indeed, the present method exhibits superior accuracy and fast convergence. For GABMM, h denotes the step size used in the simulation.

Table 2. Comparison of absolute error for Example 4.2.

x COMM [ 46 ]( N = 10 ) PSM [ 47 ]( N = 10 ) GABMM [ 48 ]( N = 6 ) h = 1 / 20 GABMM [ 48 ]( N = 6 ) h = 1 / 320 Present method( N = 2 )
0.0 6.204 × 10 − 5 - - - 0
0.1 9.450 × 10 − 6 - - - 0
0.2 8.600 × 10 − 7 - - - 0
0.3 1.261 × 10 − 5 - - - 0
0.4 1.273 × 10 − 5 - - - 0
0.5 5.170 × 10 − 6 - - - 0
0.6 1.224 × 10 − 6 - - - 0
0.7 6.27 × 10 − 6 - - - 0
0.8 8.13 × 10 − 6 - - - 0
0.9 9.30 × 10 − 6 - - - 0
1.0 2.379 × 10 − 5 6.83 × 10 − 4 3.42 × 10 − 3 5.71 × 10 − 5 0

Example4.3. Consider the FBVP

D 2 α u ( x ) = u ( x ) + x ,     0 < α ≤ 1 ,     x ∈ [ 0 , 1 ] , u ( 0 ) = 0 ,     u ( 1 ) = 0 ,

that satisfies the exact solution u ( x ) = sinh ( x ) sinh ( 1 ) − x while assuming the integer order unity.

In the same way, applying the inversion operator I 2 α from (14) to (27), and upon using properties in (7) and (10) through the conditions (28), one gets

u ( x ) = 1 Γ ( 2 + 2 α ) x 1 + 2 α − 1 Γ ( 2 + 2 α ) x + I 2 α [ u ( x ) ] ,

that yields the consequent ADM recursive scheme as follows

{ u 0 ( x ) = 1 Γ ( 2 + 2 α ) x 1 + 2 α − 1 Γ ( 2 + 2 α ) x , u n + 1 ( x ) = I 2 α [ u n ( x ) ] ,     n ≥ 0.

Subsequently, one computes some of the successive iterates as

u 1 ( x ) = I 2 α [ u 0 ( x ) ] = 1 Γ ( 2 + 4 α ) x 1 + 4 α − 1 [ Γ ( 2 + 2 α ) ] 2 x 1 + 2 α     − [ 1 Γ ( 2 + 4 α ) − 1 [ Γ ( 2 + 2 α ) ] 2 ] x , u 2 ( x ) = I 2 α [ u 1 ( x ) ] = 1 Γ ( 2 + 6 α ) x 1 + 6 α − 1 Γ ( 2 + 2 α ) Γ ( 2 + 4 α ) x 1 + 4 α     − [ 1 Γ ( 2 + 2 α ) Γ ( 2 + 4 α ) − 1 [ Γ ( 2 + 2 α ) ] 3 ] x 1 + 2 α     − [ 1 Γ ( 2 + 6 α ) − 2 Γ ( 2 + 2 α ) Γ ( 2 + 4 α ) + 1 [ Γ ( 2 + 2 α ) ] 3 ] x , u 3 ( x ) = I 2 α [ u 2 ( x ) ] = 1 Γ ( 2 + 8 α ) x 1 + 8 α − 1 Γ ( 2 + 2 α ) Γ ( 2 + 8 α ) x 1 + 6 α     − [ 1 [ Γ ( 2 + 4 α ) ] 2 − 1 [ Γ ( 2 + 2 α ) ] 2 Γ ( 2 + 4 α ) ] x 1 + 4 α     − [ 1 Γ ( 2 + 6 α ) Γ ( 2 + 2 α ) − 2 [ Γ ( 2 + 2 α ) ] 2 Γ ( 2 + 4 α ) + 1 [ Γ ( 2 + 2 α ) ] 4 ] x 1 + 2 α     − [ 1 Γ ( 2 + 8 α ) − 1 Γ ( 2 + 2 α ) Γ ( 2 + 8 α ) − 1 [ Γ ( 2 + 4 α ) ] 2     + 3 [ Γ ( 2 + 2 α ) ] 2 Γ ( 2 + 4 α ) − 1 Γ ( 2 + 6 α ) Γ ( 2 + 2 α ) + 1 [ Γ ( 2 + 2 α ) ] 4 ] x , ⋮

Accordingly, the approximate solution follows from the summation in (16). What is more, Figure 1 illustrates the nature of the solution for different fractional orders α . It is observed that when α = 1 , the reported exact solution is attained. In addition, Table 3 presents the absolute error values, where the present method consistently yields smaller errors than the method in [49], proving its superior accuracy and efficiency in solving FBVPs.

Figure 1. Exact and ADM solutions of Example 4.3 for various values of α = 0.6 , 0.8 , 1 .

Table 3. Comparison of absolute error for Example 4.3.

x FCDM [ 49 ] Present method
α = 0.6 α = 0.8 α = 1 α = 0.6 α = 0.8 α = 1
0.0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.1 2.810 × 10 − 2 9.600 × 10 − 3 2.000 × 10 − 4 1.813 × 10 − 2 7.899 × 10 − 3 1.932 × 10 − 7
0.2 4.880 × 10 − 2 1.730 × 10 − 2 1.200 × 10 − 3 3.248 × 10 − 2 1.452 × 10 − 2 3.674 × 10 − 7
0.3 6.170 × 10 − 2 2.260 × 10 − 2 3.000 × 10 − 4 4.281 × 10 − 2 1.949 × 10 − 2 5.058 × 10 − 7
0.4 6.750 × 10 − 2 2.530 × 10 − 2 3.000 × 10 − 4 4.907 × 10 − 2 2.267 × 10 − 2 5.946 × 10 − 7
0.5 6.670 × 10 − 2 2.530 × 10 − 2 2.000 × 10 − 4 5.128 × 10 − 2 2.394 × 10 − 2 6.253 × 10 − 7
0.6 6.020 × 10 − 2 2.300 × 10 − 2 6.000 × 10 − 4 4.940 × 10 − 2 2.325 × 10 − 2 5.947 × 10 − 7
0.7 4.910 × 10 − 2 1.870 × 10 − 2 5.115 × 10 − 1 4.341 × 10 − 2 2.056 × 10 − 2 5.059 × 10 − 7
0.8 3.450 × 10 − 2 1.290 × 10 − 2 3.300 × 10 − 4 3.325 × 10 − 2 1.583 × 10 − 2 3.674 × 10 − 7
0.9 1.770 × 10 − 2 6.400 × 10 − 3 4.000 × 10 − 4 1.883 × 10 − 2 8.992 × 10 − 3 1.933 × 10 − 7
1.0 0.0000 0.0000 0.0000 0.0000 9.738 × 10 − 15 0.0000

Example4.4.Consider the FBVP for the nonlinear Bratuequation

D α u ( x ) + μ e u ( x ) = 0 ,     1 < α ≤ 2 ,     x ∈ [ 0 , 1 ] , u ( 0 ) = 0 ,     u ( 1 ) = 0.

The Bratu problem, introduced by Bratu in 1914 [50] has vast applications, including the fuel ignition [50], radiative heat transfer [51], and electrospinning of polymer nanofibers [52].

In this regard, the exact solution for the governing Bratu model when α = 2 takes the expression

u ( x ) = − 2 ln [ cosh ( ( x − 1 2 ) θ 2 ) cosh ( θ 4 ) ] ,

with θ satisfying θ = 2 μ cosh ( θ 4 ) . Moreover, with the discovery of the critical value μ c as μ c = 3.513830719 , and satisfying 1 = 1 4 2 μ c sinh ( θ c 4 ) , μ is thus classified as:

1) When μ < μ c , the Bratu equation has two solutions.

2) When μ = μ c , the Bratu equation has one solution.

3) When μ > μ c , the Bratu equation has a zero solution.

Accordingly, to solve the Bratu problem in (29)-(30), one re-writes the equation using the operator notation when μ = 1 as D α u ( x ) = − e u ( x ) . Next, applying the inversion operator I α in (??) to the equation and via the property (7) and the conditions (30) to get

u ( x ) = − I α [ e u ( x ) ] .

We decompose the nonlinear term N u = e u by the Adomian polynomials from (17)-(18) as follows:

A 0 = e u 0 , A 1 = e u 0 u 1 , A 2 = 1 2 e u 0 ( u 1 2 + 2 u 2 ) , ⋮

Hence, through the classical ADM procedure, the recursive scheme is obtained as follows

{ u 0 ( x ) = 0 , u n + 1 ( x ) = − I α [ A n ] ,     n ≥ 0.

Consequently, one gets some of the solution components from the latter scheme as

u 1 ( x ) = − I α [ e u 0 ] = 1 Γ ( 1 + α ) ( x − x α ) , u 2 ( x ) = − I α [ e u 0 u 1 ] = 1 Γ ( 1 + 2 α ) ( x 2 α − x ) + 1 Γ ( 1 + α ) Γ ( 2 + α ) ( x − x 1 + α ) , u 3 ( x ) = − I α [ 1 2 e u 0 ( u 1 2 + 2 u 2 ) ] = − [ Γ ( 1 + 2 α ) 2 [ Γ ( 1 + α ) ] 2 + 1 ] x 3 α Γ ( 1 + 3 α ) + [ Γ ( 2 + α ) Γ ( 1 + α ) + 1 ] x 1 + 2 α Γ ( 1 + α ) Γ ( 2 + 2 α )     + [ 1 Γ ( 1 + 2 α ) − 1 Γ ( 1 + α ) Γ ( 2 + α ) ] x 1 + α Γ ( 2 + α )     − [ 1 [ Γ ( 1 + α ) ] 2 Γ ( 3 + α ) ] x 2 + α     − [ 2 [ Γ ( 1 + α ) ] 2 + Γ ( 1 + 3 α ) 2 [ Γ ( 1 + α ) ] 2 Γ ( 1 + 3 α ) + Γ ( 2 + α ) + Γ ( 1 + α ) [ Γ ( 1 + α ) ] 2 Γ ( 2 + 2 α )     + Γ ( 1 + 2 α ) − Γ ( 2 + α ) Γ ( 1 + α ) [ Γ ( 2 + α ) ] 2 Γ ( 1 + 2 α ) Γ ( 1 + α ) + 1 [ Γ ( 1 + α ) ] 2 Γ ( 3 + α ) ] x , ⋮

Moreover, one gets the approximate solution from (16); see Figure 2 that gives the graphical illustration of the solution while responding to variation of the fractional order α when fixing μ = 1 . The solution curves for the different values of α are very close to each other, with almost perfect overlap when α = 2 . Additionally, Table4reports absolute errors comparison, showing that the present method attains higher accuracy for all values of α using only 4 components, while the reproducing kernel method (RKM) [53] employs 20 components, yet yields less accurate results. This clearly reaffirms the superiority of the proposed method over the contending method for solving nonlinear FBVPs.

Figure 2. Exact and ADM solutions of Example 4.4 for various values of α = 1.8 , 1.9 , 2 .

Table 4. Comparison of absolute error for Example 4.4 when μ = 1 .

x RKM [ 53 ] ( N = 20 ) Present method ( N = 4 )
Error Error
α = 1.8 α = 1.9 α = 2 α = 1.8 α = 1.9 α = 2
0.0 - - - 0.0000 0.0000 0.0000
0.1 6.3195 × 10 − 3 3.1488 × 10 − 3 2.0159 × 10 − 4 6.3662 × 10 − 3 3.4466 × 10 − 3 2.9236 × 10 − 5
0.2 1.5522 × 10 − 2 7.5193 × 10 − 3 2.8370 × 10 − 4 8.2881 × 10 − 3 4.7708 × 10 − 3 5.6468 × 10 − 5
0.3 2.1884 × 10 − 2 1.0915 × 10 − 2 2.9549 × 10 − 4 7.7752 × 10 − 3 4.8318 × 10 − 3 7.9000 × 10 − 5
0.4 2.6817 × 10 − 2 1.3327 × 10 − 2 2.5664 × 10 − 4 5.8428 × 10 − 3 4.0924 × 10 − 3 9.3988 × 10 − 5
0.5 2.9159 × 10 − 2 1.4594 × 10 − 2 1.9630 × 10 − 4 3.2486 × 10 − 3 2.9089 × 10 − 3 9.9257 × 10 − 5
0.6 2.9046 × 10 − 2 1.4586 × 10 − 2 1.3779 × 10 − 4 6.2844 × 10 − 4 1.5876 × 10 − 3 9.3988 × 10 − 5
0.7 2.6148 × 10 − 2 1.3186 × 10 − 2 9.1660 × 10 − 5 1.4621 × 10 − 3 4.0252 × 10 − 4 7.9000 × 10 − 5
0.8 2.0365 × 10 − 2 1.0310 × 10 − 2 5.7160 × 10 − 5 2.5348 × 10 − 3 3.9939 × 10 − 4 5.6468 × 10 − 5
0.9 1.1642 × 10 − 2 5.9153 × 10 − 3 2.8320 × 10 − 5 2.1665 × 10 − 3 3.9939 × 10 − 4 2.9236 × 10 − 5
1.0 - - - 2.0950 × 10 − 13 9.5133 × 10 − 13 2.6886 × 10 − 13

Example4.5. Consider the system of nonhomogeneous FBVPs for x ∈ [ 0 , 1 ] as follows

{ u ″ ( x ) − v ′ ( x ) + D 0.4 v ( x ) + v ( x ) = g 1 ( x ) ,                                             ( 32 ) v ″ ( x ) − u ′ ( x ) + D 0.6 u ( x ) + u ( x ) = g 2 ( x ) ,                                             ( 33 ) u ( 0 ) = 0 = u ( 1 ) ,     v ( 0 ) = 0 = v ( 1 ) ,                                             ( 34 )

where

g 1 ( x ) = 11 x 2 + 3 x − 1 − 2 x 1.6 Γ ( 2.6 ) + 2 x 0.6 Γ ( 1.6 ) , g 2 ( x ) = x 4 − 4 x 3 − x − 1 + 24 x 3.4 Γ ( 4.4 ) − x 0.4 Γ ( 1.4 ) ,

that satisfies the exact solution v ( x ) = x − x 2 and u ( x ) = x 4 − x .

Accordingly, rewriting the fractional system in the operator form, and thereafter applying the inverses I 0.4 and I 0.6 on (32) and (33), respectively, via the application of (14) and (7), together with the conditions (34), one gets

v ( x ) = I 0.4 [ g 1 ( x ) ] + I 0.4 [ v ′ ( x ) − v ( x ) ] − I 0.4 [ u ″ ( x ) ] , u ( x ) = I 0.6 [ g 2 ( x ) ] + I 0.6 [ u ′ ( x ) − u ( x ) ] − I 0.6 [ v ″ ( x ) ] .

Next, by applying (14) and the property (10) on I 0.4 [ g 1 ( x ) ] and I 0.6 [ g 2 ( x ) ] , one obtains

v ( x ) = 22 Γ ( 3.4 ) x 2.4 + 3 Γ ( 2.4 ) x 1.4 − x 2 + x − [ 22 Γ ( 3.4 ) + 3 Γ ( 2.4 ) ] x     + I 0.4 [ v ′ ( x ) − v ( x ) ] − I 0.4 [ u ″ ( x ) ] , u ( x ) = Γ ( 5 ) Γ ( 5.6 ) x 4.6 − 4 Γ ( 4.4 ) Γ ( 4.6 ) x 3.6 − 1 Γ ( 2.6 ) x 1.6 + x 4 − x     − [ Γ ( 5 ) Γ ( 5.6 ) − 4 Γ ( 4.4 ) Γ ( 4.6 ) − 1 Γ ( 2.6 ) ] x + I 0.6 [ u ′ ( x ) − u ( x ) ] − I 0.6 [ v ″ ( x ) ] .

Lastly, through the use of the reliable modification of the ADM, the proposed new method reveals the recursive relations for v ( x ) and u ( x ) , respectively, as

{ v 0 ( x ) = x − x 2 , v 1 ( x ) = 22 Γ ( 3.4 ) x 2.4 + 3 Γ ( 2.4 ) x 1.4 − [ 22 Γ ( 3.4 ) + 3 Γ ( 2.4 ) ] x                       + I 0.4 [ v ′ 0 ( x ) − v 0 ( x ) ] − I 0.4 [ u ″ 0 ( x ) ] , v n + 1 ( x ) = I 0.4 [ v ′ n ( x ) − v n ( x ) ] − I 0.4 [ u ″ n ( x ) ] ,     n ≥ 1 ,

and

{ u 0 ( x ) = x 4 − x , u 1 ( x ) = Γ ( 5 ) Γ ( 5.6 ) x 4.6 − 4 Γ ( 4.4 ) Γ ( 4.6 ) x 3.6 − 1 Γ ( 2.6 ) x 1.6                       − [ Γ ( 5 ) Γ ( 5.6 ) − 4 Γ ( 4.4 ) Γ ( 4.6 ) − 1 Γ ( 2.6 ) ] x + I 0.6 [ u ′ 0 ( x ) ]                       − I 0.6 [ u 0 ( x ) ] − I 0.6 [ v ″ 0 ( x ) ] , u n + 1 ( x ) = I 0.6 [ u ′ n ( x ) + u n ( x ) ] − I 0.6 [ v ″ n ( x ) ] ,     n ≥ 1.

Moreover, computing the values of v 1 ( x ) and u 1 ( x ) by using I 0.4 and I 0.6 , respectively using (14) and the property (10) yields v 1 ( x ) = u 1 ( x ) = 0 . Thus, one obtains v n + 1 ( x ) = 0 ∀ n ≥ 1 and u n + 1 ( x ) = 0 ∀ n ≥ 1 , which subsequently leads to the solution of the fractional system as follows

{ v ( x ) = v 0 ( x ) = x − x 2 , u ( x ) = u 0 ( x ) = x 4 − x ,

which exactly matches the referenced analytical solution.

In the same fashion, Table 5 demonstrates the superior accuracy of the proposed iterative technique over the Sinc-collocation method (SCM) deployed in [54]. Notably, the proposed method attains the exact solution at all considered points using only two components, yielding zero absolute error for both fractional solutions v and u . In contrast, the SCM [54], despite using 40 components, produces nonzero errors of order 10−5 and 10−6. These results affirm the high correctness, competence, and swift convergence of the present techniques for treating systems of FBVPs.

Table 5. Comparison of absolute error for Example 4.5.

x SCM [ 54 ] ( N = 40 ) Present method ( N = 2 )
Error in v Error in u Error in v Error in u
0.0 0 0 0 0
0.1 1.27 × 10 − 5 7.02 × 10 − 7 0 0
0.2 2.07 × 10 − 5 2.10 × 10 − 6 0 0
0.3 1.58 × 10 − 5 3.15 × 10 − 6 0 0
0.4 1.28 × 10 − 5 3.71 × 10 − 6 0 0
0.5 1.14 × 10 − 5 3.62 × 10 − 6 0 0
0.6 1.11 × 10 − 5 3.30 × 10 − 6 0 0
0.7 1.14 × 10 − 5 3.06 × 10 − 6 0 0
0.8 1.01 × 10 − 5 2.22 × 10 − 6 0 0
0.9 8.14 × 10 − 5 1.24 × 10 − 6 0 0
1.0 0 0 0 0
5. Conclusion

In this study, the variant of modified ADM is combined with the fractional inverse operator and successfully applied to treat a wide class of FBVPs. The proposed modification allowed the direct incorporation of boundary conditions into the operator, thereby eliminating the need for unknown constants and simplifying the computational process. Through illustrative examples, the efficacy and accuracy of the devised iterative scheme were demonstrated. The numerical results showed that the new scheme not only produced exact solutions for certain problems but also provided highly accurate approximations for others, outperforming several existing numerical techniques in terms of precision and convergence speed. These results confirm the efficiency and reliability of the modified ADM as a powerful semi-analytical tool for tackling FDEs with boundary conditions. Future work may focus on extending the proposed approach to other types of boundary conditions, such as Neumann or mixed cases, to evaluate its flexibility and robustness.

 References Engheta, N. (2002) On Fractional Calculus and Fractional Multipoles in Electromagnetism. IEEETransactionsonAntennasandPropagation, 44, 554-566. https://doi.org/10.1109/8.489308 10.1109/8.489308https://doi.org/10.1109/8.489308 Engheta, N. 2002 On Fractional Calculus and Fractional Multipoles in Electromagnetism IEEE Transactions on Antennas and Propagation 44 10.1109/8.489308 Duarte, F.B.M. and Machado, J.A.T. (2002) Chaotic Phenomena and Fractional-Order Dynamics in the Trajectory Control of Redundant Manipulators. NonlinearDynamics, 29, 315-342. https://doi.org/10.1023/a:1016559314798 10.1023/a:1016559314798https://doi.org/10.1023/a:1016559314798 Duarte, F.B.M. Machado, J.A.T. 2002 Chaotic Phenomena and Fractional-Order Dynamics in the Trajectory Control of Redundant Manipulators Nonlinear Dynamics 29 101655 10.1023/a:1016559314798 Agrawal, O.P. (2004) A General Formulation and Solution Scheme for Fractional Optimal Control Problems. NonlinearDynamics, 38, 323-337. https://doi.org/10.1007/s11071-004-3764-6 10.1007/s11071-004-3764-6https://doi.org/10.1007/s11071-004-3764-6 Agrawal, O.P. 2004 A General Formulation and Solution Scheme for Fractional Optimal Control Problems Nonlinear Dynamics 38 10.1007/s11071-004-3764-6 Magin, R.L. (2010) Fractional Calculus Models of Complex Dynamics in Biological Tissues. Computers&MathematicswithApplications, 59, 1586-1593. https://doi.org/10.1016/j.camwa.2009.08.039 10.1016/j.camwa.2009.08.039https://doi.org/10.1016/j.camwa.2009.08.039 Magin, R.L. 2010 Fractional Calculus Models of Complex Dynamics in Biological Tissues Computers & Mathematics with Applications 59 10.1016/j.camwa.2009.08.039 Kulish, V.V. and Lage, J.L. (2002) Application of Fractional Calculus to Fluid Mechanics. JournalofFluidsEngineering, 124, 803-806. https://doi.org/10.1115/1.1478062 10.1115/1.1478062https://doi.org/10.1115/1.1478062 Kulish, V.V. Lage, J.L. 2002 Application of Fractional Calculus to Fluid Mechanics Journal of Fluids Engineering 124 10.1115/1.1478062 Oldham, K.B. (2010) Fractional Differential Equations in Electrochemistry. AdvancesinEngineeringSoftware, 41, 9-12. https://doi.org/10.1016/j.advengsoft.2008.12.012 10.1016/j.advengsoft.2008.12.012https://doi.org/10.1016/j.advengsoft.2008.12.012 Oldham, K.B. 2010 Fractional Differential Equations in Electrochemistry Advances in Engineering Software 41 10.1016/j.advengsoft.2008.12.012 Gafiychuk, V., Datsko, B. and Meleshko, V. (2008) Mathematical Modeling of Time Fractional Reaction-Diffusion Systems. JournalofComputationalandAppliedMathematics, 220, 215-225. https://doi.org/10.1016/j.cam.2007.08.011 10.1016/j.cam.2007.08.011https://doi.org/10.1016/j.cam.2007.08.011 Gafiychuk, V. Datsko, B. Meleshko, V. 2008 Mathematical Modeling of Time Fractional Reaction-Diffusion Systems Journal of Computational and Applied Mathematics 220 10.1016/j.cam.2007.08.011 Mainardi, F. (1997) Fractional Calculus: Some Basic Problems in Continuum and Statistical Mechanics. In: Carpinteri, A. and Mainardi, F. Eds., FractalsandFractionalCalculusinContinuumMechanics, Springer, 291-348. https://doi.org/10.1007/978-3-7091-2664-6_7 10.1007/978-3-7091-2664-6_7https://doi.org/10.1007/978-3-7091-2664-6_7 Mainardi, F. Carpinteri, A. Mainardi, F. Mechanics, S 1997 Fractional Calculus: Some Basic Problems in Continuum and Statistical Mechanics In: Carpinteri 291 10.1007/978-3-7091-2664-6_7 Meral, F.C., Royston, T.J. and Magin, R. (2010) Fractional Calculus in Viscoelasticity: An Experimental Study. CommunicationsinNonlinearScienceandNumericalSimulation, 15, 939-945. https://doi.org/10.1016/j.cnsns.2009.05.004 10.1016/j.cnsns.2009.05.004https://doi.org/10.1016/j.cnsns.2009.05.004 Meral, F.C. Royston, T.J. Magin, R. 2010 Fractional Calculus in Viscoelasticity: An Experimental Study Communications in Nonlinear Science and Numerical Simulation 15 10.1016/j.cnsns.2009.05.004 Debnath, L. (2003) Recent Applications of Fractional Calculus to Science and Engineering. International Journal of Mathematics and Mathematical Sciences, 2003, 3413-3442. https://doi.org/10.1155/s0161171203301486 10.1155/s0161171203301486https://doi.org/10.1155/s0161171203301486 Debnath, L. 2003 Recent Applications of Fractional Calculus to Science and Engineering International Journal of Mathematics and Mathematical Sciences 2003 10.1155/s0161171203301486 Sardar, T., Rana, S. and Chattopadhyay, J. (2015) A Mathematical Model of Dengue Transmission with Memory. CommunicationsinNonlinearScienceandNumericalSimulation, 22, 511-525. https://doi.org/10.1016/j.cnsns.2014.08.009 10.1016/j.cnsns.2014.08.009https://doi.org/10.1016/j.cnsns.2014.08.009 Sardar, T. Rana, S. Chattopadhyay, J. 2015 A Mathematical Model of Dengue Transmission with Memory Communications in Nonlinear Science and Numerical Simulation 22 10.1016/j.cnsns.2014.08.009 Diethelm, K. and Freed, A.D. (1999) On the Solution of Nonlinear Fractional-Order Differential Equations Used in the Modeling of Viscoplasticity. In: Keil, F., Mackens, W., Voß, H. and Werther, J. Eds., ScientificComputinginChemicalEngineeringII, Springer, 217-224. https://doi.org/10.1007/978-3-642-60185-9_24 10.1007/978-3-642-60185-9_24https://doi.org/10.1007/978-3-642-60185-9_24 Diethelm, K. Freed, A.D. Keil, F. Mackens, W. Werther, J. II, S 1999 On the Solution of Nonlinear Fractional-Order Differential Equations Used in the Modeling of Viscoplasticity In: Keil 217 10.1007/978-3-642-60185-9_24 Dalir, M. and Bashour, M. (2010) Applications of Fractional Calculus. Applied Mathematical Sciences, 4, 1021-1032. https://www.m-hikari.com/ams/ams-2010/ams-21-24-2010/bashourAMS21-24-2010.pdf Dalir, M. Bashour, M. 2010 Applications of Fractional Calculus Applied Mathematical Sciences 4 Wang, G., Yuan, H. and Baleanu, D. (2024) Stability Analysis and Solutions of Fractional Boundary Value Problem on the Cyclopentasilane Graph. Heliyon, 10, e32411. https://doi.org/10.1016/j.heliyon.2024.e32411 10.1016/j.heliyon.2024.e3241138975069https://doi.org/10.1016/j.heliyon.2024.e32411 Wang, G. Yuan, H. Baleanu, D. 2024 Stability Analysis and Solutions of Fractional Boundary Value Problem on the Cyclopentasilane Graph Heliyon 10 10.1016/j.heliyon.2024.e32411 38975069 Etemad, S. and Rezapour, S. (2020) On the Existence of Solutions for Fractional Boundary Value Problems on the Ethane Graph. AdvancesinDifferenceEquations, 2020, Article No. 276. https://doi.org/10.1186/s13662-020-02736-4 10.1186/s13662-020-02736-4https://doi.org/10.1186/s13662-020-02736-4 Etemad, S. Rezapour, S. 2020 On the Existence of Solutions for Fractional Boundary Value Problems on the Ethane Graph Advances in Difference Equations 2020 No 10.1186/s13662-020-02736-4 Rezapour, S., Deressa, C.T., Hussain, A., Etemad, S., George, R. and Ahmad, B. (2022) A Theoretical Analysis of a Fractional Multi-Dimensional System of Boundary Value Problems on the Methylpropane Graph via Fixed Point Technique. Mathematics, 10, Article 568. https://doi.org/10.3390/math10040568 10.3390/math10040568https://doi.org/10.3390/math10040568 Rezapour, S. Deressa, C.T. Hussain, A. Etemad, S. George, R. Ahmad, B. 2022 A Theoretical Analysis of a Fractional Multi-Dimensional System of Boundary Value Problems on the Methylpropane Graph via Fixed Point Technique Mathematics 10 568 10.3390/math10040568 El-Ajou, A., Abu Arqub, O. and Momani, S. (2013) Solving Fractional Two-Point Boundary Value Problems Using Continuous Analytic Method. AinShamsEngineeringJournal, 4, 539-547. https://doi.org/10.1016/j.asej.2012.11.010 10.1016/j.asej.2012.11.010https://doi.org/10.1016/j.asej.2012.11.010 El-Ajou, A. Arqub, O. Momani, S. 2013 Solving Fractional Two-Point Boundary Value Problems Using Continuous Analytic Method Ain Shams Engineering Journal 4 10.1016/j.asej.2012.11.010 Jafari, H., Bagherian, K. and Moshokoa, S.P. (2014) Homotopy Perturbation Method to Obtain Positive Solutions of Nonlinear Boundary Value Problems of Fractional Order. AbstractandAppliedAnalysis, 2014, Article ID: 919052. https://doi.org/10.1155/2014/919052 10.1155/2014/919052https://doi.org/10.1155/2014/919052 Jafari, H. Bagherian, K. Moshokoa, S.P. 2014 Homotopy Perturbation Method to Obtain Positive Solutions of Nonlinear Boundary Value Problems of Fractional Order Abstract and Applied Analysis 2014 919052 ID 10.1155/2014/919052 Wu, G.C. and Lee, E.W.M. (2010) Fractional Variational Iteration Method and Its Application. PhysicsLettersA, 374, 2506-2509. https://doi.org/10.1016/j.physleta.2010.04.034 10.1016/j.physleta.2010.04.034https://doi.org/10.1016/j.physleta.2010.04.034 Wu, G.C. Lee, E.W.M. 2010 Fractional Variational Iteration Method and Its Application Physics Letters A 374 10.1016/j.physleta.2010.04.034 Di Matteo, A. and Pirrotta, A. (2015) Generalized Differential Transform Method for Nonlinear Boundary Value Problem of Fractional Order. CommunicationsinNonlinearScienceandNumericalSimulation, 29, 88-101. https://doi.org/10.1016/j.cnsns.2015.04.017 10.1016/j.cnsns.2015.04.017https://doi.org/10.1016/j.cnsns.2015.04.017 Matteo, A. Pirrotta, A. 2015 Generalized Differential Transform Method for Nonlinear Boundary Value Problem of Fractional Order Communications in Nonlinear Science and Numerical Simulation 29 10.1016/j.cnsns.2015.04.017 Kawala, A.M. (2018) Numerical Solution for Initial and Boundary Value Problems of Fractional Order. AdvancesinPureMathematics, 8, 831-844. https://doi.org/10.4236/apm.2018.812051 10.4236/apm.2018.812051https://doi.org/10.4236/apm.2018.812051 Kawala, A.M. 2018 Numerical Solution for Initial and Boundary Value Problems of Fractional Order Advances in Pure Mathematics 8 10.4236/apm.2018.812051 Lobato, F.S. and de Oliveira Cardoso, A. (2023) Numerical Solution of Fractional Boundary Value Problem in Chemical Engineering Using the Shooting Method. RevistaCereus, 15, 149-165. https://doi.org/10.18605/2175-7275/cereus.v15n2p149-165 10.18605/2175-7275/cereus.v15n2p149-165https://doi.org/10.18605/2175-7275/cereus.v15n2p149-165 Lobato, F.S. Cardoso, A. 2023 Numerical Solution of Fractional Boundary Value Problem in Chemical Engineering Using the Shooting Method Revista Cereus 15 10.18605/2175-7275/cereus.v15n2p149-165 Doha, E.H., Bhrawy, A.H. and Ezz-Eldien, S.S. (2011) A Chebyshev Spectral Method Based on Operational Matrix for Initial and Boundary Value Problems of Fractional Order. Computers&MathematicswithApplications, 62, 2364-2373. https://doi.org/10.1016/j.camwa.2011.07.024 10.1016/j.camwa.2011.07.024https://doi.org/10.1016/j.camwa.2011.07.024 Doha, E.H. Bhrawy, A.H. Ezz-Eldien, S.S. 2011 A Chebyshev Spectral Method Based on Operational Matrix for Initial and Boundary Value Problems of Fractional Order Computers & Mathematics with Applications 62 10.1016/j.camwa.2011.07.024 Al-Mdallal, Q.M., Syam, M.I. and Anwar, M.N. (2010) A Collocation-Shooting Method for Solving Fractional Boundary Value Problems. CommunicationsinNonlinearScienceandNumericalSimulation, 15, 3814-3822. https://doi.org/10.1016/j.cnsns.2010.01.020 10.1016/j.cnsns.2010.01.020https://doi.org/10.1016/j.cnsns.2010.01.020 Al-Mdallal, Q.M. Syam, M.I. Anwar, M.N. 2010 A Collocation-Shooting Method for Solving Fractional Boundary Value Problems Communications in Nonlinear Science and Numerical Simulation 15 10.1016/j.cnsns.2010.01.020 Al-Refai, M. and Ali Hajji, M. (2011) Monotone Iterative Sequences for Nonlinear Boundary Value Problems of Fractional Order. NonlinearAnalysis: Theory, Methods&Applications, 74, 3531-3539. https://doi.org/10.1016/j.na.2011.03.006 10.1016/j.na.2011.03.006https://doi.org/10.1016/j.na.2011.03.006 Al-Refai, M. Hajji, M. Theory, M 2011 Monotone Iterative Sequences for Nonlinear Boundary Value Problems of Fractional Order Nonlinear Analysis: Theory 74 10.1016/j.na.2011.03.006 Pedas, A. and Tamme, E. (2012) Piecewise Polynomial Collocation for Linear Boundary Value Problems of Fractional Differential Equations. JournalofComputationalandAppliedMathematics, 236, 3349-3359. https://doi.org/10.1016/j.cam.2012.03.002 10.1016/j.cam.2012.03.002https://doi.org/10.1016/j.cam.2012.03.002 Pedas, A. Tamme, E. 2012 Piecewise Polynomial Collocation for Linear Boundary Value Problems of Fractional Differential Equations Journal of Computational and Applied Mathematics 236 10.1016/j.cam.2012.03.002 Adomian, G. (2013) Solving Frontier Problems of Physics: The Decomposition Method. Springer Science and Business Media. Adomian, G. 2013 Solving Frontier Problems of Physics: The Decomposition Method Wazwaz, A.-M. (2010) Partial Differential Equations and Solitary Waves Theory. Springer Science and Business Media. Wazwaz, A. 2010 Partial Differential Equations and Solitary Waves Theory Wazwaz, A.-M. (1999) A Reliable Modification of Adomian Decomposition Method. AppliedMathematicsandComputation, 102, 77-86. https://doi.org/10.1016/s0096-3003(98)10024-3 10.1016/s0096-3003(98)10024-3https://doi.org/10.1016/s0096-3003(98)10024-3 Wazwaz, A. 1999 A Reliable Modification of Adomian Decomposition Method Applied Mathematics and Computation 3003 98 10.1016/s0096-3003(98)10024-3 Wazwaz, A.-M. and El-Sayed, S.M. (2001) A New Modification of the Adomian Decomposition Method for Linear and Nonlinear Operators. AppliedMathematicsandComputation, 122, 393-405. https://doi.org/10.1016/s0096-3003(00)00060-6 10.1016/s0096-3003(00)00060-6https://doi.org/10.1016/s0096-3003(00)00060-6 Wazwaz, A. El-Sayed, S.M. 2001 A New Modification of the Adomian Decomposition Method for Linear and Nonlinear Operators Applied Mathematics and Computation 3003 00 10.1016/s0096-3003(00)00060-6 Al-Mazmumy, M. and Alsulami, M. (2024) Utilization of the Modified Adomian Decomposition Method on the Bagley-Torvik Equation Amidst Dirichlet Boundary Conditions. EuropeanJournalofPureandAppliedMathematics, 17, 546-568. https://doi.org/10.29020/nybg.ejpam.v17i1.5050 10.29020/nybg.ejpam.v17i1.5050https://doi.org/10.29020/nybg.ejpam.v17i1.5050 Al-Mazmumy, M. Alsulami, M. 2024 Utilization of the Modified Adomian Decomposition Method on the Bagley-Torvik Equation Amidst Dirichlet Boundary Conditions European Journal of Pure and Applied Mathematics 17 10.29020/nybg.ejpam.v17i1.5050 Mulenga, J. and Phiri, P.A. (2023) Solving Different Types of Differential Equations Using Modified and New Modified Adomian Decomposition Methods. Journal ofAppliedMathematicsandPhysics, 11, 1656-1676. https://doi.org/10.4236/jamp.2023.116108 10.4236/jamp.2023.116108https://doi.org/10.4236/jamp.2023.116108 Mulenga, J. Phiri, P.A. 2023 Solving Different Types of Differential Equations Using Modified and New Modified Adomian Decomposition Methods Journal of Applied Mathematics and Physics 11 10.4236/jamp.2023.116108 Alsulami, A., Al-Mazmumy, M., Bakodah, H. and Alzaid, N. (2023) On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences. AmericanJournalofComputationalMathematics, 13, 387-397. https://doi.org/10.4236/ajcm.2023.133022 10.4236/ajcm.2023.133022https://doi.org/10.4236/ajcm.2023.133022 Alsulami, A. Al-Mazmumy, M. Bakodah, H. Alzaid, N. 2023 On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences American Journal of Computational Mathematics 13 10.4236/ajcm.2023.133022 AL-Zaid, N., AL-Refaidi, A., Bakodah, H. and AL-Mazmumy, M. (2022) Solution of Second-and Higher-Order Nonlinear Two-Point Boundary-Value Problems Using Double Decomposition Method. Mathematics, 10, Article 3519. https://doi.org/10.3390/math10193519 10.3390/math10193519https://doi.org/10.3390/math10193519 AL-Zaid, N. AL-Refaidi, A. Bakodah, H. AL-Mazmumy, M. 2022 Solution of Second-and Higher-Order Nonlinear Two-Point Boundary-Value Problems Using Double Decomposition Method Mathematics 10 3519 10.3390/math10193519 Al-Hayani, W. (2015) An Efficient Advanced Adomian Decomposition Method to Solve Second-Order Boundary Value Problems with Neumann Conditions. InternationalMathematicalForum, 10, 13-23. https://doi.org/10.12988/imf.2015.410182 10.12988/imf.2015.410182https://doi.org/10.12988/imf.2015.410182 Al-Hayani, W. 2015 An Efficient Advanced Adomian Decomposition Method to Solve Second-Order Boundary Value Problems with Neumann Conditions International Mathematical Forum 10 10.12988/imf.2015.410182 Duan, J.-S. and Rach, R. (2011) A New Modification of the Adomian Decomposition Method for Solving Boundary Value Problems for Higher Order Nonlinear Differential Equations. AppliedMathematicsandComputation, 218, 4090-4118. https://doi.org/10.1016/j.amc.2011.09.037 10.1016/j.amc.2011.09.037https://doi.org/10.1016/j.amc.2011.09.037 Duan, J. Rach, R. 2011 A New Modification of the Adomian Decomposition Method for Solving Boundary Value Problems for Higher Order Nonlinear Differential Equations Applied Mathematics and Computation 218 10.1016/j.amc.2011.09.037 El-Sayed, A.M.A. and Gaber, M. (2006) The Adomian Decomposition Method for Solving Partial Differential Equations of Fractal Order in Finite Domains. PhysicsLettersA, 359, 175-182. https://doi.org/10.1016/j.physleta.2006.06.024 10.1016/j.physleta.2006.06.024https://doi.org/10.1016/j.physleta.2006.06.024 El-Sayed, A.M.A. Gaber, M. 2006 The Adomian Decomposition Method for Solving Partial Differential Equations of Fractal Order in Finite Domains Physics Letters A 359 10.1016/j.physleta.2006.06.024 Ziada, E.A.A. (2021) Solution of Nonlinear Fractional Differential Equations Using Adomain Decomposition Method. InternationalJournalofSystemsScienceandAppliedMathematics, 6, Article 111. https://doi.org/10.11648/j.ijssam.20210604.11 10.11648/j.ijssam.20210604.11https://doi.org/10.11648/j.ijssam.20210604.11 Ziada, E.A.A. 2021 Solution of Nonlinear Fractional Differential Equations Using Adomain Decomposition Method International Journal of Systems Science and Applied Mathematics 6 111 10.11648/j.ijssam.20210604.11 Wanassi, O.K., Bourguiba, R. and Torres, D.F.M. (2024) Existence and Uniqueness of Solution for Fractional Differential Equations with Integral Boundary Conditions and the Adomian Decomposition Method. MathematicalMethodsintheAppliedSciences, 47, 3582-3595. https://doi.org/10.1002/mma.8880 10.1002/mma.8880https://doi.org/10.1002/mma.8880 Wanassi, O.K. Bourguiba, R. Torres, D.F.M. 2024 Existence and Uniqueness of Solution for Fractional Differential Equations with Integral Boundary Conditions and the Adomian Decomposition Method Mathematical Methods in the Applied Sciences 47 10.1002/mma.8880 Ziada, E.A.A. (2021) Solution of Nonlinear Fractional Differential Equations Using Adomian Decomposition Method. Journal of Applied Mathematics and Physics, 9, 2544-2557. Ziada, E.A.A. 2021 Solution of Nonlinear Fractional Differential Equations Using Adomian Decomposition Method Journal of Applied Mathematics and Physics 9 Daftardar-Gejji, V. and Jafari, H. (2005) Adomian Decomposition: A Tool for Solving a System of Fractional Differential Equations. JournalofMathematicalAnalysisandApplications, 301, 508-518. https://doi.org/10.1016/j.jmaa.2004.07.039 10.1016/j.jmaa.2004.07.039https://doi.org/10.1016/j.jmaa.2004.07.039 Daftardar-Gejji, V. Jafari, H. 2005 Adomian Decomposition: A Tool for Solving a System of Fractional Differential Equations Journal of Mathematical Analysis and Applications 301 10.1016/j.jmaa.2004.07.039 Kilbas, A.A., Srivastava, H.M. and Trujillo, J.J. (2006) Theory and Applications of Fractional Differential Equations. Elsevier. Kilbas, A.A. Srivastava, H.M. Trujillo, J.J. 2006 Theory and Applications of Fractional Differential Equations Luchko, Y. and Gorenflo, R. (1999) An Operational Method for Solving Fractional Differential Equations with the Caputo Derivatives. Acta MathematicaVietnamica, 24, 207-233. https://math.ac.vn/public/uploads/files/9902207.pdf Luchko, Y. Gorenflo, R. 1999 An Operational Method for Solving Fractional Differential Equations with the Caputo Derivatives Acta Mathematica Vietnamica 24 Al-Mazmumy, M. (2023) An Efficient Maple Program for Calculating Adomian Polynomials. AmericanJournalofComputationalMathematics, 13, 512-523. https://doi.org/10.4236/ajcm.2023.134028 10.4236/ajcm.2023.134028https://doi.org/10.4236/ajcm.2023.134028 Al-Mazmumy, M. 2023 An Efficient Maple Program for Calculating Adomian Polynomials American Journal of Computational Mathematics 13 10.4236/ajcm.2023.134028 Yüzbaşı, Ş. (2013) Numerical Solution of the Bagley-Torvik Equation by the Bessel Collocation Method. MathematicalMethodsintheAppliedSciences, 36, 300-312. https://doi.org/10.1002/mma.2588 10.1002/mma.2588https://doi.org/10.1002/mma.2588 2013 Numerical Solution of the Bagley-Torvik Equation by the Bessel Collocation Method Mathematical Methods in the Applied Sciences 36 10.1002/mma.2588 Batool, A., Talib, I. and Riaz, M.B. (2025) Fractional-Order Boundary Value Problems Solutions Using Advanced Numerical Technique. PartialDifferentialEquationsinAppliedMathematics, 13, Article 101059. https://doi.org/10.1016/j.padiff.2024.101059 10.1016/j.padiff.2024.101059https://doi.org/10.1016/j.padiff.2024.101059 Batool, A. Talib, I. Riaz, M.B. 2025 Fractional-Order Boundary Value Problems Solutions Using Advanced Numerical Technique Partial Differential Equations in Applied Mathematics 13 101059 10.1016/j.padiff.2024.101059 Esmaeili, S. and Shamsi, M. (2010) A Pseudo-Spectral Scheme for the Approximate Solution of a Family of Fractional Differential Equations. CommunicationsinNonlinearScienceandNumericalSimulation, 16, 3646-3654. https://doi.org/10.1016/j.cnsns.2010.12.008 10.1016/j.cnsns.2010.12.008https://doi.org/10.1016/j.cnsns.2010.12.008 Esmaeili, S. Shamsi, M. 2010 A Pseudo-Spectral Scheme for the Approximate Solution of a Family of Fractional Differential Equations Communications in Nonlinear Science and Numerical Simulation 16 10.1016/j.cnsns.2010.12.008 Diethelm, K., Ford, N.J. and Freed, A.D. (2004) Detailed Error Analysis for a Fractional Adams Method. NumericalAlgorithms, 36, 31-52. https://doi.org/10.1023/b:numa.0000027736.85078.be 10.1023/b:numa.0000027736.85078.behttps://doi.org/10.1023/b:numa.0000027736.85078.be Diethelm, K. Ford, N.J. Freed, A.D. 2004 Detailed Error Analysis for a Fractional Adams Method Numerical Algorithms 36 10.1023/b:numa.0000027736.85078.be Al-Nana, A.A., Batiha, I.M. and Momani, S. (2023) A Numerical Approach for Dealing with Fractional Boundary Value Problems. Mathematics, 11, Article 4082. https://doi.org/10.3390/math11194082 10.3390/math11194082https://doi.org/10.3390/math11194082 Al-Nana, A.A. Batiha, I.M. Momani, S. 2023 A Numerical Approach for Dealing with Fractional Boundary Value Problems Mathematics 11 4082 10.3390/math11194082 Bratu, G. (1914) Sur les équations intégrales non linéaires. Bulletin de la SociétémathématiquedeFrance, 2, 113-142. https://doi.org/10.24033/bsmf.943 10.24033/bsmf.943https://doi.org/10.24033/bsmf.943 Bratu, G. 1914 Sur les équations intégrales non linéaires Bulletin de la Société mathématique de France 2 10.24033/bsmf.943 Jacobsen, J. and Schmitt, K. (2002) The Liouville-Bratu-Gelfand Problem for Radial Operators. JournalofDifferentialEquations, 184, 283-298. https://doi.org/10.1006/jdeq.2001.4151 10.1006/jdeq.2001.4151https://doi.org/10.1006/jdeq.2001.4151 Jacobsen, J. Schmitt, K. 2002 The Liouville-Bratu-Gelfand Problem for Radial Operators Journal of Differential Equations 184 10.1006/jdeq.2001.4151 Jafari, H. and Tajadodi, H. (2016) Electro—Spunorganic Nanofibers Elaboration Process Investigations Using BPs Operational Matrices. Iranian Journal of Mathematical Chemistry, 7, 19-27. https://ijmc.kashanu.ac.ir/article_11866_914b2e89361dc43a0182ac89574aa076.pdf Jafari, H. Tajadodi, H. 2016 Electro—Spunorganic Nanofibers Elaboration Process Investigations Using BPs Operational Matrices Iranian Journal of Mathematical Chemistry 7 Al-Smadi, M., Momani, S., Djeddi, N., El-Ajou, A. and Al-Zhour, Z. (2023) Adaptation of Reproducing Kernel Method in Solving Atangana-Baleanu Fractional Bratu Model. InternationalJournalofDynamicsandControl, 11, 136-148. https://doi.org/10.1007/s40435-022-00961-1 10.1007/s40435-022-00961-1https://doi.org/10.1007/s40435-022-00961-1 Al-Smadi, M. Momani, S. Djeddi, N. El-Ajou, A. Al-Zhour, Z. 2023 Adaptation of Reproducing Kernel Method in Solving Atangana-Baleanu Fractional Bratu Model International Journal of Dynamics and Control 11 10.1007/s40435-022-00961-1 Hatipoglu, V.F., Alkan, S. and Secer, A. (2017) An Efficient Scheme for Solving a System of Fractional Differential Equations with Boundary Conditions. AdvancesinDifferenceEquations, 2017, Article No. 204. https://doi.org/10.1186/s13662-017-1260-9 10.1186/s13662-017-1260-9https://doi.org/10.1186/s13662-017-1260-9 Hatipoglu, V.F. Alkan, S. Secer, A. 2017 An Efficient Scheme for Solving a System of Fractional Differential Equations with Boundary Conditions Advances in Difference Equations 2017 No 10.1186/s13662-017-1260-9