ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article

An Iterative Approach for Solving Fractional Differential Equations Using the α-Generalized Daftardar–Jafari Method

[version 1; peer review: awaiting peer review]
PUBLISHED 23 Dec 2025
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS AWAITING PEER REVIEW

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

Abstract

Background

Fractional calculus has gained significant attention due to its applications in modeling complex systems. This paper introduces several fundamental fractional operators, including the Riemann–Liouville, Caputo, and Hadamard operators. From these, a novel α-generalized fractional operator is derived, which plays a crucial role in solving fractional partial differential equations (FPDEs).

Methods

The α-generalized operator is analytically constructed and its properties are rigorously proven. It is then integrated with the Daftardar–Jafari iterative method (DJM) to solve both linear and nonlinear FPDEs. The convergence of the proposed method is established using the Lipschitz condition. Two benchmark equations—the linear fractional Burger’s equation and the nonlinear fractional Heat-type equation—are used to demonstrate the method’s applicability. MATLAB is employed to implement the numerical schemes and visualize the results.

Results

The DJM combined with the α-generalized operator yields superior numerical results compared to the Adomian decomposition method. This superiority is evident in both accuracy and convergence, especially for the nonlinear fractional Heat-type equation. Tabulated results and graphical comparisons confirm the effectiveness of the proposed approach across different values of α.

Conclusions

The integration of the α-generalized operator with the Daftardar–Jafari method provides a powerful tool for solving FPDEs. Its numerical advantages make it a promising alternative to existing methods. The study recommends extending this framework to other fractional operators, iterative schemes, and potential transformation techniques, with applications in quantum physics, fluid dynamics, and beyond.

Keywords

Adomian Decomposition Method, Daftardar–Jafari iterative method, Generalized fractional derivative (α-operator), Fractional partial differential equations, Existence and uniqueness, Heat-type equation, Burger’s equation.

1. Introduction

Over the past decades, fractional calculus has emerged as a powerful generalization of classical differentiation and integration. It has found applications across a wide spectrum of disciplines. Its origins trace back to a 1695 correspondence, in which Leibniz proposed the concept of a half-order derivative in his letter to L’Hôpital. Since that seminal insight, fractional operators have attracted the interest of mathematicians, physicists, biologists, engineers, and economists.

One of the earliest celebrated applications appears in Abel’s resolution of the tautochrone problem.1 In modern research, fractional models play a central role in various fields. These include biophysics, quantum mechanics, wave propagation, polymer dynamics, continuum mechanics, Lie group analysis, field theory, and spectroscopy.28 For a comprehensive treatment of these developments, see Samko et al.9

Over time, many formulations of the fractional derivative have been introduced. These include the Riemann–Liouville, Caputo, Hadamard, Erdélyi–Kober, Grünwald–Letnikov, Machado, and Riesz definitions. Each formulation addresses distinct theoretical and practical challenges.812

In fractional calculus, derivatives of non-integer order are most often introduced via fractional integrals.10,11 The Riemann–Liouville operator and its associated integral occupy a foundational role in the theory.13 The widely used Caputo derivative is itself defined through the Riemann–Liouville integral.9

In parallel, Butzer and colleagues have conducted extensive studies on the Hadamard fractional integral and derivative.10,11,1418 They derived, in particular, the Mellin transform representation of these operators.15

More recently, Pooseh et al. established expansion formulas that express Hadamard operators in terms of standard integer-order derivatives.19 For further developments and a comprehensive survey of related results, see Samko et al.,10 Srivastava et al.,11 and the references therein.

In Ref. 12, a new fractional integral was introduced. It unifies the Riemann–Liouville and Hadamard integrals within a single formulation. This α-generalized operator has been the subject of numerous recent studies and monographs.18,2024

The first part of this paper reviews the essential concepts and definitions of fractional calculus that underpin our framework. The second part presents rigorous proofs of existence and uniqueness for the generalized α-operator. In the third part, we develop a detailed analysis of the Daftardar-Jafari iterative method when combined with this operator. The fourth part illustrates the applicability of both the operator and the method through several example problems in fractional differential equations.

