Next Article in Journal
The Fourier–Legendre Series of Bessel Functions of the First Kind and the Summed Series Involving 1F2 Hypergeometric Functions That Arise from Them
Previous Article in Journal
An Algorithm Based on Compute Unified Device Architecture for Estimating Covering Functionals of Convex Bodies
Previous Article in Special Issue
Stability Analysis in a Mathematical Model for Allergic Reactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method

by
Ghaliah Alhamzi
1,
Aafrin Gouri
2,*,
Badr Saad T. Alkahtani
3 and
Ravi Shanker Dubey
2,*
1
Department of Mathematics and Statistics, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11989, Saudi Arabia
2
Department of Mathematics, Amity School of Applied Sciences, Amity University Rajasthan, Jaipur 302002, Rajasthan, India
3
Department of Mathematics, College of Science, King Saud University, Riyadh 11989, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Axioms 2024, 13(2), 133; https://doi.org/10.3390/axioms13020133
Submission received: 26 December 2023 / Revised: 8 February 2024 / Accepted: 11 February 2024 / Published: 19 February 2024

Abstract

:
In this study, we present the generalized form of the higher-order nonlinear fractional Bratu-type equation. In this generalization, we deal with a generalized fractional derivative, which is quite useful from an application point of view. Furthermore, some special cases of the generalized fractional Bratu equation are recognized and examined. To solve these nonlinear differential equations of fractional order, we employ the homotopy perturbation transform method. This work presents a useful computational method for solving these equations and advances our understanding of them. We also plot some numerical outcomes to show the efficiency of the obtained results.

1. Introduction

Numerous areas of pure and applied mathematics focus largely on fractional calculus, which works with derivatives and integrals of arbitrary order [1,2,3,4,5]. Fractional differential equations have been used more frequently in applied mathematics and physics, as well as in the modeling and interpretation of various real-world issues. Biological science, physical science, chemical engineering, and other fields of research indicate nonlinear perspectives [6,7,8,9].
In mathematical modeling, the nonlinear fractional differential equations are quite helpful. However, it becomes more challenging to solve nonlinear systems of equations. It has been seen that there are many differential equations without an exact solution. Various innovative numerical and analytical approaches have been established to solve fractional differential equations, leading to numerical and analytical solutions to these types of problems. These techniques consist of Adomian decomposition method [10,11,12], homotopy analysis method [13,14,15]), generalized differential transform approach [16], homotopy perturbation method [17,18], modified Laplace decomposition method [19], homotopy perturbation transform method [20,21], and q-homotopy analysis transform method [22,23], etc. The homotopy perturbation method is the most effective for solving the differential equations of non-linear type because of its straightforward approach and fastest convergence rate. This technique was developed by the famous mathematician Mr. He [24,25,26]. Its main benefit is that the solution determined in series form approaches its exact solution very quickly.
Many mathematicians have developed an appropriate definition and formula for fractional calculus. Major contributions to the theory of fractional calculus with singular kernels have been made by Riemann, Caputo, Liouville, Weyl, Kilbas, and others. Gorenglo, Miller–Ross, Yang—Abdel–Cattani, Atangana–Baleanu, Wiman, and many others encouraged the study of fractional integrals and derivatives having nonsingular kernels. Miller–Ross, Yang–Abdel–Cattani, Wiman, Atangana–Baleanu, and many others promoted the study of fractional integrals and derivatives with non-singular kernels. Riemann-Liouville and Caputo are the two fractional operators that are used most frequently. However, there are several problems that are not properly described by these operators. For fractional derivatives, Caputo-Fabrizio developed a new operator in 2015 with a nonsingular kernel [27]. This nonsingular kernel makes the Caputo-Fabrizio fractional derivative stronger than the Caputo derivative.
The calculus operators of fractional type have been divided into three groups by Gomez and Atangana: weak, medium, and strong [28]. They suggest that the Riemann-Liouville fractional derivative is significantly more useful in explaining physical difficulties than the Caputo version. However, the Mittag-Leffler memory pattern of the Atangana-Baleanu fractional operator provides a good clarification. Some authors have used fractional derivatives of Atangana-Baleanu-Caputo to characterize fractional differential equations [29].
A numerical slab model of a combustion problem is examined using the Bratu-type equation [30]. Several additional issues have been solved using Bratu’s problem, such as the Chandrasekhar model, the fuel ignition model, the thermal response model, and the framework for the electrospun nano-fiber creation method [31,32]. To solve the Bratu type problem, Ghazanfari and Sephvandzadeh used the homotopy perturbation method [33], the modified variational iteration method [34], and the Adomian decomposition method [35]. Recent studies have explored numerical approaches to addressing the fractional Bratu initial value problem. Investigations have been conducted using the Bézier curve [36], the fractional differential transform method [37], and a technique grounded in the CAS wavelet scheme [38] for solving fractional Bratu-type differential equations.
Applications for nonlinear differential equations, like the Bratu-type equation, can be found in many branches of science and engineering. It is used in heat transfer to investigate temperature distributions in materials, in combustion theory to understand flame propagation, and in chemical kinetics to replicate complex reaction kinetics. The modeling of complex phenomena like population dynamics, material diffusion, reactor neutron transport, electrochemical reactions, and pollutant spread is made possible by its involvement in mathematical biology, material science, nuclear physics, electrochemistry, and environmental science. In addition to having particular uses, the Bratu-type equation is a basic mathematical model that provides an understanding of the characteristics of nonlinear differential equations. Because of this, it is a flexible and essential instrument for scientific analysis and study. The “Nonlinear Local Fractional Bratu-Type Equation” is a delightful and theoretically complex differential equation that belongs to the field of fractional calculus.
As mathematicians, physicists, and engineers grapple with issues related to non-standard calculus, nonlinear dynamics, and complex systems, they are entranced by ongoing research projects centered around the “Non-linear Local Fractional Bratu-type Equation”. Having a thorough grasp of the characteristics and solutions of this equation could improve our understanding of a variety of real-world processes and possibly lead to useful applications in fields including chemical kinetics, biological modeling, and combustion theory. In the fields of science and engineering, the Bratu-type equation is quite useful. In all domains where standard differentiable calculus is relevant, this article presents some new insights into the concept of generalization. Other expansions of the differentiability concept can be found in [39,40,41]. Very few researchers have studied the generalized form of a fractional Bratu-type equation, so finding the solution of a generalized fractional Bratu-type differential equation is a new thing in itself. Realizing how important it is, we turn our attention to examining a generalized fractional Bratu-type equation.
This study aims to define the generalized Bratu-type differential equation and expand the application of the homotopy perturbation transform method for identifying specific features of the generalized Bratu-type differential equation
D t 2 n κ u ( t ) + δ E β ( u ) = 0 , 0 < κ 1 , 0 < t 1 ,
where δ is constant and E β ( u ) is the Mittag-Leffler function defined as
E β ( u ) = k = 0 u k Γ ( β k + 1 ) ,
Γ indicates the gamma function and D 2 n κ = D κ D κ D κ ( 2 n t i m e s ) w h e r e D κ u = d κ u d t κ .
The structure of this paper is divided into seven sections: We provide a brief introduction and some fundamental terminology in Section 1. The basic homotopy perturbation transform approach is defined in Section 2. In Section 3, we discuss the existence and uniqueness of the generalized fractional Bratu-type differential equation. In Section 4, we describe the main findings of this work. In Section 5, we discuss the convergence analysis. In Section 6, we describe some special cases of generalized Bratu-type equations and provide a graphic representation of the solution. This paper is concluded in the last section.

