General Analytical Solutions to Fractional Differential Equations

Abstract

This research introduces a new global technique for addressing fractional differential equations (FDEs), referred to as the N-Transform (NT). This method serves as a comprehensive integral transform, encompassing both Laplace and Sumudu transforms. A fixed point theorem is used to prove the existence and uniqueness of solutions. Finding general solutions to specific fractional differential equations and proving their uniqueness through the application of a fixed point theorem are the main goals of this study. The fractional derivatives are expressed in the Caputo framework. The findings indicate that this method is highly effective, straightforward, and applicable to differential equations across various fields of science.

Share and Cite:

Gafel, H. (2026) General Analytical Solutions to Fractional Differential Equations. Journal of Applied Mathematics and Physics, 14, 373-383. doi: 10.4236/jamp.2026.141020.

1. Introduction

The idea of fractional calculus has become notably popular and acknowledged in recent years, due to its demonstrated uses in multiple areas of science and engineering. It presents various efficient techniques for tackling fractional differential equations (FDEs) and other associated challenges that involve special functions from mathematical physics, alongside their characteristics and generalizations in one or multiple dimensions. A significant connection exists between fractional calculus and the behavior of intricate real-world issues. Owing to their non-local characteristics, fractional operators provide a more thorough and organized depiction of numeric. In many mathematical models, fractional-order differential equations provide efficient control. Because classical mathematical models are particular cases of fractional-order mathematical models, the results from the fractional mathematical model are therefore more accurate and widely applicable. FDEs, or fractional differential equations, shed light on these physical phenomena. Physical processes are represented using fractional differential equations in a variety of ways, including: [1] [2] the fractional exponential function method, [3] the modified exponential function method, [4] the fractional He-Laplace method, [5] [6] the sub-equation method. In recent years, researchers have worked to provide more useful definitions for derivatives’ connotations. The goal of requiring these activities is to prevent norm derivatives from losing their ability to portray characters. These actions resulted in the development of several concepts for derivatives of fractional phenomena, which are logical extensions of normal derivatives of integer order. These derivative operators include the Caputo derivative and others. Recent theoretical, analytical, numerical, and computational contributions to this essentially multidisciplinary subject are the main focus of this study. Fractional calculus’s connections to other branches of mathematics and physics may open up new lines of inquiry, discoveries, and uses. This study presents analytical fresh investigations of fractional differential using a revolutionary approach known as NT. Previous methods such as LT and ST suggest that NT offers new generic fractional solutions for fractional differential equations. In some cases, NT is used to receive the classical forms of fractional differential equations. They aid in our understanding and interpretation of fractional differential equation and system phenomena. Here, we’ll review a few fundamental definitions of fractional derivatives:

Definition 1.1: [7] [8] The fractional derivative of Θ( Ι ) for Φ1<Φ<Φ , ΦΝ , Φ>0 , Α,Φ,ΙR in the Caputo sense is defined by:

D I Φ Θ( Ι )= J ΡΦ D Φ Θ( Ι )= 1 Γ( ΡΦ ) Α Ι ( tτ ) ΡΦ1 Θ ( Ρ ) ( Ι )dΙ,( Φ>0 ). (1)

Definition 1.2: [9] The Mittag-Leffler function E Φ ( Η ) with Φ>0 is defined via the series representation, valid in the whole complex plane is:

E Φ ( Η )= ι=0 Η ι Γ( Φι+1 ) . (2)

A generalization of the Mittag-Leffler function is

E Φ,r ( Η )= ι=0 Η ι Γ( Φι+r ) (3)

Definition 1.3: [10] For the function with Θ( Ι ) in ( , ) the general integral transform is defined by:

Ν[ Θ( Ι ) ]( ,ƛ )= J( ,Ι )Θ( Ι )dΙ , (4)

