A $\star$-Product Approach for Analytical and Numerical Solutions of Nonautonomous Linear Fractional Differential Equations

Fabio Durastante; Pierre-Louis Giscard; Stefano Pozza

arxiv: 2507.16652 · v2 · submitted 2025-07-22 · 🧮 math.NA · cs.NA

A star-Product Approach for Analytical and Numerical Solutions of Nonautonomous Linear Fractional Differential Equations

Fabio Durastante , Pierre-Louis Giscard , Stefano Pozza This is my paper

Pith reviewed 2026-05-19 03:08 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords fractional differential equationsnonautonomous systemsstar productVolterra convolutionnumerical discretizationanalytical solutions

0 comments

The pith

A star-product reformulation yields analytical and numerical solutions for nonautonomous fractional differential equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a method to solve nonautonomous linear fractional differential equations by expressing their analytical solutions using a star-product operation that generalizes the Volterra convolution. It then discretizes this expression to obtain numerical solutions. In some special cases, the same framework produces closed-form expressions. Readers might care because many systems involving fractional derivatives with time-varying coefficients lack straightforward solution methods.

Core claim

The central claim is that reformulating the analytical solution of nonautonomous linear fractional differential equations via the star-product allows both the derivation of closed-form solutions in certain cases and an effective discretization for numerical approximation in general.

What carries the argument

The star-product, a generalization of the Volterra convolution, which reformulates the solution operator to accommodate time-dependent coefficients.

If this is right

Closed-form solutions become derivable for particular classes of time-dependent coefficients.
Numerical approximations follow directly from discretizing the star-product expression.
The approach covers linear fractional problems beyond the constant-coefficient setting.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The method could be applied to time-varying fractional models in diffusion or control problems to obtain practical solutions.
One could combine the star-product discretization with existing numerical libraries to handle higher-dimensional or nonlinear extensions.

Load-bearing premise

The star-product must correctly represent the solution operator for equations with variable coefficients, and discretization must preserve the properties required for accuracy and convergence.

What would settle it

Compare the star-product discretization output against a known exact solution for a specific nonautonomous fractional equation and verify whether the approximation error decreases with finer discretization steps.

read the original abstract

This article presents a novel solution method for nonautonomous linear ordinary fractional differential equations. The approach is based on reformulating the analytical solution using the $\star$-product, a generalization of the Volterra convolution, followed by an appropriate discretization of the resulting expression. Additionally, we demonstrate that, in certain cases, the $\star$-formalism enables the derivation of closed-form solutions, further highlighting the utility of this framework.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The star-product gives a workable reformulation for nonautonomous linear fractional DEs but the verification step for time-varying coefficients is the part that still needs tightening.

read the letter

The main takeaway is that the paper uses a generalized Volterra convolution, the star-product, to rewrite the solution operator for equations like D^α y(t) = a(t) y(t) + f(t), then discretizes the result and recovers closed forms in special cases. That framing is the actual new piece. It does a clean job of extending the convolution idea to variable coefficients and showing how the same expression can be turned into a numerical scheme without starting from scratch each time. The recovery of known solutions when a(t) is constant is a useful sanity check and gives some that the setup is consistent with existing theory. The soft spot is exactly where the stress-test note flags it: the central claim that the star-product satisfies the correct composition rule with the fractional kernel and the multiplier a(t) is stated but not shown with an explicit substitution back into the original equation or a direct comparison to the variation-of-parameters formula for non-constant a(t). Without that step the discretization inherits an unverified starting point, and referees will almost certainly ask for it. The rest of the paper looks internally coherent and the discretization appears straightforward to implement. This is aimed at people already working in numerical fractional calculus who want a convolution-based handle on the nonautonomous case. A reader in that niche could get practical value from the framework once the verification is tightened. It is solid enough on its own terms to deserve a serious referee rather than a desk reject, even if the review will focus on that missing check.

Referee Report

1 major / 1 minor

Summary. The paper proposes a ⋆-product (a generalized Volterra convolution) to reformulate the analytical solution of nonautonomous linear fractional differential equations, followed by discretization of the resulting expression for numerical solutions and derivation of closed-form solutions in special cases.

Significance. If the reformulation is rigorously verified and the discretization is shown to converge, the approach could supply a useful framework for both closed-form analysis and reliable numerics in fractional models with time-dependent coefficients, an area where standard methods often struggle.

major comments (1)

[§3] §3 (reformulation of the solution operator): the manuscript states that the ⋆-product exactly reproduces the variation-of-parameters integral operator for D^α y(t) = a(t)y(t) + f(t) with non-constant a(t), yet provides no explicit verification step such as direct substitution of the reformulated expression back into the original FDE or reduction to the constant-coefficient case. This check is load-bearing for the subsequent discretization claims.