Basic Definitions

This section provides the definitions and foundational results essential for the content presented in this paper.
A theory that deals with generalized fractional-order integration and differentiation operators is called generalized fractional calculus. Using an analogy with the concept of the Caputo and R–L fractional derivative, the analogous generalized fractional derivatives are defined as [42]
  0 C D t κ f ( t ) = 0 t f ˙ ( u ) Δ l ( t u ) d u ,   0 D t κ f ( t ) = d d t 0 t f ( u ) Δ l ( t u ) d u ,
where 0 < κ < 1 is the order of derivative, f : [ 0 , + ) R is the absolute continuous function with f ˙ L loc 1 ( 0 , + ) , Δ l is known as the general kernel function, also for f L loc 1 ( 0 , + ) ,
  0 D t κ   0 C D t κ f ( t ) = f ( t ) f ( 0 ) ,
where   0 D t κ represent the general form of the Riemann-Liouville fractional integral defined as
  0 D t κ f ( t ) = 0 t f ( u ) l ( t u ) d u .
Definition 1. 
Caputo derivative of fractional order κ is denoted by   a C D t κ and defined by [43]
  a C D t κ f ( t ) = 1 Γ ( n κ ) a t ( t u ) n κ 1 f ( n ) ( u ) d u ,
where 0 < κ < 1 a n d n 1 < κ < n . Also f ( n ) ( t ) = d n f d t n .
Definition 2. 
Caputo-Fabrizio derivative is denoted by   a C F D t κ and defined by [44]
  a C F D t κ f ( t ) = Δ ( κ ) ( 1 κ ) a t exp κ ( 1 κ ) ( t u ) f ( u ) d u ,
where Δ ( κ ) = Normalization function with Δ ( 0 ) = 1 , Δ ( 1 ) = 1 .
Definition 3. 
We have the Atangana-Baleanu derivative which is denoted by   a A B C D t κ and defined by [45]
  a A B C D t κ f ( t ) = Δ ( κ ) ( 1 κ ) a t E κ κ ( 1 κ ) ( t u ) κ f ( u ) d u ,
where Δ ( κ ) = Normalization function with Δ ( 0 ) = 1 , Δ ( 1 ) = 1 and E κ is Mittag-leffler function.
Definition 4. 
If A = { f ( t ) | M , τ 1 , τ 2 > 0 , | f ( t ) | < M exp (|t|/τj), if t ∈ (−1)j × [0,∞)} is the set of functions, then the Sumudu transform on A is defined as [46]
S [ f ( t ) ] ( u ) = 0 e t f ( u t ) d t .
Definition 5. 
In accordance with the Reference [47], we define the Sumudu transform for the Caputo derivative of order κ > 0
S   0 C D t κ f ( t ) ( s ) = S [ f ] f ( 0 ) s κ .
Definition 6. 
In accordance with the Reference [47], we define the Sumudu transform for the Caputo-Fabrizio derivative of order κ > 0
S   0 C F D t κ f ( t ) ( s ) = Δ ( κ ) 1 κ + κ s S [ f ] f ( 0 ) ,
where Δ ( κ ) = Normalization function with Δ ( 0 ) = 1 , Δ ( 1 ) = 1 .
Definition 7. 
In accordance with the Reference [47], we define the Sumudu transform for the Atangana-Baleanu fractional derivative of order κ > 0
S   0 A B C D t κ f ( t ) ( s ) = Δ ( κ ) 1 κ + κ s κ S [ f ] f ( 0 ) ,
where Δ ( κ ) = Normalization function with Δ ( 0 ) = 1 , Δ ( 1 ) = 1 .

2. Basic Homotopy Perturbation Transform Approach

