This study investigates analytical solutions of fractional differential equations (FDEs) using the MathHandbook symbolic computation platform. The objective is to evaluate how the unified operator ds (y,x,α) and the dsolve engine handle differential equations of arbitrary order, including integer-order ODEs/PDEs, fractional ODEs and PDEs, fractional integral equations, mixed differential–integral equations, systems of fractional equations, complex-order equations, and variable-order equations. Comparative analysis shows that general solutions exhibit order-dependent exponential behavior, while particular solutions remain invariant for linear equations with identical forcing terms. These results highlight the coherence of fractional calculus and demonstrate how AI-Driven symbolic computation enhances mathematical modeling, verification, and mathematics education.
Differential equations form the foundation of mathematical physics, modeling phenomena such as heat conduction, diffusion, elasticity, gravitation, and electrochemistry [1-6]. While classical differential equations involve integer-order derivatives, many real-world systems exhibit memory, nonlocality, and hereditary behavior that are better captured by fractional differential equations (FDEs) [7-9].
Fractional calculus extends differentiation and integration to arbitrary real or complex orders. FDEs appear in physics, chemistry, biology, biophysics, economics, control theory, and signal processing. Nonlinear fractional models often display rich dynamical behavior, including chaos [10-11].
This case study evaluates MathHandbook—an evolution of SymbMath—as a platform [12-15] for solving fractional and classical differential equations analytically.
MathHandbook is an online symbolic mathematics engine capable of solving differential equations of arbitrary order. Its dsolve function supports:
Integer-order odes and pdes
Fractional differential equations (Caputo definition by default)
Fractional integral equations
Mixed differential-integral equations
Systems of fractional equations
Complex-order differential equations
Variable-order differential equations
Users may input equations directly (e.g., y(1,x)-2y=exp(x)) or use symbolic operators such as ds(y,x,α) for fractional derivatives.
MathHandbook uses a unified operator
ds(y,x,α)=dαydxα,α∈ℂ.\begin{matrix} ds(y,x,\alpha) & = \frac{d^{\,\alpha}y}{dx^{\,\alpha}}, & & \alpha \in \mathbb{C}. \end{matrix}
to represent derivatives and integrals of arbitrary order:
ds(y,x,n)=dnydxn,n>0,ds(y,x,-n)=∫nydx,n>0,\begin{matrix} ds(y,x,n) & = \frac{d^{\, n}y}{dx^{\, n}}, & & n > 0, \\ ds(y,x, - n) & = \int_{}^{\, n}y\, dx, & & n > 0, \end{matrix}
Caputo derivative (default):
CDxαf(x)=1Γ(1-α)∫0x(x-t)-αf'(t)dt,0<α<1^{C}D_{x}^{\alpha}f(x) = \frac{1}{\Gamma(1 - \alpha)}\int_{0}^{x}(x - t)^{- \alpha}f'(t)\, dt,\quad\quad 0 < \alpha < 1
where Γ\Gamma is the Gamma function.
RLDxαf(x)=1Γ(1-α)∫-∞x(x-t)-αf(t)dt.^{RL}D_{x}^{\alpha}f(x) = \frac{1}{\Gamma(1 - \alpha)}\int_{- \infty}^{x}(x - t)^{- \alpha}f(t)\, dt.
dαdxαex=ex,for all α∈ℂ.\frac{d^{\alpha}}{dx^{\alpha}}e^{x} = e^{x},\quad\quad\text{for all }\alpha \in \mathbb{C}.
Since exe^{x} is an eigenfunction of both the Caputo and Riemann-Liouville operators, exponential terms retain their structure across all orders. This property ensures that exponential terms appearing in general solutions retain their structure across all differential orders, which is essential for the unified analysis presented in later sections.
Consider the first-order linear differential equation:
dydx-2y=ex.\frac{dy}{dx} - 2y = e^{x}.
dsolve (ds(y,x,1) - 2*y = exp(x) )
The system returns both the general solution and a particular solution: C1e2x-exC_{1}e^{2x} - e^{x}.
General solution: yg=C1e2xy_{g} = C_{1}e^{2x} Particular solution: yp=-exy_{p} = - e^{x}
summarizes the operator notation for a range of derivative and integral orders, including fractional, negative, complex, and variable orders. This provides a unified framework for expressing all examples used in Section 4.
| Order | Equation | Operator Form |
|---|---|---|
| 2 | d2ydx2-2y=ex\frac{d^{2}y}{dx^{2}} - 2y = e^{x} | ds(y,x,2)ds(y,x,2) |
| 1.5 | d1.5ydx1.5-2y=ex\frac{d^{1.5}y}{dx^{1.5}} - 2y = e^{x} | ds(y,x,1.5)ds(y,x,1.5) |
| 1 | dydx-2y=ex\frac{dy}{dx} - 2y = e^{x} | ds(y,x)ds(y,x) |
| 0.5 | d0.5ydx0.5-2y=ex\frac{d^{0.5}y}{dx^{0.5}} - 2y = e^{x} | ds(y,x,0.5)ds(y,x,0.5) |
| -0.5 | d-0.5ydx-0.5-2y=ex\frac{d^{- 0.5}y}{dx^{- 0.5}} - 2y = e^{x} | ds(y,x,-0.5)ds(y,x, - 0.5) |
| -1 | ∫ydx-2y=ex\int y\, dx - 2y = e^{x} | ds(y,x,-1)ds(y,x, - 1) |
| -2 | ∬y(dx)2-2y=ex\iint y\, dx^{2} - 2y = e^{x} | ds(y,x,-2)ds(y,x, - 2) |
| ii | diydxi-2y=ex\frac{d^{i}y}{dx^{i}} - 2y = e^{x} | ds(y,x,i)ds(y,x,i) |
| cos(x)\text{cos}(x) | dcos(x)ydxcos(x)-2y=ex\frac{d^{\text{cos}(x)}y}{dx^{\text{cos}(x)}} - 2y = e^{x} | ds(y,x,cos(x))ds(y,x,cos(x)) |
| Order | Equation | Operator Form |
|---|---|---|
| 2 | d2ydx2-2y=ex\frac{d^{2}y}{dx^{2}} - 2y = e^{x} | ds(y,x,2)ds(y,x,2) |
| 1.5 | d1.5ydx1.5-2y=ex\frac{d^{1.5}y}{dx^{1.5}} - 2y = e^{x} | ds(y,x,1.5)ds(y,x,1.5) |
| 1 | dydx-2y=ex\frac{dy}{dx} - 2y = e^{x} | ds(y,x)ds(y,x) |
| 0.5 | d0.5ydx0.5-2y=ex\frac{d^{0.5}y}{dx^{0.5}} - 2y = e^{x} | ds(y,x,0.5)ds(y,x,0.5) |
| -0.5 | d-0.5ydx-0.5-2y=ex\frac{d^{- 0.5}y}{dx^{- 0.5}} - 2y = e^{x} | ds(y,x,-0.5)ds(y,x, - 0.5) |
| -1 | ∫ydx-2y=ex\int y\, dx - 2y = e^{x} | ds(y,x,-1)ds(y,x, - 1) |
| -2 | ∬y(dx)2-2y=ex\iint y\, dx^{2} - 2y = e^{x} | ds(y,x,-2)ds(y,x, - 2) |
| ii | diydxi-2y=ex\frac{d^{i}y}{dx^{i}} - 2y = e^{x} | ds(y,x,i)ds(y,x,i) |
| cos(x)\text{cos}(x) | dcos(x)ydxcos(x)-2y=ex\frac{d^{\text{cos}(x)}y}{dx^{\text{cos}(x)}} - 2y = e^{x} | ds(y,x,cos(x))ds(y,x,cos(x)) |
This section presents a comprehensive set of analytical results obtained using the
MathHandbook solves both symbolically and numerically solvable ODEs.
y'=exp(y)y' = exp(y)
This first-order nonlinear ODE is solved:
y=-log(-C1-x).y = - \text{log}( - C_{1} - x).
MathHandbook supports 2D and 3D visualization of PDE solutions. Classical PDEs such as Laplace, Poisson, diffusion, and wave equations are included.
Laplace equation: Δy=0\Delta y = 0
Poisson equation: Δy=ρ\Delta y = \rho
Heat equation: yt-a2Δy=f(t,x)y_{t} - a^{2}\Delta y = f(t,x)
Wave equation: ytt-a2Δy=f(t,x,z)y_{tt} - a^{2}\Delta y = f(t,x,z)
Solutions can be visualized using plot2D and plot3D, enabling interactive exploration of time-dependent behavior.
Fractional PDEs are obtained by replacing the time derivative with a fractional order α\alpha, typically 0<α<10 < \alpha < 1.
d0.5ydt0.5=yx-ex-et\frac{d^{0.5}y}{dt^{0.5}} = y_{x} - e^{x} - e^{t}
MathHandbook solves such equations symbolically when possible, or numerically when closed-form solutions are unavailable.
Negative-order derivatives treated as integrals. Fractional integral equations extend this to non-integer negative orders.
∫ydx-2y=ex\int y\, dx - 2y = e^{x}
d-0.5ydx-0.5-2y=ex\frac{d^{- 0.5}y}{dx^{- 0.5}} - 2y = e^{x}
MathHandbook converts these to differential form by differentiation, solves them, and then verifies the solution.
Equations combining fractional derivatives and fractional integrals are also supported.
dπydxπ-∫y(x)dxπ=ex\frac{d^{\pi}y}{dx^{\pi}} - \int y(x)\, dx^{\pi} = e^{x}
Systems of equations are solved simultaneously.
d0.5xdt0.5=2x,d0.5ydt0.5=x+y\frac{d^{0.5}x}{dt^{0.5}} = 2x,\quad\quad\frac{d^{0.5}y}{dt^{0.5}} = x + y
MathHandbook extends differentiation to complex orders, producing oscillatory-exponential solutions.
diydxi-2y=ex\frac{d^{i}y}{dx^{i}} - 2y = e^{x}
The general solution involves e12ixe^{\frac{1}{2}ix}, reflecting complex-order behavior.
When the order varies with the independent variable, the solution typically involves Mittag-Leffler functions.
dcos(x)ydxcos(x)-2y=ex\frac{d^{\text{cos}(x)}y}{dx^{\text{cos}(x)}} - 2y = e^{x}
The general solution is:
yg(x)=C1Ecos(x)(2xcos(x))y_{g}(x) = C_{1}\, E_{\text{cos}(x)}(2x^{\text{cos}(x)})
where Eα(⋅)E_{\alpha}( \cdot ) is the Mittag-Leffler function.
| Order | Equation | General Solution | Particular Solution |
|---|---|---|---|
| 2 | d2ydx2-2y=ex\frac{d^{2}y}{dx^{2}} - 2y = e^{x} | C1e2xC_{1}e^{\sqrt{2}x} | -ex- e^{x} |
| 1.5 | d1.5ydx1.5-2y=ex\frac{d^{1.5}y}{dx^{1.5}} - 2y = e^{x} | C1e223xC_{1}e^{\frac{2}{\sqrt |
-ex- e^{x} |
| 1 | dydx-2y=ex\frac{dy}{dx} - 2y = e^{x} | C1e2xC_{1}e^{2x} | -ex- e^{x} |
| 0.5 | d0.5ydx0.5-2y=ex\frac{d^{0.5}y}{dx^{0.5}} - 2y = e^{x} | C1e4xC_{1}e^{4x} | -ex- e^{x} |
| -0.5 | d-0.5ydx-0.5-2y=ex\frac{d^{- 0.5}y}{dx^{- 0.5}} - 2y = e^{x} | C1ex/4C_{1}e^{x/4} | -ex- e^{x} |
| -1 | ∫ydx-2y=ex\int y\, dx - 2y = e^{x} | C1ex/2C_{1}e^{x/2} | -ex- e^{x} |
| -2 | ∬y(dx)2-2y=ex\iint y\, dx^{2} - 2y = e^{x} | C1ex/2C_{1}e^{x/\sqrt{2}} | -ex- e^{x} |
| ii | diydxi-2y=ex\frac{d^{i}y}{dx^{i}} - 2y = e^{x} | C1ei2xC_{1}e^{\frac{i}{2}x} | -ex- e^{x} |
| cos(x)\text{cos}(x) | dcos(x)ydxcos(x)-2y=ex\frac{d^{\text{cos}(x)}y}{dx^{\text{cos}(x)}} - 2y = e^{x} | C1Ecos(x)(2xcos(x))C_{1}E_{\text{cos}(x)}(2x^{\text{cos}(x)}) | -ex- e^{x} |
The results of this study reveal several structural properties that are fundamental to understanding fractional differential equations and their computational treatment. First, the preservation of linear structure across fractional orders demonstrates that fractional derivatives extend classical calculus in a mathematically coherent way.
Table 2 shows that solution structures of differential equations are the same regardless of its order:
The smooth variation of the general solution with respect to the order α\alpha highlights the continuous transition between integer-order and fractional-order dynamics, offering a unified framework for modeling systems with memory or hereditary effects.
A particularly notable finding is the invariance of the particular solution across all constant orders for the equation
dαy/dxα-2y=ex,d^{\alpha}y/dx^{\alpha} - 2y = e^{x},
which underscores the robustness of forcing-term responses in linear fractional systems. This behavior simplifies analytical interpretation and provides a stable foundation for numerical verification.
The emergence of Mittag-Leffler functions in variable-order equations reflects the richer dynamical landscape introduced by non-constant differentiation orders. These functions generalize exponential behavior and naturally encode memory effects, making them indispensable in modeling anomalous diffusion, viscoelasticity, and other nonlocal processes. Similarly, complex-order derivatives produce oscillatory-exponential hybrid solutions, suggesting potential applications in signal processing, wave propagation, and systems exhibiting combined growth and rotation.
Overall, the results demonstrate that fractional calculus offers a flexible and powerful modeling paradigm. The ability of MathHandbook to compute symbolic solutions across such a wide range of orders including fractional, integral, mixed, complex, and variable orders illustrates the maturity of modern symbolic computation and its value for both theoretical research and applied problem solving.
AI-powered symbolic computation tools such as MathHandbook significantly enhance mathematics education by:
providing instant analytical solutions
supporting experimentation with fractional orders
visualizing PDEs in 2D and 3D
enabling students to verify solutions interactively
bridging theory and computation in differential equations courses
As AI continues to evolve, such tools will play an increasingly central role in personalized and exploratory learning.
This case study provides a comprehensive examination of fractional differential equations solved using the MathHandbook symbolic computation platform. The system successfully handles an extensive spectrum of differential operators, including integer-order, fractional-order, integral, mixed-order, complex-order, and variable-order derivatives. The comparative analysis across orders reveals consistent structural patterns: general solutions exhibit order-dependent exponential behavior, while particular solutions remain invariant for linear equations with identical forcing terms.
These findings reinforce the theoretical coherence of fractional calculus and highlight its capacity to model memory-dependent and nonlocal phenomena. The integration of AI-driven symbolic computation into MathHandbook further enhances accessibility, enabling researchers, educators, and students to explore fractional models interactively and verify solutions analytically or numerically.
As fractional calculus continues to expand across physics, engineering, biology, and applied mathematics, tools such as MathHandbook will play an increasingly important role in supporting research, facilitating experimentation, and bridging the gap between theory and computation. The results presented here demonstrate that modern symbolic engines are not only capable of solving complex fractional systems but also of providing insights that deepen our understanding of their underlying mathematical structure.