minor comments (1)

[§2] The definition of the ⋆-product and its algebraic properties could be illustrated with one or two explicit low-order examples to improve readability for readers unfamiliar with the construction.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which have helped us improve the manuscript. We address the single major comment below and have incorporated the suggested verification into the revised version.

read point-by-point responses

Referee: [§3] §3 (reformulation of the solution operator): the manuscript states that the ⋆-product exactly reproduces the variation-of-parameters integral operator for D^α y(t) = a(t)y(t) + f(t) with non-constant a(t), yet provides no explicit verification step such as direct substitution of the reformulated expression back into the original FDE or reduction to the constant-coefficient case. This check is load-bearing for the subsequent discretization claims.

Authors: We agree that an explicit verification step strengthens the foundation of the reformulation. The ⋆-product was constructed precisely to reproduce the variation-of-parameters integral operator, but the original manuscript omitted a direct check. In the revised Section 3 we now include (i) direct substitution of the ⋆-product expression back into the original Caputo fractional differential equation to confirm that it satisfies the equation identically, and (ii) the reduction to the constant-coefficient case, where the expression recovers the standard Mittag-Leffler solution. These additions directly address the referee’s concern and provide the necessary justification for the subsequent discretization analysis. revision: yes

Circularity Check

0 steps flagged

No circularity: ⋆-product reformulation is a self-contained novel framework

full rationale

The manuscript introduces the ⋆-product as an explicit generalization of the Volterra convolution and applies it to recast the solution operator for nonautonomous linear fractional DEs, then discretizes the resulting expression. This construction supplies both the analytical reformulation and the numerical scheme without any quoted step reducing the claimed solution to a prior fit, self-citation, or definitional tautology. The derivation chain therefore remains independent of its own outputs and qualifies as a standard presentation of a new operator-based method.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on standard properties of fractional calculus and the definition of the ⋆-product as a generalization of Volterra convolution; no free parameters, invented entities, or ad-hoc axioms are visible in the abstract.

axioms (1)

domain assumption The ⋆-product is a well-defined associative operation that extends the Volterra convolution and correctly encodes the solution operator for linear fractional differential equations.
Invoked when the abstract states that the analytical solution can be reformulated using the ⋆-product.

pith-pipeline@v0.9.0 · 5597 in / 1217 out tokens · 36036 ms · 2026-05-19T03:08:31.834409+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Theorem 1: solution given by Θ⋆α ⋆ (δ − f̃ Θ⋆α)⋆−1 ⋆ Θ⋆(1−α) ỹs (Eq. 3.2); ⋆-product defined as integral composition (Eq. 2.1)
IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean absolute_floor_iff_bare_distinguishability unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Lemma 1: Θ⋆α(t,s) = (t−s)^{α−1}/Γ(α) Θ(t−s); Drazin inverse Dα of Θ⋆α

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages · 1 internal anchor

[1]

Aceto, L., Novati, P.: Fast and accurate approximations to fractional powers of opera- tors. IMA J. Numer. Anal.42(2), 1598–1622 (2022). DOI 10.1093/imanum/drab002. URLhttps://doi.org/10.1093/imanum/drab002

work page doi:10.1093/imanum/drab002 2022
[2]

Journal of Rheology27(3), 201–210 (1983)

Bagley, R.L., Torvik, P.J.: A Theoretical Basis for the Application of Fractional Calculus to Viscoelasticity. Journal of Rheology27(3), 201–210 (1983). DOI 10.1122/1.549724. URLhttps://doi.org/10.1122/1.549724

work page doi:10.1122/1.549724 1983
[3]

Electron

Breiten, T., Simoncini, V., Stoll, M.: Low-rank solvers for fractional differential equa- tions. Electron. Trans. Numer. Anal.45, 107–132 (2016)

work page 2016
[4]

Geophysical Journal International13(5), 529–539 (1967)

Caputo, M.: Linear models of dissipation whose q is almost frequency independent—ii. Geophysical Journal International13(5), 529–539 (1967). DOI 10.1111/j.1365-246X. 1967.tb02303.x. URLhttps://doi.org/10.1111/j.1365-246X.1967.tb02303.x

work page doi:10.1111/j.1365-246x 1967
[5]

Chen, J., Liu, F., Anh, V.: Analytical solution for the time-fractional telegraph equation by the method of separating variables. J. Math. Anal. Appl.338(2), 1364–1377 (2008). DOI 10.1016/j.jmaa.2007.06.023. URLhttps://doi.org/10.1016/j.jmaa.2007.06.023

work page doi:10.1016/j.jmaa.2007.06.023 2008
[6]

Colbrook, M.J.: Computing semigroups with error control. SIAM J. Numer. Anal. 60(1), 396–422 (2022). DOI 10.1137/21M1398616. URLhttps://doi.org/10.1137/ 21M1398616