This section covers the Homotopy Perturbation Transform Method (HPTM) procedure for solving fractional differential equations using generalized fractional operators.
The Homotopy Perturbation Technique (HPT), incorporating He’s polynomials and a transform algorithm, is formulated as the homotopy perturbation transform method. To illustrate the fundamental workings of the homotopy perturbation transform method, we apply it to a nonlinear fractional-order differential equation of the form
D t κ u ( t ) + L u ( t ) + N u ( t ) = f ( t ) , 0 < κ 1 , t > 0 .
Here, u(t) represents the probability density function, D t κ denotes the generalized fractional derivative with order κ , L stands for the linear differential operator with an order less than that of D, and N represents the general nonlinear differential operator.
Apply the Sumudu transform on both sides of (9),
S [ D t κ u ( t ) ] + S [ L u ( t ) ] + S [ N u ( t ) ] = S [ f ( t ) ] ,
as we know, Sumudu transform for generalized fractional derivatives [48],
S [ D t κ u ( t ) ] ( s ) = ς 1 S [ u ( t ) ] k = 0 n 1 s k [ u k ( t ) ] t = 0 , 0 < κ 1 ,
here,
ς = s κ ( f o r C a p u t o f r a c t i o n a l d e r i v a t i v e ) , ς = 1 κ + κ s Δ ( κ ) ( f o r C a p u t o F a b r i z i o f r a c t i o n a l d e r i v a t i v e ) , ς = 1 κ + κ s κ Δ ( κ ) ( f o r A t a n g a n a B a l e a n u C a p u t o f r a c t i o n a l d e r i v a t i v e ) ,
and Δ ( κ ) is the normalization function.
Hence
S [ u ( t ) ] = k = 0 n 1 s k [ u k ( t ) ] t = 0 ς S [ L u ( t ) ] + S [ N u ( t ) ] S [ f ( t ) ] .
Now, applying the inverse Sumudu transform, we obtain
u ( t ) = ϕ ( t ) S 1 ς S [ L u ( t ) ] + S [ N u ( t ) ] ,
where
ϕ ( t ) = S 1 k = 0 n 1 s k [ u k ( t ) ] t = 0 + ς S [ f ( t ) ] .
Now, expanding u(t) into the series of power p as u ( t ) = n = 0 p n u n ( t ) and consider N u ( t ) = n = 0 p n H n ( u ) , here H n ( u ) is He’s polynomial defined as
H n ( u 0 , u 1 , u 2 , , u n ) = 1 n ! p n N k = 0 p k u k p = 0 , n = 0 , 1 , 2 , 3 , ,
and p ( 0 , 1 ) is an embedding parameter.
Then
n = 0 p n u n ( t ) = ϕ ( t ) p S 1 ς S L n = 0 p n u n ( t ) + S n = 0 p n H n ( u ) .
This equation represents the combination of the homotopy perturbation method, He’s polynomial, and Sumudu transform.
Now, by equating the coefficients corresponding to different powers of p, we obtain
p 0 : u 0 ( t ) = ϕ ( t ) , p 1 : u 1 ( t ) = S 1 ς S L u 0 ( t ) + S H 0 ( u ) , p 2 : u 2 ( t ) = S 1 ς S L u 1 ( t ) + S H 1 ( u ) , . . .
The remaining elements of u n ( t ) , n 0 can be evaluated by applying a similar approach. As a result, the homotopy perturbation transforms method solutions are completely computed.
Thus, the approximate analytical solution of (9) is
u ( t ) = lim p 1 n = 0 p n u n ( t ) = n = 0 u n ( t ) .

3. Existence and Uniqueness of the Solution

Theorem 1. 
The solution of the generalized Bratu-type fractional differential equation defined by (1) satisfies the Lipschitz condition and has the unique solution for 0 < b < 1 , where b = ν δ Γ ( 1 + κ ) t κ is known as the Lipschitz constant.
Proof. 
Firstly, we write Equation (1) in the following form
u ( t ) = δ Γ ( 1 + κ ) 0 t G τ , u ( τ ) d τ κ
where
G t , u ( t ) = 1 Γ ( 1 + κ ) 2 n 1 0 t 0 t . . . 0 t E β u ( τ ) ( 2 n 1 ) t i m e s d τ κ d τ κ d τ κ
The non-linear part G t , u ( t ) is Lipschitz continuous with the condition
G t , u G t , v ν u v .
Lipschitz constant can be evaluated by making use of maximum norm G = max 0 t 1 G t , u ( t ) , we have
E β ( u ) E β ( v ) k = 0 u k v k Γ β k + 1 ,
k = 0 u v Γ β k + 1 u k 1 + u k 2 v + . . . + u v k 2 + v k 1 .
Since the series is convergent u k 1 + u k 2 v + . . . + u v k 2 + v k 1 is bounded for each k.
Thus,
u k 1 + u k 2 v + . . . + u v k 2 + v k 1 M , k = 1 , 2 , 3 , ,
where M is some constant.
Thus,
E β ( u ) E β ( v ) u v M k = 0 1 Γ β k + 1 ,
u v M E β ( 1 ) .
Now from (12), we obtain
G t , u G t , v 1 Γ ( 1 + κ ) 2 n 1 M E β ( 1 ) u v 0 t 0 t . . . 0 t d τ κ ( 2 n 1 ) t i m e s d τ κ d τ κ , G t , u G t , v 1 Γ ( 1 + κ ) 2 n 1 M E β ( 1 ) u v t ( 2 n 1 ) κ ( 2 n 1 ) ! , M E β ( 1 ) ( 2 n 1 ) ! Γ ( 1 + κ ) 2 n 1 u v .
Thus, we can select the Lipschitz constant as ν = M E β ( 1 ) ( 2 n 1 ) ! Γ ( 1 + κ ) 2 n 1 .
Hence, the solution of Equation (1) exists.
Now, suppose that Equation (1) has two different solutions u 1 and u 2 , so then
u 1 u 2 = δ Γ ( 1 + κ ) 0 t G τ , u 1 ( τ ) d τ κ δ Γ ( 1 + κ ) 0 t G τ , u 2 ( τ ) d τ κ , δ Γ ( 1 + κ ) ν u 1 u 2 t κ .
This produces the following inequality,
u 1 u 2 1 ν δ Γ ( 1 + κ ) t κ 0
i.e., u 1 u 2 1 b 0 , where b = ν δ Γ ( 1 + κ ) t κ .
Thus, u 1 u 2 0 for 0 < b < 1 which shows that u 1 = u 2 . This proves that the solution of (1) is unique. □

4. Main Result