where J( ,Ι ) is the kernel of the transform and the variables ƛ and are the natural transform variables. When J( ,Ι )= e st Equation (4) gives Laplace transform. When J( ,Ι )= e t Equation (4) gives Sumudu transform.

Over the set of functions

Χ={ Θ( Ι ):Υ, Τ 1 , Τ 2 >0,| Θ( Ι ) |<Υ e | Ι |/ Τ j ,ifΙ ( 1 ) j ×[ 0, ) } (5)

The N-transform of Θ( Ι ) is defined by:

Ν[ Θ( Ι ) ]=( ,ƛ )= 0 Θ( ƛΙ ) e Ι dΙ ,ƛ>0,>0 (6)

where ( ,ƛ ) is the N-transform of the time function Θ( Ι ) and the variables ƛ and are the natural transform variables. The inversion formula of the N-transform is

Ν 1 [ ( ,ƛ ) ]=Θ( Ι )= 1 2πi ci c+i e Ι ƛ ( ,ƛ )d (7)

Definition 1.4: [10], If Ν[ Θ( Ι ) ] is the N-transform of the function [ Θ( Ι ) ] , then the N-transform of fractional derivative of order Φ is defined as:

Ν[ Θ ( Φ ) ( Ι ) ]= Φ ƛ Φ ( ,ƛ ) r=0 m1 Φ( r+1 ) ƛ Φr Θ ( r ) ( 0 ) (8)

More details of N-transform found in [11] [12]. In terms of fractional differential operators, the Caputo fractional derivative appears to be the most widely used type. Its definition centers on the kernel of the power law function convolution with a function’s local derivative. For the Caputo fractional differential equation, the derivative of a constant function is zero, hence the initial conditions can be properly expressed and handled identically to those of the ordinary derivative. In this study, we choose the Caputo fractional derivative because of these factors. Problems involving non-local characteristics and phenomena that take interactions into consideration benefit from the use of the Caputo fractional derivative. The equation can be thought of as having memory in this way. A finer feature of the Caputo fractional derivative is the memory effect. It is the most effective instrument for characterizing memories. Thus, we can say that the fractional Caputo derivative’s physical meaning is an indicator of memory. Since the Caputo differential operator shares a number of traits with usual ones, it is frequently employed in fractional calculus applications to represent a variety of physical events. Lastly, to mimic real-world issues, the Caputo derivative ought to be employed as the fractional operator. This is the main motivation behind the use of the Caputo fractional-order derivative operator in different physical implementations. Using the N-transform, the study’s primary goal is to solve fractional differential equations and provide general solutions. In addition to establishing the existence and uniqueness of solutions to fractional differential equations (FDEs) using fixed-point theory as follows:

First, consider a general fractional differential equation

A general Caputo fractional differential equation of 0<Φ<1 order can be written as:

D Φ Θ( Ι )=( Ι,Θ( Ι ) ),Ι[ 0,1 ],Θ( 0 )= Θ 0 , (9)

where D Φ is the Caputo fractional derivative and ( Ι,Θ( Ι ) ) is a given function.

Using fractional integration, we rewrite Equation (8) as aVoltera integral equation:

Θ( Ι )= Θ 0 + 1 Γ( Φ ) 0 Ι ( Ι ) Φ1 ( ,Θ( ) )d (10)

Secondly, we apply fixed-point theory to prove the existence and uniqueness. A fixed-point theory relies on Banach’s fixed-point theorem. It states that a contraction mapping on a complete metric space has a unique fixed point.

Step 1: Define an Operator

Define an operator on the space of continuous functions:

( Θ )( Ι )= Θ 0 + 1 Γ( Φ ) 0 Ι ( Ι ) Φ1 ( ,Θ( ) )d (11)

Step 2: Show is a Contraction

For uniqueness, we assume ( Ι,Θ ) satisfies a Lipschitz condition:

| ( Ι, Θ 1 )( Ι, Θ 2 ) || Θ 1 Θ 2 | (12)

Then, for two functions Θ 1 , Θ 2