work page doi:10.1137/21m1398616 2022
[7]

Diethelm, K.: The analysis of fractional differential equations,Lecture Notes in Math- ematics, vol. 2004. Springer-Verlag, Berlin (2010). DOI 10.1007/978-3-642-14574-2. A⋆-Product Approach for Nonautonomous Linear FDEs 29 URLhttps://doi.org/10.1007/978-3-642-14574-2. An application-oriented exposi- tion using differential operators of Caputo type

work page doi:10.1007/978-3-642-14574-2 2004
[8]

Journal of Computational and Applied Mathematics 186(2), 482–503 (2006)

Diethelm, K., Ford, J.M., Ford, N.J., Weilbeer, M.: Pitfalls in fast numerical solvers for fractional differential equations. Journal of Computational and Applied Mathematics 186(2), 482–503 (2006). DOI https://doi.org/10.1016/j.cam.2005.03.023. URLhttps: //www.sciencedirect.com/science/article/pii/S0377042705001287

work page doi:10.1016/j.cam.2005.03.023 2006
[9]

Nonlinear Dynam.29(1-4), 3–22 (2002)

Diethelm, K., Ford, N.J., Freed, A.D.: A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dynam.29(1-4), 3–22 (2002). DOI 10.1023/A:1016592219341. URLhttps://doi.org/10.1023/A:1016592219341. Fractional order calculus and its applications 10.NIST Digital Library of Mathematical Functions. http://dlmf.n...

work page doi:10.1023/a:1016592219341 2002
[10]

Pafnuty Publications, Oxford (2014)

Driscoll, T.A., Hale, N., Trefethen, L.N.: Chebfun Guide. Pafnuty Publications, Oxford (2014)

work page 2014
[11]

Ford, N.J., Rodrigues, M.M., Vieira, N.: A numerical method for the fractional Schr¨ odinger type equation of spatial dimension two. Fract. Calc. Appl. Anal.16(2), 454–468 (2013). DOI 10.2478/s13540-013-0028-5. URLhttps://doi.org/10.2478/ s13540-013-0028-5

work page doi:10.2478/s13540-013-0028-5 2013
[12]

Garrappa, R.: Numerical evaluation of two and three parameter Mittag-Leffler functions. SIAM J. Numer. Anal.53(3), 1350–1369 (2015). DOI 10.1137/140971191. URLhttps: //doi.org/10.1137/140971191

work page doi:10.1137/140971191 2015
[13]

Garrappa, R.: Trapezoidal methods for fractional differential equations: theoretical and computational aspects. Math. Comput. Simulation110, 96–112 (2015). DOI 10.1016/ j.matcom.2013.09.012. URLhttps://doi.org/10.1016/j.matcom.2013.09.012

work page doi:10.1016/j.matcom.2013.09.012 2015
[14]

Mathematics6(2), 16 (2018)

Garrappa, R.: Numerical solution of fractional differential equations: A survey and a software tutorial. Mathematics6(2), 16 (2018)

work page 2018
[15]

Garrappa, R., Moret, I., Popolizio, M.: On the time-fractional Schr¨ odinger equation: theoretical analysis and numerical solution by matrix Mittag-Leffler functions. Comput. Math. Appl.74(5), 977–992 (2017). DOI 10.1016/j.camwa.2016.11.028. URLhttps: //doi.org/10.1016/j.camwa.2016.11.028

work page doi:10.1016/j.camwa.2016.11.028 2017
[16]

Garrappa, R., Popolizio, M.: Evaluation of generalized Mittag-Leffler functions on the real line. Adv. Comput. Math.39(1), 205–225 (2013). DOI 10.1007/s10444-012-9274-z. URLhttps://doi.org/10.1007/s10444-012-9274-z

work page doi:10.1007/s10444-012-9274-z 2013
[17]

American Mathematical Monthly377, 1026 (1967)

Gel’fand, I.M., Shilov, G.E.: Generalized functions, volume 1: Properties and oper- ations. American Mathematical Monthly377, 1026 (1967). URLhttps://api. semanticscholar.org/CorpusID:124080985

work page 1967
[18]

Gilles, M.A., Townsend, A.: Continuous analogues of Krylov subspace methods for differential operators. SIAM J. Numer. Anal.57(2), 899–924 (2019). DOI 10.1137/ 18M1177810. URLhttps://doi.org/10.1137/18M1177810

work page doi:10.1137/18m1177810 2019
[19]

Journal of Mathematical Physics56(5), 053503 (2015)

Giscard, P.L., Lui, K., Thwaite, S.J., Jaksch, D.: An exact formulation of the time- ordered exponential using path-sums. Journal of Mathematical Physics56(5), 053503 (2015). DOI 10.1063/1.4920925. URLhttps://doi.org/10.1063/1.4920925