The solution to the generalized Bratu-type equation using the homotopy perturbation transform method is given in this section.
Theorem 2. 
The solution of the generalized Bratu-type fractional differential equation
D t 2 n κ u ( t ) + δ E β ( u ) = 0 , 0 < κ 1 , 0 < t 1 ,
is given by the series u = n = 0 u n ( t ) .
Where δ is constant, E β ( u ) is the Mittag-Leffler function, and D 2 n κ is a generalized differential operator with order 2 n .
Proof. 
We have the generalized Bratu-type fractional differential equation
D t 2 n κ u ( t ) + δ E β ( u ) = 0 .
Applying Sumudu transform,
S D t 2 n κ u ( t ) + δ S E β ( u ) = 0 ,
or
ς 1 S [ u ( t ) ] k = 0 n 1 s k [ u k ( t ) ] t = 0 = δ S E β ( u ) ,
or
S u ( t ) = k = 0 n 1 s k [ u k ( t ) ] t = 0 ς δ S E β ( u ) .
Here,
ς = s 2 n κ ( f o r C a p u t o f r a c t i o n a l d e r i v a t i v e ) , ς = ( 1 2 n κ + 2 n κ s ) Δ ( 2 n κ ) ( f o r C a p u t o F a b r i z i o f r a c t i o n a l d e r i v a t i v e ) , ς = ( 1 2 n κ + 2 n κ s 2 n κ ) Δ ( 2 n κ ) ( f o r A t a n g a n a B a l e a n u C a p u t o f r a c t i o n a l d e r i v a t i v e ) ,
where Δ ( κ ) is the normalization function.
Applying the inverse Sumudu transform
u ( t ) = S 1 k = 0 n 1 s k [ u k ( t ) ] t = 0 δ S 1 ς S E β ( u ) .
Now, using the homotopy perturbation transform method, Equation (15) becomes
n = 0 p n u n ( t ) = ϕ ( t ) δ p S 1 ς S n = 0 p n H n ( u ) ,
where ϕ ( t ) = S 1 k = 0 n 1 s k [ u k ( t ) ] t = 0 .
We have the Taylor’s series expansion of E β ( u ) ,
E β ( u ) = k = 0 u k Γ ( β k + 1 ) .
Also, we have,
H 0 ( u ) = 1 + 1 Γ β + 1 u 0 + 1 Γ 2 β + 1 u 0 2 + 1 Γ 3 β + 1 u 0 3 + , H 1 ( u ) = 1 Γ β + 1 u 1 + 1 Γ 2 β + 1 ( 2 u 0 u 1 ) + 1 Γ 3 β + 1 ( 3 u 0 2 u 1 ) + , H 2 ( u ) = 1 Γ β + 1 u 2 + 1 Γ 2 β + 1 ( u 1 2 + 2 u 0 u 2 ) + 1 Γ 3 β + 1 ( 3 u 1 2 u 0 ) + , . . .
Now comparing the coefficient of the power of p in (16), we obtain
p 0 : u 0 ( t ) = ϕ ( t ) ,   p 1 : u 1 ( t ) = δ S 1 ς S H 0 ( u ) ,   p 2 : u 2 ( t ) = δ S 1 ς S H 1 ( u ) ,   p 3 : u 3 ( t ) = δ S 1 ς S H 2 ( u ) , . . .
The remaining elements of u n ( t ) , n 0 can be evaluated by applying a similar approach through using (17) wherein the series solutions are completely determined.
We have
u ( t ) = lim p 1 n = 0 p n u n ( t )
u ( t ) = u 0 ( t ) + u 1 ( t ) + u 2 ( t ) +
The generalized differential operator reduces into three well-known differential operators, namely, the Atangana-Baleanu-Caputo, Caputo, and Caputo-Fabrizio fractional operators. Here, we discuss all three cases of a generalized Bratu-type fractional differential equation. □

5. Convergence Analysis

In this section, we examine the convergence of the solution to Equation (1).
Theorem 3. 
Let u n ( t ) and u ( t ) be defined in Banach space, B = C 0 , 1 , . , then the HPTM series n = 0 u n ( t ) converges to u ( t ) of Equation (1) if there exists ζ 0 , 1 such that u n + 1 ( t ) ζ n N .
Proof. 
Let us consider how the sequence of partial sums of series (18) is S n ( t ) n 0 defined by
S l ( t ) = l = 0 u l ( t ) .
Then
S n + 1 ( t ) S n ( t ) = u n + 1 ( t ) ζ u n ( t ) ζ 2 u n 1 ( t ) ζ n + 1 u 0 ( t ) .
Now, first we have to show that S n ( t ) be a Cauchy sequence in Banach space B = C 0 , 1 , . .
S q ( t ) S r ( t ) = S q ( t ) S q 1 ( t ) + S q 1 ( t ) S q 2 ( t ) + . . . + S r + 1 ( t ) S r ( t ) , S q ( t ) S q 1 ( t ) + S q 1 ( t ) S q 2 ( t ) + S r + 1 ( t ) S r ( t ) , ζ q u 0 ( t ) + ζ q 1 u 0 ( t ) + + ζ r + 1 u 0 ( t ) , ζ r + 1 1 + ζ + ζ 2 + . . . + ζ q + . . . u 0 ( t ) , ζ r + 1 1 ζ q r 1 ζ u 0 ( t ) .
Since, ζ 0 , 1 , so 1 ζ q r > 0 , then
S q ( t ) S r ( t ) ζ r + 1 1 ζ u 0 ( t ) .
Thus, S q ( t ) S r ( t ) q , r 0 as u 0 is bounded.
This S n ( t ) is a Cauchy sequence in Banach space and hence convergent.
Therefore, there exists u B such that
l = 0 u l ( t ) = u ( t ) .
We achieve the desired result. □

6. Special Cases