Finally, we draw conclusions and compile a complete list of references.

2. Preliminaries

Definition 1.

9 The Riemann Liouville fractional integral operator of order α0, of function φ(τ) is defined as

(1)
Iταφ(τ)={1Γ(α)0τ(τs)α1φ(s)ds,α>0.τ>0φ(τ),α=0,
where Γ() is the well-known Gamma function

Definition 2.

9 The Liouville-Caputo operator (c) with order (α>0) of φ(τ) is defined as follows:

(2)
Dταcφ(τ)={1Γ(mα)0τ(τs)mα1φ(m)(s)ds,m1<αmnτnφ(τ),α=n,

Definition 3.

11 The Hadamard fractional integral introduced by J. Hadamard,11 and is given by,

(3)
Ia+α,ρφ(τ)=1Γ(α)aτ(logτt)α1φ(t)dtt.

For Re(α)>0,τ>a0 while the Hadamard fractional derivative of order α,Re(α)>0 is given by,

(4)
Da+αφ(τ)=(τd)n1Γ(nα)aτ(logτt)nα1φ(t)dtt.

Definition 4.

13 The α‐generalized fractional derivative operator of order 0<α<1,τ>0,ρ is constant of function φ(τ)ε,ε1 is defined as

(5)
Dτα,ρ0φ(τ)=ραΓ(1α)τ1ρd0τ(τρtρ)αtρ1φ(t)dt,
while the α‐generalized fractional integral operator25 of order 0<α<1,τ>0,ρ is constant of function φ(τ) is defined as
(6)
(Ia+αρφ)(τ)=ρ1αΓ(α)aτ(τρtρ)α1tρ1φ(t)dt.

Theorem 1.

26 Let αϵ,Re(α)0,n=Re(α)andρ>0.Then,forx>a,

(7)
limρ1Dτα,ρφ(τ)=(d)n1Γ(nα)0τ(τt)nα1φ(t)dt,
(8)
limρ0+Dτα,ρφ(τ)=(τd)n1Γ(nα)aτ(logτt)nα1φ(t)dtt,
(9)
limρ1Ia+αρφ(τ)=1Γ(α)aτ(τt)α1φ(t)dt,
(10)
limρ0+Ia+αρφ(τ)=1Γ(α)aτ(logτt)α1φ(t)dtt.

Proof.

This result was originally established in Ref. 26, to which we refer the reader for the full proof. This result was originally established in Ref. 26, to which we refer the reader for the full proof.

Theorem 2.

Let 0<α<1,φ(τ)Xcp(a,b),a>0,andρ>0 . Then

(11)
(Da+α,ρIa+αρ)φ(τ)=φ(τ).

Proof.

(12)
(Da+α,ρIa+αρ)φ(τ)=ραΓ(1α)τ1ρdaτ(τρtρ)αtρ1ρ1αΓ(α)at(tρxρ)α1xρ1φ(x)dxdt.

By using Fubini’s theorem, the equation will become as follows,

(13)
=ρΓ(1α)Γ(α)τ1ρdaτxρ1φ(x)dxxτ(τρtρ)α(tρxρ)α1tρ1dt.

By employing an appropriate substitution of variables,

(14)
y=(tρxρ)(τρxρ),
(15)
y1,whentτ,andy0,whentx(tρxρ)=y(τρxρ),
(16)
tρ=xρ+y(τρxρ),
(17)
ρtρ1dt=(τρxρ)dy,
(18)
tρ1dt=(τρxρ)ρdy.

Substitution Eq. (16) and Eq. (18) in Eq. (13):

(19)
=ρΓ(1α)Γ(α)τ1ρdaτxρ1φ(x)dx01(τρxρy(τρxρ))α(xρ+y(τρxρ)xρ)α1(τρxρ)ρdy.
(20)
=ρΓ(1α)Γ(α)τ1ρdaτxρ1φ(x)dx(τρtρ)α+α1+1ρ01(1y)αyα1dy.

Using the Beta function, we have

(21)
=ρΓ(1α)Γ(α)τ1ρdaτxρ1φ(x)dx.Γ(1α)Γ(α)ρ,
(22)
=τ1ρdaτxρ1φ(x)dx.

Applying Leibniz’s rule

(23)
(Da+α,ρIa+αρ)φ(τ)=φ(τ).

Theorem 3.

Let α,σ such that 0<Re(α) and Re(σ)<1. If 0<a<b< and 1ρ , then, for φ(τ)Xcp(a,b),ρ>0

(24)
[Ia+αρIa+σρ]φ(τ)=Ia+α+σρφ(τ),and(Da+α,ρDa+σ,ρ)φ(τ)=Da+α+σ,ρφ(τ).

Proof.

This result was originally established in Ref. 27, to which we refer the reader for the full proof.

Theorem 4.

Let αsuch that0<Re(α)<1.If0<a<b<and1ρ,then,forφ,ψXcp(a,b),ρ>0

(25)
Ia+αρ(φ+ψ)=Ia+αρφ+Ia+αρψ,andDa+α,ρ(φ+ψ)=Da+α,ρφ+Da+α,ρψ.

It is easy to prove this theorem, which is called the linear property theorem.

Proposition 1.

For 0<α<1 and ρ, then

(26)
Da+α,ρτn=ραΓ(nρ+1)Γ(nρ+1α)τnρα,
in which n denotes a constant of arbitrary value.

Proof.

By using definition of α-generalized fractional derivative operator, we get

(27)
Dτα,ρ0τn=ραΓ(1α)τ1ρd0τ(τρtρ)αtρ1tndt,=ραΓ(1α)τ1ρddt0τ(τρtρ)αtn+ρ1dt.

Let x=tρτρt=x1pτdt=τ1ρx1p1dx,sowe get,

(28)
Dτα,ρ0τn=ραΓ(1α)τ1ρd01(1x)αtρα(xρτ)n+ρ1.τ.1ρx1p1dx,=ραΓ(1α)τ1ρd.τρα+n+ρ1+1.1ρ01(1x)αxnρ+11ρ+1ρ1dx,=ραΓ(1α)τ1ρd.τn+ρρα.1ρ01(1x)αxnρdx,=ραΓ(1α)τ1ρ(n+ραρ).τn+ρρα1.1ρ.Γ(1α)Γ(nρ+1)Γ(nρα+2),=ραΓ(nρ+1)(nρα+1)Γ(nρα+1)(nρ+1α).τn+ρρα1+1ρ,=ραΓ(nρ+1)Γ(nρ+1α)τnρα.

Proposition 2.

For 0<α<1 and ρ, then:

(29)
Ia+αρτn=ραΓ(nρ+1)Γ(nρ+1+α)τn+ρα,
in which n denotes a constant of arbitrary value

Proof.

By using definition of α‐generalized fractional integral operator, we get

(30)
Ia+αρτn=ρ1αΓ(α)0τ(τρtρ)α1tρ1tndt,=ρ1αΓ(α)0τ(1tρτρ)α1ταρρtρ1+ndt.

Let y=tρτρt=y1ρτdt=τ1ρy1ρ1dy,sowe get,

(31)
Ia+αρτn=ρ1αΓ(α)ταρρ01(1y)α1(y1ρτ)ρ1+n1ρy1ρ1τdy,=ρ1αΓ(α)ταρρ+ρ+n1+1.1ρ01(1y)α1.y1+nρ1ρ+1ρ1dy,=ραΓ(α)ταρ+n.01(1y)α1.ynρdy,=ραΓ(α)ταρ+nΓ(α)Γ(nρ+1)Γ(nρ+α+1),=ραΓ(nρ+1)Γ(nρ+1+α)τn+ρα.

3. The α‐Generalized Fractional Daftardar-Jafari Method (α-DJM)

Consider the following form of a fractional differential equation,

(32)
Da+α,ρφ(ξ,τ)+R(φ(ξ,τ))+N(φ(ξ,τ))=g(ξ,τ),m1<α<m,τ>0,
with respect to the initial conditions
(33)
ktkφ(ξ,0)=fk(ξ),k=0,1,2,,m1,
where Da+α,ρφ(ξ,τ) α‐generalized fractional operator of order α, R linear operator, N non linear operator, and g(ξ,τ) known source function.

The method is based on applying Ia+αρ on both sides of (43), to have

(34)
φ(ξ,τ)k=0m1τkk!φ(k)(0)+Ia+αρ(R(φ))+Ia+αρ(N(φ))=Ia+αρ(g),
(35)
φ(ξ,τ)=Ia+αρ(g)+k=0m1τkk!φ(k)(0)Ia+αρ(R(φ))Ia+αρ(N(φ)).

Now, represent solution as an infinite series given below:

(36)
φ(ξ,τ)=n=0φn.

Substituting (36) into both sides of (35) gives.

(37)
m=0φm=Ia+αρ(g)+k=0m1τkk!φ(k)(0)Ia+αρ(R(n=0φn))Ia+αρ(N(n=0φn)).

The nonlinear term is decomposed as in

(38)
N(m=0φm)=N(φ0)+m=0[N(r=0mφr)N(r=0m1φr)].

Substituting (38) into (37) gives

(39)
m=0φm=Ia+αρ(g)+k=0m1τkk!φ(k)(0)Ia+αρ(R(n=0φn))Ia+αρ(N(φ0)+m=0[N(r=0mφr)N(r=0m1φr)]).

Define the recurrence relation:

(40)
F=Ia+αρ(g)+k=0m1τkk!φ(k)(0)=φ0,
(41)
L(φi)=Ia+αρ(R(φi)),i1,
(42)
Gi=Ia+αρ[N(r=0iφr)N(r=0i1φr)],i1.

Substituting Eq. (40), Eq. (41), and Eq. (42) into Eq. (39) gives

(43)
φ(x,t)=F+L(φ)+G(φ).

Where

(44)
G(φ)=i=0Gi.

And

(45)
L(φ)=L(i=0φi).

Moreover, the relation is defined with recuence, so that

(46)
φ0=F,     

And

(47)
φi+1=L(φi)+Gi.

Thus, the approximate solution of (32) is

(48)
φ(ξ,τ)=i=1φi(ξ,τ)=φ0+φ1+φ2+,
Theorem 3.1

Let us consider the fractional differential equation of the form

(49)
Da+α,ρφ(τ)=R(τ,φ(τ))+N(τ,φ(τ))0<α<1,
where Da+α,ρφ(τ) α‐generalized fractional operator of order α, R linear operator, N non linear operator, satisfying the Lipschitz condition
|N(τ,φ1)N(τ,φ2)|L|φ1φ2|,φ1,φ2,τ[0,T]
for some constant L>0,andφ(0)=φ0. Then the approximate solution obtained by the α‐Generalized Fractional Daftardar-Jafari Method (α-DJM), defined by the series
(50)
φ(τ)=n=0φn(τ),
with the recurrence relations
(51)
φ0(τ)=φ0,
(52)
φ1(τ)=Ia+αρ[R(τ,φ0)+N(τ,φ0)],
(53)
φn+1(τ)=Ia+αρ[R(τ,φ0)+(N(τ,k=0nφk(τ))N(τ,k=0n1φk(τ)))],
is convergent for all τ[0,T], provided that LTραΓ(α+1)<1.

Proof.

Assume |R(τ,φn)|Mfor someM>0. Using the boundedness of Ia+αρ we estimate:

(54)
φ1(τ)TραΓ(α+1)(M+Lφ0),
(55)
φ2(τ)TραΓ(α+1)(M+Lφ1),
(56)
φn+1(τ)TραΓ(α+1)(M+Lφn),
(57)
φn+1(τ)CM+CLφn,
where C=TραΓ(α+1).

By recursion:

(58)
φn+1(τ)CMk=0n(AL)k,

Therefore, the total norm of the series is bounded by

(59)
n=0φn(τ)n=0φn(τ)<,ifAL<1,

Hence, the method converges absolutely and uniformly for τ[0,T] under the given condition.

4. Applications

Example 1.

Assume heat-like equation with the NFD and 0<α1 ,

(60)
Dτα,ρφ(ξ,x,τ)=φξξ+φxx,φ(ξ,x,0)=sin(ξ)sin(x),
where 0ξ,x2π , τ>0.

By using the α‐Generalized Fractional Daftardar-Jafari Method (α-DJM), we get,

(61)
φ=φ(ξ,x,0)+Iταρ[φξξ+φxx],
(62)
φ=sin(ξ)sin(x)+Iταρ[φξξ+φxx],
(63)
F=sin(ξ)sin(x)=φo,
(64)
L(φ)=Iταρ[φξξ+φxx].

Let

(65)
φ=n=0φn.

Substituting Equation (38) into Equation (37) to obtain,

(66)
φn+1=Iταρ[2φnξ2+2φnx2],
(67)
φ1=Iταρ[2φoξ2+2φox2]=2sin(ξ)sin(x)ταρραΓ(α+1),
(68)
φ2=Iταρ[2φ1ξ2+2φ1x2]=4sin(ξ)sin(x)τ2αρρ2αΓ(2α+1),

As a result, this formulation serves as the approximate solution to the problem.

(69)
φ=sin(ξ)sin(x)[12ταρραΓ(α+1)+4τ2αρρ2αΓ(2α+1)+],

Therefore, the exact solution corresponding to Eq. (60) can be expressed as follows:

(70)
φ=sin(ξ)sin(x)e2τ,whenρ=α=1

Table 1 and Figure 1, show the convergence of the approximate solution to the exact solution for different α values. It seems that the closer α is to 1, the closer the approximate solution is to the exact solution. Figure 2 shows the comparison between the convergence of solutions individually for all different α values between the exact solution and the approximate solution.

Note 1: When using the Adomian method with the same operator, we will get the same approximate solution and exact solution. But the advantage of the proposed method appears in nonlinear examples, as we will explain in the second example.

Example 2.

Suppose that 0<α1 and Burger’s equation,

(71)
Da+α,ρφ(ξ,τ)+φφξ=φξξ,φ(ξ,0)=ξ,
where ξ0 , τ>0.

By using the α‐Generalized Fractional Daftardar-Jafari Method (α-DJM).

(72)
φ(ξ,τ)=φ(ξ,0)+Iταρ[φξξ]Iταρ[φφξ],
(73)
φ(ξ,τ)=ξ+Iταρ[φξξ]Iταρ[φφξ].
(74)
F=ξ=φo.

Table 1. Numerical values of the approximate solution φ(ξ,x,τ) (Eq. 69) for α=0.7,0.9,1andρ=1.5 at fixed ξandx,varyingτ .

τ φ0.7 φ0.9 φ1 φExact Error0.7 Error0.9 Error1
0.000000.042040.042040.042040.042040.000000.000000.00000
0.055560.038860.040830.041310.037610.001250.003210.00369
0.111110.035860.039040.040010.033660.002200.005380.00635
0.166670.033200.037020.038390.030120.003080.006900.00828
0.222220.030940.034900.036570.026950.003980.007950.00962
0.277780.029090.032790.034630.024120.004980.008670.01051
0.333330.027690.030750.032630.021580.006110.009170.01105
0.388890.026750.028850.030640.019310.007440.009540.01133
0.444440.026280.027140.028710.017280.009000.009860.01143
0.500000.026290.025670.026890.015460.010830.010200.01143
b584224c-70a3-4cf5-927c-4860c0d36eb4_figure1.gif

Figure 1. 2D comparison of the approximate (Eq. 69) and exact solution φ(ξ,x,τ) (Eq. 70) for α=0.7,0.9,1andρ=1.5 at fixed ξandx,varyingτ .

b584224c-70a3-4cf5-927c-4860c0d36eb4_figure2.gif

Figure 2. 3D comparison of the approximate (Eq. 69) and exact solution φ(ξ,x,τ) (Eq. 70) for α=0.7,0.9,1andρ=1.5 at fixed ξandx,varyingτ .

Let

(75)
φ=n=0φn,
(76)
φ(ξ,τ)=Iταρ[n=0φnξξ]Iταρ[n=0φnn=0φnξ].
(77)
L(φ)=Iταρ[φξξ],
(78)
G(φ)=Iταρ[φφξ].

The nonlinear terms are

(79)
Go=Iταρ[φ0ξφ0]=ξτραραΓ(α+1),
(80)
G1=Iταρ[(φ0+φ1)ξ(φ0+φ1)]+Iταρ[φ0ξφ0],=ξ(2τ2ραρ2αΓ(2α+1)τ3ραρ3α(Γ(α+1))2Γ(3α+1)),

The remaining terms can be obtained using the following recurrence relation

(81)
φn+1=L(φn)+G(φn).
(82)
φ1=L(φo)+G(φ0)=ξτραραΓ(α+1),
(83)
φ2=L(φ1)+G(φ1)=ξ(2τ2ραρ2αΓ(2α+1)τ3ραρ3α(Γ(α+1))2Γ(3α+1)).

As a result, this formulation serves as the approximate solution to the problem,

(84)
φ(ξ,τ)=ξ(1τραραΓ(α+1)+(2τ2ραρ2αΓ(2α+1)τ3ραΓ(2α+1)ρ3α(Γ(α+1))2Γ(3α+1))+).

To compare the advantage of the solution, we find the solution of Equation (71) using the Adomian decomposition method with the same operator. We will solve the following approximate solution:

(85)
φ(ξ,τ)=ξ(1τραραΓ(α+1)+2τ2ραρ2αΓ(2α+1)+).

Therefore, the exact solution corresponding to Eq. (95) can be expressed as follows:

(86)
φ=ξ1+τ,whenρ=α=1

Table 2, which presents the results using the Daftardar-Jafari Method. Table 3 shows the results using the Adomian decomposition Method. Both tables illustrate the convergence of the approximate solution toward the exact solution for different values of α. By comparing the two tables, we observe that the Daftardar-Jafari Method performs better than the Adomian Method when applied with the same operator.

Table 2. Numerical values of the approximate solution φ(ξ,τ) (Eq. 84) for α=0.7,0.9,1andρ=1 at fixed ξ, and varyingτ .

τ φ0.7 φ0.9 φ1 φExact Error0.7 Error0.9 Error1
05555000
0.055564.405244.646544.737484.736840.33160.09030.00063
0.111114.155534.390424.504754.50.344470.109580.00475
0.166674.005864.187854.300764.285710.279850.097860.01505
0.222223.914654.02744.124444.090910.176260.063510.03353
0.277783.862553.902513.974733.913040.05050.010530.06168
0.333333.838133.808563.850553.750.088130.058560.10055
0.388893.833793.741953.750853.60.233790.141950.15085
0.444443.844053.699693.674563.461540.382520.238160.21302
0.53.864753.679243.620613.333330.531420.345910.28727

Table 3. Numerical values of the approximate solution φ(ξ,τ) (Eq. 85) for α=0.7,0.9,1andρ=1 at fixed ξ, and varyingτ .

τ φ0.7 φ0.9 φ1 φExact Error0.7 Error0.9 Error1
05555000
0.055564.413154.64724.737654.736840.323690.089640.00081
0.111114.189464.394694.506174.50.310540.105310.00617
0.166674.085354.200614.305564.285710.200370.085110.01984
0.222224.060094.055144.13584.090910.030820.035760.04489
0.277784.094923.953193.996913.913040.181870.040140.08387
0.333334.17893.891463.888893.750.42890.141460.13889
0.388894.304823.867653.811733.60.704820.267650.21173
0.444444.467553.879963.765433.461541.006010.418420.30389
0.54.663213.9273.753.333331.329880.593670.41667

Figure 3 displays the convergence behavior of the approximate solution as α approaches 1, using the Daftardar-Jafari Method.

b584224c-70a3-4cf5-927c-4860c0d36eb4_figure3.gif

Figure 3. 2D comparison of the approximate (Eq.84) and exact solution φ(ξ,τ) (Eq.86) for α=0.7,0.9,1andρ=1 at fixed ξ,and varyingτ .

Figure 4 presents a comparison of the convergence between the exact and approximate solutions for the Daftardar-Jafari Method across all tested values of α.

b584224c-70a3-4cf5-927c-4860c0d36eb4_figure4.gif

Figure 4. 3D comparison of the approximate and exact solution φ(ξ,τ) (Eq.84) for α=0.7,0.9,1andρ=1 at fixed ξ,and varyingτ .

5. Discussion

This study presented two illustrative examples to evaluate the performance of the proposed α-generalized Daftardar-Jafari method both analytically and numerically.

In the first example, which involves a linear fractional partial differential equation, the numerical results obtained using the proposed method were comparable to those produced by the Adomian decomposition method. Both approaches yielded identical outcomes. However, from an analytical standpoint, the α-DJM framework offers a simpler and more direct formulation than many existing methods. Furthermore, the approximate solution was observed to converge toward the exact solution as the value of α approached 1. This behavior is clearly demonstrated in Figure 1 and Figure 2.

In contrast, the second example, which involves a nonlinear equation, revealed a distinct advantage of the proposed method. The numerical results obtained using α-DJM were more accurate than those produced by the Adomian method. This is evident from the comparative data presented in Table 2 and Table 3, respectively. Similar to the first example, the convergence of the approximate solution toward the exact solution improved as α approached 1, as illustrated in Figure 3 and Figure 4.

These findings confirm the effectiveness of the α-DJM method, particularly in handling nonlinear fractional models with improved accuracy and convergence behavior.

6. Conclusions

In this work, we introduced a new method for solving both linear and nonlinear fractional partial differential equations. The approach combines the α-generalized fractional differential operator with the Daftardar-Jafari method (DJM). The properties of the α-generalized operator, as presented in this article, demonstrate its generality. In particular cases, it reduces to well-known operators such as Caputo and Hadamard.

We applied the (α-DJM) framework to formulate a detailed solution scheme. This included deriving recurrence relations for both linear and nonlinear equations, and establishing convergence based on Lipschitz conditions through optimization techniques.

Our results show that solving a linear equation using different methods yields similar outcomes. This is clearly illustrated in the first example involving the linear heat equation. The analytical and numerical solutions are identical when substituting the parameters α=0.7,0.9,1andρ=1.5 for various values of τ. The second example addresses a nonlinear case involving the Berger equation. Upon solving and comparing with the Adomian decomposition method, we observed noticeable differences in both analytical and numerical results. In this case, we used α=0.7,0.9,1 and ρ=1 for both tables. From the comparative analysis, it is evident that the (α-DJM) method outperforms the Adomian method. It offers better accuracy and consistency.

Overall, the (α-DJM) combines analytical simplicity with numerical efficiency. It provides a unified framework for solving a wide range of fractional models with varying orders.

Future research may extend this framework to incorporate alternative methods and broader comparative studies.

Data citation

Sachit SA, Jassim HK (2025). Dataset for “An Iterative Approach for Solving Fractional Differential Equations Using the α-Generalized Daftardar–Jafari Method”. Zenodo. https://doi.org/10.5281/zenodo.17560489.

Ethics statement

This research does not involve human participants, animal subjects, or sensitive personal data. Therefore, ethical approval was not required.

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 23 Dec 2025
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Abdulkadhim Sachit S and Kamil Jassim H. An Iterative Approach for Solving Fractional Differential Equations Using the α-Generalized Daftardar–Jafari Method [version 1; peer review: awaiting peer review]. F1000Research 2025, 14:1437 (https://doi.org/10.12688/f1000research.172333.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

Current Reviewer Status:
AWAITING PEER REVIEW
AWAITING PEER REVIEW
?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 23 Dec 2025
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.