work page doi:10.1063/1.4920925 2015
[20]

Giscard, P.L., Pozza, S.: Lanczos-like algorithm for the time-ordered exponential: the∗- inverse problem. Appl. Math.65(6), 807–827 (2020). DOI 10.21136/AM.2020.0342-19. URLhttps://doi.org/10.21136/AM.2020.0342-19

work page doi:10.21136/am.2020.0342-19 2020
[21]

Giscard, P.L., Thwaite, S.J., Jaksch, D.: Evaluating matrix functions by resummations on graphs: the method of path-sums. SIAM J. Matrix Anal. Appl.34(2), 445–469 (2013). DOI 10.1137/120862880. URLhttps://doi.org/10.1137/120862880

work page doi:10.1137/120862880 2013
[22]

Walk-Sums, Continued Fractions and Unique Factorisation on Digraphs

Giscard, P.L., Thwaite, S.J., Jaksch, D.: Walk-sums, continued fractions and unique factorisation on digraphs (2015). URLhttps://arxiv.org/abs/1202.5523

work page internal anchor Pith review Pith/arXiv arXiv 2015
[23]

Theory and Computation

Higham, N.J.: Functions of Matrices. Theory and Computation. SIAM, Philadelphia (2008)

work page 2008
[24]

SIAM Journal on Matrix Analysis and Applica- tions34(3), 1341–1360 (2013)

Higham, N.J., Lin, L.: An improved Schur–Pad´ e algorithm for fractional powers of a matrix and their Fr´ echet derivatives. SIAM Journal on Matrix Analysis and Applica- tions34(3), 1341–1360 (2013). DOI 10.1137/130906118. URLhttps://doi.org/10. 1137/130906118 30 F. Durastante et al

work page doi:10.1137/130906118 2013
[25]

Pacific J

King, C.F.: A note on Drazin inverses. Pacific J. Math.73(2), 383–390 (1977). URL http://dml.mathdoc.fr/item/1102811925

work page arXiv 1977
[26]

Laskin, N.: Fractional quantum mechanics. Phys. Rev. E62, 3135–3145 (2000). DOI 10.1103/PhysRevE.62.3135. URLhttps://link.aps.org/doi/10.1103/PhysRevE.62. 3135

work page doi:10.1103/physreve.62.3135 2000
[27]

Laskin, N.: Fractional quantum mechanics and L´ evy path integrals. Phys. Lett. A 268(4-6), 298–305 (2000). DOI 10.1016/S0375-9601(00)00201-2. URLhttps://doi. org/10.1016/S0375-9601(00)00201-2

work page doi:10.1016/s0375-9601(00)00201-2 2000
[28]

Laskin, N.: Fractional Schr¨ odinger equation. Phys. Rev. E66, 056108 (2002). DOI 10. 1103/PhysRevE.66.056108. URLhttps://link.aps.org/doi/10.1103/PhysRevE.66. 056108

work page doi:10.1103/physreve.66 2002
[29]

Chaos Solitons Fractals102, 16–28 (2017)

Laskin, N.: Time fractional quantum mechanics. Chaos Solitons Fractals102, 16–28 (2017). DOI 10.1016/j.chaos.2017.04.010. URLhttps://doi.org/10.1016/j.chaos. 2017.04.010

work page doi:10.1016/j.chaos.2017.04.010 2017
[30]

Lubich, C.: Fractional linear multistep methods for Abel-Volterra integral equations of the second kind. Math. Comp.45(172), 463–469 (1985). DOI 10.2307/2008136. URL https://doi.org/10.2307/2008136

work page doi:10.2307/2008136 1985
[31]

Malinowska, A.B., Torres, D.F.M.: Towards a combined fractional mechanics and quantization. Fract. Calc. Appl. Anal.15(3), 407–417 (2012). DOI 10.2478/ s13540-012-0029-9. URLhttps://doi.org/10.2478/s13540-012-0029-9

work page doi:10.2478/s13540-012-0029-9 2012
[32]

Massei, S., Mazza, M., Robol, L.: Fast solvers for two-dimensional fractional diffusion equations using rank structured matrices. SIAM J. Sci. Comput.41(4), A2627–A2656 (2019). DOI 10.1137/18M1180803. URLhttps://doi.org/10.1137/18M1180803

work page doi:10.1137/18m1180803 2019
[33]

BIT61(1), 237–273 (2021)

Massei, S., Robol, L.: Rational Krylov for Stieltjes matrix functions: convergence and pole selection. BIT61(1), 237–273 (2021). DOI 10.1007/s10543-020-00826-z. URL https://doi.org/10.1007/s10543-020-00826-z

work page doi:10.1007/s10543-020-00826-z 2021
[34]