We have the normal generalized fractional operator. However, it may be reduced into three well-known fractional operators: Atangana-Baleanu-Caputo, Caputo, and Caputo-Fabrizio. To derive the results, we separately discuss the cases for generalized operators.
If we take n = 1, then (1) becomes
D t 2 κ u ( t ) + δ E β ( u ) = 0 , 0 < κ 1 , 0 < t 1 .
Moreover, consider the initial point u(0) = 0. This equation is known as a particular form of the generalized fractional Bratu equation.
We offer the solution to the generalized fractional Bratu equation in this section by using the process defined in Section 4.
Case 1. 
Consider generalized Bratu’s initial value problem for the Caputo fractional derivative,
  C D t 2 κ u ( t ) + δ E β ( u ) = 0 ,
with u(0) = 0.
Applying Sumudu transform on (23),
S   C D t 2 κ u ( t ) + δ S E β ( u ) = 0 ,
or
s 2 κ S u ( t ) u ( 0 ) = δ S E β ( u ) ,
or
S u ( t ) = u ( 0 ) δ s 2 κ S E β ( u ) .
Applying the inverse Sumudu transform,
u ( t ) = u ( 0 ) δ S 1 s 2 κ S E β ( u ) .
Now, using the homotopy perturbation transform method and expanding u(t) as u ( t ) = n = 0 p n u n ( t ) , then (24) becomes
n = 0 p n u n ( t ) = u ( 0 ) δ p S 1 s 2 κ S n = 0 p n H n ( u ) .
Now, we have the Taylor’s series expansion of E β ( u ) ,
E β ( u ) = k = 0 u k Γ ( β k + 1 ) .
Now, comparing the coefficients in (25), we obtain
p 0 : u 0 ( t ) = u ( 0 ) , p 1 : u 1 ( t ) = δ S 1 s 2 κ S H 0 ( u ) , p 2 : u 2 ( t ) = δ S 1 s 2 κ S H 1 ( u ) , p 3 : u 3 ( t ) = δ S 1 s 2 κ S H 2 ( u ) , . . .
Using u(0) = 0, we have
p 0 : u 0 ( t ) = 0 ,
p 1 : u 1 ( t ) = δ S 1 s 2 κ S 1 = δ t 2 κ Γ ( 2 κ + 1 ) ,
p 2 : u 2 ( t ) = δ Γ ( β + 1 ) S 1 s 2 κ S δ t 2 κ Γ ( 2 κ + 1 ) = δ 2 Γ ( β + 1 ) t 4 κ Γ ( 4 κ + 1 ) ,
p 3 : u 3 ( t ) = δ S 1 s 2 κ S u 2 Γ ( β + 1 ) + u 1 2 Γ ( 2 β + 1 ) = δ 3 1 Γ ( β + 1 ) 2 + 1 Γ ( 2 β + 1 ) Γ ( 4 κ + 1 ) Γ ( 2 κ + 1 ) 2 t 6 κ Γ ( 6 κ + 1 ) ,
. . .
Now, the solution is given by (18),
u ( t ) = ( δ ) t 2 κ Γ ( 2 κ + 1 ) + δ 2 Γ ( β + 1 ) t 4 κ Γ ( 4 κ + 1 ) + δ 3 1 Γ ( β + 1 ) 2 + 1 Γ ( 2 β + 1 ) Γ ( 4 κ + 1 ) Γ ( 2 κ + 1 ) 2 t 6 κ Γ ( 6 κ + 1 ) + . . .
Now to plot the results by considering δ = 1.0, β = 0.7 for different values of κ as illustrated in Figure 1.
Remark 1. 
If we put β = 1 in (23), then the equation becomes
  C D t 2 κ u ( t ) + δ e u = 0 ,
and the solution is
u ( t ) = ( δ ) t 2 κ Γ ( 2 κ + 1 ) + ( δ ) 2 t 4 κ Γ ( 4 κ + 1 ) + δ 3 1 + 1 2 Γ ( 4 κ + 1 ) Γ ( 2 κ + 1 ) 2 t 6 κ Γ ( 6 κ + 1 ) +
This solution for the fractional Bratu’s differential equation exactly matches the result of the analytical study of the fractional Bratu-type equation by Dubey et al. [49].
Now to plot the results by considering δ = 1.0, β = 1.0 for different values of, κ as illustrated in Figure 2.
Case 2. 
Consider generalized Bratu’s initial value problem for the Caputo-Fabrizio derivative,
  C F D t 2 κ u ( t ) + δ E β ( u ) = 0 ,
with u(0) = 0.
Applying Sumudu transform on (28),
S   C F D t 2 κ u ( t ) + δ S E β ( u ) = 0 ,
or
Δ ( 2 κ ) 1 2 κ + 2 κ s S u ( t ) u ( 0 ) = δ S E β ( u ) ,
or
S u ( t ) = u ( 0 ) δ 1 2 κ + 2 κ s Δ ( 2 κ ) S E β ( u ) .
Applying the inverse Sumudu transform,
u ( t ) = u ( 0 ) δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s S E β ( u ) .
Now, using the homotopy perturbation transform method and expanding u(t) as u ( t ) = n = 0 p n u n ( t ) , we have
n = 0 p n u n ( t ) = u ( 0 ) δ p Δ ( 2 κ ) S 1 1 2 κ + 2 κ s S n = 0 p n H n ( u ) .
Using the Taylor’s series expansion of E β ( u ) and comparing the coefficients in (28), we obtain
p 0 : u 0 ( t ) = u ( 0 ) , p 1 : u 1 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s S H 0 ( u ) , p 2 : u 2 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s S H 1 ( u ) , p 3 : u 3 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s S H 2 ( u ) , . . .
Using u(0) = 0, we have
p 0 : u 0 ( t ) = 0 , p 1 : u 1 ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t , p 2 : u 2 ( t ) = δ Δ ( 2 κ ) 2 1 Γ ( β + 1 ) 1 2 κ 2 + 4 κ 1 2 κ t + 2 κ 2 t 2 , p 3 : u 3 ( t ) = δ Δ ( 2 κ ) 3 1 Γ ( β + 1 ) 2 + 1 Γ ( 2 β + 1 ) 1 2 κ 3 + 6 κ 1 2 κ 2 t + 2 3 Γ ( β + 1 ) 2 + 4 Γ ( 2 β + 1 ) 1 2 κ κ 2 t 2 + 8 3 ! 1 Γ ( β + 1 ) 2 + 2 Γ ( 2 β + 1 ) κ 3 t 3 , . . .
Now, the solution is given by (18),
u ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t + δ Δ ( 2 κ ) 2 1 Γ ( β + 1 ) 1 2 κ 2 + 4 κ 1 2 κ t + 2 κ 2 t 2 + δ Δ ( 2 κ ) 3 1 Γ ( β + 1 ) 2 + 1 Γ ( 2 β + 1 ) 1 2 κ 3 + 6 κ 1 2 κ 2 t + 2 3 Γ ( β + 1 ) 2 + 4 Γ ( 2 β + 1 ) 1 2 κ κ 2 t 2 + 8 3 ! 1 Γ ( β + 1 ) 2 + 2 Γ ( 2 β + 1 ) κ 3 t 3 +
Now to plot the results by considering δ = 1.0, β = 0.7 for different values of, κ as illustrated in Figure 3.
Remark 2. 
If we put β = 1 in (26), then the equation becomes
  C F D t 2 κ u ( t ) + δ e u = 0 ,