| ( Θ 1 )( Ι )( Θ 2 )( Ι ) | Γ( Φ ) 0 Ι ( Ι ) Φ1 | Θ 1 ( ) Θ 2 ( ) |d (13)

Applying fractional Grönwall’s inequality, we get:

Θ 1 Θ 2 ω Θ 1 Θ 2 (14)

where ω<1 under small . This proves that is a contraction, and by Banach’s Fixed-Point Theorem, there exists a unique solution. The paper’s leftovers had been arranged as follows: The solutions to fractional differential equations and the existence of solutions to fractional differential equations and their uniqueness are presented in Section 2. The findings and discussion of the investigation are presented in Section 3.

2. Applications

In this section, we solve some fractional differential equations by NT to gain general solution of FDEs.

2.1. Problem 1

The Fractional exponential growth equation refers to an extension of the classical exponential growth model using fractional calculus, which deals with derivatives and integrals of non-integer (fractional) order. This is used when growth processes exhibit memory effects or anomalous dynamics, often seen in complex systems like biology, finance, or physics. Consider the following fractional exponential growth equation [13]:

d Φ Θ d t Φ =ΛΘ( Ι ),Θ( 0 )= Θ 0 ,0<Φ<1 , Λ is the growth rate (15)

Applying NT on both sides of Equation (14), we have

Ν[ d Φ Θ d Ι Φ ]=ΛΝ[ Θ( Ι ) ]. (16)

On simplifying 

α ƛ α ( ,ƛ ) α1 ƛ α Θ( 0 )=Λ( ,ƛ ) (17)

i.e.,

( ,ƛ )= 1 1 ƛ α s α Λ Θ 0 (18)

Taking the inverse NT of Equation (17), we obtain

Θ( Ι )= Ν 2 [ 1 ( 1 ƛ α α Λ ) 1 ] Θ 0 (19)

[14] obtained

Ν 1 [ 1 ( 1 ƛ α α Λ ) 1 ]= E Φ ( Λ Ι Φ ) (20)

where E Φ ( Λ Ι α )= ι=0 ( Λ Ι Φ ) ι Γ( Φι+1 ) Is the Mittag-Leffer function.

Then, the fractional solution of Equation (14) is

Θ( Ι )= Θ 0 [ 1+ Λ Ι Φ Φ! + ( Λ Ι Φ ) 2 2Φ! + ( Λ Ι Φ ) 3 3Φ! + ( Λ Ι Φ ) 4 4Φ! + ] (21)

Taking the case Φ=1 from Equation (21), we have the exact solution of Equation (15) obtained by [15]:

Θ( Ι )= Θ 0 e ΛΙ (22)

Figure 1. Plot of the NT solution of (15) at various values of Φ at Λ=1, Θ 0 =1 in comparing with the exact solution at Φ=1 obtained via [15].

Now, we can prove the existence and uniqueness of the solution of Equation (14) using fixed-point theory as follows:

Figure 2. Plot of the NT solution of (15) at various values of Φ at Λ=1, Θ 0 =1 in comparing with the exact solution obtained via [15].

Step 1: Consider the following fractional exponential growth equation

d Φ Θ d Ι Φ =ΛΘ( Ι ),Θ( 0 )= Θ 0 ,0<Φ<1,Ι[ 0,Ι ] (23)

Step 2: Convert Equation (22) into an integral equation

Caputo FDEs can be rewritten as a Valera-type integral equation:

Θ( Ι )= Θ 0 + Λ Γ( Φ ) 0 Ι ( Ι ) Φ1 Θ( )d (24)

Step 3: Define a fixed-point operator

Define an operator as:

( Θ )( Ι )= Θ 0 + Λ Γ( Φ ) 0 Ι ( Ι ) Φ1 Θ( )d (25)

Step 4: Prove uniqueness using Banach’s fixed-point theorem

Banach’s theorem states that if an operator is a contraction, then it has a unique fixed.

