šŸ•·ļø Crawler Inspector

URL Lookup

Direct Parameter Lookup

Raw Queries and Responses

1. Shard Calculation

Query:
Response:
Calculated Shard: 31 (from laksa130)

2. Crawled Status Check

Query:
Response:

3. Robots.txt Check

Query:
Response:

4. Spam/Ban Check

Query:
Response:

5. Seen Status Check

ā„¹ļø Skipped - page is already crawled

šŸ“„
INDEXABLE
āœ…
CRAWLED
28 days ago
šŸ¤–
ROBOTS ALLOWED

Page Info Filters

FilterStatusConditionDetails
HTTP statusPASSdownload_http_code = 200HTTP 200
Age cutoffPASSdownload_stamp > now() - 6 MONTH1 months ago
History dropPASSisNull(history_drop_reason)No drop reason
Spam/banPASSfh_dont_index != 1 AND ml_spam_score = 0ml_spam_score=0
CanonicalPASSmeta_canonical IS NULL OR = '' OR = src_unparsedNot set

Page Details

PropertyValue
URLhttps://www.mdpi.com/2227-7390/14/6/984
Last Crawled2026-03-15 02:15:21 (28 days ago)
First Indexednot set
HTTP Status Code200
Meta TitleAn Iterative Technique for Solving Time Fractional Systems with Variable Coefficients | MDPI
Meta DescriptionThis paper presents a comprehensive study on constructing exact and approximate solutions to Cauchy problems for time-fractional systems with variable coefficients.
Meta Canonicalnull
Boilerpipe Text
Article 13 March 2026 and School of Mathematics, Jilin University, Changchun 130012, China * Author to whom correspondence should be addressed. Abstract This paper presents a comprehensive study on constructing exact and approximate solutions to Cauchy problems for time-fractional systems with variable coefficients. An innovative iterative approach is developed for solving functional equations with initial conditions, combining rigorous mathematical foundations with practical computational efficiency. The proposed technique effectively handles the nonlocal nature of fractional operators through a carefully designed iterative scheme that maintains simplicity while achieving high accuracy. It demonstrates particular strength in solving nonlinear systems with well-defined conditions and variable coefficients, where traditional methods often fail. Through systematic theoretical analysis and numerical validation, we establish the method’s convergence properties and computational advantages, showing its capability to generate both exact closed-form solutions, when available, and high-precision approximations otherwise. The approach remains computationally tractable even for complex cases where variable coefficients and memory effects of fractional systems present significant challenges to conventional solution approaches. MSC: 35R11; 34A12 1. Introduction This paper presents a systematic investigation of the initial value problems arising in time-fractional partial differential equations, where the nonlocal nature of fractional derivatives introduces unique mathematical challenges that fundamentally distinguish them from their integer-order counterparts. The inherent memory effects and historical dependence of fractional operators render the initial value problem particularly intricate yet fascinating from both theoretical and computational perspectives. Our comprehensive study develops and analyzes a suite of analytical and numerical methodologies specifically designed to overcome these challenges, with the dual objectives of establishing rigorous mathematical foundations and providing practical computational tools for understanding the temporal evolution of solutions in such systems. The research builds upon and significantly extends existing work in this rapidly developing field, as documented in [ 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 ] and related references, while introducing novel approaches to address the distinctive features of fractional-order initial value problems. The historical development of fractional calculus provides crucial context for understanding its modern applications. Our exploration begins with the seminal correspondence between L’HĆ“pital and Leibniz in 1695, a foundational moment in mathematical history documented in [ 9 ]. In this historic exchange, L’HĆ“pital posed the now-famous question to Leibniz: ā€œWhat would be the result of d n y d x n if n = 1 2 ?ā€ Leibniz’s insightful reply demonstrated remarkable foresight: ā€œ d 1 / 2 x will be equal to x d x : x . This is an apparent paradox, from which, one day useful consequences will be drawn.ā€ This early conceptual breakthrough, however, remained largely undeveloped for nearly two centuries. The mathematical community of the Enlightenment era faced several barriers to advancing fractional calculus: The lack of a rigorous theoretical framework for non-integer derivatives; The absence of clear physical interpretations; Significant computational challenges; The dominance of integer-order calculus in solving contemporary problems. This historical evolution underscores how theoretical mathematics often precedes practical applications, with Leibniz’s initial insight eventually finding numerous applications across science and engineering disciplines. The fundamental nature of fractional operators distinguishes them significantly from traditional differential operators. Researchers have established that these operators are intrinsically non-local, typically expressed through integral formulations with singular kernels. This essential characteristic gives fractional calculus its unique capacity to model memory effects and hereditary properties in physical systems, where the current state depends on the complete historical evolution rather than just instantaneous conditions. A particularly noteworthy application emerges in electrochemical systems, where fractional-order operators have proven remarkably effective: They provide superior modeling of diffusion processes in electrochemical cells compared to classical approaches; They more accurately characterize electrode impedance behavior; They offer enhanced precision in describing charge transfer mechanisms. The theoretical framework of fractional calculus has consequently found extensive utility across numerous scientific and technological domains, including: modeling anomalous diffusion and viscoelastic material properties, analyzing complex reaction-diffusion processes, understanding neuronal dynamics and biomechanical tissue behavior, characterizing long-range dependent processes in financial markets and so on. These diverse applications, meticulously documented in [ 10 ] and extensively developed in subsequent research, vividly illustrate the remarkable transformation of fractional calculus from a purely theoretical construct to an indispensable analytical framework in modern science and engineering. The continually expanding scope of its applications provides compelling evidence for the fundamental importance of these non-local operators, with demonstrated successes ranging from nanoscale phenomena to large-scale economic systems. Their unique capacity to capture long-range dependencies and memory effects has enabled breakthroughs in modeling complex systems that defy description by classical approaches, particularly in situations involving power-law behaviors and anomalous dynamics. Recent advances have particularly highlighted how fractional operators offer superior fidelity in modeling multiscale phenomena with coupled temporal and spatial dependencies, while also providing natural frameworks for systems exhibiting complex hierarchical or fractal-like structures. The operational success of fractional calculus across such disparate domains not only underscores its profound theoretical significance but has also driven important innovations in computational methods and analytical techniques. This dual impact—simultaneously deepening our fundamental understanding while enhancing practical problem-solving capabilities—continues to motivate vigorous research efforts at the frontiers of applied mathematics and interdisciplinary science, promising further discoveries and applications in emerging fields. Fractional calculus, representing a significant generalization of classical differential and integral operators to non-integer orders, has emerged as a vibrant research field attracting considerable attention from both theoretical and applied perspectives. This mathematical framework extends the traditional notions of differentiation and integration to fractional orders, thereby enabling the description of complex phenomena that exhibit non-local behavior and memory effects. Throughout its historical development, numerous distinct definitions of fractional derivatives and integrals have been proposed, each offering unique advantages for specific applications and mathematical contexts. The absence of a unified definition for fractional derivatives has led to a rich diversity of operational formulations, including several prominent variants that have proven particularly useful in practice. Among these, the Grünwald–Letnikov approach provides a discrete interpretation based on finite differences, while the Riemann–Liouville formulation offers a continuous integral-based perspective. The Caputo derivative has gained prominence for its advantageous treatment of initial conditions in differential equations. Additional important formulations include the Riesz potential for symmetric fractional operations, the conformable fractional derivative for maintaining certain classical properties, and various local fractional derivatives designed for specific applications. Within this diverse landscape, the Riemann–Liouville fractional integral and Caputo fractional derivative have emerged as particularly influential, owing to their mathematical tractability and physical interpretability. These operators have found extensive application across numerous scientific and engineering disciplines, including advanced control systems where they enable more robust controller designs, the modeling of anomalous diffusion processes in complex media, and the analysis of fractal structures and multi-scale systems. The continued development and refinement of these fractional operators, as documented in [ 4 , 7 , 8 , 11 , 12 , 13 , 14 ] and related works, reflects their growing importance in addressing challenging problems that lie beyond the scope of classical calculus. Among the various formulations of fractional operators, the Grünwald–Letnikov fractional derivative (see Appendix A.2 ) occupies a special historical and conceptual position as one of the most intuitive and computationally accessible approaches. This operator, rooted in the generalization of finite difference methods, provides a natural discrete framework for fractional differentiation that closely parallels the classical limit definition of integer-order derivatives. The Grünwald–Letnikov approach is particularly valuable in numerical analysis and scientific computing, where its discrete nature allows for efficient implementation in algorithmic solutions to fractional differential equations. Its computational advantages have made it indispensable for simulations of complex systems with memory effects, as extensively documented in [ 15 ] and subsequent computational studies. In contrast, the Riemann–Liouville fractional derivative represents a continuous formulation defined through a convolution operation with a power-law kernel. This definition enjoys several theoretical advantages, including a rigorous mathematical foundation in functional analysis and a natural connection to the theory of integral equations. The Riemann–Liouville operator has proven particularly powerful in analytical approaches to fractional differential equations arising in mathematical physics, where its properties facilitate the development of existence theorems and solution representations. The complementary strengths of these two approaches—with Grünwald–Letnikov favoring computational implementations and Riemann–Liouville supporting theoretical analysis—serve to illustrate the rich diversity of tools available in fractional calculus, as discussed in detail in [ 16 ] and related mathematical physics literature. In contrast, the Caputo fractional derivative was introduced to address this issue, offering a more physically intuitive way to incorporate initial conditions in models. It is particularly popular in applied settings such as control systems. The Riesz potential, often employed in potential theory and the fractional Laplacian, generalizes the idea of a potential field to fractional dimensions and is crucial in studying processes like anomalous diffusion and LĆ©vy flights; see [ 17 ] and the references therein. In contrast to the Riemann–Liouville formulation, the Caputo fractional derivative was specifically developed to overcome certain mathematical challenges in modeling physical systems. This alternative definition provides several crucial advantages for practical applications, most notably its ability to handle initial conditions in a manner that aligns more naturally with physical interpretations. The Caputo operator’s formulation ensures that initial value problems can be specified using standard function values rather than fractional derivatives of the initial state, making it particularly valuable for modeling real-world phenomena. This property has led to its widespread adoption in applied mathematics and engineering disciplines, especially in control theory, where it enables more natural formulations of dynamical systems with memory effects. The Caputo derivative’s compatibility with conventional physical interpretations has made it a preferred choice for modeling viscoelastic materials, thermal systems, and other processes where the initial state is typically specified in terms of classical variables. Similarly important in fractional analysis is the Riesz potential, which extends classical potential theory to fractional dimensions through an elegant symmetric formulation. This operator plays a fundamental role in the theory of the fractional Laplacians and has become indispensable for studying non-local phenomena. The Riesz transform provides a powerful mathematical framework for analyzing complex diffusion processes, particularly those exhibiting anomalous behavior that deviates from classical Brownian motion. Its applications span diverse areas, including the modeling of LĆ©vy flights in statistical physics, the analysis of non-local interactions in quantum mechanics, and the study of long-range correlations in complex systems. The theoretical foundations and practical applications of the Riesz potential are thoroughly examined in [ 17 ], along with its connections to other fractional operators and its growing importance in modern mathematical physics. The rich diversity of fractional derivative definitions, while reflecting the wide-ranging applications of fractional calculus across scientific disciplines, presents significant challenges for researchers working to select, compare, and apply these various operators. This proliferation of definitions—each with distinct properties and suited for different applications—can create confusion and complicate theoretical comparisons between results obtained using different formulations. Recognizing this challenge, the mathematical community has increasingly turned to axiomatization approaches that seek to unify these diverse operators within a coherent theoretical framework. In [ 12 ], the authors adopt such an axiomatic perspective, examining several fractional derivative definitions in a systematic way. By identifying the core properties that these operators share, they developed a general framework that facilitates a deeper understanding of their similarities and differences. For example, the axioms might include the following several items: (Linearity) The operator is linear. (Identity) The zero order derivative of a function returns the function itself. (Backward compatibility) When the order is integer, the fractional derivative gives the same result as the classical derivative. (Index law) For any α < 0 and β < 0 , we have D α D β f ( t ) = D α + β f ( t ) . (Generalized Leibniz rule) D α [ f ( t ) g ( t ) ] = āˆ‘ i = 0 āˆž α i D i f ( t ) D α āˆ’ i g ( t ) . Remark Ā  1. If α ∈ N , this is exactly the classical Leibniz rule. By the way, another generalization of this rule, such as D α [ f ( t ) g ( t ) ] = āˆ‘ i = āˆ’ i n f t y āˆž α i + β D β + i f ( t ) D α āˆ’ β āˆ’ i g ( t ) . This approach not only highlights the underlying mathematical structure common to various fractional operators but also provides a way to compare and contrast their behaviors. By formalizing fractional calculus through axiomatization, it becomes easier to extend these operators to new contexts, such as multi-dimensional systems, stochastic processes, or even quantum mechanics, where fractional derivatives are increasingly finding applications; see [ 5 , 18 , 19 , 20 ] and the references therein. The axiomatization of fractional operators has profound implications for both the theoretical development of fractional calculus and its practical applications. On the theoretical side, it provides a clearer framework for exploring the connections between different types of fractional derivatives and integrals, potentially leading to the discovery of new operators or generalizations that encompass existing ones. It also helps to clarify which definitions are most suitable for specific problems, depending on the required properties of the derivative, such as memory effects, non-locality, or boundary behavior; see [ 21 , 22 ] and the references therein. From a practical standpoint, this unifying approach simplifies the application of fractional derivatives in diverse fields, from engineering to biology. For instance, in modeling materials with memory effects, the choice of a fractional derivative is crucial, and an axiomatic perspective can help guide that choice. Similarly, in complex systems where fractal geometry or anomalous diffusion is present, such as in finance, geology, or network theory, an axiomatized framework for fractional derivatives ensures that the correct operator is used to capture the necessary dynamics. In this paper, we mainly explore the construction of exact solutions or approximate solutions to the initial value problems of time fractional partial differential equations. This is a very direct problem in the study of the theory of differential equations. However, it is not easy and even rather difficult. In addition, solving fractional systems is usually more difficult than classical systems, for their operator is defined by an integral (a non-local operator with a singular kernel). Luckily, there are some different effective methods that have been developed to construct the exact solutions or the approximate solutions of fractional systems with initial conditions. Such as the iterative method [ 23 , 24 , 25 , 26 , 27 ], the homotopy analysis method [ 28 , 29 , 30 ], the Cole-Hopf transformation method [ 31 ], the residual power series method [ 32 , 33 , 34 , 35 , 36 ], the Laplace transform method [ 37 ], and so on. Outline of Paper This manuscript presents a systematic investigation of iterative techniques for solving initial value problems in time-fractional partial differential equations. The paper is organized as follows: In Section 2 , we establish the fundamental mathematical framework by introducing essential concepts from fractional calculus, including the various definitions of fractional derivatives and their key properties. This section provides the necessary theoretical foundation by presenting important lemmas regarding fractional operators, particularly focusing on their non-local characteristics and memory effects. We carefully define all mathematical notations to ensure clarity throughout subsequent developments, paying special attention to the function spaces and operator norms relevant to our analysis of fractional differential equations. Section 3 develops an iterative approach for solving functional equations with initial conditions, with particular emphasis on applications to fractional-order systems. We present a comprehensive convergence analysis of the proposed iterative scheme, establishing sufficient conditions under which the method guarantees convergence to the true solution. The section includes detailed proofs of convergence theorems and examines the rate of convergence under various regularity assumptions. Furthermore, we discuss the computational advantages of this approach compared to traditional methods, highlighting its efficiency in handling the non-local nature of fractional operators. Section 4 demonstrates the practical implementation of our iterative technique through an in-depth study of the one-dimensional time-fractional Tricomi equation. We derive both exact solutions (where obtainable) and high-precision approximate solutions, comparing our results with known analytical solutions to validate the method’s accuracy. The section includes detailed case studies with different fractional orders and initial conditions, accompanied by rigorous error analysis and convergence rate verification. Special attention is given to the physical interpretation of solutions and their relationship to the fractional order parameter. The analysis is extended in Section 5 to multidimensional time-fractional Tricomi equations, where we address the additional complexities introduced by higher spatial dimensions. We develop dimensionally-adaptive versions of our iterative scheme and examine their performance for various initial boundary value problems. The section presents numerical results for two- and three-dimensional cases, including visualization of solution surfaces and comparison with alternative numerical methods. We particularly focus on the method’s scalability and its ability to handle different types of initial conditions, from smooth analytic functions to discontinuous data. Finally, Section 6 synthesizes our findings into a comprehensive conclusion, summarizing the key advantages and limitations of the proposed iterative technique. We discuss broader implications for the field of fractional calculus and suggest several directions for future research, including potential applications to nonlinear fractional systems and coupled fractional-order equations. The section concludes with practical recommendations for implementing the method and its possible extensions to other classes of fractional differential equations. 2. Fractional Calculus This subsection introduces the fundamental theoretical framework necessary for our investigation of fractional systems. We begin with the essential special functions that form the backbone of fractional calculus. The Gamma function, which extends the factorial concept to continuous values, plays an indispensable role in defining fractional operators due to its appearance in their integral kernels. Equally crucial is the Mittag–Leffler function, a sophisticated generalization of the exponential function that naturally describes processes with memory effects and power-law relaxation. The mathematical landscape of fractional calculus features a diverse collection of definitions for fractional integrals and derivatives, each developed to serve specific theoretical or applied purposes. Among these various formulations, several have emerged as particularly significant. The Riemann–Liouville fractional integral represents the most direct extension of classical integration theory to non-integer orders, while the Caputo fractional derivative has become preferred for physical modeling due to its more intuitive treatment of initial conditions. Other important variants include the Riesz potential, which provides symmetric fractional operations valuable for certain boundary value problems, and the Weyl derivative, especially useful for analyzing periodic phenomena. The comparative advantages and specific applications of these different definitions are thoroughly examined in the comprehensive references [ 7 , 8 , 11 , 12 , 13 , 14 ]. Within this spectrum of fractional operators, the Riemann–Liouville integral and Caputo derivative have gained particular prominence in both theoretical and applied contexts. The Riemann–Liouville approach provides the most natural mathematical foundation for fractional integration, while the Caputo formulation offers distinct practical advantages for solving initial value problems that model physical systems. This preference stems from the Caputo derivative’s compatibility with conventional initial conditions specified in terms of standard derivatives, making it more suitable for applications in engineering and physics. Our subsequent analysis will build upon several key lemmas concerning the operational properties of these fractional operators and the characteristic behaviors of Mittag–Leffler functions. These fundamental results enable rigorous investigation of existence, uniqueness, and qualitative properties for solutions to fractional differential equations. The interplay between these various concepts creates a robust framework for examining the nonlinear fractional systems that constitute the primary focus of this work, while also providing the necessary tools for developing effective solution methods. 2.1. The Gamma Function and the Mittag–Leffler Function This subsection introduces the definitions for the Gamma function and the Mittag–Leffler function, which will be used in the definitions of the fractional integral and derivative. Definition Ā  1 ([ 38 ]) . The Gamma function is defined by Ī“ ( x ) = ∫ 0 āˆž t x āˆ’ 1 e āˆ’ t d t Re ( x ) > 0 , where Re ( x ) denotes the real part of the complex number x. Proposition Ā  1 ([ 38 ]) . Some main properties of the Gamma function: Ī“ ( x + 1 ) = x Ī“ ( x ) for all x with Re x > 0 . Ī“ ( 1 ) = 1 . Ī“ ( 1 2 ) = Ļ€ . Definition Ā  2 ([ 39 ]) . The two-parameters Mittag–Leffler function is defined by letting E α , β ( x ) = āˆ‘ k = 0 āˆž x k Ī“ ( α k + β ) α , β > 0 . And the one-parameter Mittag–Leffler function is defined by letting E α ( x ) = āˆ‘ k = 0 āˆž x k Ī“ ( α k + 1 ) α > 0 . Remark Ā  2. The two-parameter Mittag–Leffler function is a generalization of the one-parameter Mittag–Leffler function, that is, E α , 1 ( x ) = E α ( x ) . The one-parameter Mittag–Leffler function is a generalization of the exponential function, or rather, E 1 ( x ) = e x . 2.2. The Riemann–Liouville Fractional Integral and the Caputo Derivative This subsection recalls the definition of the Riemann–Liouville fractional integral and Caputon derivative, see [ 4 , 5 , 6 , 13 , 40 ] and the references therein. Definition Ā  3. Let Īø : [ 0 , T ] → R is a continuous function. Then, the fractional-order integral of Īø in the Riemann–Liouville sense is defined by J α Īø ( t ) = 1 Ī“ ( α ) ∫ 0 t ( t āˆ’ s ) α āˆ’ 1 Īø ( s ) d s , α > 0 , t ≄ 0 , (1) where Ī“ ( α ) = ∫ 0 āˆž e āˆ’ u u α āˆ’ 1 d u . Definition Ā  4. Let Īø : [ 0 , T ] → R is a continuous function. Then, the Caputo fractional-order derivative of Īø is given by D α c Īø ( t ) = 1 Ī“ ( n āˆ’ α ) ∫ 0 t ( t āˆ’ s ) n āˆ’ α āˆ’ 1 Īø ( m ) ( s ) d s , where m āˆ’ 1 < α ≤ m ; m ∈ N . Lemma Ā  1 ([ 27 ]) . For any m āˆ’ 1 < β ≤ m ∈ N āˆ— , one has I β D β u ( x , t ) = u ( x , t ) āˆ’ āˆ‘ k = 0 m āˆ’ 1 āˆ‚ k u ( x , 0 ) t k k ! . Lemma Ā  2. For any β , γ , one has I β t γ = Ī“ ( γ + 1 ) Ī“ ( γ + β + 1 ) t γ + β . 3. An Iterative Technique for Solving a Class of Nonlinear FunctionalĀ Equations In this section, we present a fundamental lemma that plays a pivotal role in our subsequent mathematical analysis. It provides the essential theoretical foundation for establishing the key propositions and theorems developed throughout this paper. It has been chosen for its particular relevance to our investigation of fractional differential systems, offering crucial insights into the behavior of solutions and the properties of fractional operators. Its proper application will be carefully demonstrated in subsequent sections where we develop our main results. Lemma Ā  3. Consider a class of abstract functional equations of the following form u ( x , t ) = f ( x , t ) + L u ( x , t ) + N u ( x , t ) , (2) where ( t , x ) ∈ D = { ( t , x ) : t ∈ R , x ∈ R n , n ∈ N } , u ( x , t ) is the unknown function, f is a given known function, L is a linear operator and N is a nonlinear operator with respect to u from a Banach space Dom ( L ) ∩ Dom ( N ) to itself (Dom ( L ) ∩ Dom ( N ) stands for the domain of L and N respectively). For abstract function Equation ( 2 ), consider the series solutions of the following form u ( t , x ) = āˆ‘ n = 0 āˆž u n ( t , x ) = u 0 + u 1 + u 2 + ⋯ , where the iterative sequence { u n } can be obtained by the following recurrence equations u 0 = f , (3) u 1 = L ( u 0 ) + N ( u 0 ) , (4) u m + 1 = L ( u m ) + N āˆ‘ n = 0 m u n āˆ’ N āˆ‘ n = 0 m āˆ’ 1 u n , m = 1 , 2 ⋯ . (5) Furthermore, if the operators L and N are contractive, then the series solution āˆ‘ n = 0 āˆž u n ( t , x ) converges absolutely and uniformly. Proof. Obviously, the nonlinear operator N can be decomposed as N ( u ( t , x ) ) = N āˆ‘ n = 0 āˆž u n ( t , x ) = N ( u 0 ) + āˆ‘ n = 1 āˆž N āˆ‘ j = 0 n u j āˆ’ N āˆ‘ j = 0 n āˆ’ 1 u j . Similarly, the linear operator L can also be decomposed as L ( u ( t , x ) ) = L āˆ‘ n = 0 āˆž u n ( t , x ) = L ( u 0 ) + āˆ‘ n = 1 āˆž L āˆ‘ j = 0 n u j āˆ’ L āˆ‘ j = 0 n āˆ’ 1 u j . Denote by M u ( t , x ) = L u ( t , x ) + N u ( t , x ) . then one has M ( u ( t , x ) ) = M ( u 0 ) + āˆ‘ n = 1 āˆž M āˆ‘ j = 0 n u j āˆ’ M āˆ‘ j = 0 n āˆ’ 1 u j . Define the following recurrence equations u 0 = f , u 1 = M ( u 0 ) , u m + 1 = M āˆ‘ n = 0 m u n āˆ’ M āˆ‘ n = 0 m āˆ’ 1 u n , m = 1 , 2 ⋯ . Since operators L and N are contractive, then M is also contractive, i.e., there exists a constant 0 < K < 1 , such that ∄ M ( v i ) āˆ’ M ( v j ) ∄ ≤ K ∄ v i āˆ’ v j ∄ , āˆ€ v i , v j ∈ Dom ( L ) ∩ D o m ( N ) , where ∄ Ā· ∄ denotes the norm on Banach space Dom ( L ) ∩ Dom ( N ) . Furthermore, for u m + 1 , one can obtain ∄ u m + 1 ∄ = M āˆ‘ n = 0 m u n āˆ’ M āˆ‘ n = 0 m āˆ’ 1 u n ≤ K ∄ u m ∄ ≤ K m + 1 ∄ u 0 ∄ . Since 0 < K < 1 , the series āˆ‘ m = 1 āˆž K m + 1 ‖ u 0 ‖ converges absolutely as well as uniformly. According to the Weierstrass M-test, one can obtain that the series āˆ‘ i = 0 āˆž u i converges absolutely as well as uniformly. ā–” Some more specific details about this type of iterative method could be found in [ 23 , 24 , 25 , 26 , 27 ] and the references therein. 4. One-Dimensional Time Fractional Tricomi Equation with VariableĀ Coefficients Consider the following one-dimensional time-fractional Tricomi equation of order α , as studied in [ 41 ]: D t α u + ξ ( t ) āˆ‚ x 2 u = f ( x , t ) , (6) subject to the initial conditions u ( x , 0 ) = ω 1 ( x ) , D t u ( x , 0 ) = ω 2 ( x ) , (7) where 1 < α ≤ 2 represents the fractional order of the time derivative, t ∈ [ 0 , T ] denotes the temporal variable, and x ∈ R is the spatial variable. The special case when α = 2 reduces ( 6 ) to the classical Tricomi equation, first introduced by Francesco Tricomi in his seminal 1923 work [ 42 ]. This equation occupies a particularly important position in the theory of partial differential equations due to its remarkable property of exhibiting mixed-type characteristics. Specifically, the nature of the equation undergoes a fundamental change depending on the sign of the coefficient function ξ ( t ) : when ξ ( t ) > 0 , the equation demonstrates elliptic behavior, while for ξ ( t ) < 0 , it becomes hyperbolic. This dual character makes the Tricomi equation a paradigmatic example of mixed-type partial differential equations. The physical significance of this equation was further elucidated by Frankl [ 43 ], who established its crucial connection to transonic gas dynamics. More precisely, the Tricomi equation provides a mathematical model describing the transition between subsonic and supersonic flow regimes in gas dynamics. The elliptic region ( ξ ( t ) > 0 ) corresponds to subsonic flow, where disturbances propagate in all directions, while the hyperbolic region ( ξ ( t ) < 0 ) represents supersonic flow, characterized by directional wave propagation and the formation of shock waves. This interpretation underscores the equation’s importance in aerodynamics and fluid mechanics. In the following analysis, we will develop an iterative solution technique to solve the initial value problem ( 6 ) and ( 7 ) for the time-fractional generalization of this important equation. Our approach will carefully account for both the fractional temporal derivative and the mixed-type spatial operator, while maintaining the essential physical interpretation of the solutions. By integrating both sides of Equation ( 6 ) with respect to the temporal variable from 0 to t , we derive the equivalent integral formulation: u ( x , t ) = u ( x , 0 ) + D t u ( x , 0 ) t āˆ’ I t α ξ ( t ) āˆ‚ x 2 u + I t α f ( x , t ) . (8) Incorporating the initial conditions ( 7 ) into the integral Equation ( 8 ) yields: u ( x , t ) = ω 1 ( x ) + ω 2 ( x ) t + I t α f ( x , t ) āˆ’ I t α ξ ( t ) āˆ‚ x 2 u . From the perspective of iterative methods, we can interpret this equation as follows: The terms ω 1 ( x ) + ω 2 ( x ) + I t α f ( x , t ) constitute the initial iterative data, incorporating both the given initial conditions and the source term f ( x , t ) . The remaining term āˆ’ I t α ξ ( t ) āˆ‚ x 2 u acts as a nonlinear feedback operator that propagates the solution through successive iterations. We propose a series solution of the form: u ( x , t ) = āˆ‘ n = 0 āˆž u n ( x , t ) = u 0 + u 1 + u 2 + ⋯ , where the iterative sequence is determined by the following recursive scheme: u 0 = ω 1 ( x ) + ω 2 ( x ) t + I t α f ( x , t ) , u 1 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 , u 2 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 , ā‹® u m + 1 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u m . To demonstrate the effectiveness of this iterative approach, we now examine several specific cases with carefully chosen initial conditions, source terms, and coefficient functions. Proposition Ā  2 (Trivial Solution Case) . For the simplest case where ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , and f ( x , t ) ≔ 0 , we observe that for any well-behaved coefficient function ξ ( t ) , the solution identically vanishes: u ( x , t ) ≔ 0 . This trivial solution serves as an important consistency check for our method. Proposition Ā  3 (Exponential Source Term with Constant Coefficient) . Consider the case with zero initial conditions ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , an exponential source term f ( x , t ) ≔ e x , and constant coefficient ξ ( t ) ≔ āˆ’ 1 . The iterative scheme yields: u n = e x t ( n + 1 ) α Ī“ ( ( n + 1 ) α + 1 ) āˆ€ n ∈ N . The complete solution then takes the elegant form: u ( x , t ) = āˆ‘ n = 0 āˆž u n = e x āˆ‘ n = 0 āˆž t ( n + 1 ) α Ī“ ( ( n + 1 ) α + 1 ) = e x E α ( t α ) āˆ’ 1 , where E α ( Ā· ) denotes the Mittag–Leffler function, which naturally emerges in solutions to fractional differential equations. Proposition Ā  4 (Modified Exponential Source Term) . For the case with ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x t α Ī“ ( āˆ’ α + 1 ) , and ξ ( t ) ≔ āˆ’ 1 , we obtain: u n = e x t n α Ī“ ( n α + 1 ) āˆ€ n ∈ N . The solution series converges to: u ( x , t ) = e x E α ( t α ) , demonstrating how different source terms can lead to distinct Mittag–Leffler function representations. Proposition Ā  5 (Non-constant Coefficient Case) . Consider the problem with ω 1 ( x ) = e x , ω 2 ( x ) = 0 , f ( x , t ) = 0 , ξ ( t ) = t α . Then the iterative scheme generates the sequence u 0 = e x , u 1 = āˆ’ e x Ī“ ( α + 1 ) Ī“ ( 2 α + 1 ) t 2 α , u 2 = e x Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) t 4 α . By the induction of mathematics, we establish the general term: u n = ( āˆ’ 1 ) n e x āˆ k = 1 n Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) āˆ k = 1 n Ī“ ( 2 k α + 1 ) t 2 n α , n ≄ 0 . The complete solution is expressed as the infinite series: u ( x , t ) = e x āˆ‘ n = 0 āˆž ( āˆ’ 1 ) n āˆ k = 1 n Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) āˆ k = 1 n Ī“ ( 2 k α + 1 ) t 2 n α , revealing the intricate relationship between the solution’s temporal behavior and the fractional order α. Proof. Step 1: Iterative structure. Under the given assumptions, the iteration takes the form u n + 1 ( x , t ) = āˆ’ I t α t α u n ( x , t ) , where I t α denotes the fractional integral operator. Because the initial data is separable and u 0 = e x , we observe that each iterate preserves the spatial factor e x . Hence we may write u n ( x , t ) = e x a n ( t ) , reducing the problem to determining the scalar sequence a n ( t ) . Step 2: First iterates. For n = 0 , a 0 ( t ) = 1 . Using the fractional integral formula I t α t β = Ī“ ( β + 1 ) Ī“ ( β + α + 1 ) t β + α , we obtain a 1 ( t ) = āˆ’ Ī“ ( α + 1 ) Ī“ ( 2 α + 1 ) t 2 α . Repeating the same computation yields a 2 ( t ) = Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) t 4 α . Step 3: Inductive hypothesis. Assume that for some n ≄ 1 , u n = ( āˆ’ 1 ) n e x āˆ k = 1 n Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) āˆ k = 1 n Ī“ ( 2 k α + 1 ) t 2 n α . Step 4: Inductive step. Applying the iteration formula, u n + 1 = āˆ’ I t α t α u n . Substituting the inductive expression gives t α u n = ( āˆ’ 1 ) n e x āˆ k = 1 n Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) āˆ k = 1 n Ī“ ( 2 k α + 1 ) t ( 2 n + 1 ) α . Applying the fractional integral formula once more, I t α t ( 2 n + 1 ) α = Ī“ ( ( 2 n + 1 ) α + 1 ) Ī“ ( ( 2 n + 2 ) α + 1 ) t ( 2 n + 2 ) α . Therefore, u n + 1 = ( āˆ’ 1 ) n + 1 e x āˆ k = 1 n + 1 Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) āˆ k = 1 n + 1 Ī“ ( 2 k α + 1 ) t 2 ( n + 1 ) α . This completes the induction. Step 5: Series representation. Summing the iterates yields the solution as the convergent power-type series u ( x , t ) = āˆ‘ n = 0 āˆž u n ( x , t ) , which gives the stated expansion. ā–” 5. Higher-Dimensional Time Fractional Tricomi Equation with VariableĀ Coefficients Consider the following higher-dimensional time fractional Tricomi equation D t α u + ξ ( t ) Ī” x u = f ( x , t ) (9) with initial conditions u ( x , 0 ) = ω 1 ( x ) , D t u ( x , 0 ) = ω 2 ( x ) . (10) where 1 < α ≤ 2 , t ∈ [ 0 , T ] , x ∈ R n , and Ī” x denotes the standard Laplace operator and it is defined by Ī” x u ( x , t ) = āˆ‘ i = 1 n āˆ‚ x i 2 u ( x , t ) . Similarly, one can obtain the equivalent integral Equations ( 9 ) and ( 10 ) as follows u ( x , t ) = ω 1 ( x ) + ω 2 ( x ) t + I t α f ( x , t ) āˆ’ I t α ξ ( t ) Ī” x u ( x , t ) . Suppose that the unknown function u has the series form u ( x , t ) = āˆ‘ n = 0 āˆž u n ( x , t ) = u 0 + u 1 + u 2 + ⋯ , where the sequence { u n } can be obtained by the following iterative process u 0 = ω 1 ( x ) + ω 2 ( x ) t + I t α f ( x , t ) u 1 = āˆ’ I t α ξ ( t ) Ī” x u 0 u 2 = āˆ’ I t α ξ ( t ) Ī” x u 1 ā‹® u m + 1 = āˆ’ I t α ξ ( t ) Ī” x u m . Proposition Ā  6. If ω 1 ( x ) ≔ 0 , ω 2 ≔ 0 , f ( x , t ) ≔ 0 , then for any reasonable ξ ( t ) , u ( x , t ) ≔ 0 . Proposition Ā  7. Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 + x 2 + ⋯ + x n and ξ ( t ) ≔ āˆ’ 1 . Iterative construction. Since the initial data vanish, the first term of the iteration is obtained directly from the source term: u 0 = e x 1 + ⋯ + x n t α Ī“ ( α + 1 ) . Applying the iterative scheme once more and using I t α t β = Ī“ ( β + 1 ) Ī“ ( β + α + 1 ) t β + α , we obtain u 1 = e x 1 + ⋯ + x n t 2 α Ī“ ( 2 α + 1 ) . Proceeding inductively, one easily verifies that u n = e x 1 + ⋯ + x n t ( n + 1 ) α Ī“ ( ( n + 1 ) α + 1 ) , āˆ€ n ∈ N . Series solution. Summing all iterates, u ( x , t ) = āˆ‘ n = 0 āˆž u n = e x 1 + ⋯ + x n āˆ‘ n = 0 āˆž t ( n + 1 ) α Ī“ ( ( n + 1 ) α + 1 ) = e x 1 + ⋯ + x n āˆ‘ n = 1 āˆž t n α Ī“ ( n α + 1 ) . Rewriting the series, u ( x , t ) = e x 1 + ⋯ + x n āˆ‘ n = 0 āˆž t n α Ī“ ( n α + 1 ) āˆ’ 1 . Hence, u ( x , t ) = e x 1 + ⋯ + x n E α ( t α ) āˆ’ 1 . Proposition Ā  8. Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 + x 2 + ⋯ + x n t α Ī“ ( 1 āˆ’ α ) , ξ ( t ) ≔ āˆ’ 1 . Iterative construction. Since the initial data vanish and the source term already contains t α , the first iterate reduces to u 0 = e x 1 + ⋯ + x n . Applying the iterative scheme once, u 1 = e x 1 + ⋯ + x n t α Ī“ ( α + 1 ) . Continuing in the same manner, one verifies inductively that u n = e x 1 + ⋯ + x n t n α Ī“ ( n α + 1 ) , āˆ€ n ∈ N . Series solution. Summing the iterates gives u ( x , t ) = āˆ‘ n = 0 āˆž u n = e x 1 + ⋯ + x n āˆ‘ n = 0 āˆž t n α Ī“ ( n α + 1 ) . Recognizing the Mittag–Leffler series, E α ( z ) = āˆ‘ n = 0 āˆž z n Ī“ ( n α + 1 ) , we conclude that u ( x , t ) = e x 1 + ⋯ + x n E α ( t α ) . Proposition Ā  9. Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ e x , f ( x , t ) ≔ 0 , and ξ ( t ) ≔ āˆ’ t α Ī“ ( α + 1 ) . Iterative construction. Since the source term vanishes and the initial velocity is e x , the first term of the iteration is u 0 = e x . Proceeding with the iterative scheme, one obtains successively terms of increasing order in t α . By induction it follows that u n = e x t n α Ī“ ( n α + 1 ) , āˆ€ n ∈ N . Series solution. Summing the iterates, u ( x , t ) = āˆ‘ n = 0 āˆž u n = e x āˆ‘ n = 0 āˆž t n α Ī“ ( n α + 1 ) . Recognizing the Mittag–Leffler expansion, we conclude that u ( x , t ) = e x E α ( t α ) . Proposition Ā  10. Assume ω 1 ( x ) ≔ e x 1 + ⋯ + x n , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ 0 , ξ ( t ) ≔ t α . First iterates. The initial term is u 0 = e x 1 + ⋯ + x n . Since āˆ‚ x 2 e x 1 + ⋯ + x n = e x 1 + ⋯ + x n , the first iteration gives u 1 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 = āˆ’ e x 1 + ⋯ + x n I t α t α = āˆ’ e x 1 + ⋯ + x n Ī“ ( α + 1 ) Ī“ ( 2 α + 1 ) t 2 α . Applying the scheme once more, u 2 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 = e x 1 + ⋯ + x n Ī“ ( α + 1 ) Ī“ ( 2 α + 1 ) I t α t 3 α = e x 1 + ⋯ + x n Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) t 4 α . Inductive step. Assume that for some k ≄ 1 , u k = ( āˆ’ 1 ) k + 1 e x 1 + ⋯ + x n Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) ⋯ Ī“ ( 2 k α + 1 ) t 2 k α . Then u k + 1 = āˆ’ I t α ξ ( t ) āˆ‚ x 2 u k = ( āˆ’ 1 ) k + 2 e x 1 + ⋯ + x n Ī“ ( α + 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α + 1 ) Ī“ ( 2 α + 1 ) ⋯ Ī“ ( 2 k α + 1 ) I t α t ( 2 k + 1 ) α . Using I t α t β = Ī“ ( β + 1 ) Ī“ ( β + α + 1 ) t β + α , we obtain u k + 1 = ( āˆ’ 1 ) k + 2 e x 1 + ⋯ + x n Ī“ ( α + 1 ) ⋯ Ī“ ( ( 2 k + 1 ) α + 1 ) Ī“ ( 2 α + 1 ) ⋯ Ī“ ( ( 2 k + 2 ) α + 1 ) t ( 2 k + 2 ) α . Hence, by induction, for all n ∈ N , u n = ( āˆ’ 1 ) n + 1 e x 1 + ⋯ + x n Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) ⋯ Ī“ ( 2 n α + 1 ) t 2 n α . Series representation. Summing the iterates, u ( x , t ) = āˆ‘ n = 0 āˆž u n , that is, u ( x , t ) = e x 1 + ⋯ + x n āˆ‘ n = 0 āˆž ( āˆ’ 1 ) n + 1 Ī“ ( α + 1 ) Ī“ ( 3 α + 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α + 1 ) Ī“ ( 2 α + 1 ) Ī“ ( 4 α + 1 ) ⋯ Ī“ ( 2 n α + 1 ) t 2 n α . Author Contributions Methodology, W.L. and J.H.; Validation, J.H.; Writing—original draft, J.H.; Writing—review and editing, W.L.; Visualization, W.L. and J.H. All authors have read and agreed to the published version of the manuscript. Funding This research received no external funding. Data Availability Statement No new data were created or analyzed in this study. Conflicts of Interest The authors declare no conflicts of interest. Appendix A Appendix A.1. The Conformable Fractional Derivative This subsection recalls the definition of the conformable fractional derivative, which looks like the limit of the difference quotient (the classical Newtonian derivative). The α -th conformable fractional derivative of a function f is defined by letting ( T α f ) ( t ) = lim ϵ → 0 f ( t + ϵ t 1 āˆ’ α ) āˆ’ f ( t ) ϵ . Here 0 ≤ α ≤ 1 . It is clear that the case of α = 1 coincides with the classical Newtonian derivative. Furthermore, one can verify that The zero order derivative of a function does not return the function itself. In fact, we have ( T 0 f ) ( t ) = lim ϵ → 0 f ( t + ϵ t ) āˆ’ f ( t ) ϵ = lim ϵ → 0 f ( ( 1 + ϵ ) t ) āˆ’ f ( t ) 1 + ϵ Ā· lim ϵ → 0 1 + ϵ ϵ = ( D f ) ( t ) Ā· āˆž ≠ f ( t ) , where D denotes the quantum derivative. The index law does not hold, that is, ( T α T β f ) ( t ) ≠ ( T α + β f ) ( t ) . In fact, we have LHS = lim ϵ → 0 1 ϵ 2 [ f ( t + ϵ t 1 āˆ’ α + ϵ ( t + ϵ t 1 āˆ’ α ) 1 āˆ’ β ) āˆ’ f ( t + ϵ t 1 āˆ’ α ) āˆ’ f ( t + ϵ t 1 āˆ’ α ) āˆ’ f ( t + ϵ t 1 āˆ’ β ) + f ( t ) ] and RHS = lim ϵ → 0 f ( t + ϵ t 1 āˆ’ ( α + β ) ) āˆ’ f ( t ) ϵ . It is not difficult to see that ( T α T β f ) ( t ) ≠ ( T α + β f ) ( t ) without any restriction. In this perspective of the criteria, the conformable fractional derivative is not a fractional derivative; see [ 13 ]. By the way, regarding the definition of the higher order conformable fractional derivative, please refer to [ 44 ] (Definition 2.1) and the references therein. Appendix A.2. The Grünwald–Letnikov Derivative This subsection recalls the definition of the Grünwald–Letnikov derivative. It is a type of fractional derivative that extends the concept of differentiation to non-integer orders. It is defined using a limit process similar to the traditional definition of the derivative but incorporates fractional differences. Specifically, it is expressed as a limit of a sum involving the values of a function at discrete points, weighted by binomial coefficients. This derivative is useful in various applications, particularly in fields like physics and engineering, where systems exhibit non-local or memory effects. Its flexibility allows for modeling behaviors that standard derivatives cannot capture. More precisely, it is defined by letting ( GL D Īø α f ) ( z ) = e āˆ’ i α Īø lim | h | → 0 āˆ‘ k = 0 āˆž ( āˆ’ 1 ) k C α k f ( z āˆ’ k h ) | h | α , where z ∈ C and h = | h | e i Īø , āˆ’ Ļ€ < Īø ≤ Ļ€ . When the complex independent variable z becomes the real number t , we have the following two types of fractional derivative: Forward derivative ( h = | h | ): ( D F α GL f ) ( t ) = lim | h | → 0 āˆ‘ k = 0 āˆž ( āˆ’ 1 ) k C α k f ( t āˆ’ k h ) | h | α Backward derivative ( h = āˆ’ | h | ): ( D B α GL f ) ( t ) = lim | h | → 0 āˆ‘ k = 0 āˆž ( āˆ’ 1 ) k C α k f ( t + k | h | ) | h | α We introduce only the forward derivative. The operator is linear, that is, D F α GL ( Ī» f + ν g ) ( z ) = Ī» ( GL D F α f ) ( z ) + μ ( GL D F α g ) ( z ) , āˆ€ Ī» , ν ∈ R . Additivity and commutativity: D F α GL D F β GL ( f ) ( t ) = D F β GL D F α GL ( f ) ( t ) = ( D F α + β GL ( f ) ) ( t ) . Neutral element. Letting β = āˆ’ α yields that D F α GL D F āˆ’ α GL ( f ) ( t ) = f ( t ) and D F āˆ’ α GL D F α GL ( f ) ( t ) = f ( t ) . Backward compatibility. The generalized Leibniz rule. That is, for the derivative of the product φ ψ , we have D F α G L [ φ ψ ] ( t ) = āˆ‘ n = 0 āˆž C n α D n φ ( t ) D F α āˆ’ n G L ψ ( t ) . References Nguyen, A.T.; Caraballo, T.; Tuan, N.H. On the initial value problem for a class of nonlinear biharmonic equation with time-fractional derivative. Proc. R. Soc. Edinb. Sect. A Math. 2022 , 152 , 989–1031. [ Google Scholar ] [ CrossRef ] 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 ] Luchko, Y. Initial-boundary-value problems for the one-dimensional time-fractional diffusion equation. Fract. Calc. Appl. Anal. 2012 , 15 , 141–160. [ Google Scholar ] [ CrossRef ] Karthikeyan, K.; Karthikeyan, P.; Patel, D.K.; Al-Askar, F.M.; Cesarano, C. Solving fractional integro-differential equations with delay and relaxation impulsive terms by fixed point techniques. Bound. Value Probl. 2024 , 2024 , 146. [ Google Scholar ] [ CrossRef ] Abbas, M.; Ciobanescu, C.; Asghar, M.W.; Omame, A. Solution approximation of fractional boundary value problems and convergence analysis using AA-iterative scheme. AIMS Math. 2024 , 9 , 13129–13158. [ Google Scholar ] [ CrossRef ] Alsaedi, A.; Ahmad, B.; Alharbi, N.; Alharbi, N.; Alharbi, N. Fixed point and stability analysis of a tripled system of nonlinear fractional differential equations with n-nonlinear terms. Fractal Fract. 2024 , 8 , 697. [ Google Scholar ] [ CrossRef ] Mali, A.D.; JĆørgensen, J.B.; KĆøstner, N.; LĆønborg, J.; Dutta, H. Numerical investigation of generalized tempered-type integrodifferential equations with respect to another function. Fract. Calc. Appl. Anal. 2023 , 26 , 2580–2601. [ Google Scholar ] [ CrossRef ] Gunasekar, T.; Raghavendran, P.; Santra, S.S.; Majumder, D.; Baleanu, D.; Balasundaram, H. Application of laplace transform to solve fractional integro-differential equations. J. Math. Comput. Sci. 2024 , 33 , 225–237. [ Google Scholar ] [ CrossRef ] Ross, B. The development of fractional calculus 1695–1900. Historia Math. 1977 , 4 , 75–89. [ Google Scholar ] [ CrossRef ] El-Ajou, A.; Arqub, O.A.; Zhour, Z.A.; Momani, S. New results on fractional power series: Theories and applications. Entropy 2013 , 15 , 5305–5323. [ Google Scholar ] [ CrossRef ] Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations ; Volume 204 of North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [ Google Scholar ] Ortigueira, M.D.; Machado, J.A.T. What is a fractional derivative? J. Comput. Phys. 2015 , 293 , 4–13. [ Google Scholar ] [ CrossRef ] Katugampola, U.N. Correction to ā€œWhat is a fractional derivative?ā€ by Ortigueira and Machado [Journal of Computational Physics, Volume 293, 15 July 2015, Pages 4–13. Special issue on fractional PDEs]. J. Comput. Phys. 2016 , 321 , 1255–1257. [ Google Scholar ] [ CrossRef ] Kuroda, L.K.B.; Gomes, A.V.; Tavoni, R.; Mancera, P.F.d.; Varalta, N.; Camargo, R.d.F. Unexpected behavior of Caputo fractional derivative. Comput. Appl. Math. 2017 , 36 , 1173–1183. [ Google Scholar ] [ CrossRef ] Murio, D.A. Stable numerical evaluation of grünwald–letnikov fractional derivatives applied to a fractional ihcp. Inverse Probl. Sci. Eng. 2009 , 17 , 229–243. [ Google Scholar ] [ CrossRef ] Chikrii, A.; Matychyn, I. Riemann–liouville, caputo, and sequential fractional derivatives in differential games. In Advances in Dynamic Games: Theory, Applications, and Numerical Methods for Differential and Stochastic Games ; Springer: Berlin/Heidelberg, Germany, 2011; pp. 61–81. [ Google Scholar ] Agrawal, O.P. Fractional variational calculus in terms of riesz fractional derivatives. J. Phys. Math. Theor. 2007 , 40 , 6287. [ Google Scholar ] [ CrossRef ] Laskin, N. Principles of fractional quantum mechanics. arXiv 2010 , arXiv:1009.5533. [ Google Scholar ] [ CrossRef ] Laskin, N. Fractional quantum mechanics. Phys. Rev. E 2000 , 62 , 3135. [ Google Scholar ] [ CrossRef ] Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E 1997 , 55 , 3581. [ Google Scholar ] [ CrossRef ] Falcini, F.; Garra, R.; Voller, V. Modeling anomalous heat diffusion: Comparing fractional derivative and non-linear diffusivity treatments. Int. J. Therm. Sci. 2019 , 137 , 584–588. [ Google Scholar ] [ CrossRef ] Garra, R.; Falcini, F.; Voller, V.R.; Pagnini, G. A generalized stefan model accounting for system memory and non-locality. Int. Commun. Heat Mass Transf. 2020 , 114 , 104584. [ Google Scholar ] [ CrossRef ] Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. 2006 , 316 , 753–763. [ Google Scholar ] [ CrossRef ] Noor, M.A.; Noor, K.I.; Mohyud-Din, S.T.; Shabbir, A. An iterative method with cubic convergence for nonlinear equations. Appl. Math. Comput. 2006 , 183 , 1249–1255. [ Google Scholar ] [ CrossRef ] Bhalekar, S.; Daftardar-Gejji, V. New iterative method: Application to partial differential equations. Appl. Math. Comput. 2008 , 203 , 778–783. [ Google Scholar ] [ CrossRef ] Daftardar-Gejji, V.; Bhalekar, S. Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method. Comput. Math. Appl. 2010 , 59 , 1801–1809. [ Google Scholar ] [ CrossRef ] Dhaigude, C.D.; Nikam, V.R. Solution of fractional partial differential equations using iterative method. Fract. Calc. Appl. Anal. 2012 , 15 , 684–699. [ Google Scholar ] [ CrossRef ] Liao, S. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. 2004 , 147 , 499–513. [ Google Scholar ] [ CrossRef ] Das, S.; Vishal, K. Homotopy analysis method for fractional Swift-Hohenberg equation. In Advances in the Homotopy Analysis Method ; World Scientific Publishing: Hackensack, NJ, USA, 2014; pp. 291–308. [ Google Scholar ] [ CrossRef ] Demir, A.; Bayrak, M.A.; Ozbilge, E. A new approach for the approximate analytical solution of space-time fractional differential equations by the homotopy analysis method. Adv. Math. Phys. 2019 , 2019 , 5602565. [ Google Scholar ] [ CrossRef ] Li, L.; Li, D. Exact solutions and numerical study of time fractional Burgers’ equations. Appl. Math. Lett. 2020 , 100 , 106011. [ Google Scholar ] [ CrossRef ] Arqub, O.A. Series solution of fuzzy differential equations under strongly generalized differentiability. J. Adv. Res. Appl. Math. 2013 , 5 , 31–52. [ Google Scholar ] [ CrossRef ] El-Ajou, A.; Arqub, O.A.; Momani, S.; Baleanu, D.; Alsaedi, A. A novel expansion iterative method for solving linear partial differential equations of fractional order. Appl. Math. Comput. 2015 , 257 , 119–133. [ Google Scholar ] [ CrossRef ] Dubey, V.P.; Kumar, R.; Kumar, D. A reliable treatment of residual power series method for time-fractional Black–Scholes European option pricing equations. Physical A 2019 , 533 , 122040. [ Google Scholar ] [ CrossRef ] Syam, M.I. On approximate solutions of Fokker Planck equation by the modified residual power series method. J. Fract. Calc. Appl. 2019 , 10 , 225–232. [ Google Scholar ] Az-Zo’bi, E.A.; Yıldırım, A.; AlZoubi, W.A. The residual power series method for the one-dimensional unsteady flow of a van der Waals gas. Physical A 2019 , 517 , 188–196. [ Google Scholar ] [ CrossRef ] Mosaleheh, K.; Vakilzadeh, A. The Laplace transform method for linear ordinary differential equations of fractional order. J. Math. Ext. 2007 , 2 , 93–102. [ Google Scholar ] Artin, E. The Gamma Function ; Butler, M., Translator; Athena Series: Selected Topics in Mathematics; Holt, Rinehart and Winston: New York, NY, USA; Toronto, ON, Canada; London, UK, 1964. [ Google Scholar ] Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications ; Springer Monographs in Mathematics; Springer: Heidelberg, Germany, 2014. [ Google Scholar ] El-Ajou, A.; Arqub, O.A.; Momani, S. Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: A new iterative algorithm. J. Comput. Phys. 2015 , 293 , 81–95. [ Google Scholar ] [ CrossRef ] Arqub, O.A. Solutions of time-fractional Tricomi and Keldysh equations of Dirichlet functions types in Hilbert space. Numer. Methods Partial. Differ. Equ. 2018 , 34 , 1759–1780. [ Google Scholar ] [ CrossRef ] Zhang, X.; Liu, J.; Wen, J.; Tang, B.; He, Y. Analysis for one-dimensional time-fractional Tricomi-type equations by LDG methods. Numer. Algorithms 2013 , 63 , 143–164. [ Google Scholar ] [ CrossRef ] Frankl, F.L. On the Problems of Chaplygin for Mixed Sub- and Supersonic Flows. 1947. Available online: https://ntrs.nasa.gov/search.jsp?R=20050019646 (accessed on 10 March 2026). El-Ajou, A.; Oqielat, M.N.; Al-Zhour, Z.; Kumar, S.; Momani, S. Solitary solutions for time-fractional nonlinear dispersive PDEs in the sense of conformable fractional derivative. Chaos 2019 , 29 , 093102. [ Google Scholar ] [ CrossRef ] 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. Ā© 2026 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license . Article Metrics Citations Multiple requests from the same IP address are counted as one view.
Markdown
[Skip to Content](https://www.mdpi.com/2227-7390/14/6/984#main) [![MDPI](https://mdpi-res.com/cdn-cgi/image/w=56,h=30/https://mdpi-res.com/data/mdpi-logo-black.svg)![MDPI](https://mdpi-res.com/cdn-cgi/image/w=70,h=56/https://mdpi-res.com/data/mdpi-logo-black.svg)](https://www.mdpi.com/) - Journals - [All Journals](https://www.mdpi.com/about/journals) - [Journal Finder](https://www.mdpi.com/about/journalfinder) - [Proceedings Series](https://www.mdpi.com/about/proceedings) - [Propose a Journal](https://www.mdpi.com/about/journals/proposal) - Topics Our Topics - [All Topics](https://www.mdpi.com/topics) - [About Topics](https://www.mdpi.com/topics/about) - [Topics Awards](https://www.mdpi.com/topics/award) - [Subject Editors](https://www.mdpi.com/topics/editors) - [Propose a Topic](https://www.mdpi.com/topics/proposal) By Subjects - [Biology & Life Sciences](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wBGQsc0lDcmgIQgHtyuBzAJ4ACAGRCAMhABmGIQGUsEDADssGfAF8AukA) - [Chemistry & Materials](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wBmQsc0lDcmgYQAsMBbCBHACcAngAIAZKICypHBkEQKCUQGUsEDADss1AL4BdIA) - [Computer Science & Mathematics](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wAWQsc0lDcmgYQHsBbAByRwYATgAIAylggYAdlgwiAZCICypHAAsM3NRCz4AvgF0gA) - [Engineering](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wBWQsc0lDcmgUQDsBzCTjDACc%2B3MAF8AukA) - [Environmental & Earth Sciences](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wA2Qsc0lDcmgUQDtiIAnAe04BbDJxwUABADIJ7UrxwALCQGUsEUVgz4AvgF0gA) - [Medicine & Pharmacology](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wB2Qsc0lDcmgWQygiwgDsMAAgBkggAoALUgCcAtqSwB7cooDmATzABfALpA) - [Public Health & Healthcare](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wBGABkLHNJQ3JoAVVyIsACABIYKOABb8AZEJHlxWUgCdqAXwC6QA) - [Social Sciences, Arts & Humanities](https://www.mdpi.com/topics?facets=NobwRAlgJmBcYGcCuAjAVgUwMYBcFgBowA3AQwBskM4wBOQsc0lDcmgZQHssIKACdjwwA7LBgQE%2BAQQBOePqWFQ%2BACSQBbRRBwRxYAL4BdIA) - [Author Services](https://www.mdpi.com/authors/english) - Information General Content - [Publishing Services](https://www.mdpi.com/publishing_services) - [Open Access Policy](https://www.mdpi.com/openaccess) - [Editorial Process](https://www.mdpi.com/editorial_process) - [Publication Ethics](https://www.mdpi.com/ethics) - [Special Issues Guidelines](https://www.mdpi.com/special_issues_guidelines) - [Article Processing Charge](https://www.mdpi.com/apc) Guidelines - [For Authors](https://www.mdpi.com/authors) - [For Reviewers](https://www.mdpi.com/reviewers) - [For Editors](https://www.mdpi.com/editors) - [For Librarians](https://www.mdpi.com/librarians) Partnerships - [Societies](https://www.mdpi.com/societies) - [Conferences](https://www.mdpi.com/partnerships/conferences) - [Institutional Open Access Program](https://www.mdpi.com/ioap) - About Company - [About Us](https://www.mdpi.com/about) - [Mission](https://www.mdpi.com/about/mission) - [Impact](https://www.mdpi.com/about/impact) - [History](https://www.mdpi.com/about/history) - [Leadership](https://www.mdpi.com/about/leadership) - [Office Locations](https://www.mdpi.com/about/office-locations) - [Careers](https://careers.mdpi.com/) - [Awards](https://www.mdpi.com/awards) - [Products](https://www.mdpi.com/about/products) - [Accessibility](https://www.mdpi.com/accessibility) Media - [News](https://www.mdpi.com/news) - [Blog](https://blog.mdpi.com/) Contact - [Support](https://www.mdpi.com/support) - [Send Feedback](https://www.mdpi.com/support/feedback) Open mobile navigation 1. [Home](https://www.mdpi.com/) 2. [Journals](https://www.mdpi.com/about/journals) 3. [Mathematics](https://www.mdpi.com/journal/mathematics) 4. [Volume 14](https://www.mdpi.com/2227-7390/14) 5. [Issue 6](https://www.mdpi.com/2227-7390/14/6) 6. [10\.3390/math14060984](https://www.mdpi.com/2227-7390/14/6/984) [![](https://mdpi-res.com/cdn-cgi/image/w=32,h=32/https://mdpi-res.com/img/journals/mathematics-logo-sq.png?72e5113f8e55c761)MathematicsMathematics](https://www.mdpi.com/journal/mathematics) Get Alerted Get Alerted Open Cite Cite Open Share Share [Download PDF PDF](https://www.mdpi.com/2227-7390/14/6/984/pdf) [Abstract](https://www.mdpi.com/2227-7390/14/6/984#Abstract)[Introduction](https://www.mdpi.com/2227-7390/14/6/984#Introduction)[Fractional Calculus](https://www.mdpi.com/2227-7390/14/6/984#Fractional_Calculus)[An Iterative Technique for Solving a Class of Nonlinear Functional Equations](https://www.mdpi.com/2227-7390/14/6/984#An_Iterative_Technique_for_Solving_a_Class_of_Nonlinear_Functional%C2%A0Equations)[One-Dimensional Time Fractional Tricomi Equation with Variable Coefficients](https://www.mdpi.com/2227-7390/14/6/984#OneDimensional_Time_Fractional_Tricomi_Equation_with_Variable%C2%A0Coefficients)[Higher-Dimensional Time Fractional Tricomi Equation with Variable Coefficients](https://www.mdpi.com/2227-7390/14/6/984#HigherDimensional_Time_Fractional_Tricomi_Equation_with_Variable%C2%A0Coefficients)[Concluding Remarks](https://www.mdpi.com/2227-7390/14/6/984#Concluding_Remarks)[Author Contributions](https://www.mdpi.com/2227-7390/14/6/984#Author_Contributions)[Funding](https://www.mdpi.com/2227-7390/14/6/984#Funding)[Data Availability Statement](https://www.mdpi.com/2227-7390/14/6/984#Data_Availability_Statement)[Conflicts of Interest](https://www.mdpi.com/2227-7390/14/6/984#Conflicts_of_Interest)[Appendix A](https://www.mdpi.com/2227-7390/14/6/984#Appendix_A)[References](https://www.mdpi.com/2227-7390/14/6/984#References)[Article Metrics](https://www.mdpi.com/2227-7390/14/6/984#Article_Metrics) - Article - ![Open Access](https://mdpi-res.com/cdn-cgi/image/w=14,h=14/https://mdpi-res.com/data/open-access.svg) 13 March 2026 # An Iterative Technique for Solving Time Fractional Systems with Variable Coefficients Wenlei Li and Juanying Huo\* School of Mathematics, Jilin University, Changchun 130012, China \* Author to whom correspondence should be addressed. [Mathematics](https://www.mdpi.com/journal/mathematics)**2026**, *14*(6), 984;<https://doi.org/10.3390/math14060984> This article belongs to the Section [C1: Difference and Differential Equations](https://www.mdpi.com/journal/mathematics/sections/difference_differential_equations) [Version Notes](https://www.mdpi.com/2227-7390/14/6/984/notes) [Order Reprints](https://www.mdpi.com/2227-7390/14/6/984/reprints) ## Abstract This paper presents a comprehensive study on constructing exact and approximate solutions to Cauchy problems for time-fractional systems with variable coefficients. An innovative iterative approach is developed for solving functional equations with initial conditions, combining rigorous mathematical foundations with practical computational efficiency. The proposed technique effectively handles the nonlocal nature of fractional operators through a carefully designed iterative scheme that maintains simplicity while achieving high accuracy. It demonstrates particular strength in solving nonlinear systems with well-defined conditions and variable coefficients, where traditional methods often fail. Through systematic theoretical analysis and numerical validation, we establish the method’s convergence properties and computational advantages, showing its capability to generate both exact closed-form solutions, when available, and high-precision approximations otherwise. The approach remains computationally tractable even for complex cases where variable coefficients and memory effects of fractional systems present significant challenges to conventional solution approaches. Keywords: [time fractional tricomi equation with variable coefficients](https://www.mdpi.com/search?q=time+fractional+tricomi+equation+with+variable+coefficients); [initial value problem](https://www.mdpi.com/search?q=initial+value+problem); [iterative technique](https://www.mdpi.com/search?q=iterative+technique) MSC: 35R11; 34A12 ## 1\. Introduction This paper presents a systematic investigation of the initial value problems arising in time-fractional partial differential equations, where the nonlocal nature of fractional derivatives introduces unique mathematical challenges that fundamentally distinguish them from their integer-order counterparts. The inherent memory effects and historical dependence of fractional operators render the initial value problem particularly intricate yet fascinating from both theoretical and computational perspectives. Our comprehensive study develops and analyzes a suite of analytical and numerical methodologies specifically designed to overcome these challenges, with the dual objectives of establishing rigorous mathematical foundations and providing practical computational tools for understanding the temporal evolution of solutions in such systems. The research builds upon and significantly extends existing work in this rapidly developing field, as documented in \[1,2,3,4,5,6,7,8\] and related references, while introducing novel approaches to address the distinctive features of fractional-order initial value problems. The historical development of fractional calculus provides crucial context for understanding its modern applications. Our exploration begins with the seminal correspondence between L’HĆ“pital and Leibniz in 1695, a foundational moment in mathematical history documented in \[9\]. In this historic exchange, L’HĆ“pital posed the now-famous question to Leibniz: ā€œWhat would be the result of d n y d x n if n \= 1 2 ?ā€ Leibniz’s insightful reply demonstrated remarkable foresight: ā€œ d 1 / 2 x will be equal to x d x : x . This is an apparent paradox, from which, one day useful consequences will be drawn.ā€ This early conceptual breakthrough, however, remained largely undeveloped for nearly two centuries. The mathematical community of the Enlightenment era faced several barriers to advancing fractional calculus: The lack of a rigorous theoretical framework for non-integer derivatives; The absence of clear physical interpretations; Significant computational challenges; The dominance of integer-order calculus in solving contemporary problems. This historical evolution underscores how theoretical mathematics often precedes practical applications, with Leibniz’s initial insight eventually finding numerous applications across science and engineering disciplines. The fundamental nature of fractional operators distinguishes them significantly from traditional differential operators. Researchers have established that these operators are intrinsically non-local, typically expressed through integral formulations with singular kernels. This essential characteristic gives fractional calculus its unique capacity to model memory effects and hereditary properties in physical systems, where the current state depends on the complete historical evolution rather than just instantaneous conditions. A particularly noteworthy application emerges in electrochemical systems, where fractional-order operators have proven remarkably effective: They provide superior modeling of diffusion processes in electrochemical cells compared to classical approaches; They more accurately characterize electrode impedance behavior; They offer enhanced precision in describing charge transfer mechanisms. The theoretical framework of fractional calculus has consequently found extensive utility across numerous scientific and technological domains, including: modeling anomalous diffusion and viscoelastic material properties, analyzing complex reaction-diffusion processes, understanding neuronal dynamics and biomechanical tissue behavior, characterizing long-range dependent processes in financial markets and so on. These diverse applications, meticulously documented in \[10\] and extensively developed in subsequent research, vividly illustrate the remarkable transformation of fractional calculus from a purely theoretical construct to an indispensable analytical framework in modern science and engineering. The continually expanding scope of its applications provides compelling evidence for the fundamental importance of these non-local operators, with demonstrated successes ranging from nanoscale phenomena to large-scale economic systems. Their unique capacity to capture long-range dependencies and memory effects has enabled breakthroughs in modeling complex systems that defy description by classical approaches, particularly in situations involving power-law behaviors and anomalous dynamics. Recent advances have particularly highlighted how fractional operators offer superior fidelity in modeling multiscale phenomena with coupled temporal and spatial dependencies, while also providing natural frameworks for systems exhibiting complex hierarchical or fractal-like structures. The operational success of fractional calculus across such disparate domains not only underscores its profound theoretical significance but has also driven important innovations in computational methods and analytical techniques. This dual impact—simultaneously deepening our fundamental understanding while enhancing practical problem-solving capabilities—continues to motivate vigorous research efforts at the frontiers of applied mathematics and interdisciplinary science, promising further discoveries and applications in emerging fields. Fractional calculus, representing a significant generalization of classical differential and integral operators to non-integer orders, has emerged as a vibrant research field attracting considerable attention from both theoretical and applied perspectives. This mathematical framework extends the traditional notions of differentiation and integration to fractional orders, thereby enabling the description of complex phenomena that exhibit non-local behavior and memory effects. Throughout its historical development, numerous distinct definitions of fractional derivatives and integrals have been proposed, each offering unique advantages for specific applications and mathematical contexts. The absence of a unified definition for fractional derivatives has led to a rich diversity of operational formulations, including several prominent variants that have proven particularly useful in practice. Among these, the Grünwald–Letnikov approach provides a discrete interpretation based on finite differences, while the Riemann–Liouville formulation offers a continuous integral-based perspective. The Caputo derivative has gained prominence for its advantageous treatment of initial conditions in differential equations. Additional important formulations include the Riesz potential for symmetric fractional operations, the conformable fractional derivative for maintaining certain classical properties, and various local fractional derivatives designed for specific applications. Within this diverse landscape, the Riemann–Liouville fractional integral and Caputo fractional derivative have emerged as particularly influential, owing to their mathematical tractability and physical interpretability. These operators have found extensive application across numerous scientific and engineering disciplines, including advanced control systems where they enable more robust controller designs, the modeling of anomalous diffusion processes in complex media, and the analysis of fractal structures and multi-scale systems. The continued development and refinement of these fractional operators, as documented in \[4,7,8,11,12,13,14\] and related works, reflects their growing importance in addressing challenging problems that lie beyond the scope of classical calculus. Among the various formulations of fractional operators, the Grünwald–Letnikov fractional derivative (see [Appendix A.2](https://www.mdpi.com/2227-7390/14/6/984#secAdot2-mathematics-14-00984)) occupies a special historical and conceptual position as one of the most intuitive and computationally accessible approaches. This operator, rooted in the generalization of finite difference methods, provides a natural discrete framework for fractional differentiation that closely parallels the classical limit definition of integer-order derivatives. The Grünwald–Letnikov approach is particularly valuable in numerical analysis and scientific computing, where its discrete nature allows for efficient implementation in algorithmic solutions to fractional differential equations. Its computational advantages have made it indispensable for simulations of complex systems with memory effects, as extensively documented in \[15\] and subsequent computational studies. In contrast, the Riemann–Liouville fractional derivative represents a continuous formulation defined through a convolution operation with a power-law kernel. This definition enjoys several theoretical advantages, including a rigorous mathematical foundation in functional analysis and a natural connection to the theory of integral equations. The Riemann–Liouville operator has proven particularly powerful in analytical approaches to fractional differential equations arising in mathematical physics, where its properties facilitate the development of existence theorems and solution representations. The complementary strengths of these two approaches—with Grünwald–Letnikov favoring computational implementations and Riemann–Liouville supporting theoretical analysis—serve to illustrate the rich diversity of tools available in fractional calculus, as discussed in detail in \[16\] and related mathematical physics literature. In contrast, the Caputo fractional derivative was introduced to address this issue, offering a more physically intuitive way to incorporate initial conditions in models. It is particularly popular in applied settings such as control systems. The Riesz potential, often employed in potential theory and the fractional Laplacian, generalizes the idea of a potential field to fractional dimensions and is crucial in studying processes like anomalous diffusion and LĆ©vy flights; see \[17\] and the references therein. In contrast to the Riemann–Liouville formulation, the Caputo fractional derivative was specifically developed to overcome certain mathematical challenges in modeling physical systems. This alternative definition provides several crucial advantages for practical applications, most notably its ability to handle initial conditions in a manner that aligns more naturally with physical interpretations. The Caputo operator’s formulation ensures that initial value problems can be specified using standard function values rather than fractional derivatives of the initial state, making it particularly valuable for modeling real-world phenomena. This property has led to its widespread adoption in applied mathematics and engineering disciplines, especially in control theory, where it enables more natural formulations of dynamical systems with memory effects. The Caputo derivative’s compatibility with conventional physical interpretations has made it a preferred choice for modeling viscoelastic materials, thermal systems, and other processes where the initial state is typically specified in terms of classical variables. Similarly important in fractional analysis is the Riesz potential, which extends classical potential theory to fractional dimensions through an elegant symmetric formulation. This operator plays a fundamental role in the theory of the fractional Laplacians and has become indispensable for studying non-local phenomena. The Riesz transform provides a powerful mathematical framework for analyzing complex diffusion processes, particularly those exhibiting anomalous behavior that deviates from classical Brownian motion. Its applications span diverse areas, including the modeling of LĆ©vy flights in statistical physics, the analysis of non-local interactions in quantum mechanics, and the study of long-range correlations in complex systems. The theoretical foundations and practical applications of the Riesz potential are thoroughly examined in \[17\], along with its connections to other fractional operators and its growing importance in modern mathematical physics. The rich diversity of fractional derivative definitions, while reflecting the wide-ranging applications of fractional calculus across scientific disciplines, presents significant challenges for researchers working to select, compare, and apply these various operators. This proliferation of definitions—each with distinct properties and suited for different applications—can create confusion and complicate theoretical comparisons between results obtained using different formulations. Recognizing this challenge, the mathematical community has increasingly turned to axiomatization approaches that seek to unify these diverse operators within a coherent theoretical framework. In \[12\], the authors adopt such an axiomatic perspective, examining several fractional derivative definitions in a systematic way. By identifying the core properties that these operators share, they developed a general framework that facilitates a deeper understanding of their similarities and differences. For example, the axioms might include the following several items: - (Linearity) The operator is linear. - (Identity) The zero order derivative of a function returns the function itself. - (Backward compatibility) When the order is integer, the fractional derivative gives the same result as the classical derivative. - (Index law) For any α \< 0 and β \< 0 , we have D α D β f ( t ) \= D α \+ β f ( t ) . - (Generalized Leibniz rule) D α \[ f ( t ) g ( t ) \] \= āˆ‘ i \= 0 āˆž α i D i f ( t ) D α āˆ’ i g ( t ) . **Remark** **1\.** If α ∈ N , this is exactly the classical Leibniz rule. By the way, another generalization of this rule, such as D α \[ f ( t ) g ( t ) \] \= āˆ‘ i \= āˆ’ i n f t y āˆž α i \+ β D β \+ i f ( t ) D α āˆ’ β āˆ’ i g ( t ) . This approach not only highlights the underlying mathematical structure common to various fractional operators but also provides a way to compare and contrast their behaviors. By formalizing fractional calculus through axiomatization, it becomes easier to extend these operators to new contexts, such as multi-dimensional systems, stochastic processes, or even quantum mechanics, where fractional derivatives are increasingly finding applications; see \[5,18,19,20\] and the references therein. The axiomatization of fractional operators has profound implications for both the theoretical development of fractional calculus and its practical applications. On the theoretical side, it provides a clearer framework for exploring the connections between different types of fractional derivatives and integrals, potentially leading to the discovery of new operators or generalizations that encompass existing ones. It also helps to clarify which definitions are most suitable for specific problems, depending on the required properties of the derivative, such as memory effects, non-locality, or boundary behavior; see \[21,22\] and the references therein. From a practical standpoint, this unifying approach simplifies the application of fractional derivatives in diverse fields, from engineering to biology. For instance, in modeling materials with memory effects, the choice of a fractional derivative is crucial, and an axiomatic perspective can help guide that choice. Similarly, in complex systems where fractal geometry or anomalous diffusion is present, such as in finance, geology, or network theory, an axiomatized framework for fractional derivatives ensures that the correct operator is used to capture the necessary dynamics. In this paper, we mainly explore the construction of exact solutions or approximate solutions to the initial value problems of time fractional partial differential equations. This is a very direct problem in the study of the theory of differential equations. However, it is not easy and even rather difficult. In addition, solving fractional systems is usually more difficult than classical systems, for their operator is defined by an integral (a non-local operator with a singular kernel). Luckily, there are some different effective methods that have been developed to construct the exact solutions or the approximate solutions of fractional systems with initial conditions. Such as the iterative method \[23,24,25,26,27\], the homotopy analysis method \[28,29,30\], the Cole-Hopf transformation method \[31\], the residual power series method \[32,33,34,35,36\], the Laplace transform method \[37\], and so on. ### Outline of Paper This manuscript presents a systematic investigation of iterative techniques for solving initial value problems in time-fractional partial differential equations. The paper is organized as follows: In [Section 2](https://www.mdpi.com/2227-7390/14/6/984#sec2-mathematics-14-00984), we establish the fundamental mathematical framework by introducing essential concepts from fractional calculus, including the various definitions of fractional derivatives and their key properties. This section provides the necessary theoretical foundation by presenting important lemmas regarding fractional operators, particularly focusing on their non-local characteristics and memory effects. We carefully define all mathematical notations to ensure clarity throughout subsequent developments, paying special attention to the function spaces and operator norms relevant to our analysis of fractional differential equations. [Section 3](https://www.mdpi.com/2227-7390/14/6/984#sec3-mathematics-14-00984) develops an iterative approach for solving functional equations with initial conditions, with particular emphasis on applications to fractional-order systems. We present a comprehensive convergence analysis of the proposed iterative scheme, establishing sufficient conditions under which the method guarantees convergence to the true solution. The section includes detailed proofs of convergence theorems and examines the rate of convergence under various regularity assumptions. Furthermore, we discuss the computational advantages of this approach compared to traditional methods, highlighting its efficiency in handling the non-local nature of fractional operators. [Section 4](https://www.mdpi.com/2227-7390/14/6/984#sec4-mathematics-14-00984) demonstrates the practical implementation of our iterative technique through an in-depth study of the one-dimensional time-fractional Tricomi equation. We derive both exact solutions (where obtainable) and high-precision approximate solutions, comparing our results with known analytical solutions to validate the method’s accuracy. The section includes detailed case studies with different fractional orders and initial conditions, accompanied by rigorous error analysis and convergence rate verification. Special attention is given to the physical interpretation of solutions and their relationship to the fractional order parameter. The analysis is extended in [Section 5](https://www.mdpi.com/2227-7390/14/6/984#sec5-mathematics-14-00984) to multidimensional time-fractional Tricomi equations, where we address the additional complexities introduced by higher spatial dimensions. We develop dimensionally-adaptive versions of our iterative scheme and examine their performance for various initial boundary value problems. The section presents numerical results for two- and three-dimensional cases, including visualization of solution surfaces and comparison with alternative numerical methods. We particularly focus on the method’s scalability and its ability to handle different types of initial conditions, from smooth analytic functions to discontinuous data. Finally, [Section 6](https://www.mdpi.com/2227-7390/14/6/984#sec6-mathematics-14-00984) synthesizes our findings into a comprehensive conclusion, summarizing the key advantages and limitations of the proposed iterative technique. We discuss broader implications for the field of fractional calculus and suggest several directions for future research, including potential applications to nonlinear fractional systems and coupled fractional-order equations. The section concludes with practical recommendations for implementing the method and its possible extensions to other classes of fractional differential equations. ## 2\. Fractional Calculus This subsection introduces the fundamental theoretical framework necessary for our investigation of fractional systems. We begin with the essential special functions that form the backbone of fractional calculus. The Gamma function, which extends the factorial concept to continuous values, plays an indispensable role in defining fractional operators due to its appearance in their integral kernels. Equally crucial is the Mittag–Leffler function, a sophisticated generalization of the exponential function that naturally describes processes with memory effects and power-law relaxation. The mathematical landscape of fractional calculus features a diverse collection of definitions for fractional integrals and derivatives, each developed to serve specific theoretical or applied purposes. Among these various formulations, several have emerged as particularly significant. The Riemann–Liouville fractional integral represents the most direct extension of classical integration theory to non-integer orders, while the Caputo fractional derivative has become preferred for physical modeling due to its more intuitive treatment of initial conditions. Other important variants include the Riesz potential, which provides symmetric fractional operations valuable for certain boundary value problems, and the Weyl derivative, especially useful for analyzing periodic phenomena. The comparative advantages and specific applications of these different definitions are thoroughly examined in the comprehensive references \[7,8,11,12,13,14\]. Within this spectrum of fractional operators, the Riemann–Liouville integral and Caputo derivative have gained particular prominence in both theoretical and applied contexts. The Riemann–Liouville approach provides the most natural mathematical foundation for fractional integration, while the Caputo formulation offers distinct practical advantages for solving initial value problems that model physical systems. This preference stems from the Caputo derivative’s compatibility with conventional initial conditions specified in terms of standard derivatives, making it more suitable for applications in engineering and physics. Our subsequent analysis will build upon several key lemmas concerning the operational properties of these fractional operators and the characteristic behaviors of Mittag–Leffler functions. These fundamental results enable rigorous investigation of existence, uniqueness, and qualitative properties for solutions to fractional differential equations. The interplay between these various concepts creates a robust framework for examining the nonlinear fractional systems that constitute the primary focus of this work, while also providing the necessary tools for developing effective solution methods. ### 2\.1. The Gamma Function and the Mittag–Leffler Function This subsection introduces the definitions for the Gamma function and the Mittag–Leffler function, which will be used in the definitions of the fractional integral and derivative. **Definition** **1** (\[38\])**.** The Gamma function is defined by Ī“ ( x ) \= ∫ 0 āˆž t x āˆ’ 1 e āˆ’ t d t Re ( x ) \> 0 , where Re ( x ) denotes the real part of the complex number x. **Proposition** **1** (\[38\])**.** Some main properties of the Gamma function: - Ī“ ( x \+ 1 ) \= x Ī“ ( x ) for all x with Re x \> 0 . - Ī“ ( 1 ) \= 1 . - Ī“ ( 1 2 ) \= Ļ€ . **Definition** **2** (\[39\])**.** The two-parameters Mittag–Leffler function is defined by letting E α , β ( x ) \= āˆ‘ k \= 0 āˆž x k Ī“ ( α k \+ β ) α , β \> 0 . And the one-parameter Mittag–Leffler function is defined by letting E α ( x ) \= āˆ‘ k \= 0 āˆž x k Ī“ ( α k \+ 1 ) α \> 0 . **Remark** **2\.** The two-parameter Mittag–Leffler function is a generalization of the one-parameter Mittag–Leffler function, that is, E α , 1 ( x ) \= E α ( x ) . The one-parameter Mittag–Leffler function is a generalization of the exponential function, or rather, E 1 ( x ) \= e x . ### 2\.2. The Riemann–Liouville Fractional Integral and the Caputo Derivative This subsection recalls the definition of the Riemann–Liouville fractional integral and Caputon derivative, see \[4,5,6,13,40\] and the references therein. **Definition** **3\.** Let Īø : \[ 0 , T \] → R is a continuous function. Then, the fractional-order integral of Īø in the Riemann–Liouville sense is defined by J α Īø ( t ) \= 1 Ī“ ( α ) ∫ 0 t ( t āˆ’ s ) α āˆ’ 1 Īø ( s ) d s , α \> 0 , t ≄ 0 , (1) where Ī“ ( α ) \= ∫ 0 āˆž e āˆ’ u u α āˆ’ 1 d u . **Definition** **4\.** Let Īø : \[ 0 , T \] → R is a continuous function. Then, the Caputo fractional-order derivative of Īø is given by D α c Īø ( t ) \= 1 Ī“ ( n āˆ’ α ) ∫ 0 t ( t āˆ’ s ) n āˆ’ α āˆ’ 1 Īø ( m ) ( s ) d s , where m āˆ’ 1 \< α ≤ m ; m ∈ N . **Lemma** **1** (\[27\])**.** For any m āˆ’ 1 \< β ≤ m ∈ N āˆ— , one has I β D β u ( x , t ) \= u ( x , t ) āˆ’ āˆ‘ k \= 0 m āˆ’ 1 āˆ‚ k u ( x , 0 ) t k k \! . **Lemma** **2\.** For any β , γ , one has I β t γ \= Ī“ ( γ \+ 1 ) Ī“ ( γ \+ β \+ 1 ) t γ \+ β . ## 3\. An Iterative Technique for Solving a Class of Nonlinear Functional Equations In this section, we present a fundamental lemma that plays a pivotal role in our subsequent mathematical analysis. It provides the essential theoretical foundation for establishing the key propositions and theorems developed throughout this paper. It has been chosen for its particular relevance to our investigation of fractional differential systems, offering crucial insights into the behavior of solutions and the properties of fractional operators. Its proper application will be carefully demonstrated in subsequent sections where we develop our main results. **Lemma** **3\.** Consider a class of abstract functional equations of the following form u ( x , t ) \= f ( x , t ) \+ L u ( x , t ) \+ N u ( x , t ) , (2) where ( t , x ) ∈ D \= { ( t , x ) : t ∈ R , x ∈ R n , n ∈ N } , u ( x , t ) is the unknown function, f is a given known function, L is a linear operator and N is a nonlinear operator with respect to u from a Banach space Dom ( L ) ∩ Dom ( N ) to itself (Dom ( L ) ∩ Dom ( N ) stands for the domain of L and N respectively). For abstract function Equation ([2](https://www.mdpi.com/2227-7390/14/6/984#FD2-mathematics-14-00984)), consider the series solutions of the following form u ( t , x ) \= āˆ‘ n \= 0 āˆž u n ( t , x ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the iterative sequence { u n } can be obtained by the following recurrence equations u 0 \= f , (3) u 1 \= L ( u 0 ) \+ N ( u 0 ) , (4) u m \+ 1 \= L ( u m ) \+ N āˆ‘ n \= 0 m u n āˆ’ N āˆ‘ n \= 0 m āˆ’ 1 u n , m \= 1 , 2 ⋯ . (5) Furthermore, if the operators L and N are contractive, then the series solution āˆ‘ n \= 0 āˆž u n ( t , x ) converges absolutely and uniformly. **Proof.** Obviously, the nonlinear operator N can be decomposed as N ( u ( t , x ) ) \= N āˆ‘ n \= 0 āˆž u n ( t , x ) \= N ( u 0 ) \+ āˆ‘ n \= 1 āˆž N āˆ‘ j \= 0 n u j āˆ’ N āˆ‘ j \= 0 n āˆ’ 1 u j . Similarly, the linear operator L can also be decomposed as L ( u ( t , x ) ) \= L āˆ‘ n \= 0 āˆž u n ( t , x ) \= L ( u 0 ) \+ āˆ‘ n \= 1 āˆž L āˆ‘ j \= 0 n u j āˆ’ L āˆ‘ j \= 0 n āˆ’ 1 u j . Denote by M u ( t , x ) \= L u ( t , x ) \+ N u ( t , x ) . then one has M ( u ( t , x ) ) \= M ( u 0 ) \+ āˆ‘ n \= 1 āˆž M āˆ‘ j \= 0 n u j āˆ’ M āˆ‘ j \= 0 n āˆ’ 1 u j . Define the following recurrence equations u 0 \= f , u 1 \= M ( u 0 ) , u m \+ 1 \= M āˆ‘ n \= 0 m u n āˆ’ M āˆ‘ n \= 0 m āˆ’ 1 u n , m \= 1 , 2 ⋯ . Since operators L and N are contractive, then M is also contractive, i.e., there exists a constant 0 \< K \< 1 , such that ∄ M ( v i ) āˆ’ M ( v j ) ∄ ≤ K ∄ v i āˆ’ v j ∄ , āˆ€ v i , v j ∈ Dom ( L ) ∩ D o m ( N ) , where ∄ Ā· ∄ denotes the norm on Banach space Dom ( L ) ∩ Dom ( N ) . Furthermore, for u m \+ 1 , one can obtain ∄ u m \+ 1 ∄ \= M āˆ‘ n \= 0 m u n āˆ’ M āˆ‘ n \= 0 m āˆ’ 1 u n ≤ K ∄ u m ∄ ≤ K m \+ 1 ∄ u 0 ∄ . Since 0 \< K \< 1 , the series āˆ‘ m \= 1 āˆž K m \+ 1 ‖ u 0 ‖ converges absolutely as well as uniformly. According to the Weierstrass M-test, one can obtain that the series āˆ‘ i \= 0 āˆž u i converges absolutely as well as uniformly. ā–” Some more specific details about this type of iterative method could be found in \[23,24,25,26,27\] and the references therein. ## 4\. One-Dimensional Time Fractional Tricomi Equation with Variable Coefficients Consider the following one-dimensional time-fractional Tricomi equation of order α , as studied in \[41\]: D t α u \+ ξ ( t ) āˆ‚ x 2 u \= f ( x , t ) , (6) subject to the initial conditions u ( x , 0 ) \= ω 1 ( x ) , D t u ( x , 0 ) \= ω 2 ( x ) , (7) where 1 \< α ≤ 2 represents the fractional order of the time derivative, t ∈ \[ 0 , T \] denotes the temporal variable, and x ∈ R is the spatial variable. The special case when α \= 2 reduces ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) to the classical Tricomi equation, first introduced by Francesco Tricomi in his seminal 1923 work \[42\]. This equation occupies a particularly important position in the theory of partial differential equations due to its remarkable property of exhibiting mixed-type characteristics. Specifically, the nature of the equation undergoes a fundamental change depending on the sign of the coefficient function ξ ( t ) : when ξ ( t ) \> 0 , the equation demonstrates elliptic behavior, while for ξ ( t ) \< 0 , it becomes hyperbolic. This dual character makes the Tricomi equation a paradigmatic example of mixed-type partial differential equations. The physical significance of this equation was further elucidated by Frankl \[43\], who established its crucial connection to transonic gas dynamics. More precisely, the Tricomi equation provides a mathematical model describing the transition between subsonic and supersonic flow regimes in gas dynamics. The elliptic region ( ξ ( t ) \> 0 ) corresponds to subsonic flow, where disturbances propagate in all directions, while the hyperbolic region ( ξ ( t ) \< 0 ) represents supersonic flow, characterized by directional wave propagation and the formation of shock waves. This interpretation underscores the equation’s importance in aerodynamics and fluid mechanics. In the following analysis, we will develop an iterative solution technique to solve the initial value problem ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) and ([7](https://www.mdpi.com/2227-7390/14/6/984#FD7-mathematics-14-00984)) for the time-fractional generalization of this important equation. Our approach will carefully account for both the fractional temporal derivative and the mixed-type spatial operator, while maintaining the essential physical interpretation of the solutions. By integrating both sides of Equation ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) with respect to the temporal variable from 0 to t, we derive the equivalent integral formulation: u ( x , t ) \= u ( x , 0 ) \+ D t u ( x , 0 ) t āˆ’ I t α ξ ( t ) āˆ‚ x 2 u \+ I t α f ( x , t ) . (8) Incorporating the initial conditions ([7](https://www.mdpi.com/2227-7390/14/6/984#FD7-mathematics-14-00984)) into the integral Equation ([8](https://www.mdpi.com/2227-7390/14/6/984#FD8-mathematics-14-00984)) yields: u ( x , t ) \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) āˆ’ I t α ξ ( t ) āˆ‚ x 2 u . From the perspective of iterative methods, we can interpret this equation as follows: The terms ω 1 ( x ) \+ ω 2 ( x ) \+ I t α f ( x , t ) constitute the initial iterative data, incorporating both the given initial conditions and the source term f ( x , t ) . The remaining term āˆ’ I t α ξ ( t ) āˆ‚ x 2 u acts as a nonlinear feedback operator that propagates the solution through successive iterations. We propose a series solution of the form: u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the iterative sequence is determined by the following recursive scheme: u 0 \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) , u 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 , u 2 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 , ā‹® u m \+ 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u m . To demonstrate the effectiveness of this iterative approach, we now examine several specific cases with carefully chosen initial conditions, source terms, and coefficient functions. **Proposition** **2** (Trivial Solution Case)**.** For the simplest case where ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , and f ( x , t ) ≔ 0 , we observe that for any well-behaved coefficient function ξ ( t ) , the solution identically vanishes: u ( x , t ) ≔ 0 . This trivial solution serves as an important consistency check for our method. **Proposition** **3** (Exponential Source Term with Constant Coefficient)**.** Consider the case with zero initial conditions ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , an exponential source term f ( x , t ) ≔ e x , and constant coefficient ξ ( t ) ≔ āˆ’ 1 . The iterative scheme yields: u n \= e x t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) āˆ€ n ∈ N . The complete solution then takes the elegant form: u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x āˆ‘ n \= 0 āˆž t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) \= e x E α ( t α ) āˆ’ 1 , where E α ( Ā· ) denotes the Mittag–Leffler function, which naturally emerges in solutions to fractional differential equations. **Proposition** **4** (Modified Exponential Source Term)**.** For the case with ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x t α Ī“ ( āˆ’ α \+ 1 ) , and ξ ( t ) ≔ āˆ’ 1 , we obtain: u n \= e x t n α Ī“ ( n α \+ 1 ) āˆ€ n ∈ N . The solution series converges to: u ( x , t ) \= e x E α ( t α ) , demonstrating how different source terms can lead to distinct Mittag–Leffler function representations. **Proposition** **5** (Non-constant Coefficient Case)**.** Consider the problem with ω 1 ( x ) \= e x , ω 2 ( x ) \= 0 , f ( x , t ) \= 0 , ξ ( t ) \= t α . Then the iterative scheme generates the sequence u 0 \= e x , u 1 \= āˆ’ e x Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α , u 2 \= e x Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . By the induction of mathematics, we establish the general term: u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α , n ≄ 0 . The complete solution is expressed as the infinite series: u ( x , t ) \= e x āˆ‘ n \= 0 āˆž ( āˆ’ 1 ) n āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α , revealing the intricate relationship between the solution’s temporal behavior and the fractional order α. **Proof.** **Step 1: Iterative structure.** Under the given assumptions, the iteration takes the form u n \+ 1 ( x , t ) \= āˆ’ I t α t α u n ( x , t ) , where I t α denotes the fractional integral operator. Because the initial data is separable and u 0 \= e x , we observe that each iterate preserves the spatial factor e x . Hence we may write u n ( x , t ) \= e x a n ( t ) , reducing the problem to determining the scalar sequence a n ( t ) . **Step 2: First iterates.** For n \= 0 , a 0 ( t ) \= 1 . Using the fractional integral formula I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain a 1 ( t ) \= āˆ’ Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α . Repeating the same computation yields a 2 ( t ) \= Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . **Step 3: Inductive hypothesis.** Assume that for some n ≄ 1 , u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α . **Step 4: Inductive step.** Applying the iteration formula, u n \+ 1 \= āˆ’ I t α t α u n . Substituting the inductive expression gives t α u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t ( 2 n \+ 1 ) α . Applying the fractional integral formula once more, I t α t ( 2 n \+ 1 ) α \= Ī“ ( ( 2 n \+ 1 ) α \+ 1 ) Ī“ ( ( 2 n \+ 2 ) α \+ 1 ) t ( 2 n \+ 2 ) α . Therefore, u n \+ 1 \= ( āˆ’ 1 ) n \+ 1 e x āˆ k \= 1 n \+ 1 Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n \+ 1 Ī“ ( 2 k α \+ 1 ) t 2 ( n \+ 1 ) α . This completes the induction. **Step 5: Series representation.** Summing the iterates yields the solution as the convergent power-type series u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) , which gives the stated expansion. ā–” ## 5\. Higher-Dimensional Time Fractional Tricomi Equation with Variable Coefficients Consider the following higher-dimensional time fractional Tricomi equation D t α u \+ ξ ( t ) Ī” x u \= f ( x , t ) (9) with initial conditions u ( x , 0 ) \= ω 1 ( x ) , D t u ( x , 0 ) \= ω 2 ( x ) . (10) where 1 \< α ≤ 2 , t ∈ \[ 0 , T \] , x ∈ R n , and Ī” x denotes the standard Laplace operator and it is defined by Ī” x u ( x , t ) \= āˆ‘ i \= 1 n āˆ‚ x i 2 u ( x , t ) . Similarly, one can obtain the equivalent integral Equations ([9](https://www.mdpi.com/2227-7390/14/6/984#FD9-mathematics-14-00984)) and ([10](https://www.mdpi.com/2227-7390/14/6/984#FD10-mathematics-14-00984)) as follows u ( x , t ) \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) āˆ’ I t α ξ ( t ) Ī” x u ( x , t ) . Suppose that the unknown function u has the series form u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the sequence { u n } can be obtained by the following iterative process u 0 \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) u 1 \= āˆ’ I t α ξ ( t ) Ī” x u 0 u 2 \= āˆ’ I t α ξ ( t ) Ī” x u 1 ā‹® u m \+ 1 \= āˆ’ I t α ξ ( t ) Ī” x u m . **Proposition** **6\.** If ω 1 ( x ) ≔ 0 , ω 2 ≔ 0 , f ( x , t ) ≔ 0 , then for any reasonable ξ ( t ) , u ( x , t ) ≔ 0 . **Proposition** **7\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 \+ x 2 \+ ⋯ \+ x n and ξ ( t ) ≔ āˆ’ 1 . - **Iterative construction.** Since the initial data vanish, the first term of the iteration is obtained directly from the source term: u 0 \= e x 1 \+ ⋯ \+ x n t α Ī“ ( α \+ 1 ) . Applying the iterative scheme once more and using I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain u 1 \= e x 1 \+ ⋯ \+ x n t 2 α Ī“ ( 2 α \+ 1 ) . Proceeding inductively, one easily verifies that u n \= e x 1 \+ ⋯ \+ x n t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing all iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 1 āˆž t n α Ī“ ( n α \+ 1 ) . Rewriting the series, u ( x , t ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) āˆ’ 1 . Hence, u ( x , t ) \= e x 1 \+ ⋯ \+ x n E α ( t α ) āˆ’ 1 . **Proposition** **8\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 \+ x 2 \+ ⋯ \+ x n t α Ī“ ( 1 āˆ’ α ) , ξ ( t ) ≔ āˆ’ 1 . **Iterative construction.** Since the initial data vanish and the source term already contains t α , the first iterate reduces to u 0 \= e x 1 \+ ⋯ \+ x n . Applying the iterative scheme once, u 1 \= e x 1 \+ ⋯ \+ x n t α Ī“ ( α \+ 1 ) . Continuing in the same manner, one verifies inductively that u n \= e x 1 \+ ⋯ \+ x n t n α Ī“ ( n α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing the iterates gives u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) . Recognizing the Mittag–Leffler series, E α ( z ) \= āˆ‘ n \= 0 āˆž z n Ī“ ( n α \+ 1 ) , we conclude that u ( x , t ) \= e x 1 \+ ⋯ \+ x n E α ( t α ) . **Proposition** **9\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ e x , f ( x , t ) ≔ 0 , and ξ ( t ) ≔ āˆ’ t α Ī“ ( α \+ 1 ) . **Iterative construction.** Since the source term vanishes and the initial velocity is e x , the first term of the iteration is u 0 \= e x . Proceeding with the iterative scheme, one obtains successively terms of increasing order in t α . By induction it follows that u n \= e x t n α Ī“ ( n α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing the iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) . Recognizing the Mittag–Leffler expansion, we conclude that u ( x , t ) \= e x E α ( t α ) . **Proposition** **10\.** Assume ω 1 ( x ) ≔ e x 1 \+ ⋯ \+ x n , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ 0 , ξ ( t ) ≔ t α . **First iterates.** The initial term is u 0 \= e x 1 \+ ⋯ \+ x n . Since āˆ‚ x 2 e x 1 \+ ⋯ \+ x n \= e x 1 \+ ⋯ \+ x n , the first iteration gives u 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 \= āˆ’ e x 1 \+ ⋯ \+ x n I t α t α \= āˆ’ e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α . Applying the scheme once more, u 2 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 \= e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) I t α t 3 α \= e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . **Inductive step.** Assume that for some k ≄ 1 , u k \= ( āˆ’ 1 ) k \+ 1 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 k α \+ 1 ) t 2 k α . Then u k \+ 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u k \= ( āˆ’ 1 ) k \+ 2 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) ⋯ Ī“ ( 2 k α \+ 1 ) I t α t ( 2 k \+ 1 ) α . Using I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain u k \+ 1 \= ( āˆ’ 1 ) k \+ 2 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) ⋯ Ī“ ( ( 2 k \+ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) ⋯ Ī“ ( ( 2 k \+ 2 ) α \+ 1 ) t ( 2 k \+ 2 ) α . Hence, by induction, for all n ∈ N , u n \= ( āˆ’ 1 ) n \+ 1 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 n α \+ 1 ) t 2 n α . **Series representation.** Summing the iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n , that is, u ( x , t ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž ( āˆ’ 1 ) n \+ 1 Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 n α \+ 1 ) t 2 n α . ## 6\. Concluding Remarks In this paper, we systematically investigate an innovative iterative technique for solving a broad class of functional equations that frequently arise in mathematical analysis and applied mathematics. Building upon established theoretical foundations, we develop and rigorously analyze this computational method, demonstrating its effectiveness through both theoretical considerations and practical applications. The primary focus of our study lies in the application of this methodology to initial value problems associated with time fractional partial differential equations. Through careful implementation of our approach, we are able to derive either exact analytical solutions or highly accurate approximate solutions, depending on the specific characteristics of each problem under consideration. A particularly noteworthy aspect of our research reveals that this iterative technique offers several significant advantages: (1) conceptual simplicity that facilitates straightforward implementation, (2) computational efficiency that enables practical solutions to complex problems, and (3) remarkable adaptability to various types of nonlinear systems. Most importantly, we establish that our method maintains these advantageous properties even when applied to challenging nonlinear systems of differential equations with initial conditions, where many traditional approaches often encounter difficulties. The robustness and versatility of our proposed method are thoroughly examined through multiple carefully selected examples, which serve to illustrate both the theoretical framework and practical applications. Our results suggest that this approach may provide a valuable addition to the existing toolkit for analyzing fractional differential equations and related functional equations in applied mathematics. ## Author Contributions Methodology, W.L. and J.H.; Validation, J.H.; Writing—original draft, J.H.; Writing—review and editing, W.L.; Visualization, W.L. and J.H. All authors have read and agreed to the published version of the manuscript. ## Funding This research received no external funding. ## Data Availability Statement No new data were created or analyzed in this study. ## Conflicts of Interest The authors declare no conflicts of interest. ## Appendix A ### Appendix A.1. The Conformable Fractional Derivative This subsection recalls the definition of the conformable fractional derivative, which looks like the limit of the difference quotient (the classical Newtonian derivative). The α \-th conformable fractional derivative of a function f is defined by letting ( T α f ) ( t ) \= lim ϵ → 0 f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t ) ϵ . Here 0 ≤ α ≤ 1 . It is clear that the case of α \= 1 coincides with the classical Newtonian derivative. Furthermore, one can verify that - The zero order derivative of a function does not return the function itself. In fact, we have ( T 0 f ) ( t ) \= lim ϵ → 0 f ( t \+ ϵ t ) āˆ’ f ( t ) ϵ \= lim ϵ → 0 f ( ( 1 \+ ϵ ) t ) āˆ’ f ( t ) 1 \+ ϵ Ā· lim ϵ → 0 1 \+ ϵ ϵ \= ( D f ) ( t ) Ā· āˆž ≠ f ( t ) , where D denotes the quantum derivative. - The index law does not hold, that is, ( T α T β f ) ( t ) ≠ ( T α \+ β f ) ( t ) . In fact, we have LHS \= lim ϵ → 0 1 ϵ 2 \[ f ( t \+ ϵ t 1 āˆ’ α \+ ϵ ( t \+ ϵ t 1 āˆ’ α ) 1 āˆ’ β ) āˆ’ f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t \+ ϵ t 1 āˆ’ β ) \+ f ( t ) \] and RHS \= lim ϵ → 0 f ( t \+ ϵ t 1 āˆ’ ( α \+ β ) ) āˆ’ f ( t ) ϵ . It is not difficult to see that ( T α T β f ) ( t ) ≠ ( T α \+ β f ) ( t ) without any restriction. - In this perspective of the criteria, the conformable fractional derivative is not a fractional derivative; see \[13\]. By the way, regarding the definition of the higher order conformable fractional derivative, please refer to \[44\] (Definition 2.1) and the references therein. ### Appendix A.2. The Grünwald–Letnikov Derivative This subsection recalls the definition of the Grünwald–Letnikov derivative. It is a type of fractional derivative that extends the concept of differentiation to non-integer orders. It is defined using a limit process similar to the traditional definition of the derivative but incorporates fractional differences. Specifically, it is expressed as a limit of a sum involving the values of a function at discrete points, weighted by binomial coefficients. This derivative is useful in various applications, particularly in fields like physics and engineering, where systems exhibit non-local or memory effects. Its flexibility allows for modeling behaviors that standard derivatives cannot capture. More precisely, it is defined by letting ( GL D Īø α f ) ( z ) \= e āˆ’ i α Īø lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( z āˆ’ k h ) \| h \| α , where z ∈ C and h \= \| h \| e i Īø , āˆ’ Ļ€ \< Īø ≤ Ļ€ . When the complex independent variable z becomes the real number t, we have the following two types of fractional derivative: - Forward derivative ( h \= \| h \| ): ( D F α GL f ) ( t ) \= lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( t āˆ’ k h ) \| h \| α - Backward derivative ( h \= āˆ’ \| h \| ): ( D B α GL f ) ( t ) \= lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( t \+ k \| h \| ) \| h \| α We introduce only the forward derivative. - The operator is linear, that is, D F α GL ( Ī» f \+ ν g ) ( z ) \= Ī» ( GL D F α f ) ( z ) \+ μ ( GL D F α g ) ( z ) , āˆ€ Ī» , ν ∈ R . - Additivity and commutativity: D F α GL D F β GL ( f ) ( t ) \= D F β GL D F α GL ( f ) ( t ) \= ( D F α \+ β GL ( f ) ) ( t ) . - Neutral element. Letting β \= āˆ’ α yields that D F α GL D F āˆ’ α GL ( f ) ( t ) \= f ( t ) and D F āˆ’ α GL D F α GL ( f ) ( t ) \= f ( t ) . - Backward compatibility. - The generalized Leibniz rule. That is, for the derivative of the product φ ψ , we have D F α G L \[ φ ψ \] ( t ) \= āˆ‘ n \= 0 āˆž C n α D n φ ( t ) D F α āˆ’ n G L ψ ( t ) . ## References 1. Nguyen, A.T.; Caraballo, T.; Tuan, N.H. On the initial value problem for a class of nonlinear biharmonic equation with time-fractional derivative. Proc. R. Soc. Edinb. Sect. A Math. **2022**, 152, 989–1031. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+the+initial+value+problem+for+a+class+of+nonlinear+biharmonic+equation+with+time-fractional+derivative&author=Nguyen,+A.T.&author=Caraballo,+T.&author=Tuan,+N.H.&publication_year=2022&journal=Proc.+R.+Soc.+Edinb.+Sect.+A+Math.&volume=152&pages=989%E2%80%931031&doi=10.1017/prm.2021.44)\] \[[CrossRef](https://doi.org/10.1017/prm.2021.44)\] 2. 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](https://scholar.google.com/scholar_lookup?title=General+time-fractional+diffusion+equation:+Some+uniqueness+and+existence+results+for+the+initial-boundary-value+problems&author=Luchko,+Y.&author=Yamamoto,+M.&publication_year=2016&journal=Fract.+Calc.+Appl.+Anal.&volume=19&pages=676%E2%80%93695&doi=10.1515/fca-2016-0036)\] \[[CrossRef](https://doi.org/10.1515/fca-2016-0036)\] 3. Luchko, Y. Initial-boundary-value problems for the one-dimensional time-fractional diffusion equation. Fract. Calc. Appl. Anal. **2012**, 15, 141–160. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Initial-boundary-value+problems+for+the+one-dimensional+time-fractional+diffusion+equation&author=Luchko,+Y.&publication_year=2012&journal=Fract.+Calc.+Appl.+Anal.&volume=15&pages=141%E2%80%93160&doi=10.2478/s13540-012-0010-7)\] \[[CrossRef](https://doi.org/10.2478/s13540-012-0010-7)\] 4. Karthikeyan, K.; Karthikeyan, P.; Patel, D.K.; Al-Askar, F.M.; Cesarano, C. Solving fractional integro-differential equations with delay and relaxation impulsive terms by fixed point techniques. Bound. Value Probl. **2024**, 2024, 146. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solving+fractional+integro-differential+equations+with+delay+and+relaxation+impulsive+terms+by+fixed+point+techniques&author=Karthikeyan,+K.&author=Karthikeyan,+P.&author=Patel,+D.K.&author=Al-Askar,+F.M.&author=Cesarano,+C.&publication_year=2024&journal=Bound.+Value+Probl.&volume=2024&pages=146&doi=10.1186/s13661-024-01957-w)\] \[[CrossRef](https://doi.org/10.1186/s13661-024-01957-w)\] 5. Abbas, M.; Ciobanescu, C.; Asghar, M.W.; Omame, A. Solution approximation of fractional boundary value problems and convergence analysis using AA-iterative scheme. AIMS Math. **2024**, 9, 13129–13158. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solution+approximation+of+fractional+boundary+value+problems+and+convergence+analysis+using+AA-iterative+scheme&author=Abbas,+M.&author=Ciobanescu,+C.&author=Asghar,+M.W.&author=Omame,+A.&publication_year=2024&journal=AIMS+Math.&volume=9&pages=13129%E2%80%9313158&doi=10.3934/math.2024641)\] \[[CrossRef](https://doi.org/10.3934/math.2024641)\] 6. Alsaedi, A.; Ahmad, B.; Alharbi, N.; Alharbi, N.; Alharbi, N. Fixed point and stability analysis of a tripled system of nonlinear fractional differential equations with n-nonlinear terms. Fractal Fract. **2024**, 8, 697. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fixed+point+and+stability+analysis+of+a+tripled+system+of+nonlinear+fractional+differential+equations+with+n-nonlinear+terms&author=Alsaedi,+A.&author=Ahmad,+B.&author=Alharbi,+N.&author=Alharbi,+N.&author=Alharbi,+N.&publication_year=2024&journal=Fractal+Fract.&volume=8&pages=697&doi=10.3390/fractalfract8120697)\] \[[CrossRef](https://doi.org/10.3390/fractalfract8120697)\] 7. Mali, A.D.; JĆørgensen, J.B.; KĆøstner, N.; LĆønborg, J.; Dutta, H. Numerical investigation of generalized tempered-type integrodifferential equations with respect to another function. Fract. Calc. Appl. Anal. **2023**, 26, 2580–2601. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Numerical+investigation+of+generalized+tempered-type+integrodifferential+equations+with+respect+to+another+function&author=Mali,+A.D.&author=J%C3%B8rgensen,+J.B.&author=K%C3%B8stner,+N.&author=L%C3%B8nborg,+J.&author=Dutta,+H.&publication_year=2023&journal=Fract.+Calc.+Appl.+Anal.&volume=26&pages=2580%E2%80%932601&doi=10.1007/s13540-023-00198-5)\] \[[CrossRef](https://doi.org/10.1007/s13540-023-00198-5)\] 8. Gunasekar, T.; Raghavendran, P.; Santra, S.S.; Majumder, D.; Baleanu, D.; Balasundaram, H. Application of laplace transform to solve fractional integro-differential equations. J. Math. Comput. Sci. **2024**, 33, 225–237. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Application+of+laplace+transform+to+solve+fractional+integro-differential+equations&author=Gunasekar,+T.&author=Raghavendran,+P.&author=Santra,+S.S.&author=Majumder,+D.&author=Baleanu,+D.&author=Balasundaram,+H.&publication_year=2024&journal=J.+Math.+Comput.+Sci.&volume=33&pages=225%E2%80%93237&doi=10.22436/jmcs.033.03.02)\] \[[CrossRef](https://doi.org/10.22436/jmcs.033.03.02)\] 9. Ross, B. The development of fractional calculus 1695–1900. Historia Math. **1977**, 4, 75–89. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+development+of+fractional+calculus+1695%E2%80%931900&author=Ross,+B.&publication_year=1977&journal=Historia+Math.&volume=4&pages=75%E2%80%9389&doi=10.1016/0315-0860\(77\)90039-8)\] \[[CrossRef](https://doi.org/10.1016/0315-0860\(77\)90039-8)\] 10. El-Ajou, A.; Arqub, O.A.; Zhour, Z.A.; Momani, S. New results on fractional power series: Theories and applications. Entropy **2013**, 15, 5305–5323. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=New+results+on+fractional+power+series:+Theories+and+applications&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Zhour,+Z.A.&author=Momani,+S.&publication_year=2013&journal=Entropy&volume=15&pages=5305%E2%80%935323&doi=10.3390/e15125305)\] \[[CrossRef](https://doi.org/10.3390/e15125305)\] 11. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Volume 204 of North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Theory+and+Applications+of+Fractional+Differential+Equations&author=Kilbas,+A.A.&author=Srivastava,+H.M.&author=Trujillo,+J.J.&publication_year=2006)\] 12. Ortigueira, M.D.; Machado, J.A.T. What is a fractional derivative? J. Comput. Phys. **2015**, 293, 4–13. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=What+is+a+fractional+derivative?&author=Ortigueira,+M.D.&author=Machado,+J.A.T.&publication_year=2015&journal=J.+Comput.+Phys.&volume=293&pages=4%E2%80%9313&doi=10.1016/j.jcp.2014.07.019)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2014.07.019)\] 13. Katugampola, U.N. Correction to ā€œWhat is a fractional derivative?ā€ by Ortigueira and Machado \[Journal of Computational Physics, Volume 293, 15 July 2015, Pages 4–13. Special issue on fractional PDEs\]. J. Comput. Phys. **2016**, 321, 1255–1257. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Correction+to+%E2%80%9CWhat+is+a+fractional+derivative?%E2%80%9D+by+Ortigueira+and+Machado+%5BJournal+of+Computational+Physics,+Volume+293,+15+July+2015,+Pages+4%E2%80%9313.+Special+issue+on+fractional+PDEs%5D&author=Katugampola,+U.N.&publication_year=2016&journal=J.+Comput.+Phys.&volume=321&pages=1255%E2%80%931257&doi=10.1016/j.jcp.2016.05.052)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2016.05.052)\] 14. Kuroda, L.K.B.; Gomes, A.V.; Tavoni, R.; Mancera, P.F.d.; Varalta, N.; Camargo, R.d.F. Unexpected behavior of Caputo fractional derivative. Comput. Appl. Math. **2017**, 36, 1173–1183. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Unexpected+behavior+of+Caputo+fractional+derivative&author=Kuroda,+L.K.B.&author=Gomes,+A.V.&author=Tavoni,+R.&author=Mancera,+P.F.d.&author=Varalta,+N.&author=Camargo,+R.d.F.&publication_year=2017&journal=Comput.+Appl.+Math.&volume=36&pages=1173%E2%80%931183&doi=10.1007/s40314-015-0301-9)\] \[[CrossRef](https://doi.org/10.1007/s40314-015-0301-9)\] 15. Murio, D.A. Stable numerical evaluation of grünwald–letnikov fractional derivatives applied to a fractional ihcp. Inverse Probl. Sci. Eng. **2009**, 17, 229–243. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Stable+numerical+evaluation+of+gr%C3%BCnwald%E2%80%93letnikov+fractional+derivatives+applied+to+a+fractional+ihcp&author=Murio,+D.A.&publication_year=2009&journal=Inverse+Probl.+Sci.+Eng.&volume=17&pages=229%E2%80%93243&doi=10.1080/17415970802082872)\] \[[CrossRef](https://doi.org/10.1080/17415970802082872)\] 16. Chikrii, A.; Matychyn, I. Riemann–liouville, caputo, and sequential fractional derivatives in differential games. In Advances in Dynamic Games: Theory, Applications, and Numerical Methods for Differential and Stochastic Games; Springer: Berlin/Heidelberg, Germany, 2011; pp. 61–81. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Riemann%E2%80%93liouville,+caputo,+and+sequential+fractional+derivatives+in+differential+games&author=Chikrii,+A.&author=Matychyn,+I.&publication_year=2011&pages=61%E2%80%9381)\] 17. Agrawal, O.P. Fractional variational calculus in terms of riesz fractional derivatives. J. Phys. Math. Theor. **2007**, 40, 6287. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fractional+variational+calculus+in+terms+of+riesz+fractional+derivatives&author=Agrawal,+O.P.&publication_year=2007&journal=J.+Phys.+Math.+Theor.&volume=40&pages=6287&doi=10.1088/1751-8113/40/24/003)\] \[[CrossRef](https://doi.org/10.1088/1751-8113/40/24/003)\] 18. Laskin, N. Principles of fractional quantum mechanics. arXiv **2010**, arXiv:1009.5533. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Principles+of+fractional+quantum+mechanics&author=Laskin,+N.&publication_year=2010&journal=arXiv&doi=10.48550/arXiv.1009.5533)\] \[[CrossRef](https://doi.org/10.48550/arXiv.1009.5533)\] 19. Laskin, N. Fractional quantum mechanics. Phys. Rev. E **2000**, 62, 3135. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fractional+quantum+mechanics&author=Laskin,+N.&publication_year=2000&journal=Phys.+Rev.+E&volume=62&pages=3135&doi=10.1103/PhysRevE.62.3135)\] \[[CrossRef](https://doi.org/10.1103/PhysRevE.62.3135)\] 20. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E **1997**, 55, 3581. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Mechanics+with+fractional+derivatives&author=Riewe,+F.&publication_year=1997&journal=Phys.+Rev.+E&volume=55&pages=3581&doi=10.1103/PhysRevE.55.3581)\] \[[CrossRef](https://doi.org/10.1103/PhysRevE.55.3581)\] 21. Falcini, F.; Garra, R.; Voller, V. Modeling anomalous heat diffusion: Comparing fractional derivative and non-linear diffusivity treatments. Int. J. Therm. Sci. **2019**, 137, 584–588. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Modeling+anomalous+heat+diffusion:+Comparing+fractional+derivative+and+non-linear+diffusivity+treatments&author=Falcini,+F.&author=Garra,+R.&author=Voller,+V.&publication_year=2019&journal=Int.+J.+Therm.+Sci.&volume=137&pages=584%E2%80%93588&doi=10.1016/j.ijthermalsci.2018.11.016)\] \[[CrossRef](https://doi.org/10.1016/j.ijthermalsci.2018.11.016)\] 22. Garra, R.; Falcini, F.; Voller, V.R.; Pagnini, G. A generalized stefan model accounting for system memory and non-locality. Int. Commun. Heat Mass Transf. **2020**, 114, 104584. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+generalized+stefan+model+accounting+for+system+memory+and+non-locality&author=Garra,+R.&author=Falcini,+F.&author=Voller,+V.R.&author=Pagnini,+G.&publication_year=2020&journal=Int.+Commun.+Heat+Mass+Transf.&volume=114&pages=104584&doi=10.1016/j.icheatmasstransfer.2020.104584)\] \[[CrossRef](https://doi.org/10.1016/j.icheatmasstransfer.2020.104584)\] 23. Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. **2006**, 316, 753–763. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=An+iterative+method+for+solving+nonlinear+functional+equations&author=Daftardar-Gejji,+V.&author=Jafari,+H.&publication_year=2006&journal=J.+Math.+Anal.+Appl.&volume=316&pages=753%E2%80%93763&doi=10.1016/j.jmaa.2005.05.009)\] \[[CrossRef](https://doi.org/10.1016/j.jmaa.2005.05.009)\] 24. Noor, M.A.; Noor, K.I.; Mohyud-Din, S.T.; Shabbir, A. An iterative method with cubic convergence for nonlinear equations. Appl. Math. Comput. **2006**, 183, 1249–1255. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=An+iterative+method+with+cubic+convergence+for+nonlinear+equations&author=Noor,+M.A.&author=Noor,+K.I.&author=Mohyud-Din,+S.T.&author=Shabbir,+A.&publication_year=2006&journal=Appl.+Math.+Comput.&volume=183&pages=1249%E2%80%931255&doi=10.1016/j.amc.2006.05.133)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2006.05.133)\] 25. Bhalekar, S.; Daftardar-Gejji, V. New iterative method: Application to partial differential equations. Appl. Math. Comput. **2008**, 203, 778–783. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=New+iterative+method:+Application+to+partial+differential+equations&author=Bhalekar,+S.&author=Daftardar-Gejji,+V.&publication_year=2008&journal=Appl.+Math.+Comput.&volume=203&pages=778%E2%80%93783&doi=10.1016/j.amc.2008.05.071)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2008.05.071)\] 26. Daftardar-Gejji, V.; Bhalekar, S. Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method. Comput. Math. Appl. **2010**, 59, 1801–1809. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solving+fractional+boundary+value+problems+with+Dirichlet+boundary+conditions+using+a+new+iterative+method&author=Daftardar-Gejji,+V.&author=Bhalekar,+S.&publication_year=2010&journal=Comput.+Math.+Appl.&volume=59&pages=1801%E2%80%931809&doi=10.1016/j.camwa.2009.08.018)\] \[[CrossRef](https://doi.org/10.1016/j.camwa.2009.08.018)\] 27. Dhaigude, C.D.; Nikam, V.R. Solution of fractional partial differential equations using iterative method. Fract. Calc. Appl. Anal. **2012**, 15, 684–699. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solution+of+fractional+partial+differential+equations+using+iterative+method&author=Dhaigude,+C.D.&author=Nikam,+V.R.&publication_year=2012&journal=Fract.+Calc.+Appl.+Anal.&volume=15&pages=684%E2%80%93699&doi=10.2478/s13540-012-0046-8)\] \[[CrossRef](https://doi.org/10.2478/s13540-012-0046-8)\] 28. Liao, S. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. **2004**, 147, 499–513. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+the+homotopy+analysis+method+for+nonlinear+problems&author=Liao,+S.&publication_year=2004&journal=Appl.+Math.+Comput.&volume=147&pages=499%E2%80%93513&doi=10.1016/S0096-3003\(02\)00790-7)\] \[[CrossRef](https://doi.org/10.1016/S0096-3003\(02\)00790-7)\] 29. Das, S.; Vishal, K. Homotopy analysis method for fractional Swift-Hohenberg equation. In Advances in the Homotopy Analysis Method; World Scientific Publishing: Hackensack, NJ, USA, 2014; pp. 291–308. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Homotopy+analysis+method+for+fractional+Swift-Hohenberg+equation&author=Das,+S.&author=Vishal,+K.&publication_year=2014&pages=291%E2%80%93308)\] \[[CrossRef](https://doi.org/10.1142/9789814551250_0007)\] 30. Demir, A.; Bayrak, M.A.; Ozbilge, E. A new approach for the approximate analytical solution of space-time fractional differential equations by the homotopy analysis method. Adv. Math. Phys. **2019**, 2019, 5602565. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+new+approach+for+the+approximate+analytical+solution+of+space-time+fractional+differential+equations+by+the+homotopy+analysis+method&author=Demir,+A.&author=Bayrak,+M.A.&author=Ozbilge,+E.&publication_year=2019&journal=Adv.+Math.+Phys.&volume=2019&pages=5602565&doi=10.1155/2019/5602565)\] \[[CrossRef](https://doi.org/10.1155/2019/5602565)\] 31. Li, L.; Li, D. Exact solutions and numerical study of time fractional Burgers’ equations. Appl. Math. Lett. **2020**, 100, 106011. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Exact+solutions+and+numerical+study+of+time+fractional+Burgers%E2%80%99+equations&author=Li,+L.&author=Li,+D.&publication_year=2020&journal=Appl.+Math.+Lett.&volume=100&pages=106011&doi=10.1016/j.aml.2019.106011)\] \[[CrossRef](https://doi.org/10.1016/j.aml.2019.106011)\] 32. Arqub, O.A. Series solution of fuzzy differential equations under strongly generalized differentiability. J. Adv. Res. Appl. Math. **2013**, 5, 31–52. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Series+solution+of+fuzzy+differential+equations+under+strongly+generalized+differentiability&author=Arqub,+O.A.&publication_year=2013&journal=J.+Adv.+Res.+Appl.+Math.&volume=5&pages=31%E2%80%9352&doi=10.5373/jaram.1447.051912)\] \[[CrossRef](https://doi.org/10.5373/jaram.1447.051912)\] 33. El-Ajou, A.; Arqub, O.A.; Momani, S.; Baleanu, D.; Alsaedi, A. A novel expansion iterative method for solving linear partial differential equations of fractional order. Appl. Math. Comput. **2015**, 257, 119–133. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+novel+expansion+iterative+method+for+solving+linear+partial+differential+equations+of+fractional+order&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Momani,+S.&author=Baleanu,+D.&author=Alsaedi,+A.&publication_year=2015&journal=Appl.+Math.+Comput.&volume=257&pages=119%E2%80%93133&doi=10.1016/j.amc.2014.12.121)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2014.12.121)\] 34. Dubey, V.P.; Kumar, R.; Kumar, D. A reliable treatment of residual power series method for time-fractional Black–Scholes European option pricing equations. Physical A **2019**, 533, 122040. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+reliable+treatment+of+residual+power+series+method+for+time-fractional+Black%E2%80%93Scholes+European+option+pricing+equations&author=Dubey,+V.P.&author=Kumar,+R.&author=Kumar,+D.&publication_year=2019&journal=Physical+A&volume=533&pages=122040&doi=10.1016/j.physa.2019.122040)\] \[[CrossRef](https://doi.org/10.1016/j.physa.2019.122040)\] 35. Syam, M.I. On approximate solutions of Fokker Planck equation by the modified residual power series method. J. Fract. Calc. Appl. **2019**, 10, 225–232. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+approximate+solutions+of+Fokker+Planck+equation+by+the+modified+residual+power+series+method&author=Syam,+M.I.&publication_year=2019&journal=J.+Fract.+Calc.+Appl.&volume=10&pages=225%E2%80%93232)\] 36. Az-Zo’bi, E.A.; Yıldırım, A.; AlZoubi, W.A. The residual power series method for the one-dimensional unsteady flow of a van der Waals gas. Physical A **2019**, 517, 188–196. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+residual+power+series+method+for+the+one-dimensional+unsteady+flow+of+a+van+der+Waals+gas&author=Az-Zo%E2%80%99bi,+E.A.&author=Y%C4%B1ld%C4%B1r%C4%B1m,+A.&author=AlZoubi,+W.A.&publication_year=2019&journal=Physical+A&volume=517&pages=188%E2%80%93196&doi=10.1016/j.physa.2018.11.030)\] \[[CrossRef](https://doi.org/10.1016/j.physa.2018.11.030)\] 37. Mosaleheh, K.; Vakilzadeh, A. The Laplace transform method for linear ordinary differential equations of fractional order. J. Math. Ext. **2007**, 2, 93–102. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+Laplace+transform+method+for+linear+ordinary+differential+equations+of+fractional+order&author=Mosaleheh,+K.&author=Vakilzadeh,+A.&publication_year=2007&journal=J.+Math.+Ext.&volume=2&pages=93%E2%80%93102)\] 38. Artin, E. The Gamma Function; Butler, M., Translator; Athena Series: Selected Topics in Mathematics; Holt, Rinehart and Winston: New York, NY, USA; Toronto, ON, Canada; London, UK, 1964. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+Gamma+Function&author=Artin,+E.&publication_year=1964)\] 39. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer Monographs in Mathematics; Springer: Heidelberg, Germany, 2014. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Mittag-Leffler+Functions,+Related+Topics+and+Applications&author=Gorenflo,+R.&author=Kilbas,+A.A.&author=Mainardi,+F.&author=Rogosin,+S.V.&publication_year=2014)\] 40. El-Ajou, A.; Arqub, O.A.; Momani, S. Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: A new iterative algorithm. J. Comput. Phys. **2015**, 293, 81–95. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Approximate+analytical+solution+of+the+nonlinear+fractional+KdV-Burgers+equation:+A+new+iterative+algorithm&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Momani,+S.&publication_year=2015&journal=J.+Comput.+Phys.&volume=293&pages=81%E2%80%9395&doi=10.1016/j.jcp.2014.08.004)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2014.08.004)\] 41. Arqub, O.A. Solutions of time-fractional Tricomi and Keldysh equations of Dirichlet functions types in Hilbert space. Numer. Methods Partial. Differ. Equ. **2018**, 34, 1759–1780. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solutions+of+time-fractional+Tricomi+and+Keldysh+equations+of+Dirichlet+functions+types+in+Hilbert+space&author=Arqub,+O.A.&publication_year=2018&journal=Numer.+Methods+Partial.+Differ.+Equ.&volume=34&pages=1759%E2%80%931780&doi=10.1002/num.22236)\] \[[CrossRef](https://doi.org/10.1002/num.22236)\] 42. Zhang, X.; Liu, J.; Wen, J.; Tang, B.; He, Y. Analysis for one-dimensional time-fractional Tricomi-type equations by LDG methods. Numer. Algorithms **2013**, 63, 143–164. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Analysis+for+one-dimensional+time-fractional+Tricomi-type+equations+by+LDG+methods&author=Zhang,+X.&author=Liu,+J.&author=Wen,+J.&author=Tang,+B.&author=He,+Y.&publication_year=2013&journal=Numer.+Algorithms&volume=63&pages=143%E2%80%93164&doi=10.1007/s11075-012-9617-3)\] \[[CrossRef](https://doi.org/10.1007/s11075-012-9617-3)\] 43. Frankl, F.L. On the Problems of Chaplygin for Mixed Sub- and Supersonic Flows. 1947. Available online: <https://ntrs.nasa.gov/search.jsp?R=20050019646> (accessed on 10 March 2026). 44. El-Ajou, A.; Oqielat, M.N.; Al-Zhour, Z.; Kumar, S.; Momani, S. Solitary solutions for time-fractional nonlinear dispersive PDEs in the sense of conformable fractional derivative. Chaos **2019**, 29, 093102. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solitary+solutions+for+time-fractional+nonlinear+dispersive+PDEs+in+the+sense+of+conformable+fractional+derivative&author=El-Ajou,+A.&author=Oqielat,+M.N.&author=Al-Zhour,+Z.&author=Kumar,+S.&author=Momani,+S.&publication_year=2019&journal=Chaos&volume=29&pages=093102&doi=10.1063/1.5100234)\] \[[CrossRef](https://doi.org/10.1063/1.5100234)\] | | | |---|---| | | **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. | Ā© 2026 by the authors. Licensee MDPI, Basel, Switzerland. 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/). ## Article Metrics Citations ### Article Access Statistics [Journal Statistics](https://www.mdpi.com/journal/mathematics/stats) Multiple requests from the same IP address are counted as one view. [Reciprocal Neural State–Disturbance Observer for Model-Free Trajectory Tracking of Robotic ManipulatorsPrevious](https://www.mdpi.com/2227-7390/14/6/983) [What Makes a Transformer Solve the TSP? A Component-Wise AnalysisNext](https://www.mdpi.com/2227-7390/14/6/985) ### Discover - [Articles](https://www.mdpi.com/search) - [Journals](https://www.mdpi.com/about/journals) - [Research Awards](https://www.mdpi.com/awards) - [Open Access Policy](https://www.mdpi.com/openaccess) ### Services - [Author Services](https://www.mdpi.com/authors/english) - [Conferences](https://www.mdpi.com/partnerships/conferences) - [Societies](https://www.mdpi.com/societies) - [Products & Solutions](https://www.mdpi.com/about/products) ### Guidelines - [For Authors](https://www.mdpi.com/authors) - [For Reviewers](https://www.mdpi.com/reviewers) - [For Editors](https://www.mdpi.com/editors) - [For Librarians](https://www.mdpi.com/librarians) ### Company - [About Us](https://www.mdpi.com/about) - [Careers](https://careers.mdpi.com/) - [Blog](https://blog.mdpi.com/) ### Support - [Pay an Invoice](https://www.mdpi.com/about/payment) - [Processing Charges](https://www.mdpi.com/apc) - [Discount Policy](https://www.mdpi.com/authors#Discounts_on_APCs) - [Contact Us](https://www.mdpi.com/support) [![MDPI](https://mdpi-res.com/cdn-cgi/image/w=56,h=30/https://mdpi-res.com/data/mdpi-logo-black.svg)![MDPI](https://mdpi-res.com/cdn-cgi/image/w=70,h=56/https://mdpi-res.com/data/mdpi-logo-black.svg)](https://www.mdpi.com/) [![Facebook](https://mdpi-res.com/data/facebook.svg)](https://www.facebook.com/MDPIOpenAccessPublishing)[![LinkedIn](https://mdpi-res.com/data/linkedin.svg)](https://www.linkedin.com/company/mdpi/)[![Twitter](https://mdpi-res.com/data/x.svg)](https://twitter.com/MDPIOpenAccess)[![YouTube](https://mdpi-res.com/data/youtube.svg)](https://www.youtube.com/@MDPIOA) Ā© 2026 MDPI AG (Basel, Switzerland) [Terms and Conditions](https://www.mdpi.com/about/terms-and-conditions)[Terms of Use](https://www.mdpi.com/about/termsofuse)[Privacy Policy](https://www.mdpi.com/about/privacy)[Accessibility](https://www.mdpi.com/accessibility) Privacy Settings Disclaimer
Readable Markdown
- Article - ![Open Access](https://mdpi-res.com/cdn-cgi/image/w=14,h=14/https://mdpi-res.com/data/open-access.svg) 13 March 2026 and School of Mathematics, Jilin University, Changchun 130012, China \* Author to whom correspondence should be addressed. ## Abstract This paper presents a comprehensive study on constructing exact and approximate solutions to Cauchy problems for time-fractional systems with variable coefficients. An innovative iterative approach is developed for solving functional equations with initial conditions, combining rigorous mathematical foundations with practical computational efficiency. The proposed technique effectively handles the nonlocal nature of fractional operators through a carefully designed iterative scheme that maintains simplicity while achieving high accuracy. It demonstrates particular strength in solving nonlinear systems with well-defined conditions and variable coefficients, where traditional methods often fail. Through systematic theoretical analysis and numerical validation, we establish the method’s convergence properties and computational advantages, showing its capability to generate both exact closed-form solutions, when available, and high-precision approximations otherwise. The approach remains computationally tractable even for complex cases where variable coefficients and memory effects of fractional systems present significant challenges to conventional solution approaches. MSC: 35R11; 34A12 ## 1\. Introduction This paper presents a systematic investigation of the initial value problems arising in time-fractional partial differential equations, where the nonlocal nature of fractional derivatives introduces unique mathematical challenges that fundamentally distinguish them from their integer-order counterparts. The inherent memory effects and historical dependence of fractional operators render the initial value problem particularly intricate yet fascinating from both theoretical and computational perspectives. Our comprehensive study develops and analyzes a suite of analytical and numerical methodologies specifically designed to overcome these challenges, with the dual objectives of establishing rigorous mathematical foundations and providing practical computational tools for understanding the temporal evolution of solutions in such systems. The research builds upon and significantly extends existing work in this rapidly developing field, as documented in \[1,2,3,4,5,6,7,8\] and related references, while introducing novel approaches to address the distinctive features of fractional-order initial value problems. The historical development of fractional calculus provides crucial context for understanding its modern applications. Our exploration begins with the seminal correspondence between L’HĆ“pital and Leibniz in 1695, a foundational moment in mathematical history documented in \[9\]. In this historic exchange, L’HĆ“pital posed the now-famous question to Leibniz: ā€œWhat would be the result of d n y d x n if n \= 1 2?ā€ Leibniz’s insightful reply demonstrated remarkable foresight: ā€œd 1 / 2 x will be equal to x d x : x. This is an apparent paradox, from which, one day useful consequences will be drawn.ā€ This early conceptual breakthrough, however, remained largely undeveloped for nearly two centuries. The mathematical community of the Enlightenment era faced several barriers to advancing fractional calculus: The lack of a rigorous theoretical framework for non-integer derivatives; The absence of clear physical interpretations; Significant computational challenges; The dominance of integer-order calculus in solving contemporary problems. This historical evolution underscores how theoretical mathematics often precedes practical applications, with Leibniz’s initial insight eventually finding numerous applications across science and engineering disciplines. The fundamental nature of fractional operators distinguishes them significantly from traditional differential operators. Researchers have established that these operators are intrinsically non-local, typically expressed through integral formulations with singular kernels. This essential characteristic gives fractional calculus its unique capacity to model memory effects and hereditary properties in physical systems, where the current state depends on the complete historical evolution rather than just instantaneous conditions. A particularly noteworthy application emerges in electrochemical systems, where fractional-order operators have proven remarkably effective: They provide superior modeling of diffusion processes in electrochemical cells compared to classical approaches; They more accurately characterize electrode impedance behavior; They offer enhanced precision in describing charge transfer mechanisms. The theoretical framework of fractional calculus has consequently found extensive utility across numerous scientific and technological domains, including: modeling anomalous diffusion and viscoelastic material properties, analyzing complex reaction-diffusion processes, understanding neuronal dynamics and biomechanical tissue behavior, characterizing long-range dependent processes in financial markets and so on. These diverse applications, meticulously documented in \[10\] and extensively developed in subsequent research, vividly illustrate the remarkable transformation of fractional calculus from a purely theoretical construct to an indispensable analytical framework in modern science and engineering. The continually expanding scope of its applications provides compelling evidence for the fundamental importance of these non-local operators, with demonstrated successes ranging from nanoscale phenomena to large-scale economic systems. Their unique capacity to capture long-range dependencies and memory effects has enabled breakthroughs in modeling complex systems that defy description by classical approaches, particularly in situations involving power-law behaviors and anomalous dynamics. Recent advances have particularly highlighted how fractional operators offer superior fidelity in modeling multiscale phenomena with coupled temporal and spatial dependencies, while also providing natural frameworks for systems exhibiting complex hierarchical or fractal-like structures. The operational success of fractional calculus across such disparate domains not only underscores its profound theoretical significance but has also driven important innovations in computational methods and analytical techniques. This dual impact—simultaneously deepening our fundamental understanding while enhancing practical problem-solving capabilities—continues to motivate vigorous research efforts at the frontiers of applied mathematics and interdisciplinary science, promising further discoveries and applications in emerging fields. Fractional calculus, representing a significant generalization of classical differential and integral operators to non-integer orders, has emerged as a vibrant research field attracting considerable attention from both theoretical and applied perspectives. This mathematical framework extends the traditional notions of differentiation and integration to fractional orders, thereby enabling the description of complex phenomena that exhibit non-local behavior and memory effects. Throughout its historical development, numerous distinct definitions of fractional derivatives and integrals have been proposed, each offering unique advantages for specific applications and mathematical contexts. The absence of a unified definition for fractional derivatives has led to a rich diversity of operational formulations, including several prominent variants that have proven particularly useful in practice. Among these, the Grünwald–Letnikov approach provides a discrete interpretation based on finite differences, while the Riemann–Liouville formulation offers a continuous integral-based perspective. The Caputo derivative has gained prominence for its advantageous treatment of initial conditions in differential equations. Additional important formulations include the Riesz potential for symmetric fractional operations, the conformable fractional derivative for maintaining certain classical properties, and various local fractional derivatives designed for specific applications. Within this diverse landscape, the Riemann–Liouville fractional integral and Caputo fractional derivative have emerged as particularly influential, owing to their mathematical tractability and physical interpretability. These operators have found extensive application across numerous scientific and engineering disciplines, including advanced control systems where they enable more robust controller designs, the modeling of anomalous diffusion processes in complex media, and the analysis of fractal structures and multi-scale systems. The continued development and refinement of these fractional operators, as documented in \[4,7,8,11,12,13,14\] and related works, reflects their growing importance in addressing challenging problems that lie beyond the scope of classical calculus. Among the various formulations of fractional operators, the Grünwald–Letnikov fractional derivative (see [Appendix A.2](https://www.mdpi.com/2227-7390/14/6/984#secAdot2-mathematics-14-00984)) occupies a special historical and conceptual position as one of the most intuitive and computationally accessible approaches. This operator, rooted in the generalization of finite difference methods, provides a natural discrete framework for fractional differentiation that closely parallels the classical limit definition of integer-order derivatives. The Grünwald–Letnikov approach is particularly valuable in numerical analysis and scientific computing, where its discrete nature allows for efficient implementation in algorithmic solutions to fractional differential equations. Its computational advantages have made it indispensable for simulations of complex systems with memory effects, as extensively documented in \[15\] and subsequent computational studies. In contrast, the Riemann–Liouville fractional derivative represents a continuous formulation defined through a convolution operation with a power-law kernel. This definition enjoys several theoretical advantages, including a rigorous mathematical foundation in functional analysis and a natural connection to the theory of integral equations. The Riemann–Liouville operator has proven particularly powerful in analytical approaches to fractional differential equations arising in mathematical physics, where its properties facilitate the development of existence theorems and solution representations. The complementary strengths of these two approaches—with Grünwald–Letnikov favoring computational implementations and Riemann–Liouville supporting theoretical analysis—serve to illustrate the rich diversity of tools available in fractional calculus, as discussed in detail in \[16\] and related mathematical physics literature. In contrast, the Caputo fractional derivative was introduced to address this issue, offering a more physically intuitive way to incorporate initial conditions in models. It is particularly popular in applied settings such as control systems. The Riesz potential, often employed in potential theory and the fractional Laplacian, generalizes the idea of a potential field to fractional dimensions and is crucial in studying processes like anomalous diffusion and LĆ©vy flights; see \[17\] and the references therein. In contrast to the Riemann–Liouville formulation, the Caputo fractional derivative was specifically developed to overcome certain mathematical challenges in modeling physical systems. This alternative definition provides several crucial advantages for practical applications, most notably its ability to handle initial conditions in a manner that aligns more naturally with physical interpretations. The Caputo operator’s formulation ensures that initial value problems can be specified using standard function values rather than fractional derivatives of the initial state, making it particularly valuable for modeling real-world phenomena. This property has led to its widespread adoption in applied mathematics and engineering disciplines, especially in control theory, where it enables more natural formulations of dynamical systems with memory effects. The Caputo derivative’s compatibility with conventional physical interpretations has made it a preferred choice for modeling viscoelastic materials, thermal systems, and other processes where the initial state is typically specified in terms of classical variables. Similarly important in fractional analysis is the Riesz potential, which extends classical potential theory to fractional dimensions through an elegant symmetric formulation. This operator plays a fundamental role in the theory of the fractional Laplacians and has become indispensable for studying non-local phenomena. The Riesz transform provides a powerful mathematical framework for analyzing complex diffusion processes, particularly those exhibiting anomalous behavior that deviates from classical Brownian motion. Its applications span diverse areas, including the modeling of LĆ©vy flights in statistical physics, the analysis of non-local interactions in quantum mechanics, and the study of long-range correlations in complex systems. The theoretical foundations and practical applications of the Riesz potential are thoroughly examined in \[17\], along with its connections to other fractional operators and its growing importance in modern mathematical physics. The rich diversity of fractional derivative definitions, while reflecting the wide-ranging applications of fractional calculus across scientific disciplines, presents significant challenges for researchers working to select, compare, and apply these various operators. This proliferation of definitions—each with distinct properties and suited for different applications—can create confusion and complicate theoretical comparisons between results obtained using different formulations. Recognizing this challenge, the mathematical community has increasingly turned to axiomatization approaches that seek to unify these diverse operators within a coherent theoretical framework. In \[12\], the authors adopt such an axiomatic perspective, examining several fractional derivative definitions in a systematic way. By identifying the core properties that these operators share, they developed a general framework that facilitates a deeper understanding of their similarities and differences. For example, the axioms might include the following several items: - (Linearity) The operator is linear. - (Identity) The zero order derivative of a function returns the function itself. - (Backward compatibility) When the order is integer, the fractional derivative gives the same result as the classical derivative. - (Index law) For any α \< 0 and β \< 0, we have D α D β f ( t ) \= D α \+ β f ( t ) . - (Generalized Leibniz rule) D α \[ f ( t ) g ( t ) \] \= āˆ‘ i \= 0 āˆž α i D i f ( t ) D α āˆ’ i g ( t ) . **Remark** **1\.** If α ∈ N , this is exactly the classical Leibniz rule. By the way, another generalization of this rule, such as D α \[ f ( t ) g ( t ) \] \= āˆ‘ i \= āˆ’ i n f t y āˆž α i \+ β D β \+ i f ( t ) D α āˆ’ β āˆ’ i g ( t ) . This approach not only highlights the underlying mathematical structure common to various fractional operators but also provides a way to compare and contrast their behaviors. By formalizing fractional calculus through axiomatization, it becomes easier to extend these operators to new contexts, such as multi-dimensional systems, stochastic processes, or even quantum mechanics, where fractional derivatives are increasingly finding applications; see \[5,18,19,20\] and the references therein. The axiomatization of fractional operators has profound implications for both the theoretical development of fractional calculus and its practical applications. On the theoretical side, it provides a clearer framework for exploring the connections between different types of fractional derivatives and integrals, potentially leading to the discovery of new operators or generalizations that encompass existing ones. It also helps to clarify which definitions are most suitable for specific problems, depending on the required properties of the derivative, such as memory effects, non-locality, or boundary behavior; see \[21,22\] and the references therein. From a practical standpoint, this unifying approach simplifies the application of fractional derivatives in diverse fields, from engineering to biology. For instance, in modeling materials with memory effects, the choice of a fractional derivative is crucial, and an axiomatic perspective can help guide that choice. Similarly, in complex systems where fractal geometry or anomalous diffusion is present, such as in finance, geology, or network theory, an axiomatized framework for fractional derivatives ensures that the correct operator is used to capture the necessary dynamics. In this paper, we mainly explore the construction of exact solutions or approximate solutions to the initial value problems of time fractional partial differential equations. This is a very direct problem in the study of the theory of differential equations. However, it is not easy and even rather difficult. In addition, solving fractional systems is usually more difficult than classical systems, for their operator is defined by an integral (a non-local operator with a singular kernel). Luckily, there are some different effective methods that have been developed to construct the exact solutions or the approximate solutions of fractional systems with initial conditions. Such as the iterative method \[23,24,25,26,27\], the homotopy analysis method \[28,29,30\], the Cole-Hopf transformation method \[31\], the residual power series method \[32,33,34,35,36\], the Laplace transform method \[37\], and so on. ### Outline of Paper This manuscript presents a systematic investigation of iterative techniques for solving initial value problems in time-fractional partial differential equations. The paper is organized as follows: In [Section 2](https://www.mdpi.com/2227-7390/14/6/984#sec2-mathematics-14-00984), we establish the fundamental mathematical framework by introducing essential concepts from fractional calculus, including the various definitions of fractional derivatives and their key properties. This section provides the necessary theoretical foundation by presenting important lemmas regarding fractional operators, particularly focusing on their non-local characteristics and memory effects. We carefully define all mathematical notations to ensure clarity throughout subsequent developments, paying special attention to the function spaces and operator norms relevant to our analysis of fractional differential equations. [Section 3](https://www.mdpi.com/2227-7390/14/6/984#sec3-mathematics-14-00984) develops an iterative approach for solving functional equations with initial conditions, with particular emphasis on applications to fractional-order systems. We present a comprehensive convergence analysis of the proposed iterative scheme, establishing sufficient conditions under which the method guarantees convergence to the true solution. The section includes detailed proofs of convergence theorems and examines the rate of convergence under various regularity assumptions. Furthermore, we discuss the computational advantages of this approach compared to traditional methods, highlighting its efficiency in handling the non-local nature of fractional operators. [Section 4](https://www.mdpi.com/2227-7390/14/6/984#sec4-mathematics-14-00984) demonstrates the practical implementation of our iterative technique through an in-depth study of the one-dimensional time-fractional Tricomi equation. We derive both exact solutions (where obtainable) and high-precision approximate solutions, comparing our results with known analytical solutions to validate the method’s accuracy. The section includes detailed case studies with different fractional orders and initial conditions, accompanied by rigorous error analysis and convergence rate verification. Special attention is given to the physical interpretation of solutions and their relationship to the fractional order parameter. The analysis is extended in [Section 5](https://www.mdpi.com/2227-7390/14/6/984#sec5-mathematics-14-00984) to multidimensional time-fractional Tricomi equations, where we address the additional complexities introduced by higher spatial dimensions. We develop dimensionally-adaptive versions of our iterative scheme and examine their performance for various initial boundary value problems. The section presents numerical results for two- and three-dimensional cases, including visualization of solution surfaces and comparison with alternative numerical methods. We particularly focus on the method’s scalability and its ability to handle different types of initial conditions, from smooth analytic functions to discontinuous data. Finally, [Section 6](https://www.mdpi.com/2227-7390/14/6/984#sec6-mathematics-14-00984) synthesizes our findings into a comprehensive conclusion, summarizing the key advantages and limitations of the proposed iterative technique. We discuss broader implications for the field of fractional calculus and suggest several directions for future research, including potential applications to nonlinear fractional systems and coupled fractional-order equations. The section concludes with practical recommendations for implementing the method and its possible extensions to other classes of fractional differential equations. ## 2\. Fractional Calculus This subsection introduces the fundamental theoretical framework necessary for our investigation of fractional systems. We begin with the essential special functions that form the backbone of fractional calculus. The Gamma function, which extends the factorial concept to continuous values, plays an indispensable role in defining fractional operators due to its appearance in their integral kernels. Equally crucial is the Mittag–Leffler function, a sophisticated generalization of the exponential function that naturally describes processes with memory effects and power-law relaxation. The mathematical landscape of fractional calculus features a diverse collection of definitions for fractional integrals and derivatives, each developed to serve specific theoretical or applied purposes. Among these various formulations, several have emerged as particularly significant. The Riemann–Liouville fractional integral represents the most direct extension of classical integration theory to non-integer orders, while the Caputo fractional derivative has become preferred for physical modeling due to its more intuitive treatment of initial conditions. Other important variants include the Riesz potential, which provides symmetric fractional operations valuable for certain boundary value problems, and the Weyl derivative, especially useful for analyzing periodic phenomena. The comparative advantages and specific applications of these different definitions are thoroughly examined in the comprehensive references \[7,8,11,12,13,14\]. Within this spectrum of fractional operators, the Riemann–Liouville integral and Caputo derivative have gained particular prominence in both theoretical and applied contexts. The Riemann–Liouville approach provides the most natural mathematical foundation for fractional integration, while the Caputo formulation offers distinct practical advantages for solving initial value problems that model physical systems. This preference stems from the Caputo derivative’s compatibility with conventional initial conditions specified in terms of standard derivatives, making it more suitable for applications in engineering and physics. Our subsequent analysis will build upon several key lemmas concerning the operational properties of these fractional operators and the characteristic behaviors of Mittag–Leffler functions. These fundamental results enable rigorous investigation of existence, uniqueness, and qualitative properties for solutions to fractional differential equations. The interplay between these various concepts creates a robust framework for examining the nonlinear fractional systems that constitute the primary focus of this work, while also providing the necessary tools for developing effective solution methods. ### 2\.1. The Gamma Function and the Mittag–Leffler Function This subsection introduces the definitions for the Gamma function and the Mittag–Leffler function, which will be used in the definitions of the fractional integral and derivative. **Definition** **1** (\[38\])**.** The Gamma function is defined by Ī“ ( x ) \= ∫ 0 āˆž t x āˆ’ 1 e āˆ’ t d t Re ( x ) \> 0 , where Re( x ) denotes the real part of the complex number x. **Proposition** **1** (\[38\])**.** Some main properties of the Gamma function: - Ī“ ( x \+ 1 ) \= x Ī“ ( x ) for all x with Rex \> 0 . - Ī“ ( 1 ) \= 1 . - Ī“ ( 1 2 ) \= Ļ€ . **Definition** **2** (\[39\])**.** The two-parameters Mittag–Leffler function is defined by letting E α , β ( x ) \= āˆ‘ k \= 0 āˆž x k Ī“ ( α k \+ β ) α , β \> 0 . And the one-parameter Mittag–Leffler function is defined by letting E α ( x ) \= āˆ‘ k \= 0 āˆž x k Ī“ ( α k \+ 1 ) α \> 0 . **Remark** **2\.** The two-parameter Mittag–Leffler function is a generalization of the one-parameter Mittag–Leffler function, that is, E α , 1 ( x ) \= E α ( x ) . The one-parameter Mittag–Leffler function is a generalization of the exponential function, or rather, E 1 ( x ) \= e x . ### 2\.2. The Riemann–Liouville Fractional Integral and the Caputo Derivative This subsection recalls the definition of the Riemann–Liouville fractional integral and Caputon derivative, see \[4,5,6,13,40\] and the references therein. **Definition** **3\.** Let Īø : \[ 0 , T \] → R is a continuous function. Then, the fractional-order integral of Īø in the Riemann–Liouville sense is defined by J α Īø ( t ) \= 1 Ī“ ( α ) ∫ 0 t ( t āˆ’ s ) α āˆ’ 1 Īø ( s ) d s , α \> 0 , t ≄ 0 , (1) where Ī“ ( α ) \= ∫ 0 āˆž e āˆ’ u u α āˆ’ 1 d u . **Definition** **4\.** Let Īø : \[ 0 , T \] → R is a continuous function. Then, the Caputo fractional-order derivative of Īø is given by D α c Īø ( t ) \= 1 Ī“ ( n āˆ’ α ) ∫ 0 t ( t āˆ’ s ) n āˆ’ α āˆ’ 1 Īø ( m ) ( s ) d s , where m āˆ’ 1 \< α ≤ m ; m ∈ N . **Lemma** **1** (\[27\])**.** For any m āˆ’ 1 \< β ≤ m ∈ N āˆ— , one has I β D β u ( x , t ) \= u ( x , t ) āˆ’ āˆ‘ k \= 0 m āˆ’ 1 āˆ‚ k u ( x , 0 ) t k k \! . **Lemma** **2\.** For any β , γ , one has I β t γ \= Ī“ ( γ \+ 1 ) Ī“ ( γ \+ β \+ 1 ) t γ \+ β . ## 3\. An Iterative Technique for Solving a Class of Nonlinear Functional Equations In this section, we present a fundamental lemma that plays a pivotal role in our subsequent mathematical analysis. It provides the essential theoretical foundation for establishing the key propositions and theorems developed throughout this paper. It has been chosen for its particular relevance to our investigation of fractional differential systems, offering crucial insights into the behavior of solutions and the properties of fractional operators. Its proper application will be carefully demonstrated in subsequent sections where we develop our main results. **Lemma** **3\.** Consider a class of abstract functional equations of the following form u ( x , t ) \= f ( x , t ) \+ L u ( x , t ) \+ N u ( x , t ) , (2) where ( t , x ) ∈ D \= { ( t , x ) : t ∈ R , x ∈ R n , n ∈ N } , u ( x , t ) is the unknown function, f is a given known function, L is a linear operator and N is a nonlinear operator with respect to u from a Banach space Dom ( L ) ∩ Dom ( N ) to itself (Dom ( L ) ∩ Dom ( N ) stands for the domain of L and N respectively). For abstract function Equation ([2](https://www.mdpi.com/2227-7390/14/6/984#FD2-mathematics-14-00984)), consider the series solutions of the following form u ( t , x ) \= āˆ‘ n \= 0 āˆž u n ( t , x ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the iterative sequence { u n } can be obtained by the following recurrence equations u 0 \= f , (3) u 1 \= L ( u 0 ) \+ N ( u 0 ) , (4) u m \+ 1 \= L ( u m ) \+ N āˆ‘ n \= 0 m u n āˆ’ N āˆ‘ n \= 0 m āˆ’ 1 u n , m \= 1 , 2 ⋯ . (5) Furthermore, if the operators L and N are contractive, then the series solution āˆ‘ n \= 0 āˆž u n ( t , x ) converges absolutely and uniformly. **Proof.** Obviously, the nonlinear operator N can be decomposed as N ( u ( t , x ) ) \= N āˆ‘ n \= 0 āˆž u n ( t , x ) \= N ( u 0 ) \+ āˆ‘ n \= 1 āˆž N āˆ‘ j \= 0 n u j āˆ’ N āˆ‘ j \= 0 n āˆ’ 1 u j . Similarly, the linear operator L can also be decomposed as L ( u ( t , x ) ) \= L āˆ‘ n \= 0 āˆž u n ( t , x ) \= L ( u 0 ) \+ āˆ‘ n \= 1 āˆž L āˆ‘ j \= 0 n u j āˆ’ L āˆ‘ j \= 0 n āˆ’ 1 u j . Denote by M u ( t , x ) \= L u ( t , x ) \+ N u ( t , x ) . then one has M ( u ( t , x ) ) \= M ( u 0 ) \+ āˆ‘ n \= 1 āˆž M āˆ‘ j \= 0 n u j āˆ’ M āˆ‘ j \= 0 n āˆ’ 1 u j . Define the following recurrence equations u 0 \= f , u 1 \= M ( u 0 ) , u m \+ 1 \= M āˆ‘ n \= 0 m u n āˆ’ M āˆ‘ n \= 0 m āˆ’ 1 u n , m \= 1 , 2 ⋯ . Since operators L and N are contractive, then M is also contractive, i.e., there exists a constant 0 \< K \< 1, such that ∄ M ( v i ) āˆ’ M ( v j ) ∄ ≤ K ∄ v i āˆ’ v j ∄ , āˆ€ v i , v j ∈ Dom ( L ) ∩ D o m ( N ) , where ∄ Ā· ∄ denotes the norm on Banach space Dom( L ) ∩Dom( N ). Furthermore, for u m \+ 1, one can obtain ∄ u m \+ 1 ∄ \= M āˆ‘ n \= 0 m u n āˆ’ M āˆ‘ n \= 0 m āˆ’ 1 u n ≤ K ∄ u m ∄ ≤ K m \+ 1 ∄ u 0 ∄ . Since 0 \< K \< 1, the series āˆ‘ m \= 1 āˆž K m \+ 1 ‖ u 0 ‖ converges absolutely as well as uniformly. According to the Weierstrass M-test, one can obtain that the series āˆ‘ i \= 0 āˆž u i converges absolutely as well as uniformly. ā–” Some more specific details about this type of iterative method could be found in \[23,24,25,26,27\] and the references therein. ## 4\. One-Dimensional Time Fractional Tricomi Equation with Variable Coefficients Consider the following one-dimensional time-fractional Tricomi equation of order α, as studied in \[41\]: D t α u \+ ξ ( t ) āˆ‚ x 2 u \= f ( x , t ) , (6) subject to the initial conditions u ( x , 0 ) \= ω 1 ( x ) , D t u ( x , 0 ) \= ω 2 ( x ) , (7) where 1 \< α ≤ 2 represents the fractional order of the time derivative, t ∈ \[ 0 , T \] denotes the temporal variable, and x ∈ R is the spatial variable. The special case when α \= 2 reduces ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) to the classical Tricomi equation, first introduced by Francesco Tricomi in his seminal 1923 work \[42\]. This equation occupies a particularly important position in the theory of partial differential equations due to its remarkable property of exhibiting mixed-type characteristics. Specifically, the nature of the equation undergoes a fundamental change depending on the sign of the coefficient function ξ ( t ): when ξ ( t ) \> 0, the equation demonstrates elliptic behavior, while for ξ ( t ) \< 0, it becomes hyperbolic. This dual character makes the Tricomi equation a paradigmatic example of mixed-type partial differential equations. The physical significance of this equation was further elucidated by Frankl \[43\], who established its crucial connection to transonic gas dynamics. More precisely, the Tricomi equation provides a mathematical model describing the transition between subsonic and supersonic flow regimes in gas dynamics. The elliptic region (ξ ( t ) \> 0) corresponds to subsonic flow, where disturbances propagate in all directions, while the hyperbolic region (ξ ( t ) \< 0) represents supersonic flow, characterized by directional wave propagation and the formation of shock waves. This interpretation underscores the equation’s importance in aerodynamics and fluid mechanics. In the following analysis, we will develop an iterative solution technique to solve the initial value problem ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) and ([7](https://www.mdpi.com/2227-7390/14/6/984#FD7-mathematics-14-00984)) for the time-fractional generalization of this important equation. Our approach will carefully account for both the fractional temporal derivative and the mixed-type spatial operator, while maintaining the essential physical interpretation of the solutions. By integrating both sides of Equation ([6](https://www.mdpi.com/2227-7390/14/6/984#FD6-mathematics-14-00984)) with respect to the temporal variable from 0 to t, we derive the equivalent integral formulation: u ( x , t ) \= u ( x , 0 ) \+ D t u ( x , 0 ) t āˆ’ I t α ξ ( t ) āˆ‚ x 2 u \+ I t α f ( x , t ) . (8) Incorporating the initial conditions ([7](https://www.mdpi.com/2227-7390/14/6/984#FD7-mathematics-14-00984)) into the integral Equation ([8](https://www.mdpi.com/2227-7390/14/6/984#FD8-mathematics-14-00984)) yields: u ( x , t ) \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) āˆ’ I t α ξ ( t ) āˆ‚ x 2 u . From the perspective of iterative methods, we can interpret this equation as follows: The terms ω 1 ( x ) \+ ω 2 ( x ) \+ I t α f ( x , t ) constitute the initial iterative data, incorporating both the given initial conditions and the source term f ( x , t ). The remaining term āˆ’ I t α ξ ( t ) āˆ‚ x 2 u acts as a nonlinear feedback operator that propagates the solution through successive iterations. We propose a series solution of the form: u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the iterative sequence is determined by the following recursive scheme: u 0 \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) , u 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 , u 2 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 , ā‹® u m \+ 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u m . To demonstrate the effectiveness of this iterative approach, we now examine several specific cases with carefully chosen initial conditions, source terms, and coefficient functions. **Proposition** **2** (Trivial Solution Case)**.** For the simplest case where ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , and f ( x , t ) ≔ 0 , we observe that for any well-behaved coefficient function ξ ( t ) , the solution identically vanishes: u ( x , t ) ≔ 0 . This trivial solution serves as an important consistency check for our method. **Proposition** **3** (Exponential Source Term with Constant Coefficient)**.** Consider the case with zero initial conditions ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , an exponential source term f ( x , t ) ≔ e x , and constant coefficient ξ ( t ) ≔ āˆ’ 1 . The iterative scheme yields: u n \= e x t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) āˆ€ n ∈ N . The complete solution then takes the elegant form: u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x āˆ‘ n \= 0 āˆž t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) \= e x E α ( t α ) āˆ’ 1 , where E α ( Ā· ) denotes the Mittag–Leffler function, which naturally emerges in solutions to fractional differential equations. **Proposition** **4** (Modified Exponential Source Term)**.** For the case with ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x t α Ī“ ( āˆ’ α \+ 1 ) , and ξ ( t ) ≔ āˆ’ 1 , we obtain: u n \= e x t n α Ī“ ( n α \+ 1 ) āˆ€ n ∈ N . The solution series converges to: u ( x , t ) \= e x E α ( t α ) , demonstrating how different source terms can lead to distinct Mittag–Leffler function representations. **Proposition** **5** (Non-constant Coefficient Case)**.** Consider the problem with ω 1 ( x ) \= e x , ω 2 ( x ) \= 0 , f ( x , t ) \= 0 , ξ ( t ) \= t α . Then the iterative scheme generates the sequence u 0 \= e x , u 1 \= āˆ’ e x Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α , u 2 \= e x Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . By the induction of mathematics, we establish the general term: u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α , n ≄ 0 . The complete solution is expressed as the infinite series: u ( x , t ) \= e x āˆ‘ n \= 0 āˆž ( āˆ’ 1 ) n āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α , revealing the intricate relationship between the solution’s temporal behavior and the fractional order α. **Proof.** **Step 1: Iterative structure.** Under the given assumptions, the iteration takes the form u n \+ 1 ( x , t ) \= āˆ’ I t α t α u n ( x , t ) , where I t α denotes the fractional integral operator. Because the initial data is separable and u 0 \= e x, we observe that each iterate preserves the spatial factor e x. Hence we may write u n ( x , t ) \= e x a n ( t ) , reducing the problem to determining the scalar sequence a n ( t ). **Step 2: First iterates.** For n \= 0, a 0 ( t ) \= 1 . Using the fractional integral formula I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain a 1 ( t ) \= āˆ’ Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α . Repeating the same computation yields a 2 ( t ) \= Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . **Step 3: Inductive hypothesis.** Assume that for some n ≄ 1, u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t 2 n α . **Step 4: Inductive step.** Applying the iteration formula, u n \+ 1 \= āˆ’ I t α t α u n . Substituting the inductive expression gives t α u n \= ( āˆ’ 1 ) n e x āˆ k \= 1 n Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n Ī“ ( 2 k α \+ 1 ) t ( 2 n \+ 1 ) α . Applying the fractional integral formula once more, I t α t ( 2 n \+ 1 ) α \= Ī“ ( ( 2 n \+ 1 ) α \+ 1 ) Ī“ ( ( 2 n \+ 2 ) α \+ 1 ) t ( 2 n \+ 2 ) α . Therefore, u n \+ 1 \= ( āˆ’ 1 ) n \+ 1 e x āˆ k \= 1 n \+ 1 Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) āˆ k \= 1 n \+ 1 Ī“ ( 2 k α \+ 1 ) t 2 ( n \+ 1 ) α . This completes the induction. **Step 5: Series representation.** Summing the iterates yields the solution as the convergent power-type series u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) , which gives the stated expansion. ā–” ## 5\. Higher-Dimensional Time Fractional Tricomi Equation with Variable Coefficients Consider the following higher-dimensional time fractional Tricomi equation D t α u \+ ξ ( t ) Ī” x u \= f ( x , t ) (9) with initial conditions u ( x , 0 ) \= ω 1 ( x ) , D t u ( x , 0 ) \= ω 2 ( x ) . (10) where 1 \< α ≤ 2 , t ∈ \[ 0 , T \] , x ∈ R n, and Ī” x denotes the standard Laplace operator and it is defined by Ī” x u ( x , t ) \= āˆ‘ i \= 1 n āˆ‚ x i 2 u ( x , t ). Similarly, one can obtain the equivalent integral Equations ([9](https://www.mdpi.com/2227-7390/14/6/984#FD9-mathematics-14-00984)) and ([10](https://www.mdpi.com/2227-7390/14/6/984#FD10-mathematics-14-00984)) as follows u ( x , t ) \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) āˆ’ I t α ξ ( t ) Ī” x u ( x , t ) . Suppose that the unknown function u has the series form u ( x , t ) \= āˆ‘ n \= 0 āˆž u n ( x , t ) \= u 0 \+ u 1 \+ u 2 \+ ⋯ , where the sequence { u n } can be obtained by the following iterative process u 0 \= ω 1 ( x ) \+ ω 2 ( x ) t \+ I t α f ( x , t ) u 1 \= āˆ’ I t α ξ ( t ) Ī” x u 0 u 2 \= āˆ’ I t α ξ ( t ) Ī” x u 1 ā‹® u m \+ 1 \= āˆ’ I t α ξ ( t ) Ī” x u m . **Proposition** **6\.** If ω 1 ( x ) ≔ 0 , ω 2 ≔ 0 , f ( x , t ) ≔ 0 , then for any reasonable ξ ( t ) , u ( x , t ) ≔ 0 . **Proposition** **7\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 \+ x 2 \+ ⋯ \+ x n and ξ ( t ) ≔ āˆ’ 1 . - **Iterative construction.** Since the initial data vanish, the first term of the iteration is obtained directly from the source term: u 0 \= e x 1 \+ ⋯ \+ x n t α Ī“ ( α \+ 1 ) . Applying the iterative scheme once more and using I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain u 1 \= e x 1 \+ ⋯ \+ x n t 2 α Ī“ ( 2 α \+ 1 ) . Proceeding inductively, one easily verifies that u n \= e x 1 \+ ⋯ \+ x n t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing all iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t ( n \+ 1 ) α Ī“ ( ( n \+ 1 ) α \+ 1 ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 1 āˆž t n α Ī“ ( n α \+ 1 ) . Rewriting the series, u ( x , t ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) āˆ’ 1 . Hence, u ( x , t ) \= e x 1 \+ ⋯ \+ x n E α ( t α ) āˆ’ 1 . **Proposition** **8\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ e x 1 \+ x 2 \+ ⋯ \+ x n t α Ī“ ( 1 āˆ’ α ) , ξ ( t ) ≔ āˆ’ 1 . **Iterative construction.** Since the initial data vanish and the source term already contains t α , the first iterate reduces to u 0 \= e x 1 \+ ⋯ \+ x n . Applying the iterative scheme once, u 1 \= e x 1 \+ ⋯ \+ x n t α Ī“ ( α \+ 1 ) . Continuing in the same manner, one verifies inductively that u n \= e x 1 \+ ⋯ \+ x n t n α Ī“ ( n α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing the iterates gives u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) . Recognizing the Mittag–Leffler series, E α ( z ) \= āˆ‘ n \= 0 āˆž z n Ī“ ( n α \+ 1 ) , we conclude that u ( x , t ) \= e x 1 \+ ⋯ \+ x n E α ( t α ) . **Proposition** **9\.** Assume ω 1 ( x ) ≔ 0 , ω 2 ( x ) ≔ e x , f ( x , t ) ≔ 0 , and ξ ( t ) ≔ āˆ’ t α Ī“ ( α \+ 1 ) . **Iterative construction.** Since the source term vanishes and the initial velocity is e x , the first term of the iteration is u 0 \= e x . Proceeding with the iterative scheme, one obtains successively terms of increasing order in t α . By induction it follows that u n \= e x t n α Ī“ ( n α \+ 1 ) , āˆ€ n ∈ N . **Series solution.** Summing the iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n \= e x āˆ‘ n \= 0 āˆž t n α Ī“ ( n α \+ 1 ) . Recognizing the Mittag–Leffler expansion, we conclude that u ( x , t ) \= e x E α ( t α ) . **Proposition** **10\.** Assume ω 1 ( x ) ≔ e x 1 \+ ⋯ \+ x n , ω 2 ( x ) ≔ 0 , f ( x , t ) ≔ 0 , ξ ( t ) ≔ t α . **First iterates.** The initial term is u 0 \= e x 1 \+ ⋯ \+ x n . Since āˆ‚ x 2 e x 1 \+ ⋯ \+ x n \= e x 1 \+ ⋯ \+ x n , the first iteration gives u 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 0 \= āˆ’ e x 1 \+ ⋯ \+ x n I t α t α \= āˆ’ e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) t 2 α . Applying the scheme once more, u 2 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u 1 \= e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 2 α \+ 1 ) I t α t 3 α \= e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) t 4 α . **Inductive step.** Assume that for some k ≄ 1 , u k \= ( āˆ’ 1 ) k \+ 1 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 k α \+ 1 ) t 2 k α . Then u k \+ 1 \= āˆ’ I t α ξ ( t ) āˆ‚ x 2 u k \= ( āˆ’ 1 ) k \+ 2 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) ⋯ Ī“ ( ( 2 k āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) ⋯ Ī“ ( 2 k α \+ 1 ) I t α t ( 2 k \+ 1 ) α . Using I t α t β \= Ī“ ( β \+ 1 ) Ī“ ( β \+ α \+ 1 ) t β \+ α , we obtain u k \+ 1 \= ( āˆ’ 1 ) k \+ 2 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) ⋯ Ī“ ( ( 2 k \+ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) ⋯ Ī“ ( ( 2 k \+ 2 ) α \+ 1 ) t ( 2 k \+ 2 ) α . Hence, by induction, for all n ∈ N , u n \= ( āˆ’ 1 ) n \+ 1 e x 1 \+ ⋯ \+ x n Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 n α \+ 1 ) t 2 n α . **Series representation.** Summing the iterates, u ( x , t ) \= āˆ‘ n \= 0 āˆž u n , that is, u ( x , t ) \= e x 1 \+ ⋯ \+ x n āˆ‘ n \= 0 āˆž ( āˆ’ 1 ) n \+ 1 Ī“ ( α \+ 1 ) Ī“ ( 3 α \+ 1 ) ⋯ Ī“ ( ( 2 n āˆ’ 1 ) α \+ 1 ) Ī“ ( 2 α \+ 1 ) Ī“ ( 4 α \+ 1 ) ⋯ Ī“ ( 2 n α \+ 1 ) t 2 n α . ## Author Contributions Methodology, W.L. and J.H.; Validation, J.H.; Writing—original draft, J.H.; Writing—review and editing, W.L.; Visualization, W.L. and J.H. All authors have read and agreed to the published version of the manuscript. ## Funding This research received no external funding. ## Data Availability Statement No new data were created or analyzed in this study. ## Conflicts of Interest The authors declare no conflicts of interest. ## Appendix A ### Appendix A.1. The Conformable Fractional Derivative This subsection recalls the definition of the conformable fractional derivative, which looks like the limit of the difference quotient (the classical Newtonian derivative). The α\-th conformable fractional derivative of a function f is defined by letting ( T α f ) ( t ) \= lim ϵ → 0 f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t ) ϵ . Here 0 ≤ α ≤ 1. It is clear that the case of α \= 1 coincides with the classical Newtonian derivative. Furthermore, one can verify that - The zero order derivative of a function does not return the function itself. In fact, we have ( T 0 f ) ( t ) \= lim ϵ → 0 f ( t \+ ϵ t ) āˆ’ f ( t ) ϵ \= lim ϵ → 0 f ( ( 1 \+ ϵ ) t ) āˆ’ f ( t ) 1 \+ ϵ Ā· lim ϵ → 0 1 \+ ϵ ϵ \= ( D f ) ( t ) Ā· āˆž ≠ f ( t ) , where D denotes the quantum derivative. - The index law does not hold, that is, ( T α T β f ) ( t ) ≠ ( T α \+ β f ) ( t ). In fact, we have LHS \= lim ϵ → 0 1 ϵ 2 \[ f ( t \+ ϵ t 1 āˆ’ α \+ ϵ ( t \+ ϵ t 1 āˆ’ α ) 1 āˆ’ β ) āˆ’ f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t \+ ϵ t 1 āˆ’ α ) āˆ’ f ( t \+ ϵ t 1 āˆ’ β ) \+ f ( t ) \] and RHS \= lim ϵ → 0 f ( t \+ ϵ t 1 āˆ’ ( α \+ β ) ) āˆ’ f ( t ) ϵ . It is not difficult to see that ( T α T β f ) ( t ) ≠ ( T α \+ β f ) ( t ) without any restriction. - In this perspective of the criteria, the conformable fractional derivative is not a fractional derivative; see \[13\]. By the way, regarding the definition of the higher order conformable fractional derivative, please refer to \[44\] (Definition 2.1) and the references therein. ### Appendix A.2. The Grünwald–Letnikov Derivative This subsection recalls the definition of the Grünwald–Letnikov derivative. It is a type of fractional derivative that extends the concept of differentiation to non-integer orders. It is defined using a limit process similar to the traditional definition of the derivative but incorporates fractional differences. Specifically, it is expressed as a limit of a sum involving the values of a function at discrete points, weighted by binomial coefficients. This derivative is useful in various applications, particularly in fields like physics and engineering, where systems exhibit non-local or memory effects. Its flexibility allows for modeling behaviors that standard derivatives cannot capture. More precisely, it is defined by letting ( GL D Īø α f ) ( z ) \= e āˆ’ i α Īø lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( z āˆ’ k h ) \| h \| α , where z ∈ C and h \= \| h \| e i Īø , āˆ’ Ļ€ \< Īø ≤ Ļ€. When the complex independent variable z becomes the real number t, we have the following two types of fractional derivative: - Forward derivative (h \= \| h \|): ( D F α GL f ) ( t ) \= lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( t āˆ’ k h ) \| h \| α - Backward derivative (h \= āˆ’ \| h \|): ( D B α GL f ) ( t ) \= lim \| h \| → 0 āˆ‘ k \= 0 āˆž ( āˆ’ 1 ) k C α k f ( t \+ k \| h \| ) \| h \| α We introduce only the forward derivative. - The operator is linear, that is, D F α GL ( Ī» f \+ ν g ) ( z ) \= Ī» ( GL D F α f ) ( z ) \+ μ ( GL D F α g ) ( z ) , āˆ€ Ī» , ν ∈ R . - Additivity and commutativity: D F α GL D F β GL ( f ) ( t ) \= D F β GL D F α GL ( f ) ( t ) \= ( D F α \+ β GL ( f ) ) ( t ) . - Neutral element. Letting β \= āˆ’ α yields that D F α GL D F āˆ’ α GL ( f ) ( t ) \= f ( t ) and D F āˆ’ α GL D F α GL ( f ) ( t ) \= f ( t ) . - Backward compatibility. - The generalized Leibniz rule. That is, for the derivative of the product φ ψ, we have D F α G L \[ φ ψ \] ( t ) \= āˆ‘ n \= 0 āˆž C n α D n φ ( t ) D F α āˆ’ n G L ψ ( t ) . ## References 1. Nguyen, A.T.; Caraballo, T.; Tuan, N.H. On the initial value problem for a class of nonlinear biharmonic equation with time-fractional derivative. Proc. R. Soc. Edinb. Sect. A Math. **2022**, 152, 989–1031. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+the+initial+value+problem+for+a+class+of+nonlinear+biharmonic+equation+with+time-fractional+derivative&author=Nguyen,+A.T.&author=Caraballo,+T.&author=Tuan,+N.H.&publication_year=2022&journal=Proc.+R.+Soc.+Edinb.+Sect.+A+Math.&volume=152&pages=989%E2%80%931031&doi=10.1017/prm.2021.44)\] \[[CrossRef](https://doi.org/10.1017/prm.2021.44)\] 2. 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](https://scholar.google.com/scholar_lookup?title=General+time-fractional+diffusion+equation:+Some+uniqueness+and+existence+results+for+the+initial-boundary-value+problems&author=Luchko,+Y.&author=Yamamoto,+M.&publication_year=2016&journal=Fract.+Calc.+Appl.+Anal.&volume=19&pages=676%E2%80%93695&doi=10.1515/fca-2016-0036)\] \[[CrossRef](https://doi.org/10.1515/fca-2016-0036)\] 3. Luchko, Y. Initial-boundary-value problems for the one-dimensional time-fractional diffusion equation. Fract. Calc. Appl. Anal. **2012**, 15, 141–160. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Initial-boundary-value+problems+for+the+one-dimensional+time-fractional+diffusion+equation&author=Luchko,+Y.&publication_year=2012&journal=Fract.+Calc.+Appl.+Anal.&volume=15&pages=141%E2%80%93160&doi=10.2478/s13540-012-0010-7)\] \[[CrossRef](https://doi.org/10.2478/s13540-012-0010-7)\] 4. Karthikeyan, K.; Karthikeyan, P.; Patel, D.K.; Al-Askar, F.M.; Cesarano, C. Solving fractional integro-differential equations with delay and relaxation impulsive terms by fixed point techniques. Bound. Value Probl. **2024**, 2024, 146. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solving+fractional+integro-differential+equations+with+delay+and+relaxation+impulsive+terms+by+fixed+point+techniques&author=Karthikeyan,+K.&author=Karthikeyan,+P.&author=Patel,+D.K.&author=Al-Askar,+F.M.&author=Cesarano,+C.&publication_year=2024&journal=Bound.+Value+Probl.&volume=2024&pages=146&doi=10.1186/s13661-024-01957-w)\] \[[CrossRef](https://doi.org/10.1186/s13661-024-01957-w)\] 5. Abbas, M.; Ciobanescu, C.; Asghar, M.W.; Omame, A. Solution approximation of fractional boundary value problems and convergence analysis using AA-iterative scheme. AIMS Math. **2024**, 9, 13129–13158. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solution+approximation+of+fractional+boundary+value+problems+and+convergence+analysis+using+AA-iterative+scheme&author=Abbas,+M.&author=Ciobanescu,+C.&author=Asghar,+M.W.&author=Omame,+A.&publication_year=2024&journal=AIMS+Math.&volume=9&pages=13129%E2%80%9313158&doi=10.3934/math.2024641)\] \[[CrossRef](https://doi.org/10.3934/math.2024641)\] 6. Alsaedi, A.; Ahmad, B.; Alharbi, N.; Alharbi, N.; Alharbi, N. Fixed point and stability analysis of a tripled system of nonlinear fractional differential equations with n-nonlinear terms. Fractal Fract. **2024**, 8, 697. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fixed+point+and+stability+analysis+of+a+tripled+system+of+nonlinear+fractional+differential+equations+with+n-nonlinear+terms&author=Alsaedi,+A.&author=Ahmad,+B.&author=Alharbi,+N.&author=Alharbi,+N.&author=Alharbi,+N.&publication_year=2024&journal=Fractal+Fract.&volume=8&pages=697&doi=10.3390/fractalfract8120697)\] \[[CrossRef](https://doi.org/10.3390/fractalfract8120697)\] 7. Mali, A.D.; JĆørgensen, J.B.; KĆøstner, N.; LĆønborg, J.; Dutta, H. Numerical investigation of generalized tempered-type integrodifferential equations with respect to another function. Fract. Calc. Appl. Anal. **2023**, 26, 2580–2601. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Numerical+investigation+of+generalized+tempered-type+integrodifferential+equations+with+respect+to+another+function&author=Mali,+A.D.&author=J%C3%B8rgensen,+J.B.&author=K%C3%B8stner,+N.&author=L%C3%B8nborg,+J.&author=Dutta,+H.&publication_year=2023&journal=Fract.+Calc.+Appl.+Anal.&volume=26&pages=2580%E2%80%932601&doi=10.1007/s13540-023-00198-5)\] \[[CrossRef](https://doi.org/10.1007/s13540-023-00198-5)\] 8. Gunasekar, T.; Raghavendran, P.; Santra, S.S.; Majumder, D.; Baleanu, D.; Balasundaram, H. Application of laplace transform to solve fractional integro-differential equations. J. Math. Comput. Sci. **2024**, 33, 225–237. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Application+of+laplace+transform+to+solve+fractional+integro-differential+equations&author=Gunasekar,+T.&author=Raghavendran,+P.&author=Santra,+S.S.&author=Majumder,+D.&author=Baleanu,+D.&author=Balasundaram,+H.&publication_year=2024&journal=J.+Math.+Comput.+Sci.&volume=33&pages=225%E2%80%93237&doi=10.22436/jmcs.033.03.02)\] \[[CrossRef](https://doi.org/10.22436/jmcs.033.03.02)\] 9. Ross, B. The development of fractional calculus 1695–1900. Historia Math. **1977**, 4, 75–89. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+development+of+fractional+calculus+1695%E2%80%931900&author=Ross,+B.&publication_year=1977&journal=Historia+Math.&volume=4&pages=75%E2%80%9389&doi=10.1016/0315-0860\(77\)90039-8)\] \[[CrossRef](https://doi.org/10.1016/0315-0860\(77\)90039-8)\] 10. El-Ajou, A.; Arqub, O.A.; Zhour, Z.A.; Momani, S. New results on fractional power series: Theories and applications. Entropy **2013**, 15, 5305–5323. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=New+results+on+fractional+power+series:+Theories+and+applications&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Zhour,+Z.A.&author=Momani,+S.&publication_year=2013&journal=Entropy&volume=15&pages=5305%E2%80%935323&doi=10.3390/e15125305)\] \[[CrossRef](https://doi.org/10.3390/e15125305)\] 11. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Volume 204 of North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Theory+and+Applications+of+Fractional+Differential+Equations&author=Kilbas,+A.A.&author=Srivastava,+H.M.&author=Trujillo,+J.J.&publication_year=2006)\] 12. Ortigueira, M.D.; Machado, J.A.T. What is a fractional derivative? J. Comput. Phys. **2015**, 293, 4–13. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=What+is+a+fractional+derivative?&author=Ortigueira,+M.D.&author=Machado,+J.A.T.&publication_year=2015&journal=J.+Comput.+Phys.&volume=293&pages=4%E2%80%9313&doi=10.1016/j.jcp.2014.07.019)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2014.07.019)\] 13. Katugampola, U.N. Correction to ā€œWhat is a fractional derivative?ā€ by Ortigueira and Machado \[Journal of Computational Physics, Volume 293, 15 July 2015, Pages 4–13. Special issue on fractional PDEs\]. J. Comput. Phys. **2016**, 321, 1255–1257. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Correction+to+%E2%80%9CWhat+is+a+fractional+derivative?%E2%80%9D+by+Ortigueira+and+Machado+%5BJournal+of+Computational+Physics,+Volume+293,+15+July+2015,+Pages+4%E2%80%9313.+Special+issue+on+fractional+PDEs%5D&author=Katugampola,+U.N.&publication_year=2016&journal=J.+Comput.+Phys.&volume=321&pages=1255%E2%80%931257&doi=10.1016/j.jcp.2016.05.052)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2016.05.052)\] 14. Kuroda, L.K.B.; Gomes, A.V.; Tavoni, R.; Mancera, P.F.d.; Varalta, N.; Camargo, R.d.F. Unexpected behavior of Caputo fractional derivative. Comput. Appl. Math. **2017**, 36, 1173–1183. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Unexpected+behavior+of+Caputo+fractional+derivative&author=Kuroda,+L.K.B.&author=Gomes,+A.V.&author=Tavoni,+R.&author=Mancera,+P.F.d.&author=Varalta,+N.&author=Camargo,+R.d.F.&publication_year=2017&journal=Comput.+Appl.+Math.&volume=36&pages=1173%E2%80%931183&doi=10.1007/s40314-015-0301-9)\] \[[CrossRef](https://doi.org/10.1007/s40314-015-0301-9)\] 15. Murio, D.A. Stable numerical evaluation of grünwald–letnikov fractional derivatives applied to a fractional ihcp. Inverse Probl. Sci. Eng. **2009**, 17, 229–243. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Stable+numerical+evaluation+of+gr%C3%BCnwald%E2%80%93letnikov+fractional+derivatives+applied+to+a+fractional+ihcp&author=Murio,+D.A.&publication_year=2009&journal=Inverse+Probl.+Sci.+Eng.&volume=17&pages=229%E2%80%93243&doi=10.1080/17415970802082872)\] \[[CrossRef](https://doi.org/10.1080/17415970802082872)\] 16. Chikrii, A.; Matychyn, I. Riemann–liouville, caputo, and sequential fractional derivatives in differential games. In Advances in Dynamic Games: Theory, Applications, and Numerical Methods for Differential and Stochastic Games; Springer: Berlin/Heidelberg, Germany, 2011; pp. 61–81. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Riemann%E2%80%93liouville,+caputo,+and+sequential+fractional+derivatives+in+differential+games&author=Chikrii,+A.&author=Matychyn,+I.&publication_year=2011&pages=61%E2%80%9381)\] 17. Agrawal, O.P. Fractional variational calculus in terms of riesz fractional derivatives. J. Phys. Math. Theor. **2007**, 40, 6287. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fractional+variational+calculus+in+terms+of+riesz+fractional+derivatives&author=Agrawal,+O.P.&publication_year=2007&journal=J.+Phys.+Math.+Theor.&volume=40&pages=6287&doi=10.1088/1751-8113/40/24/003)\] \[[CrossRef](https://doi.org/10.1088/1751-8113/40/24/003)\] 18. Laskin, N. Principles of fractional quantum mechanics. arXiv **2010**, arXiv:1009.5533. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Principles+of+fractional+quantum+mechanics&author=Laskin,+N.&publication_year=2010&journal=arXiv&doi=10.48550/arXiv.1009.5533)\] \[[CrossRef](https://doi.org/10.48550/arXiv.1009.5533)\] 19. Laskin, N. Fractional quantum mechanics. Phys. Rev. E **2000**, 62, 3135. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Fractional+quantum+mechanics&author=Laskin,+N.&publication_year=2000&journal=Phys.+Rev.+E&volume=62&pages=3135&doi=10.1103/PhysRevE.62.3135)\] \[[CrossRef](https://doi.org/10.1103/PhysRevE.62.3135)\] 20. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E **1997**, 55, 3581. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Mechanics+with+fractional+derivatives&author=Riewe,+F.&publication_year=1997&journal=Phys.+Rev.+E&volume=55&pages=3581&doi=10.1103/PhysRevE.55.3581)\] \[[CrossRef](https://doi.org/10.1103/PhysRevE.55.3581)\] 21. Falcini, F.; Garra, R.; Voller, V. Modeling anomalous heat diffusion: Comparing fractional derivative and non-linear diffusivity treatments. Int. J. Therm. Sci. **2019**, 137, 584–588. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Modeling+anomalous+heat+diffusion:+Comparing+fractional+derivative+and+non-linear+diffusivity+treatments&author=Falcini,+F.&author=Garra,+R.&author=Voller,+V.&publication_year=2019&journal=Int.+J.+Therm.+Sci.&volume=137&pages=584%E2%80%93588&doi=10.1016/j.ijthermalsci.2018.11.016)\] \[[CrossRef](https://doi.org/10.1016/j.ijthermalsci.2018.11.016)\] 22. Garra, R.; Falcini, F.; Voller, V.R.; Pagnini, G. A generalized stefan model accounting for system memory and non-locality. Int. Commun. Heat Mass Transf. **2020**, 114, 104584. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+generalized+stefan+model+accounting+for+system+memory+and+non-locality&author=Garra,+R.&author=Falcini,+F.&author=Voller,+V.R.&author=Pagnini,+G.&publication_year=2020&journal=Int.+Commun.+Heat+Mass+Transf.&volume=114&pages=104584&doi=10.1016/j.icheatmasstransfer.2020.104584)\] \[[CrossRef](https://doi.org/10.1016/j.icheatmasstransfer.2020.104584)\] 23. Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. **2006**, 316, 753–763. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=An+iterative+method+for+solving+nonlinear+functional+equations&author=Daftardar-Gejji,+V.&author=Jafari,+H.&publication_year=2006&journal=J.+Math.+Anal.+Appl.&volume=316&pages=753%E2%80%93763&doi=10.1016/j.jmaa.2005.05.009)\] \[[CrossRef](https://doi.org/10.1016/j.jmaa.2005.05.009)\] 24. Noor, M.A.; Noor, K.I.; Mohyud-Din, S.T.; Shabbir, A. An iterative method with cubic convergence for nonlinear equations. Appl. Math. Comput. **2006**, 183, 1249–1255. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=An+iterative+method+with+cubic+convergence+for+nonlinear+equations&author=Noor,+M.A.&author=Noor,+K.I.&author=Mohyud-Din,+S.T.&author=Shabbir,+A.&publication_year=2006&journal=Appl.+Math.+Comput.&volume=183&pages=1249%E2%80%931255&doi=10.1016/j.amc.2006.05.133)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2006.05.133)\] 25. Bhalekar, S.; Daftardar-Gejji, V. New iterative method: Application to partial differential equations. Appl. Math. Comput. **2008**, 203, 778–783. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=New+iterative+method:+Application+to+partial+differential+equations&author=Bhalekar,+S.&author=Daftardar-Gejji,+V.&publication_year=2008&journal=Appl.+Math.+Comput.&volume=203&pages=778%E2%80%93783&doi=10.1016/j.amc.2008.05.071)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2008.05.071)\] 26. Daftardar-Gejji, V.; Bhalekar, S. Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method. Comput. Math. Appl. **2010**, 59, 1801–1809. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solving+fractional+boundary+value+problems+with+Dirichlet+boundary+conditions+using+a+new+iterative+method&author=Daftardar-Gejji,+V.&author=Bhalekar,+S.&publication_year=2010&journal=Comput.+Math.+Appl.&volume=59&pages=1801%E2%80%931809&doi=10.1016/j.camwa.2009.08.018)\] \[[CrossRef](https://doi.org/10.1016/j.camwa.2009.08.018)\] 27. Dhaigude, C.D.; Nikam, V.R. Solution of fractional partial differential equations using iterative method. Fract. Calc. Appl. Anal. **2012**, 15, 684–699. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solution+of+fractional+partial+differential+equations+using+iterative+method&author=Dhaigude,+C.D.&author=Nikam,+V.R.&publication_year=2012&journal=Fract.+Calc.+Appl.+Anal.&volume=15&pages=684%E2%80%93699&doi=10.2478/s13540-012-0046-8)\] \[[CrossRef](https://doi.org/10.2478/s13540-012-0046-8)\] 28. Liao, S. On the homotopy analysis method for nonlinear problems. Appl. Math. Comput. **2004**, 147, 499–513. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+the+homotopy+analysis+method+for+nonlinear+problems&author=Liao,+S.&publication_year=2004&journal=Appl.+Math.+Comput.&volume=147&pages=499%E2%80%93513&doi=10.1016/S0096-3003\(02\)00790-7)\] \[[CrossRef](https://doi.org/10.1016/S0096-3003\(02\)00790-7)\] 29. Das, S.; Vishal, K. Homotopy analysis method for fractional Swift-Hohenberg equation. In Advances in the Homotopy Analysis Method; World Scientific Publishing: Hackensack, NJ, USA, 2014; pp. 291–308. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Homotopy+analysis+method+for+fractional+Swift-Hohenberg+equation&author=Das,+S.&author=Vishal,+K.&publication_year=2014&pages=291%E2%80%93308)\] \[[CrossRef](https://doi.org/10.1142/9789814551250_0007)\] 30. Demir, A.; Bayrak, M.A.; Ozbilge, E. A new approach for the approximate analytical solution of space-time fractional differential equations by the homotopy analysis method. Adv. Math. Phys. **2019**, 2019, 5602565. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+new+approach+for+the+approximate+analytical+solution+of+space-time+fractional+differential+equations+by+the+homotopy+analysis+method&author=Demir,+A.&author=Bayrak,+M.A.&author=Ozbilge,+E.&publication_year=2019&journal=Adv.+Math.+Phys.&volume=2019&pages=5602565&doi=10.1155/2019/5602565)\] \[[CrossRef](https://doi.org/10.1155/2019/5602565)\] 31. Li, L.; Li, D. Exact solutions and numerical study of time fractional Burgers’ equations. Appl. Math. Lett. **2020**, 100, 106011. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Exact+solutions+and+numerical+study+of+time+fractional+Burgers%E2%80%99+equations&author=Li,+L.&author=Li,+D.&publication_year=2020&journal=Appl.+Math.+Lett.&volume=100&pages=106011&doi=10.1016/j.aml.2019.106011)\] \[[CrossRef](https://doi.org/10.1016/j.aml.2019.106011)\] 32. Arqub, O.A. Series solution of fuzzy differential equations under strongly generalized differentiability. J. Adv. Res. Appl. Math. **2013**, 5, 31–52. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Series+solution+of+fuzzy+differential+equations+under+strongly+generalized+differentiability&author=Arqub,+O.A.&publication_year=2013&journal=J.+Adv.+Res.+Appl.+Math.&volume=5&pages=31%E2%80%9352&doi=10.5373/jaram.1447.051912)\] \[[CrossRef](https://doi.org/10.5373/jaram.1447.051912)\] 33. El-Ajou, A.; Arqub, O.A.; Momani, S.; Baleanu, D.; Alsaedi, A. A novel expansion iterative method for solving linear partial differential equations of fractional order. Appl. Math. Comput. **2015**, 257, 119–133. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+novel+expansion+iterative+method+for+solving+linear+partial+differential+equations+of+fractional+order&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Momani,+S.&author=Baleanu,+D.&author=Alsaedi,+A.&publication_year=2015&journal=Appl.+Math.+Comput.&volume=257&pages=119%E2%80%93133&doi=10.1016/j.amc.2014.12.121)\] \[[CrossRef](https://doi.org/10.1016/j.amc.2014.12.121)\] 34. Dubey, V.P.; Kumar, R.; Kumar, D. A reliable treatment of residual power series method for time-fractional Black–Scholes European option pricing equations. Physical A **2019**, 533, 122040. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=A+reliable+treatment+of+residual+power+series+method+for+time-fractional+Black%E2%80%93Scholes+European+option+pricing+equations&author=Dubey,+V.P.&author=Kumar,+R.&author=Kumar,+D.&publication_year=2019&journal=Physical+A&volume=533&pages=122040&doi=10.1016/j.physa.2019.122040)\] \[[CrossRef](https://doi.org/10.1016/j.physa.2019.122040)\] 35. Syam, M.I. On approximate solutions of Fokker Planck equation by the modified residual power series method. J. Fract. Calc. Appl. **2019**, 10, 225–232. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=On+approximate+solutions+of+Fokker+Planck+equation+by+the+modified+residual+power+series+method&author=Syam,+M.I.&publication_year=2019&journal=J.+Fract.+Calc.+Appl.&volume=10&pages=225%E2%80%93232)\] 36. Az-Zo’bi, E.A.; Yıldırım, A.; AlZoubi, W.A. The residual power series method for the one-dimensional unsteady flow of a van der Waals gas. Physical A **2019**, 517, 188–196. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+residual+power+series+method+for+the+one-dimensional+unsteady+flow+of+a+van+der+Waals+gas&author=Az-Zo%E2%80%99bi,+E.A.&author=Y%C4%B1ld%C4%B1r%C4%B1m,+A.&author=AlZoubi,+W.A.&publication_year=2019&journal=Physical+A&volume=517&pages=188%E2%80%93196&doi=10.1016/j.physa.2018.11.030)\] \[[CrossRef](https://doi.org/10.1016/j.physa.2018.11.030)\] 37. Mosaleheh, K.; Vakilzadeh, A. The Laplace transform method for linear ordinary differential equations of fractional order. J. Math. Ext. **2007**, 2, 93–102. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+Laplace+transform+method+for+linear+ordinary+differential+equations+of+fractional+order&author=Mosaleheh,+K.&author=Vakilzadeh,+A.&publication_year=2007&journal=J.+Math.+Ext.&volume=2&pages=93%E2%80%93102)\] 38. Artin, E. The Gamma Function; Butler, M., Translator; Athena Series: Selected Topics in Mathematics; Holt, Rinehart and Winston: New York, NY, USA; Toronto, ON, Canada; London, UK, 1964. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=The+Gamma+Function&author=Artin,+E.&publication_year=1964)\] 39. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer Monographs in Mathematics; Springer: Heidelberg, Germany, 2014. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Mittag-Leffler+Functions,+Related+Topics+and+Applications&author=Gorenflo,+R.&author=Kilbas,+A.A.&author=Mainardi,+F.&author=Rogosin,+S.V.&publication_year=2014)\] 40. El-Ajou, A.; Arqub, O.A.; Momani, S. Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: A new iterative algorithm. J. Comput. Phys. **2015**, 293, 81–95. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Approximate+analytical+solution+of+the+nonlinear+fractional+KdV-Burgers+equation:+A+new+iterative+algorithm&author=El-Ajou,+A.&author=Arqub,+O.A.&author=Momani,+S.&publication_year=2015&journal=J.+Comput.+Phys.&volume=293&pages=81%E2%80%9395&doi=10.1016/j.jcp.2014.08.004)\] \[[CrossRef](https://doi.org/10.1016/j.jcp.2014.08.004)\] 41. Arqub, O.A. Solutions of time-fractional Tricomi and Keldysh equations of Dirichlet functions types in Hilbert space. Numer. Methods Partial. Differ. Equ. **2018**, 34, 1759–1780. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solutions+of+time-fractional+Tricomi+and+Keldysh+equations+of+Dirichlet+functions+types+in+Hilbert+space&author=Arqub,+O.A.&publication_year=2018&journal=Numer.+Methods+Partial.+Differ.+Equ.&volume=34&pages=1759%E2%80%931780&doi=10.1002/num.22236)\] \[[CrossRef](https://doi.org/10.1002/num.22236)\] 42. Zhang, X.; Liu, J.; Wen, J.; Tang, B.; He, Y. Analysis for one-dimensional time-fractional Tricomi-type equations by LDG methods. Numer. Algorithms **2013**, 63, 143–164. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Analysis+for+one-dimensional+time-fractional+Tricomi-type+equations+by+LDG+methods&author=Zhang,+X.&author=Liu,+J.&author=Wen,+J.&author=Tang,+B.&author=He,+Y.&publication_year=2013&journal=Numer.+Algorithms&volume=63&pages=143%E2%80%93164&doi=10.1007/s11075-012-9617-3)\] \[[CrossRef](https://doi.org/10.1007/s11075-012-9617-3)\] 43. Frankl, F.L. On the Problems of Chaplygin for Mixed Sub- and Supersonic Flows. 1947. Available online: <https://ntrs.nasa.gov/search.jsp?R=20050019646> (accessed on 10 March 2026). 44. El-Ajou, A.; Oqielat, M.N.; Al-Zhour, Z.; Kumar, S.; Momani, S. Solitary solutions for time-fractional nonlinear dispersive PDEs in the sense of conformable fractional derivative. Chaos **2019**, 29, 093102. \[[Google Scholar](https://scholar.google.com/scholar_lookup?title=Solitary+solutions+for+time-fractional+nonlinear+dispersive+PDEs+in+the+sense+of+conformable+fractional+derivative&author=El-Ajou,+A.&author=Oqielat,+M.N.&author=Al-Zhour,+Z.&author=Kumar,+S.&author=Momani,+S.&publication_year=2019&journal=Chaos&volume=29&pages=093102&doi=10.1063/1.5100234)\] \[[CrossRef](https://doi.org/10.1063/1.5100234)\] | | | |---|---| | | **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. | Ā© 2026 by the authors. Licensee MDPI, Basel, Switzerland. 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/). ## Article Metrics Citations Multiple requests from the same IP address are counted as one view.
Shard31 (laksa)
Root Hash11974975279136771031
Unparsed URLcom,mdpi!www,/2227-7390/14/6/984 s443