and the solution is given by (18),
u ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t + δ Δ ( 2 κ ) 2 1 2 κ 2 + 4 κ 1 2 κ t + 2 κ 2 t 2 + δ Δ ( 2 κ ) 3 3 2 1 2 κ 3 + 9 κ 1 2 κ 2 t + 10 1 2 κ κ 2 t 2 + 8 3 κ 3 t 3 +
Now to plot the results by considering δ = 1.0, β = 1.0 for different values of, κ as illustrated in Figure 4.
Case 3. 
Consider generalized Bratu’s initial value problem for the Atangana-Baleanu-Caputo derivative,
  A B C D t 2 κ u ( t ) + δ E β ( u ) = 0 ,
with u(0) = 0.
Applying the Sumudu transform on (29),
S   A B C D t 2 κ u ( t ) + δ S E β ( u ) = 0 ,
or
Δ ( 2 κ ) 1 2 κ + 2 κ s 2 κ S u ( t ) u ( 0 ) = δ S E β ( u ) ,
or
S u ( t ) = u ( 0 ) δ 1 2 κ + 2 κ s 2 κ Δ ( 2 κ ) S E β ( u ) .
Applying the inverse Sumudu transform,
u ( t ) = u ( 0 ) δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s 2 κ S E β ( u ) .
Now, using the homotopy perturbation transform method, we have
n = 0 p n u n ( t ) = u ( 0 ) δ p Δ ( 2 κ ) S 1 1 2 κ + 2 κ s 2 κ S n = 0 p n H n ( u ) .
Using the Taylor’s series expansion of E β ( u ) and comparing the coefficients in (30), we obtain
p 0 : u 0 ( t ) = u ( 0 ) , p 1 : u 1 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s 2 κ S H 0 ( u ) , p 2 : u 2 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s 2 κ S H 1 ( u ) , p 3 : u 3 ( t ) = δ Δ ( 2 κ ) S 1 1 2 κ + 2 κ s 2 κ S H 2 ( u ) , . . .
Using u(0) = 0, we have
p 0 : u 0 ( t ) = 0 , p 1 : u 1 ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t 2 κ Γ ( 2 κ + 1 ) , p 2 : u 2 ( t ) = δ Δ ( 2 κ ) 2 1 Γ ( β + 1 ) 1 2 κ 2 + 4 κ 1 2 κ t 2 κ Γ ( 2 κ + 1 ) + 4 κ 2 t 4 κ Γ ( 4 κ + 1 ) , . . .
Now, the solution is given by (18),
u ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t 2 κ Γ ( 2 κ + 1 ) + δ Δ ( 2 κ ) 2 1 Γ ( β + 1 ) 1 2 κ 2 + 4 κ 1 2 κ t 2 κ Γ ( 2 κ + 1 ) + 4 κ 2 t 4 κ Γ ( 4 κ + 1 ) +
Now to plot the results by considering δ = 1.0, β = 0.7 for different values of, κ as illustrated in Figure 5.
Remark 3. 
If we put β = 1 in (29), then the equation becomes
  A B C D t 2 κ u ( t ) + δ e u = 0 ,
and the solution is given by (18),
u ( t ) = δ Δ ( 2 κ ) ( 1 2 κ ) + 2 κ t 2 κ Γ ( 2 κ + 1 ) + δ Δ ( 2 κ ) 1 2 κ 2 + 4 κ 1 2 κ t 2 κ Γ ( 2 κ + 1 ) + 4 κ 2 t 4 κ Γ ( 4 κ + 1 ) + . . .
Now to plot the results by considering δ = 1.0, β = 1.0 for different values of, κ as illustrated in Figure 6.

7. Conclusions

In this investigation, we have introduced and addressed the generalized nonlinear Bratu-type equation incorporating the generalized fractional derivatives operator. The solution is derived by applying an innovative analytical technique known as the homotopy perturbation transformation method. We have demonstrated the effectiveness of the generalized fractional derivatives operator and presented reduced results expressed in well-known fractional operators, including Caputo, Caputo-Fabrizio, and Atangana-Baleanu. Furthermore, various special cases of the generalized fractional Bratu’s equation have been identified and thoroughly examined. We also plot some numerical results to show the behavior of the obtained results.

Author Contributions

G.A. organised the necessary research materials, project administration, and Funding acquisition. The methodology, was directed by R.S.D., who also constructed the validation, and software. A.G. did the formal analysis, and investigation. While B.S.T.A. prepared the Resources and carried out data curation. He also did the project administration and supervision. The draft was read, corrected, and polished by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-RP23066).

Data Availability Statement

Availability of data is already cited in the article.

Acknowledgments