Step 4.1: Define a function space

Consider the Banach space C ([0, T]) with the norm:

Θ sup Ι[ 0,Ι ] | Θ( Ι ) | (26)

Applying fractional Grönwall’s inequality, we get:

Θ sup Ι[ 0,Ι ] | Θ( Ι ) | (27)

Step 4.2: Show is a contraction

Take two functions Θ 1 and Θ 2 compute:

| ( Θ 1 )( Ι )( Θ 2 )( Ι ) |=| Λ 0 Ι ( Ι ) Φ Γ( Φ ) ( Θ 2 ( ) Θ 1 ( ) )d | (28)

Applying the norm:

| ( Θ 1 )( Ι )( Θ 2 )( Ι ) | Λ Γ( Φ ) sup [ 0,Ι ] | Θ 1 ( ) Θ 2 ( ) | 0 Ι ( Ι ) Φ1 d (29)

Since

0 Ι ( Ι ) Φ1 d= Ι Φ Φ (30)

we get:

| ( Θ 1 )( Ι )( Θ 2 )( Ι ) | Λ Ι Φ ΦΓ( Φ ) Θ 1 Θ 2 (31)

If

ω= Λ Ι Φ ΦΓ( Φ ) <1 (32)

Then is a contraction, and by Banach’s fixed-point. So, there exists a unique solution of the fractional exponential growth problem.

2.2. Problem 2

Temperature difference in any situation results from energy flow into a system or energy flow from a system to surroundings. The former leads to heating, whereas the latter leads to cooling. Newton’s law of cooling states that the rate of change of temperature of the body is proportional to the difference between the temperature of the body and that of the surrounding medium. [16] Let Θ( τ ) be the temperature of the object at time τ , and η be the ambient temperature (constant). The fractional Newton’s Law of Cooling is given by:

D t α Θ( τ )=ε( Θ( τ )η ),0<Φ1, (33)

D t α is the Caputo fractional derivative of order Φ , Θ( τ ) temperature of the object, ε>0 is the cooling rate constant, and Θ( 0 )= Θ 0 is the initial temperature of the object. Applying NT on both sides of Equation (32), we have

Φ ƛ Φ ( ,ƛ ) Φ1 ƛ Φ Θ( 0 )=ε( ,ƛ )+ εη (34)

On simplifying

( Φ ƛ Φ +ε )( ,ƛ )= Φ1 ƛ Φ Θ 0 + εη (35)

i.e.,

( ,ƛ )= 1 ƛ 1 ƛ ( ƛ ) Φ ( ƛ ) Φ + ƛ Φ Φ ε Θ 0 + 1 ƛ 1 ƛ ε ( ƛ ) Φ +ε η (36)

Taking the inverse NT of Equation (35), we obtain

Θ( τ )= Ν 1 [ 1 ƛ 1 ƛ ( ƛ ) Φ ( ƛ ) Φ + ƛ Φ Φ ε Θ 0 + 1 ƛ 1 ƛ ε ( ƛ ) Φ +ε η ] (37)

[14] obtained

Ν 1 [ 1 ƛ 1 ƛ ( ƛ ) Φ ( ƛ ) Φ + ƛ Φ Φ ε ]= Ε Φ ( ε τ Φ ), (38)

Ν 1 [ 1 ƛ 1 ƛ ε ( ƛ ) Φ +ε η ]=1 Ε α ( Κ t α ) (39)

Then, the fractional solution of Equation (32) is

Θ( τ )= Ε α ( ε τ Φ ) Θ 0 +( 1 Ε α ( ε τ α ) )η (40)

Then,

Θ( τ )=η+( Θ 0 η ) Ε α ( ε τ α ) (41)

Θ( τ )=η+( Θ 0 η ) e ετ (42)

This is the exact solution obtained by [16].