Journal of Mathematical Physics 45(8), 3339–3352 (2004)

Naber, M.: Time fractional Schr¨ odinger equation. Journal of Mathematical Physics 45(8), 3339–3352 (2004). DOI 10.1063/1.1769611. URLhttps://doi.org/10.1063/1. 1769611

work page doi:10.1063/1.1769611 2004
[35]

Mathematics in science and engineering

Podlubny, I.: Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applica- tions. Mathematics in science and engineering. Academic Press, London (1999)

work page 1999
[36]

Linear and Multilinear Algebra0(0), 1–11 (2024)

Pozza, S.: A new closed-form expression for the solution of ODEs in a ring of distribu- tions and its connection with the matrix algebra. Linear and Multilinear Algebra0(0), 1–11 (2024). DOI 10.1080/03081087.2024.2303058. URLhttps://doi.org/10.1080/ 03081087.2024.2303058

work page doi:10.1080/03081087.2024.2303058 2024
[37]

PAMM22, e202200117 (2023)

Pozza, S., Van Buggenhout, N.: A new matrix equation expression for the solution of non-autonomous linear systems of ODEs. PAMM22, e202200117 (2023). DOI 10.1002/PAMM.202200117. URLhttps://onlinelibrary.wiley.com/doi/full/10. 1002/pamm.202200117https://onlinelibrary.wiley.com/doi/abs/10.1002/pamm. 202200117https://onlinelibrary.wiley.com/doi/10.1002/pamm....

work page doi:10.1002/pamm.202200117 2023
[38]

ETNA - Electronic Transactions on Numerical Analysis pp

Pozza, S., Van Buggenhout, N.: A new Legendre polynomial-based approach for non- autonomous linear ODEs. ETNA - Electronic Transactions on Numerical Analysis pp. 292–326 (2024). URLhttps://epub.oeaw.ac.at/?arp=0x003f234e

work page 2024
[39]

Ryckebusch, M., Bouhamidi, A., Giscard, P.L.: A Fr´ echet Lie group on distribu- tions. J. Math. Anal. Appl.546(1), 129195 (2025). DOI https://doi.org/10. 1016/j.jmaa.2024.129195. URLhttps://www.sciencedirect.com/science/article/ pii/S0022247X2401117X

work page arXiv 2025
[40]

Society for Industrial and Applied Mathematics, Philadelphia, PA (2003)

Saad, Y.: Iterative Methods for Sparse Linear Systems, 2nd edn. Society for Industrial and Applied Mathematics, Philadelphia, PA (2003)

work page 2003
[41]

SIAM Rev.58, 377–441 (2016)

Simoncini, V.: Computational methods for linear matrix equations. SIAM Rev.58, 377–441 (2016). DOI 10.1137/130912839. URLhttps://epubs.siam.org/doi/10.1137/ 130912839

work page doi:10.1137/130912839 2016
[42]

Volterra, V., P´ er` es, J.: Le¸ cons sur la composition et les fonctions permutables.´Editions Jacques Gabay, Paris (1928)

work page 1928
[43]

Chaos Solitons Fractals102, 29–46 (2017)

Zhang, Y., Sun, H., Stowell, H.H., Zayernouri, M., Hansen, S.E.: A review of applications of fractional calculus in Earth system dynamics. Chaos Solitons Fractals102, 29–46 (2017). DOI 10.1016/j.chaos.2017.03.051. URLhttps://doi.org/10.1016/j.chaos. 2017.03.051 A⋆-Product Approach for Nonautonomous Linear FDEs 31 Publisher’s NoteSpringer Nature remains ne...

work page doi:10.1016/j.chaos.2017.03.051 2017

[1] [1]

Aceto, L., Novati, P.: Fast and accurate approximations to fractional powers of opera- tors. IMA J. Numer. Anal.42(2), 1598–1622 (2022). DOI 10.1093/imanum/drab002. URLhttps://doi.org/10.1093/imanum/drab002

work page doi:10.1093/imanum/drab002 2022

[2] [2]

Journal of Rheology27(3), 201–210 (1983)

Bagley, R.L., Torvik, P.J.: A Theoretical Basis for the Application of Fractional Calculus to Viscoelasticity. Journal of Rheology27(3), 201–210 (1983). DOI 10.1122/1.549724. URLhttps://doi.org/10.1122/1.549724

work page doi:10.1122/1.549724 1983

[3] [3]

Electron

Breiten, T., Simoncini, V., Stoll, M.: Low-rank solvers for fractional differential equa- tions. Electron. Trans. Numer. Anal.45, 107–132 (2016)

work page 2016

[4] [4]

Geophysical Journal International13(5), 529–539 (1967)

Caputo, M.: Linear models of dissipation whose q is almost frequency independent—ii. Geophysical Journal International13(5), 529–539 (1967). DOI 10.1111/j.1365-246X. 1967.tb02303.x. URLhttps://doi.org/10.1111/j.1365-246X.1967.tb02303.x

work page doi:10.1111/j.1365-246x 1967

[5] [5]

Chen, J., Liu, F., Anh, V.: Analytical solution for the time-fractional telegraph equation by the method of separating variables. J. Math. Anal. Appl.338(2), 1364–1377 (2008). DOI 10.1016/j.jmaa.2007.06.023. URLhttps://doi.org/10.1016/j.jmaa.2007.06.023

work page doi:10.1016/j.jmaa.2007.06.023 2008

[6] [6]

Colbrook, M.J.: Computing semigroups with error control. SIAM J. Numer. Anal. 60(1), 396–422 (2022). DOI 10.1137/21M1398616. URLhttps://doi.org/10.1137/ 21M1398616

work page doi:10.1137/21m1398616 2022

[7] [7]

Diethelm, K.: The analysis of fractional differential equations,Lecture Notes in Math- ematics, vol. 2004. Springer-Verlag, Berlin (2010). DOI 10.1007/978-3-642-14574-2. A⋆-Product Approach for Nonautonomous Linear FDEs 29 URLhttps://doi.org/10.1007/978-3-642-14574-2. An application-oriented exposi- tion using differential operators of Caputo type

work page doi:10.1007/978-3-642-14574-2 2004

[8] [8]

Journal of Computational and Applied Mathematics 186(2), 482–503 (2006)

Diethelm, K., Ford, J.M., Ford, N.J., Weilbeer, M.: Pitfalls in fast numerical solvers for fractional differential equations. Journal of Computational and Applied Mathematics 186(2), 482–503 (2006). DOI https://doi.org/10.1016/j.cam.2005.03.023. URLhttps: //www.sciencedirect.com/science/article/pii/S0377042705001287

work page doi:10.1016/j.cam.2005.03.023 2006

[9] [9]

Nonlinear Dynam.29(1-4), 3–22 (2002)

Diethelm, K., Ford, N.J., Freed, A.D.: A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dynam.29(1-4), 3–22 (2002). DOI 10.1023/A:1016592219341. URLhttps://doi.org/10.1023/A:1016592219341. Fractional order calculus and its applications 10.NIST Digital Library of Mathematical Functions. http://dlmf.n...

work page doi:10.1023/a:1016592219341 2002

[10] [10]

Pafnuty Publications, Oxford (2014)

Driscoll, T.A., Hale, N., Trefethen, L.N.: Chebfun Guide. Pafnuty Publications, Oxford (2014)

work page 2014

[11] [11]

Ford, N.J., Rodrigues, M.M., Vieira, N.: A numerical method for the fractional Schr¨ odinger type equation of spatial dimension two. Fract. Calc. Appl. Anal.16(2), 454–468 (2013). DOI 10.2478/s13540-013-0028-5. URLhttps://doi.org/10.2478/ s13540-013-0028-5

work page doi:10.2478/s13540-013-0028-5 2013

[12] [12]

Garrappa, R.: Numerical evaluation of two and three parameter Mittag-Leffler functions. SIAM J. Numer. Anal.53(3), 1350–1369 (2015). DOI 10.1137/140971191. URLhttps: //doi.org/10.1137/140971191

work page doi:10.1137/140971191 2015

[13] [13]

Garrappa, R.: Trapezoidal methods for fractional differential equations: theoretical and computational aspects. Math. Comput. Simulation110, 96–112 (2015). DOI 10.1016/ j.matcom.2013.09.012. URLhttps://doi.org/10.1016/j.matcom.2013.09.012

work page doi:10.1016/j.matcom.2013.09.012 2015

[14] [14]

Mathematics6(2), 16 (2018)

Garrappa, R.: Numerical solution of fractional differential equations: A survey and a software tutorial. Mathematics6(2), 16 (2018)

work page 2018

[15] [15]

Garrappa, R., Moret, I., Popolizio, M.: On the time-fractional Schr¨ odinger equation: theoretical analysis and numerical solution by matrix Mittag-Leffler functions. Comput. Math. Appl.74(5), 977–992 (2017). DOI 10.1016/j.camwa.2016.11.028. URLhttps: //doi.org/10.1016/j.camwa.2016.11.028

work page doi:10.1016/j.camwa.2016.11.028 2017

[16] [16]

Garrappa, R., Popolizio, M.: Evaluation of generalized Mittag-Leffler functions on the real line. Adv. Comput. Math.39(1), 205–225 (2013). DOI 10.1007/s10444-012-9274-z. URLhttps://doi.org/10.1007/s10444-012-9274-z

work page doi:10.1007/s10444-012-9274-z 2013

[17] [17]

American Mathematical Monthly377, 1026 (1967)

Gel’fand, I.M., Shilov, G.E.: Generalized functions, volume 1: Properties and oper- ations. American Mathematical Monthly377, 1026 (1967). URLhttps://api. semanticscholar.org/CorpusID:124080985

work page 1967

[18] [18]

Gilles, M.A., Townsend, A.: Continuous analogues of Krylov subspace methods for differential operators. SIAM J. Numer. Anal.57(2), 899–924 (2019). DOI 10.1137/ 18M1177810. URLhttps://doi.org/10.1137/18M1177810

work page doi:10.1137/18m1177810 2019

[19] [19]

Journal of Mathematical Physics56(5), 053503 (2015)

Giscard, P.L., Lui, K., Thwaite, S.J., Jaksch, D.: An exact formulation of the time- ordered exponential using path-sums. Journal of Mathematical Physics56(5), 053503 (2015). DOI 10.1063/1.4920925. URLhttps://doi.org/10.1063/1.4920925

work page doi:10.1063/1.4920925 2015

[20] [20]

Giscard, P.L., Pozza, S.: Lanczos-like algorithm for the time-ordered exponential: the∗- inverse problem. Appl. Math.65(6), 807–827 (2020). DOI 10.21136/AM.2020.0342-19. URLhttps://doi.org/10.21136/AM.2020.0342-19

work page doi:10.21136/am.2020.0342-19 2020

[21] [21]

Giscard, P.L., Thwaite, S.J., Jaksch, D.: Evaluating matrix functions by resummations on graphs: the method of path-sums. SIAM J. Matrix Anal. Appl.34(2), 445–469 (2013). DOI 10.1137/120862880. URLhttps://doi.org/10.1137/120862880

work page doi:10.1137/120862880 2013

[22] [22]

Walk-Sums, Continued Fractions and Unique Factorisation on Digraphs

Giscard, P.L., Thwaite, S.J., Jaksch, D.: Walk-sums, continued fractions and unique factorisation on digraphs (2015). URLhttps://arxiv.org/abs/1202.5523

work page internal anchor Pith review Pith/arXiv arXiv 2015

[23] [23]

Theory and Computation

Higham, N.J.: Functions of Matrices. Theory and Computation. SIAM, Philadelphia (2008)

work page 2008

[24] [24]

SIAM Journal on Matrix Analysis and Applica- tions34(3), 1341–1360 (2013)

Higham, N.J., Lin, L.: An improved Schur–Pad´ e algorithm for fractional powers of a matrix and their Fr´ echet derivatives. SIAM Journal on Matrix Analysis and Applica- tions34(3), 1341–1360 (2013). DOI 10.1137/130906118. URLhttps://doi.org/10. 1137/130906118 30 F. Durastante et al

work page doi:10.1137/130906118 2013

[25] [25]

Pacific J

King, C.F.: A note on Drazin inverses. Pacific J. Math.73(2), 383–390 (1977). URL http://dml.mathdoc.fr/item/1102811925

work page arXiv 1977

[26] [26]

Laskin, N.: Fractional quantum mechanics. Phys. Rev. E62, 3135–3145 (2000). DOI 10.1103/PhysRevE.62.3135. URLhttps://link.aps.org/doi/10.1103/PhysRevE.62. 3135

work page doi:10.1103/physreve.62.3135 2000

[27] [27]

Laskin, N.: Fractional quantum mechanics and L´ evy path integrals. Phys. Lett. A 268(4-6), 298–305 (2000). DOI 10.1016/S0375-9601(00)00201-2. URLhttps://doi. org/10.1016/S0375-9601(00)00201-2

work page doi:10.1016/s0375-9601(00)00201-2 2000

[28] [28]

Laskin, N.: Fractional Schr¨ odinger equation. Phys. Rev. E66, 056108 (2002). DOI 10. 1103/PhysRevE.66.056108. URLhttps://link.aps.org/doi/10.1103/PhysRevE.66. 056108

work page doi:10.1103/physreve.66 2002

[29] [29]

Chaos Solitons Fractals102, 16–28 (2017)

Laskin, N.: Time fractional quantum mechanics. Chaos Solitons Fractals102, 16–28 (2017). DOI 10.1016/j.chaos.2017.04.010. URLhttps://doi.org/10.1016/j.chaos. 2017.04.010

work page doi:10.1016/j.chaos.2017.04.010 2017

[30] [30]

Lubich, C.: Fractional linear multistep methods for Abel-Volterra integral equations of the second kind. Math. Comp.45(172), 463–469 (1985). DOI 10.2307/2008136. URL https://doi.org/10.2307/2008136

work page doi:10.2307/2008136 1985

[31] [31]

Malinowska, A.B., Torres, D.F.M.: Towards a combined fractional mechanics and quantization. Fract. Calc. Appl. Anal.15(3), 407–417 (2012). DOI 10.2478/ s13540-012-0029-9. URLhttps://doi.org/10.2478/s13540-012-0029-9

work page doi:10.2478/s13540-012-0029-9 2012

[32] [32]

Massei, S., Mazza, M., Robol, L.: Fast solvers for two-dimensional fractional diffusion equations using rank structured matrices. SIAM J. Sci. Comput.41(4), A2627–A2656 (2019). DOI 10.1137/18M1180803. URLhttps://doi.org/10.1137/18M1180803

work page doi:10.1137/18m1180803 2019

[33] [33]

BIT61(1), 237–273 (2021)

Massei, S., Robol, L.: Rational Krylov for Stieltjes matrix functions: convergence and pole selection. BIT61(1), 237–273 (2021). DOI 10.1007/s10543-020-00826-z. URL https://doi.org/10.1007/s10543-020-00826-z

work page doi:10.1007/s10543-020-00826-z 2021

[34] [34]

Journal of Mathematical Physics 45(8), 3339–3352 (2004)

Naber, M.: Time fractional Schr¨ odinger equation. Journal of Mathematical Physics 45(8), 3339–3352 (2004). DOI 10.1063/1.1769611. URLhttps://doi.org/10.1063/1. 1769611

work page doi:10.1063/1.1769611 2004

[35] [35]

Mathematics in science and engineering

Podlubny, I.: Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applica- tions. Mathematics in science and engineering. Academic Press, London (1999)

work page 1999

[36] [36]

Linear and Multilinear Algebra0(0), 1–11 (2024)

Pozza, S.: A new closed-form expression for the solution of ODEs in a ring of distribu- tions and its connection with the matrix algebra. Linear and Multilinear Algebra0(0), 1–11 (2024). DOI 10.1080/03081087.2024.2303058. URLhttps://doi.org/10.1080/ 03081087.2024.2303058

work page doi:10.1080/03081087.2024.2303058 2024

[37] [37]

PAMM22, e202200117 (2023)

Pozza, S., Van Buggenhout, N.: A new matrix equation expression for the solution of non-autonomous linear systems of ODEs. PAMM22, e202200117 (2023). DOI 10.1002/PAMM.202200117. URLhttps://onlinelibrary.wiley.com/doi/full/10. 1002/pamm.202200117https://onlinelibrary.wiley.com/doi/abs/10.1002/pamm. 202200117https://onlinelibrary.wiley.com/doi/10.1002/pamm....

work page doi:10.1002/pamm.202200117 2023

[38] [38]

ETNA - Electronic Transactions on Numerical Analysis pp

Pozza, S., Van Buggenhout, N.: A new Legendre polynomial-based approach for non- autonomous linear ODEs. ETNA - Electronic Transactions on Numerical Analysis pp. 292–326 (2024). URLhttps://epub.oeaw.ac.at/?arp=0x003f234e

work page 2024

[39] [39]

Ryckebusch, M., Bouhamidi, A., Giscard, P.L.: A Fr´ echet Lie group on distribu- tions. J. Math. Anal. Appl.546(1), 129195 (2025). DOI https://doi.org/10. 1016/j.jmaa.2024.129195. URLhttps://www.sciencedirect.com/science/article/ pii/S0022247X2401117X

work page arXiv 2025

[40] [40]

Society for Industrial and Applied Mathematics, Philadelphia, PA (2003)

Saad, Y.: Iterative Methods for Sparse Linear Systems, 2nd edn. Society for Industrial and Applied Mathematics, Philadelphia, PA (2003)

work page 2003

[41] [41]

SIAM Rev.58, 377–441 (2016)

Simoncini, V.: Computational methods for linear matrix equations. SIAM Rev.58, 377–441 (2016). DOI 10.1137/130912839. URLhttps://epubs.siam.org/doi/10.1137/ 130912839

work page doi:10.1137/130912839 2016

[42] [42]

Volterra, V., P´ er` es, J.: Le¸ cons sur la composition et les fonctions permutables.´Editions Jacques Gabay, Paris (1928)

work page 1928

[43] [43]

Chaos Solitons Fractals102, 29–46 (2017)

Zhang, Y., Sun, H., Stowell, H.H., Zayernouri, M., Hansen, S.E.: A review of applications of fractional calculus in Earth system dynamics. Chaos Solitons Fractals102, 29–46 (2017). DOI 10.1016/j.chaos.2017.03.051. URLhttps://doi.org/10.1016/j.chaos. 2017.03.051 A⋆-Product Approach for Nonautonomous Linear FDEs 31 Publisher’s NoteSpringer Nature remains ne...

work page doi:10.1016/j.chaos.2017.03.051 2017