We would like to express our sincere thanks to the editor and reviewers for their valuable comments and suggestions. The authors thank the Dean of Scientific Research at Imam Mohammad Ibn Saud Islamic University for their support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; A Wiley-Interscience Publication: New York, NY, USA, 1993; ISBN 0-471-58884-9. [Google Scholar]
  2. Oldham, K.B.; Spanier, J. The Fractional Calculus-Theory and Applications of Differentiation and Integration to Arbitrary Order, 1st ed.; Mathematics in Science and Engineering: New York, NY, USA; London, UK, 1974; ISBN 978-012-525-550-9. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; North-Holland Mathematics Studies: Amsterdam, The Netherlands, 2006; ISBN 978-0-444-51832-3. [Google Scholar]
  4. Podlubny, I. Fractional Differential Equations—An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, 1st ed.; Mathematics in Science and Engineering: New York, NY, USA; London, UK, 1999; ISBN 978-0-12-558840-9. [Google Scholar]
  5. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Beijing, China, 2014; ISBN 978-981-4579-89-6. [Google Scholar]
  6. Yang, X.J. Advanced Local Fractional Calculus and Its Applications; World Science Publisher: New York, NY, USA, 2012; ISBN 978-1-938576-01-0. [Google Scholar]
  7. Yang, X.J.; Gao, F.; Srivastava, H.M. Exact travelling wave solutions for the local fractional two-dimensional Burger’s type equations. J. Comput. Appl. Math. 2017, 73, 203–210. [Google Scholar] [CrossRef]
  8. Atangana, A. On the new fractional derivative and application to non-linear Fisher’s reaction-diffusion equation. Appl. Math. Comput. 2016, 273, 948–956. [Google Scholar] [CrossRef]
  9. Saad, K.M.; Atangana, A.; Baleanu, D. New fractional derivatives with non-singular kernel applied to the Burgers equation. Chaos 2018, 28, 63–109. [Google Scholar] [CrossRef]
  10. Daftardar-Gejji, V.; Jafari, H. Solving a multi-order fractional differential equation using adomian decomposition. Appl. Math. Comput. 2007, 189, 541–548. [Google Scholar] [CrossRef]
  11. Wazwaz, A.M. A reliable modification of Adomian decomposition method. Appl. Math. Comput. 1999, 102, 77–86. [Google Scholar] [CrossRef]
  12. Biazar, J.; Babolian, E.; Islam, R. Solution of the system of ordinary differential equations by Adomian decomposition method. Appl. Math. Comput. 2004, 147, 713–719. [Google Scholar] [CrossRef]
  13. Hashim, I.; Abdulaziz, O.; Momani, S. Homotopy Analysis Method for Fractional IVPs. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 674–684. [Google Scholar] [CrossRef]
  14. Liao, S. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. 2014, 147, 499–513. [Google Scholar] [CrossRef]
  15. Liao, S. Notes on the homotopy analysis method: Some definitions and theorems. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 983–997. [Google Scholar] [CrossRef]
  16. Odibat, Z.; Momani, S.; Erturk, V.S. Generalized differential transform method for solving a space- and time-fractional diffusion-wave equation. Phys. Lett. A 2007, 370, 379–387. [Google Scholar] [CrossRef]
  17. Abdulaziz, O.; Hashim, I.; Momani, S. Solving systems of fractional differential equations by homotopy-perturbation method. Phys. Lett. A 2008, 372, 451–459. [Google Scholar] [CrossRef]
  18. He, J.H. Homotopy perturbation method: A new nonlinear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar] [CrossRef]
  19. Khan, Y. An effective modification of the laplace decomposition method for nonlinear equations. Int. J. Nonlinear Sci. Numer. Simul. 2009, 10, 1373–1376. [Google Scholar] [CrossRef]
  20. Khan, Y. Homotopy perturbation transform method for nonlinear equations using He’s polynomials. Comput. Math. Appl. 2011, 61, 1963–1967. [Google Scholar] [CrossRef]
  21. Khan, M.; Gondal, M.A.; Hussain, I.; Vanani, S.K. A new comparative study between homotopy analysis transform method and homotopy perturbation transform method on a semi infinite domain. Math. Comput. Model. 2012, 55, 1143–1150. [Google Scholar] [CrossRef]
  22. Singh, J.; Kumar, D.; Swroop, R. Numerical solution of time- and space-fractional coupled Burger’s equations via homotopy algorithm. Alex. Eng. J. 2016, 55, 1753–1763. [Google Scholar] [CrossRef]
  23. Mohamed, M.S.; Hamed, Y.S. Solving the convection–diffusion equation by means of the optimal q-homotopy analysis method (Oq-HAM). Results Phys. 2016, 6, 20–25. [Google Scholar] [CrossRef]
  24. He, J.H. Homotopy perturbation technique. Comput. Methods Appl. Mech. Eng. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  25. He, J.H. Homotopy perturbation method for solving boundary value problems. Phys. Lett. A 2006, 350, 87–88. [Google Scholar] [CrossRef]
  26. He, J.H. The homotopy perturbation method for nonlinear oscillators with discontinuities. Appl. Math. Comput. 2004, 151, 287–292. [Google Scholar] [CrossRef]
  27. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  28. Gomez-Aguilar, J.F.; Atangana, A. New insight in fractional differentiation: Power, exponential decay and Mittag–Leffler laws and applications. Eur. Phys. J. Plus 2017, 132, 1–21. [Google Scholar] [CrossRef]
  29. Gomez-Aguilar, J.F.; Atangana, A. Decolonisation of fractional calculus rules: Breaking commutativity and associativity to capture more natural phenomena. Eur. Phys. J. Plus 2018, 133, 166. [Google Scholar] [CrossRef]
  30. Singh, J.; Kumar, D.; Swroop, R.; Kumar, S. An efficient computational approach for time-fractional Rosenau–Hyman equation. Neural Comput. Appl. 2017, 30, 3063–3070. [Google Scholar] [CrossRef]
  31. Ascher, U.M.; Matheij, R.; Russell, R.D. Numerical Solution of Boundary Value Problems for Ordinary Differential Equations; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1995; ISBN 978-089-871-354-1. [Google Scholar]
  32. Chandrasekhar, S. An Introduction to the Study of Stellar Structure; Dover: Chicago, IL, USA, 1957; ISBN 978-0-486-60413-8. [Google Scholar]
  33. Ghazanfari, B.; Sepahvandzadeh, A. Homotopy perturbation method for solving fractional Bratu-type equation. J. Math. Model. 2015, 2, 143–155. [Google Scholar]
  34. Ghazanfari, B.; Sepahvandzadeh, A. Solving fractional Bratu-type equations by modified variational iteration method. Selcuk. J. Appl. Math. 2014, 39, 23–29. [Google Scholar]
  35. Ghazanfari, B.; Sepahvandzadeh, A. Adomian decomposition method for solving fractional Bratu-type equations. J. Math. Comput. Sci. 2014, 8, 236–244. [Google Scholar] [CrossRef]
  36. Ghomanjani, F.; Shateyi, S. Numerical solution for fractional Bratu’s initial value problem. Open Phys. 2017, 15, 1045–1048. [Google Scholar] [CrossRef]
  37. Demir, D.D.; Zeybek, A. The numerical solution of fractional Bratu-type differential equations. ITM Web Conf. 2017, 13, 1–11. [Google Scholar] [CrossRef]
  38. Yi, M.; Sun, K.; Huang, J.; Wang, L. Numerical solution of fractional integrodifferential equations of bratu type by using CAS wavelets. J. Appl. Math. 2013, 2013, 801395. [Google Scholar] [CrossRef]
  39. Kukushkin, M.V. Abstract Fractional Calculus for m-accretive Operators. Int. J. Appl. Math. 2021, 34, 1–41. [Google Scholar] [CrossRef]
  40. Kukushkin, M.V. On Solvability of the Sonin–Abel Equation in the Weighted Lebesgue Space. Fractal Fract. 2021, 5, 77. [Google Scholar] [CrossRef]
  41. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications, 1st ed.; Gordon and Breach, Yverdon; CRC Press: Boca Raton, FL, USA, 1993; ISBN 2881248640. [Google Scholar]
  42. Luchko, Y.; Yamamoto, M. General time-fractional diffusion equation: Some uniqueness and existence results for the initial boundary value problems. Fract. Calc. Appl. Anal. 2016, 19, 676–695. [Google Scholar] [CrossRef]
  43. Oliveira, D.S.; Oliveira, E.C.D. On a Caputo-type fractional derivative. Adv. Pure Appl. Math. 2018, 10, 81–91. [Google Scholar] [CrossRef]
  44. Qureshi, S.; Rangaig, N.A.; Baleanu, D. New Numerical Aspects of Caputo-Fabrizio Fractional Derivative Operator. Mathematics 2019, 7, 374. [Google Scholar] [CrossRef]
  45. Yadav, S.; Pandey, R.K.; Shukla, A.K. Numerical approximations of Atangana–Baleanu Caputo derivative and its application. Chaos Solitons Fractals 2019, 118, 58–64. [Google Scholar] [CrossRef]
  46. Demiray, S.T.; Bulut, H.; Belgacem, F.B.M. Sumudu Transform Method for Analytical Solutions of Fractional Type Ordinary Differential Equations. Math. Probl. Eng. 2014, 2015, 131690. [Google Scholar] [CrossRef]
  47. Atangana, A.; Akgul, A. Can transfer function and Bode diagram be obtained from Sumudu transform. Alex. Eng. J. 2019, 59, 1971–1984. [Google Scholar] [CrossRef]
  48. Bodkhe, D.S.; Panchal, S.K. On Sumudu Transform Of Fractional Derivatives and Its Applications to Fractional Differential Equations. Asian J. Math. Comput. 2016, 11, 69–77. [Google Scholar]
  49. Dubey, V.P.; Kumar, R.; Kumar, D. Analytical study of fractional Bratu-type equation arising in electro-spun organic nanofibers elaboration. Phys. A 2019, 521, 762–772. [Google Scholar] [CrossRef]