Figure 3. Plot of the NT solution of (32) at various values of Φ at ε=0.1, Θ 0 =100,η=25 in comparing with the exact solution at Φ=1 obtained via [16].

Figure 4. Plot of the NT solution of (32) at various values of Φ at ε=0.1, Θ 0 =100,η=25 in comparing with the solution via STVIM and the exact solution obtained via [16].

Now, we can prove the existence and uniqueness of the solution of Equation (32) using fixed-point theory as follows:

Step 1: Consider the following fractional Newton’s law of cooling is given by:

D t α Θ( τ )=ε( Θ( τ )η ),0<Φ1, (43)

Step 2: Convert Equation (41) into an integral equation

Caputo FDEs can be rewritten as a Valera-type integral equation:

Θ( τ )= Θ 0 ε Γ( Φ ) 0 τ ( Θ( )η ) ( τ ) Φ1 d (44)

This is a good starting point for applying fixed-point theory.

Step 3: Define a fixed-point operator

Define an operator as:

( Θ )( τ )= Θ 0 ε Γ( Φ ) 0 τ ( Θ( )η ) ( τ ) Φ1 d (45)

Step 4: Prove uniqueness using Banach’s fixed-point. It states that if an operator is a contraction, then it has a unique fixed point, as follows:

Step 4.1: Define a function space

Consider the Banach space C ([0, R]) with the norm:

Θ sup τ[ 0,Ι ] | Θ( τ ) | (46)

Step 4.2: Show is a Contraction

Take two functions Θ 1 and Θ 2 compute:

| ( Θ 1 )( τ )( Θ 2 )( τ ) |=| ε 0 τ ( τ ) Φ Γ( Φ ) ( Θ 2 ( ) Θ 1 ( ) )d | (47)

Applying the norm:

| ( Θ 1 )( τ )( Θ 2 )( τ ) | ε Γ( Φ ) sup [ 0,Ι ] | Θ 1 ( ) Θ 2 ( ) | 0 Ι ( Ι ) Φ1 d (48)

Since

0 I ( I ) Φ1 d= I Φ Φ (49)

we get:

| ( Θ 1 )( τ )( Θ 2 )( τ ) | ε Ι Φ ΦΓ( Φ ) Θ 1 Θ 2 (50)

If

ω= ε τ Φ ΦΓ( Φ ) <1 (51)

Then is a contraction, and by Banach’s Fixed-Point, there exists a unique solution of the fractional exponential growth problem.

3. Conclusion

This approach displayed a novel technique named the N-Transform (NT). It is a general algorithm of the Laplace transform and the Sumudu transform. It gives novel general solutions to fractional differential equations in different fields of science. The results give no vel generic fractional solutions to exponential growth equation and Newton’s law of cooling. To clarify the generalization of the algorithm, we compared the results of the current methods with those of other methods like the Sumudu transform method ( Τ ), the Laplace transform method ( Τ ), and the exact solutions. At =1 , in Equation (19) and Equation (37),we have the solutions via Τ . At ƛ=1 , in Equation (19) and Equation (37), we have solutions via Τ . At Φ1 we have the classical form solutions. Fractional exponential growth equation and Newton’s law of cooling were studied analytically and numerically and novel general solutions were gained via the proposed algorithm. The existence and uniqueness of the solutions were proved using Banach’s fixed-point theorem. Figures 1-4 clarify the simplicity, quality, and generalization of the proposed algorithm (NT). The outcomes showed the significance of the proposed method (NT) for solving various differential equations in different aspects of science. The NT method has many advantages compared with other methods. The NT method is more effective, efficient, accurate and can easily handle a wide class of linear or nonlinear differential equations. The NT method demonstrates fast convergence of the solution. It will be successfully used to handle most types of differential equations that appear in several physical models and scientific applications. Its method attacks the problem in a direct way and in a straightforward fashion without using linearization, perturbation, or any other restrictive assumption. Lastly, the NT method is a novel generic method that can be employed to solve problems in different aspects of science.

Acknowledgements

The researcher would like to acknowledge the Deanship of Graduate Studies and Scientific Research at Taif University for funding this work.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Alotaibi, H., Abo-Dahab, S.M., Abdlrahim, H.R. and Kilany, A.A. (2020) Fractional Calculus of Thermoelastic P-Waves Reflection under Influence of Gravity and Electromagnetic Fields. Fractals, 28, Article ID: 2040037.[CrossRef]
[2] Gurefe, Y., Pandir, Y. and Akturk, T. (2022) On the Nonlinear Mathematical Model Representing the Coriolis Effect. Mathematical Problems in Engineering, 2022, Article ID: 2504907.[CrossRef]
[3] Pandir, Y., Akturk, T., Gurefe, Y. and Juya, H. (2023) The Modified Exponential Function Method for Beta Time Fractional Biswas-Arshed Equation. Advances in Mathematical Physics, 2023, Article ID: 1091355.[CrossRef]
[4] Habib, S., Batool, A., Islam, A., Nadeem, M., Gepreel, K.A. and He, J. (2021) Study of Nonlinear Hirota-Satsuma Coupled KdV and Coupled mKdV System with Time Fractional Derivative. Fractals, 29, Article ID: 2150108.[CrossRef]
[5] Nadeem, M., Feng, G., Islam, A. and He, C. (2022) Chemical Reaction and Radiation on Boundary-Layer Flow of Electrically Conduction Micropolar Fluid through a Porous Shrinking Sheet. Thermal Science, 26, 2593-2598.[CrossRef]
[6] Yépez-Martínez, H. and Gómez-Aguilar, J.F. (2019) Optical Solitons Solution of Resonance Nonlinear Schrödinger Type Equation with Atangana’s-Conformable Derivative Using Sub-Equation Method. Waves in Random and Complex Media, 31, 573-596.[CrossRef]
[7] Podlubny, I. (1999) Fractional Differential Equations, Mathematics in SCIENCE and Engineering. Academic Press.
[8] Caputo, M. and Fabrizio, M. (2015) A New Definition of Fractional Derivative without Singular Kernel Fractional Differential Equations. Progress in Fractional Differentiation and Applications, 2, 73-85.
[9] Mittag-Leffler, G. (1903) Sur la Nouvelle Fonction Eα(x). Comptes Rendus de lAcademie des Sciences Paris, 137, 554-558.
[10] Khan, Z.H. and Khan, W.A. (2008) N-Transform-Properties and Applications. NUST Journal of Engineering Sciences, 1, 127-133.
[11] Abdl-Rahim, H.R., Zayed, M. and Ismail, G.M. (2022) Analytical Study of Fractional Epidemic Model via Natural Transform Homotopy Analysis Method. Symmetry, 14, Article 1695.[CrossRef]
[12] Ismail, G.M., Abdl-Rahim, H.R., Abdel-Aty, A., Kharabsheh, R., Alharbi, W. and Abdel-Aty, M. (2020) An Analytical Solution for Fractional Oscillator in a Resisting Medium. Chaos, Solitons & Fractals, 130, Article ID: 109395.[CrossRef]
[13] Omran, M. and Kiliçman, A. (2016) Natural Transform of Fractional Order and Some Properties. Cogent Mathematics, 3, Article ID: 1251874.[CrossRef]
[14] Kilbas, A.A., Srivastava, H.M. and Trujillo, J.J. (2006) Theory and Applications of Fractional Differential Equations. Elsevier.
[15] Han, H. and Kim, H. (2020) The Solution of Exponential Growth and Exponential Decay by Using Laplace Transform. International Journal of Difference Equations, 15, 191-195.[CrossRef]
[16] Bhangale, N. and Kachhia, K.B. (2022) A Fractional Calculus Approach to Study Newton’s Law of Cooling. Progress in Fractional Differentiation and Applications, 8, 275-287.

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.