Abstract
Compartmental models based on coupled differential equations of fractional order have been widely employed in the literature for modeling. An abstraction of these models is given by a system with polynomial vector field. We investigate the use of power series for solving generic polynomial differential equations in any dimension, with Caputo fractional derivative. As is well known, power series convert a continuous formulation into a discrete system of difference equations, which are easily solved by recursion. The novelty of this paper is that we rigorously prove that the series converge on a neighborhood of the initial instant, which is an analogue of the Cauchy–Kovalevskaya theorem. Besides, these series are proved to be continuous with respect to the fractional index. For applications, a general-purpose symbolic implementation of truncated power series is developed, and its execution is illustrated for the fractional SIR epidemiological model.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
When modeling the dynamics of a phenomenon in a population, the interaction between the elements must be considered. For example, when an epidemic is happening in a region (flu, COVID-19, etc.), contacts between susceptible and infected citizens are crucial, because they cause the transmission of the virus (Lotfi et al. 2020). For social behaviors, such as market choices, ideology, criminality, or health habits, contacts with peers are also important, due to attitudes of imitation, stimulation, pressure, persuasion, etc. Harkins et al. (2017), Esiri (2016), and Blanchower et al. (2009). Elements in the population could be other than persons; for example, for an ill patient undergoing a clinical follow-up, the dynamics of the disease in the body and the interaction between the virus and uninfected cells are of relevance (Masenga et al. 2023).
Compartmental models consider that a population is divided into different subgroups, according to certain conditions. Elements interact and move between the compartments. Since rates of flow need to be described, differential equations play a key role for compartmental modeling. Interactions are formulated by nonlinear terms, specifically products (multiplications). Simulation of the differential equations, for certain input values, permits understanding the dynamics of the phenomenon under investigation. Some examples of ordinary differential equation models with division into compartments are Brauer et al. (2008), Cooper et al. (2020), Santonja et al. (2010), Acedo et al. (2010a), Cervello et al. (2014), Sánchez et al. (2011), and Calatayud and Jornet (2020). These works deal with Epidemiology, COVID-19, alcohol consumption, respiratory syncytial virus, mobile telecommunications, cocaine use, and excess weight, respectively.
The word “ordinary” refers to the standard derivative. There are operators of fractional order, that extend the integer-order derivative. The literature is vast and many definitions of fractional derivatives have been given (De Oliveira and Tenreiro Machado 2014; Ortigueira and Machado 2015). In general, the operators are based on a convolution with respect to a kernel, so that history effects appear in the model. Some sort of continuous past delay is incorporated. There are a lot of contributions in the literature on compartmental models based on fractional differential equations; some examples are Carvalho et al. (2020), Carvalho et al. (2018), Martinez et al. (2021), Pinto and Carvalho (2017), Biala and Khaliq (2021), Alshomrani et al. (2021), Ndaïrou et al. (2021), and Area et al. (2015), for HIV, COVID-19 and Ebola dynamics, respectively.
Thus, this exposition motivates the study of fractional-order compartmental models, to extend the integer-order counterpart. Specifically, in this paper, we aim at investigating systems of fractional differential equations with polynomial vector field, in the Caputo-derivative sense. The polynomial is multivariate and depends on the state variable (autonomous equation). Polynomial expressions arise in many compartmental models, because nonlinearities are products (i.e., monomials) to account for contacts that occur homogeneously in the region. Thus, the problem proposed in this article is an abstraction of usual compartmental systems.
A great deal of research work in applied mathematics consists in the investigation of explicit, semi-explicit, and numerical solutions to physical models. Our research in this paper is focused on power-series solutions (i.e., analytic solutions). Power series possess appealing properties in terms of algebraic manipulations. When a power series is put into a fractional or ordinary differential equation, one derives an associated discrete model for the coefficients of the expansion. The difference equations can easily be solved by recursion, explicitly or in the computer, so that a simple functional representation of the model’s response is obtained. Although in practice the method works, a difficult issue is to establish convergence theoretically, at least around the initial condition.
Of course, the use of power series for differential equations is not new (Teschl 2012, Chapter 4). Many linear ordinary differential equations, such as Airy, Hermite, Legendre, etc. in mathematical physics, have solutions expressed in terms of power series (Fröbenius method) (Koekoek et al. 2010). Even when stochastic effects are incorporated (Calatayud Gregori et al. 2020; Jornet 2021). Nonlinear ordinary differential equations may be solved with power series too Srivastava et al. (2021). The old Cauchy–Kovalevskaya theorem (Himonas and Petronilho 2020) ensures, essentially, that analytic inputs imply analytic output. This theorem applies for a certain class of partial differential equations as well. For fractional differential equations, specific applications of power series are available in the literature; for example, Area and Nieto (2021, 2023), Jornet (2023) for fractional logistic equations, and Villafuerte (2023) for random linear equations.
To the best of my knowledge, the generic fractional model and the mathematical methodology of this paper have not previously been considered in the literature. Only particular cases have been addressed before (Srivastava et al. 2021; Area and Nieto 2021, 2023; Jornet 2023), which are the bases for this submitted contribution. For systems of fractional differential equations with polynomial vector field, we rigorously investigate the convergence of the power series. This provides a foundation for applications. Besides, the continuity of the model with respect to the fractional order is studied. Finally, a general-purpose symbolic implementation of the method is tackled.
The novelties of this work are the following:
-
The mathematical proof of the convergence of the power series for generic Caputo fractional differential equations with polynomial input function. This is an analogue of the Cauchy–Kovalevskaya theorem. In the literature, only applications without a rigorous analysis of convergence (Area and Nieto 2023), or specific models (such as logistic) with an analysis of convergence (Area and Nieto 2021), have been addressed. Hence, the present paper is an advance in the mathematical analysis of fractional systems.
-
The generic implementation of the power-series solutions to compartmental models. Hence, the present paper is an advance in the computational analysis of fractional systems.
The paper is organized as follows. In Sect. 2, the fractional system of study is introduced. The notation is given and some examples provide the link with compartmental models. In Sect. 3, the formal power-series solution of the system is derived, by obtaining the discrete model for the expansion coefficients. These parts are a preparation (concepts, notation, etc.) for Sect. 4, where convergence of the power series is proved. This is the main novelty of the paper. Later, in Sect. 5, continuity with respect to the fractional index is proved. In Sect. 6, a computer code for the methodology is implemented, with some examples of execution. Finally, in Sect. 7, a summary with the main aspects of the article and possible extensions are presented.
2 Context
We work with a fractional-order differential equation model of the form
The state variable of the system (1) is x, which is a function of time \(t\ge 0\) over the temporal range \({\mathcal {J}}=[0,T]\), where \(T>0\) is a time horizon. Its image takes vector values, on \({\mathbb {R}}^d\), for dimension \(d\ge 1\). By components, it is expressed as \(x=(x_1,\ldots ,x_d)\). The scalar case corresponds to \(d=1\).
The differential present in (1) is of fractional order, through the Caputo operator \({\mathcal {D}}^\alpha \) (Gerasimov 1948; De Oliveira and Tenreiro Machado 2014; Ortigueira and Machado 2015; Caputo 1967). For order \(0<\alpha <1\), it is defined as
where f is an absolutely continuous function with the first-order derivative \(f'\) and \(\varGamma \) is the gamma function. It extends classical differentiation, which is retrieved for \(\alpha \rightarrow 1^{-}\). The convolution in (2) includes memory effects somehow, because \(f'\) depends on \(\tau <t\).
The autonomous vector field in (1), p, is of polynomial type, with multivariate evaluation, from \({\mathbb {R}}^d\) to \({\mathbb {R}}^d\). Its coordinates are \(p=(p_1,\ldots ,p_d)\). Componentwise, model (1) reads as
The data in (1) are the initial state, x(0), and the real coefficients of p. For Caputo fractional differential equations, initial conditions can be handled analogously to the classical integer-order counterpart.
To express p, we use multi-index notation (Dunkl and Xu 2014). A multi-index is a vector of the form \(I=(i_1,\ldots ,i_d)\), formed by integer elements \(i_1,\ldots ,i_d\ge 0\). Associated notation is
to compactly denote monomials, and
to compute its degree. In consequence, each coordinate of p may be expressed as
where \(a_{j,I}\in {\mathbb {R}}\) are the coefficients, \(x^I\) are the monomials, \(\delta _j\) is the degree, and the sum runs over all multi-indices I with associated degree from 0 to \(\delta _j\). The number of multi-indices I is \(\left( {\begin{array}{c}d+\delta _j\\ d\end{array}}\right) =(d+\delta _j)!/(d!\delta _j!)\). In vector notation, if \(a_I=(a_{1,I},\ldots ,a_{d,I})\), then
Although clear in the formulation, the multi-index notation is cumbersome in practice. For clarity, one may use the graded lexicographic order: given two multi-indices I and L of length d, \(I > L\) if and only if \(|I| \ge |L|\) and the first nonzero entry in the difference, \(I - L\), is positive. For example, when \(d=4\), the multi-indices in ascending order are (0, 0, 0, 0), (0, 0, 0, 1), (0, 0, 1, 0), (0, 1, 0, 0), (1, 0, 0, 0), (0, 0, 0, 2), (0, 0, 1, 1), (0, 0, 2, 0), (0, 1, 0, 1), (0, 1, 1, 0), (0, 2, 0, 0), (1, 0, 0, 1), (1, 0, 1, 0), (1, 1, 0, 0), (2, 0, 0, 0), etc. In terms of monomials, these multi-indices translate into 1, \(x_4\), \(x_3\), \(x_2\), \(x_1\), \(x_4^2\), \(x_3x_4\), \(x_3^2\), \(x_2x_4\), \(x_2x_3\), \(x_2^2\), \(x_1x_4\), \(x_1x_3\), \(x_1x_2\), \(x_1^2\), etc.
The investigation of (1) relies on its widespread use in the literature for modeling. Many models have polynomial vector fields, to account for interactions between a finite number of dynamic subgroups.
Example 1
The fractional logistic equation (Area and Nieto 2021)
is a particular case of (1), for dimension \(d=1\) and \(p(x)=x-x^2\). The degree is \(\delta =2\), and the coefficients are \(a_0=0\), \(a_1=1\) and \(a_2=-1\). Logistic growth can be generalized by incorporating the Allee effect (Area and Nieto 2023)
where \(\theta \in (0,1)\). The parameters are \(d=1\), \(\delta =3\), \(a_0=0\), \(a_1=-\theta \), \(a_2=\theta +1\), and \(a_3=-1\).
Example 2
The fractional SIR (susceptible-infected-recovered) model (Alshomrani et al. 2021; Area et al. 2015)
The parameters are positive and are raised to the power of \(\alpha \) for dimensional consistency (the units are \(\text {time}^{-\alpha }\)). Another approach to be consistent with dimensionality would consist in introducing an artificial constant, \(\textrm{d}/\textrm{d}t \rightarrow (1/\sigma ^{1-\alpha }){\mathcal {D}}^\alpha \) (Gómez-Aguilar et al. 2014; Popović et al. 2010). The written SIR model is a particular case of (1) for dimension \(d=3\)
and
The degrees are \(\delta _1=2\), \(\delta _2=2\) and \(\delta _3=1\), and the nonzero coefficients are \(a_{1,(1,1,0)}=-\beta ^\alpha \), \(a_{2,(1,1,0)}=\beta ^\alpha \), \(a_{2,(0,0,1)}=-\gamma ^\alpha \) and \(a_{3,(0,0,1)}=\gamma ^\alpha \). The first subscript denotes \(j\in \{1,2,3\}\) and the second subscript is the multi-index \(I=(i_1,i_2,i_3)\).
Example 3
A fractional model for HIV (human immunodeficiency virus) is Pinto and Carvalho (2017)
Except \(\eta \in (0,1)\), the parameters are positive and are raised to the power of \(\alpha \) for dimensional consistency (the units are \(\text {time}^{-\alpha }\)). For the details on the formulation, there reader is referred to Pinto and Carvalho (2017). Essentially, there are seven compartments for a patient under clinical follow-up: the uninfected CD 4\(^{+}\) T cells, the latently infected CD 4\(^{+}\) T cells, the productively infected CD 4\(^{+}\) T cells, the uninfected macrophages, the infected macrophages, the virus, and the CTLs. In our notation
and
The degrees are \(\delta _1=2\), \(\delta _2=2\), \(\delta _3=2\), \(\delta _4=2\), \(\delta _5=2\), \(\delta _6=1\), and \(\delta _7=2\). The nonzero coefficients are \(a_{1,(0,0,0,0,0,0,0)}=s^\alpha \), \(a_{1,(1,0,0,0,0,0,0)}=-\mu _U^\alpha \), \(a_{1,(1,0,0,0,0,1,0)}=-k_1^\alpha \), \(a_{2,(1,0,0,0,0,1,0)}=k_1^\alpha \eta \), \(a_{2,(0,1,0,0,0,0,0)}=-a_L^\alpha -\mu _L^\alpha \), \(a_{3,(1,0,0,0,0,1,0)}=k_1^\alpha (1-\eta )\), \(a_{3,(0,1,0,0,0,0,0)}=a_L^\alpha \), \(a_{3,(0,0,1,0,0,0,0)}=-\delta _F^\alpha \), etc. The first subscript denotes \(j\in \{1,2,3,4,5,6,7\}\) and the second subscript is the multi-index \(I=(i_1,i_2,i_3,i_4,i_5,i_6,i_7)\).
3 Formal power series: from continuous to discrete
As widely studied in the literature, to solve (1), a power series at \(t^\alpha \) is proposed
where \(t\in {\mathcal {J}}\) and \(b_n\in {\mathbb {R}}^d\). If we denote \(b_n=(b_{1,n},\ldots ,b_{d,n})\), then (5) is given by components as
For now, the series (5) is formal, meaning that one uses the operations associated with power series, but ignoring whether convergence holds on an interval \([0,\epsilon )\subseteq {\mathcal {J}}\), for some \(\epsilon >0\). See Area and Nieto (2021, 2023) for the algebra associated with fractional power series. The aim of this paper is to show that (5), indeed, converges on a relative neighborhood of 0 in \({\mathcal {J}}\). To solve (1) and analyze the convergence of (5), we need to see how the coefficients \(b_n\) behave.
With the fractional differentiation of powers
the power series is differentiated, term by term, as
When (5) converges on \([0,\epsilon )\), the identity (7) is not merely formal, but holds with convergence on \([0,\epsilon )\) too.
On the other hand, the power series of p(x) is required. Taking the form (4) into account, we calculate \(x_j^{i_j}(t)\) and \(x^I(t)\). Recall the expression of the Cauchy product, defined by a discrete convolution. Some examples are the following:
etc. For an arbitrary power \(i_j\), we have
We employ the shorthand notation
so that (8) becomes
For a multi-index power, by applying (10)
If
we can express (11) as
If (5) converges on \([0,\epsilon )\), it is well known that the Cauchy product (13) converges on \([0,\epsilon )\), as well.
Then, the coefficients of the power series must be equal. After isolating \(b_{j,n+1}\), we arrive at the following scalar recursive equations:
for \(j=1,\ldots ,d\) and \(n=0,1,2,\ldots \), with initial data \(b_{j,0}=x_j(0)\). For each n, one runs (14) for all \(j=1,\ldots ,d\), before moving to \(n+1\). That is, the loop over j is nested in the loop over n. Recall that \(C_I\) is defined by (12) and (9). Given (14), the main question is whether (5) converges around 0. Under convergence, the previous series manipulations are rigorous and we obtain the proper solution to (1).
4 Convergence of the power series
The aim is to prove that (5) converges around 0. We first need to upper bound the absolute value of the coefficients, \(|b_{j,n+1}|\), given by (14). These depend on (12) and (9).
By the triangular inequality on (9)
As a consequence, by the triangular inequality on (12)
On the other hand, from a certain \(n\ge 0\) (no need to specify it for easiness of notation), the monotone condition
holds (Kershaw 1983). In fact
because it is well known that \(\varGamma (y+\alpha )\sim \varGamma (y)y^\alpha \) when \(y\rightarrow \infty \).
By (15) and (16), the coefficients in (14) are upper bounded as follows:
From (17), we build a “majorizing” series for \(\sum _{n=0}^\infty |b_{j,n}| t^{\alpha n}\). Let
and
be new coefficients, for \(j=1,\ldots ,d\) and \(n\ge 0\). Again, for each n, one runs (18) for all \(j=1,\ldots ,d\), before advancing to \(n+1\). By induction on n, it is trivially justified that
for all subscripts. Let
and with vector notation
where \(\psi =(\psi _1,\ldots ,\psi _d)\), \(z\in {\mathbb {R}}\) and \(h_n=(h_{1,n},\ldots ,h_{d,n})\in {\mathbb {R}}^d\). For now, the power series (20) is formal. We need to check its convergence on a neighborhood of 0. In such a case, the power series (5), with coefficients (14), will be well defined around 0, by (19), as wanted.
Notice that, by (20), (18) and the Cauchy product for multi-indices [as in (13)], it holds
Viewing this identity (22) for \(\psi \) as a functional root, we consider the maps
for \(j=1,\ldots ,d\), which are the coordinates of a map
By (22), \(\psi \) in (21) and its coefficients (18) are characterized by \(\psi (0)=h_0\) and \(\phi (z,\psi (z))=0\). We have
and the Jacobian
with
Here, \(*\) denotes any expression and \(\textrm{Id}_d\) is the identity matrix of size \(d\times d\). By the analytic version of the implicit-function theorem (Kaup and Kaup 2011, Section 8, Chapter 0), there exists a unique analytic function \(w\equiv w(z)\) on a neighborhood \((-\mu ,\mu )\) of zero, \(\mu >0\), \(w:(-\mu ,\mu )\rightarrow {\mathbb {R}}^d\), such that \(w(0)=d_0\) and \(\phi (z,w(z))=0\). Then, \(\psi =w\) is convergent on \((-\mu ,\mu )\). By (19)
for \(t\in [0,\mu ^{1/\alpha })\). In conclusion, the power series (5), with coefficients (14), is well defined around 0 and defines a suitable solution to (1).
5 Consistency with the integer-order model
When \(\alpha =1\), problem (1) is no more fractional, since \({\mathcal {D}}^1 x(t)=x'(t)\). Let us denote the coefficients in (5) and (14) as \(b_{j,n}|_{\alpha }\), that is, we make the fractional-order explicit. The aim is to prove consistency, namely
This last series is the classical integer-order power-series solution.
We distinguish two situations, taking into account the role of \(\alpha \):
- Case 1:
-
The coefficients \(a_{j,I}\) of \(p_j\), see (4) for \(j=1,\ldots ,d\), do not depend on \(\alpha \); that is, they do not change with variations of \(\alpha \).
The coefficients \(b_{j,n}|_{\alpha }\) satisfy (19), i.e., \(|b_{j,n}|_{\alpha }|\le h_{j,n}\), where the terms \(h_{j,n}\) are defined by (18) and are independent of \(\alpha \). Independence holds due to (16).
On the other hand, \(t^{\alpha n}\le t^n\) if \(t\ge 1\), and \(t^{\alpha n}\le (\sqrt{t})^n\) if \(t<1\) and \(\alpha >1/2\).
In consequence, if (21) converges on \((-\mu ,\mu )\) and \(t,\sqrt{t}\in (-\mu ,\mu )\), then \(|b_n||t|^{\alpha n}\) is upper bounded by summable \(\alpha \)-independent terms. By the dominated convergence theorem (Rudin 1976, result 11.32, page 321), which permits interchanging \(\lim _{\alpha \rightarrow 1^-}\) and \(\sum _{n=0}^\infty \), (23) holds.
- Case 2:
-
The coefficients \(a_{j,I}\) of \(p_j\), see (4) for \(j=1,\ldots ,d\), depend on \(\alpha \), as \(a_{j,I}=\overline{a}_{j,I}^\alpha \); this is often the case when dimensional consistency is sought, where \(\overline{a}_{j,I}>0\) is independent of \(\alpha \).
The proof is very similar to Case 1. Actually, we proceed to reduce Case 2 to Case 1.
We start by \(|b_{j,n}|_{\alpha }|\le h_{j,n}\), see (19), but now the terms \(h_{j,n}\) defined by (18) are functions of \(\alpha \). We bound \(|a_{j,I}|\le \tilde{a}_{j,I}\) for all subindices, where the coefficients \(\tilde{a}_{j,I}\) are independent of \(\alpha \). This can be done by setting \(\tilde{a}_{j,I}=\max \{\overline{a}_{j,I},\sqrt{\overline{a}_{j,I}}\}\), considering \(\alpha >1/2\) and whether \(\overline{a}_{j,I}\le 1\) or \(\overline{a}_{j,I}>1\). Define new \(\tilde{h}_{j,n}\) (independent of \(\alpha \)), as (18) but with \(\tilde{a}_{j,I}\) instead
The inequalities \(|b_{j,n}|_{\alpha }|\le h_{j,n}\le \tilde{h}_{j,n}\) hold, by induction on n. The previous case and calculations then apply, with the “majorizing” series \(\tilde{\psi }_j(z)=\sum _{n=0}^\infty \tilde{h}_{j,n}z^n\) (i.e., with tildes).
6 Symbolic implementation
We implement the computation of a truncated sum of (5)
We use computer algebra to obtain the function \(t\mapsto x^{[N]}(t)\), where the variable t is symbolic. If a numerical value is required, then t is substituted accordingly.
Symbolic or algebraic computation is concerned with exact operations and expressions, with variables that have no given value (symbols) or are fractions (no decimal numbers).
A general-purpose algorithm is developed in Mathematica® (Wolfram Research, Inc. 2020) and illustrated here. The function is expansionFractionalCompartmentODE, with inputs:
-
d: The dimension \(d\ge 1\) of the system. It has to be a number.
-
NN: The truncation index N, which is repeated in writing to distinguish from a built-in command N of the software. It has to be a number.
-
p: The polynomial p, as a function of vectors in \({\mathbb {R}}^d\), defined through p[v_]:=...,...,.... In the dots, one writes the equations of the compartments, in terms of the components v[[1]], v[[2]], etc. The coefficients in p can be symbolic or numeric.
-
alpha: The fractional order \(\alpha \in (0,1]\). It can be symbolic or numeric.
-
x0: The initial condition \(x(0)\in {\mathbb {R}}^d\). It can be symbolic or numeric.
The output is \(x^{[N]}(t)\), in terms of the symbolic variable t.
The code is the following:
We show some executions for the SIR model of Example 2. This is a very important compartmental system, originated in the work (Kermack and McKendrick 1927), and closed-form solutions are being investigated for it and its extensions (Srivastava et al. 2021; Harko et al. 2014; Heng and Althaus 2020; Acedo et al. 2010b).
Example 4
Set the following variables:
We have
This example for integer order was conducted in Srivastava et al. (2021) too.
Example 5
Set the following variables:
We have
We report some timings. The computer is an Intel® Core\(^{\textrm{TM}}\) i7 CPU 2.9 GHz, in Windows 10. For order \(N=20\), if beta = 1/2, gamma = 1/9, s0 = 99/100, i0 = 1/100, r0 = 0 and alpha = 1/2 are specified with symbolic fractions, then 4.5 s are required to calculate the three-coordinates function \(x^{[20]}(t)\). However, if those parameters are specified with decimals, numerically, then the time is 0.02 s. With decimal numbers, other timings are 0.08 s for \(N=50\), 0.35 s for \(N=100\), 2.2 s for \(N=200\), and 4.5 s for \(N=250\). Once the symbolic expression \(x^{[250]}(t)\) is computed, the evaluation at a specific numerical value of t lasts around 0.0005 s.
Example 6
Set the following variables:
We have
We report some timings, with the previous computer’s specifications. Compared to the earlier example, the calculations are more costly, because all of the variables are symbolic. For \(N=8\), the code runs in 25 s; for \(N=7\), in 2.6 s; and for \(N=6\), in 0.5 s.
7 Conclusions
Compartmental models, based on coupled differential equations, have been used in a lot of contributions for dynamical modeling. The fields of application include Epidemiology, Ecology and Sociology, for example, because individuals or species interact (nonlinear terms) and transmit diseases, behaviors, etc. This paper is a contribution to the study of this type of models, with fractional derivatives incorporated. We used the Caputo definition.
We investigated the use of power-series expansions to find closed-form expressions for the solution. Of course, the applications of power series for differential equations, even nonlinear, are well known. However, to my knowledge, there is no study in the literature that conducts the presented mathematical treatment for a generic fractional model of the form (1), in any state dimensionality.
The coefficients of the formal power series, for the candidate solution, satisfy difference equations. This process may be viewed as a conversion from a continuous into a discrete model, as is well known in the literature. After introducing the necessary concepts and notation, we rigorously proved that the power series converge on a neighborhood of the initial instant zero. We based on “majorizing” series and the analytic implicit-function theorem, dealing with the cumbersome notation from the Cauchy products. Key facts were the monotone condition of the gamma function and the absolute-value inequality met by the Cauchy products. On the other hand, continuity with respect to the fractional order was established, by applying the dominated convergence theorem. Two scenarios were distinguished, depending on whether the coefficients of the polynomial were related to the fractional parameter or not. This result showed consistency between the fractional- and the integer-order formulations.
The series, which must be truncated sums for applications, were implemented in a computer-algebra software. The code was explicitly given in the article. The output of our developed function depends on the symbolic time variable. Some executions were exhibited and commented for the SIR epidemiological model of fractional order, with different symbolic and numeric parameters (force of infection, recovery rate, initial conditions, and fractional index). Explicit sums and CPU times were reported. Other examples, with other inputs, would be analogous.
Some extensions of this paper could be the following:
-
The estimation of the radius of convergence for the power-series solution (5). The analytic version of the implicit-function theorem only provides the existence of a neighborhood of analyticity. The estimation would likely require obtaining more information about the “majorizing” series \(\psi \), through the identity (22). For example, the scalar case \(d=1\) of polynomial degree \(\delta =2\), which corresponds to a logistic equation, gives an explicit root \(\psi (z)\) in (22), and the largest interval of convergence for \(\psi \) could be calculated or estimated.
-
The estimation of the rate of convergence for the power-series solution (5). Also, the rate of convergence as the fractional order tends to integer-order one. Both the implicit-function theorem and the dominated convergence theorem do not give that information. Probably, a proof with specific inequalities would be required.
-
A more efficient implementation of the symbolic code presented for the truncated power series (24).
-
The analysis of other problems (1) with a general fractional order \(\alpha \in (0,\infty )\) or with alternative fractional derivatives.
-
The analysis of problems \({\mathcal {D}}^\alpha x=f(x)\), where f is any analytic function. If f is not a polynomial, then one could use a polynomial approximation \(f\approx p\) and solve \({\mathcal {D}}^\alpha \tilde{x}=p(\tilde{x})\), by truncating the Taylor series, but the new analytic solution \(\tilde{x}\) would be biased. The methodology of the paper does not directly apply to \({\mathcal {D}}^\alpha x=f(x)\).
-
The investigation of the corresponding stochastic problem, where x(0) and/or coefficients of the polynomial p are random variables. The power series may be of use for uncertainty-propagation computations, but mean-square convergence should be established to ensure convergence for the second-order statistics.
Data Availability
No data were used for this study.
References
Acedo L, Díez-Domingo J, Moraño JA, Villanueva RJ (2010) Mathematical modelling of respiratory syncytial virus (RSV): vaccination strategies and budget applications. Epidemiol Infect 138(6):853–860
Acedo L, González-Parra G, Arenas AJ (2010) An exact global solution for the classical SIRS epidemic model. Nonlinear Anal Real World Appl 11(3):1819–1825
Alshomrani AS, Ullah MZ, Baleanu D (2021) Caputo SIR model for COVID-19 under optimized fractional order. Adv Differ Equ 1:185
Area I, Nieto JJ (2021) Power series solution of the fractional logistic equation. Phys A 573:125947
Area I, Nieto JJ (2023) On the fractional Allee logistic equation in the Caputo sense. Ex Counterex 4:100121
Area I, Batarfi H, Losada J, Nieto JJ, Shammakh W, Torres Á (2015) On a fractional order Ebola epidemic model. Adv Differ Equ 1:1–12
Biala TA, Khaliq AQM (2021) A fractional-order compartmental model for the spread of the COVID-19 pandemic. Commun Nonlinear Sci Numer Simul 98:105764
Blanchower DG, Oswald AJ, Landeghem BV (2009) Imitative obesity and relative utility. J Eur Econ Assoc 7:528–538
Brauer F, van den Driessche P, Wu J (2008) Mathematical epidemiology. Springer, Berlin, pp 19–79
Calatayud J, Jornet M (2020) Mathematical modeling of adulthood obesity epidemic in Spain using deterministic, frequentist and Bayesian approaches. Chaos Soliton Fract 140:110179
Calatayud Gregori J, Cortés JC, Jornet Sanz M (2020) Beyond the hypothesis of boundedness for the random coefficient of Airy, Hermite and Laguerre differential equations with uncertainties. Stoch Anal Appl 38(5):875–885
Caputo M (1967) Linear model of dissipation whose Q is almost frequency independent-II. Geophys. J R Astron Soc 13:529–539
Carvalho AR, Pinto C, Baleanu D (2018) HIV/HCV coinfection model: a fractional-order perspective for the effect of the HIV viral load. Adv Differ Equ 1:1–22
Carvalho AR, Pinto CM, de Carvalho JM (2020) Fractional model for type 1 diabetes. Math Model Optim Eng Probl 175–185
Cervello R, Cortés JC, Santonja FJ, Villanueva RJ (2014) The dynamics over the next few years of the Spanish mobile telecommunications market share: a mathematical modelling approach. Math Comput Model Dyn 20(6):557–565
Cooper I, Mondal A, Antonopoulos CG (2020) A SIR model assumption for the spread of COVID-19 in different communities. Chaos Solitons Fractals 139:110057
De Oliveira EC, Tenreiro Machado JA (2014) A review of definitions for fractional derivatives and integral. Math Probl Eng 2014:238459
Dunkl C, Xu Y (2014) Orthogonal polynomials of several variables. Encyclopedia of Mathematics and its Applications, 2nd edn, vol 155, Cambridge University Press, Cambridge
Esiri MO (2016) The influence of peer pressure on criminal behaviour. J Humanit Soc Sci 21(1):08–14
Gerasimov AN (1948) A generalization of linear laws of deformation and its application to problems of internal friction. Akad Nauk SSSR Prikladnaya Matematika i Mekhanika 12:251–259 ((in Russian))
Gómez-Aguilar JF, Razo-Hernández R, Granados-Lieberman D (2014) A physical interpretation of fractional calculus in observables terms: analysis of the fractional time constant and the transitory response. Rev Mex Fís 60:32–38
Harkins SG, Williams KD, Burger J (2017) The Oxford handbook of social influence. Oxford University Press, Oxford
Harko T, Lobo FSN, Mak MK (2014) Exact analytical solutions of the Susceptible-Infected-Recovered (SIR) epidemic model and of the SIR model with equal death and birth rates. Appl Math Comput 236:184–194
Heng K, Althaus CL (2020) The approximately universal shapes of epidemic curves in the Susceptible-Exposed-Infectious-Recovered (SEIR) model. Sci Rep 10:19365
Himonas AA, Petronilho G (2020) Analyticity in partial differential equations. Complex Anal Synerg 6(15):1–16
Jornet M (2021) Beyond the hypothesis of boundedness for the random coefficient of the Legendre differential equation with uncertainties. Appl Math Comput 391:125638
Jornet M (2023) On the convergence of the power-series method for the fractional Allee logistic equation. Submitted
Kaup L, Kaup B (2011) Holomorphic functions of several variables: an introduction to the fundamental theory, vol 3. Walter de Gruyter, Berlin
Kermack WO, McKendrick AG (1927) Contribution to the mathematical theory of epidemics. Proc R Soc Lond A 115:700–721
Kershaw D (1983) Some extensions of W Gautschi’s inequalities for the Gamma function. Math Comput 41:607–611
Koekoek R, Lesky PA, Swarttouw RF (2010) Hypergeometric orthogonal polynomials and their q-analogues. Springer, Berlin
Lotfi M, Hamblin MR, Rezaei N (2020) COVID-19: Transmission, prevention, and potential therapeutic opportunities. Clin Chim Acta 508:254–266
Martinez VM, Barbosa AN, Mancera PFA, Rodrigues DS, Camargo RF (2021) A fractional calculus model for HIV dynamics: real data, parameter estimation and computational strategies. Chaos Solitons Fractals 152:111398
Masenga SK, Mweene BC, Luwaya E, Muchaili L, Chona M, Kirabo A (2023) HIV-host cell interactions. Cells 12(10):1351
Ndaïrou F, Area I, Nieto JJ, Silva CJ, Torres DF (2021) Fractional model of COVID-19 applied to Galicia. Spain and Portugal. Chaos Solitons Fractals 144:110652
Ortigueira MD, Machado JT (2015) What is a fractional derivative? J Comput Phys 293:4–13
Pinto CM, Carvalho AR (2017) A latency fractional order model for HIV dynamics. J Comput Appl Math 312:240–256
Popović JK, Atanacković MT, Pilipović AS, Rapaić MR, Pilipović S, Atanacković TM (2010) A new approach to the compartmental analysis in pharmacokinetics: fractional time evolution of diclofenac. J Pharmacokinet Pharmacodyn 37:119–134
Rudin W (1976) Principles of mathematical analysis, 3rd edn, International Series in Pure & Applied Mathematics
Sánchez E, Villanueva RJ, Santonja FJ, Rubio M (2011) Predicting cocaine consumption in Spain: a mathematical modelling approach. Drug-Educ Prev Polic 18(2):108–115
Santonja FJ, Sánchez E, Rubio M, Morera JL (2010) Alcohol consumption in Spain and its economic cost: a mathematical modeling approach. Math Comput Model. 52(7–8):999–1003
Srivastava HM, Area Carracedo IC, Nieto JJ (2021) Power-series solution of compartmental epidemiological models. Math Biosci Eng 18(4):3274–3290
Teschl G (2012) Ordinary differential equations and dynamical systems. American Mathematical Society, Providence
Villafuerte L (2023) Solution processes for second-order linear fractional differential equations with random inhomogeneous parts. Math Comput Simul 210:17–48
Wolfram Research, Inc. (2020) Mathematica, Version 12.1, Champaign
Funding
Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The author declares that there is no conflict of interest regarding the publication of this article.
Additional information
Communicated by Carla M. A. Pinto.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Jornet, M. Power-series solutions of fractional-order compartmental models. Comp. Appl. Math. 43, 67 (2024). https://doi.org/10.1007/s40314-023-02579-1
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-023-02579-1
Keywords
- Compartmental differential equation system
- Caputo fractional derivative
- Convergent power-series solution
- Discrete equation
- Symbolic computation
- Fractional SIR model