Figure 1. Change in u for Case 1 with respect to time, for different values of κ .
Figure 1. Change in u for Case 1 with respect to time, for different values of κ .
Axioms 13 00133 g001
Figure 2. Change in u for remark of Case 1 with respect to time, for different values of κ .
Figure 2. Change in u for remark of Case 1 with respect to time, for different values of κ .
Axioms 13 00133 g002
Figure 3. Change in u for Case 2 with respect to time, for different values of κ .
Figure 3. Change in u for Case 2 with respect to time, for different values of κ .
Axioms 13 00133 g003
Figure 4. Change in u for remark of Case 2 with respect to time, for different values of κ .
Figure 4. Change in u for remark of Case 2 with respect to time, for different values of κ .
Axioms 13 00133 g004
Figure 5. Change in u for Case 3 with respect to time, for different values of κ .
Figure 5. Change in u for Case 3 with respect to time, for different values of κ .
Axioms 13 00133 g005
Figure 6. Change in u for remark of Case 3 with respect to time, for different values of κ .
Figure 6. Change in u for remark of Case 3 with respect to time, for different values of κ .
Axioms 13 00133 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alhamzi, G.; Gouri, A.; Alkahtani, B.S.T.; Dubey, R.S. Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method. Axioms 2024, 13, 133. https://doi.org/10.3390/axioms13020133

AMA Style

Alhamzi G, Gouri A, Alkahtani BST, Dubey RS. Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method. Axioms. 2024; 13(2):133. https://doi.org/10.3390/axioms13020133

Chicago/Turabian Style

Alhamzi, Ghaliah, Aafrin Gouri, Badr Saad T. Alkahtani, and Ravi Shanker Dubey. 2024. "Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method" Axioms 13, no. 2: 133. https://doi.org/10.3390/axioms13020133

APA Style

Alhamzi, G., Gouri, A., Alkahtani, B. S. T., & Dubey, R. S. (2024). Analytical Solution of Generalized Bratu-Type Fractional Differential Equations Using the Homotopy Perturbation Transform Method. Axioms, 13(2), 133. https://doi.org/10.3390/axioms13020133

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop