Next Article in Journal / Special Issue
Solution of Inhomogeneous Differential Equations with Polynomial Coefficients in Terms of the Green’s Function
Previous Article in Journal
Graph Structures in Bipolar Neutrosophic Environment
Previous Article in Special Issue
An Investigation of Radial Basis Function-Finite Difference (RBF-FD) Method for Numerical Solution of Elliptic Partial Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mixed Order Fractional Differential Equations

1
Department of Mathematical Analysis and Numerical Mathematics, Comenius University in Bratislava, Mlynská dolina, 842 48 Bratislava, Slovakia
2
Mathematical Institute of Slovak Academy of Sciences, Štefánikova 49, 814 73 Bratislava, Slovakia
3
Department of Mathematics, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Mathematics 2017, 5(4), 61; https://doi.org/10.3390/math5040061
Submission received: 8 September 2017 / Revised: 30 October 2017 / Accepted: 31 October 2017 / Published: 7 November 2017
(This article belongs to the Special Issue Operators of Fractional Calculus and Their Applications)

Abstract

:
This paper studies fractional differential equations (FDEs) with mixed fractional derivatives. Existence, uniqueness, stability, and asymptotic results are derived.

1. Introduction

Recently, fractional differential equations (FDEs) arise naturally in various fields, such as economics, engineering, and physics. For some existence results of FDEs we refer the reader to [1,2,3,4,5,6] and the references cited therein.
Bonilla et al. [1] studied linear systems of the same order linear FDEs and obtained an explicit representation of the solution. However, there are very few works on the study of mixed order nonlinear fractional differential equations (MOFDEs), which is a natural extension of [1].
This paper is devoted to the study of MOFDEs of the form
D 0 + q i x i = f i ( t , x 1 , , x n ) , x i ( 0 ) = u i , i = 1 , , n
where D 0 + q i denotes the Caputo derivative with the lower limit at 0, q i ( 0 , 1 ] , f : R + × R n R n are continuous, specified below and u i R n . We may suppose q 1 q n . Here R + = [ 0 , ) . We are interested in the existence of solutions of (1), then their stability and asymptotic properties under reasonable conditions on f i . FDEs with equal order (i.e., q 1 = = q n ) are widely studied, and we refer the reader to the basic books describing FDEs, such as [7,8]. On the other hand, there are many MOFDEs with interesting applications—for example, to economic systems in [9]. In fact, (1) formulates a model of the national economies in a case of the study of n commonwealth countries, which cannot be simply divided into clear groups of independent and dependent variables. The purpose of this paper is to set a rigorous theoretical background for (1).
The main contributions are stated as follows:
We give some existence and uniqueness results for solutions of (1) when the nonlinear term satisfies global and local Lipschitz conditions.
We analyze the upper bound for Lyapunov exponents of solutions of (1).
We show that the zero solution of an autonomous version system of (1) is asymptotically stable.

2. Existence Results

First we prove an existence and uniqueness result for globally Lipschitzian f = ( f 1 , , f n ) . Let x = max i = 1 , , n | x i | for x = ( x 1 , , x n ) . By C ( J , R n ) we denote the Banach space of all continuous functions from a compact interval J R to R n with the uniform convergence topology on J.
Theorem 1.
Let T > 0 and suppose the existence of L > 0 such that
f ( t , x ) f ( t , y ) L x y ( t , x ) , ( t , y ) [ 0 , T ] × R n .
Then (1) has a unique solution x C ( I , R n ) , I = [ 0 , T ] .
Proof of Theorem 1.
Note that (2) implies
f ( t , x ) L x + M f ( t , x ) I × R n
for M f = max t I f ( t , 0 ) . Next, (1) is equivalent to the fixed point problem
x = F ( x , u ) , x ( t ) = ( x 1 ( t ) , , x n ( t ) ) , u = ( u 1 , , u n ) , F ( x , u ) = ( F 1 ( x , u ) , , F n ( x , u ) ) , F i ( x , u ) ( t ) = u i + 1 Γ ( q i ) 0 t ( t s ) q i 1 f i ( s , x 1 ( s ) , , x n ( s ) ) d s , i = 1 , , n .
Fix α 0 and set x α = max t [ 0 , T ] x ( t ) e α t for any x C ( I , R n ) . Let x C ( I , R n ) be a solution of (1). For α > 0 , (3) and (4) give
| x i ( t ) | | u i | + 1 Γ ( q i ) 0 t ( t s ) q i 1 L x ( s ) + M f d s u + M f T q i Γ ( q i + 1 ) + L x α Γ ( q i ) 0 t ( t s ) q i 1 e α s d s u + M f T q i Γ ( q i + 1 ) + L x α α q i e α t ,
which implies
x α u + M f Θ + L x α α q m
for α > 1 and
Θ = max i = 1 , , m T q i Γ ( q i + 1 ) .
Hence
x α u + M f Θ 1 L α q m ,
for
α > max { 1 , L 1 / q m } .
So
| x i ( t ) | u + M f Θ 1 L α q m e α T , t I , i = 1 , , n .
Similarly, we derive
F ( x , u ) F ( y , u ) α L α q m x y α , x , y C ( I , R n ) .
Consequently, assuming (5), we can apply the Banach fixed point theorem to get a unique solution x C ( I , R n ) of (1), which also satisfies (6). The proof is finished. ☐
Now we prove an existence and uniqueness result for locally Lipschitzian f.
Theorem 2.
Suppose that for any r > 0 there is an L r > 0 such that
f ( t , x ) f ( t , y ) L r x y ( t , x ) , ( t , y ) [ 0 , T ] × R n , max { x , y } r .
Then (1) has a unique solution x C ( I 0 , R n ) , I 0 = [ 0 , T 0 ] for some 0 < T 0 T .
Proof of Theorem 2.
Set r 0 = u + M f Θ + 1 and B ( r 0 ) = { x R n x r 0 } . Then we extend f from the set I × B ( r 0 ) to f ˜ on I × R n such that f ˜ satisfies (2) for some L > 0 . This extension is given by
f ˜ ( t , x ) = χ ( x / ( r 0 + 1 ) ) f ( t , x ) for t I , | | x | | 2 ( r 0 + 1 ) 0 for t I , | | x | | 2 ( r 0 + 1 )
for a Lipschitz function χ : R + [ 0 , 1 ] with χ ( r ) = 1 for r [ 0 , 1 ] and χ ( r ) = 0 for r 2 . Applying Theorem 1 to
D 0 + q i x i = f ˜ i ( t , x 1 , , x n ) , x i ( 0 ) = u i , i = 1 , , n ,
there is a unique solution x C ( I , R n ) of (8) which also satisfies (6) for α satisfying (5). Note M f = M f ˜ . Let us take T T 0 > 0 , α > 0 satisfying (5) and
u + M f Θ 1 L α q m e α T 0 < r 0 .
Then the unique solution x C ( I , R n ) of (8) satisfies
x ( t ) r 0 t I 0 .
However, this is also a unique solution of (1) on I 0 . The proof is finished. ☐
Remark 1.
Let us denote by x u the solution from Theorem 1. Then, following the proof of Theorem 1, for any u , v R n , we derive
x u x v α = F ( x u , u ) F ( x v , v ) α L α q m x u x v α + F ( x v , u ) F ( x v , v ) α L α q m x u x v α + u v ,
which implies
x u x v α u v 1 L α q m ,
i.e.,
x u ( t ) x v ( t ) u v 1 L α q m t I ,
provided that (5) holds. So, the continuous dependence of the solution of (1) is shown on the initial value u under conditions of Theorem 1 or Theorem 2.

3. Asymptotic Results

We find the upper bound for Lyapunov exponents of solutions of (1).
Theorem 3.
Suppose assumptions of Theorem 2 are satisfied. Moreover, we suppose the existence of a nonnegative n × n -matrix M = { m i j } i , j 1 n and a nonegative vector v = ( v 1 , , v n ) such that
| f i ( t , x ) | j = 1 n m i j | x j | + v i ( t , x ) R + × R n , i = 1 , , n .
Then the Lyapunov exponent
Λ ( u ) = lim sup t ln x ( t ) t
of the unique solution x C ( R + , R n ) of (1) satisfies
Λ ( u ) ρ ( M ) max i = 1 , , n 1 / q i ,
where ρ ( M ) is the spectral radius of M. Note x ( 0 ) = u , so we consider as usually Λ : R n R + (see [10]).
Proof. 
Clearly (9) implies
f ( t , x ) max i = 1 , , n j = 1 n m i j x + v ( t , x ) R + × R n .
Then, like in the proof of Theorem 1, for any T > 0 there is a unique solution of (1) on I. However, since T > 0 is arbitrary, we get a unique solution x ( t ) on R + . Then, we compute
| x i ( t ) | | u i | + 1 Γ ( q i ) 0 t ( t s ) q i 1 j = 1 n m i j | x j ( s ) | + v j d s = | u i | + t q i Γ ( q i + 1 ) j = 1 n m i j v j + 1 Γ ( q i ) j = 1 n m i j 0 t ( t s ) q i 1 e α s | x j ( s ) | e α s d s | u i | + t q i Γ ( q i + 1 ) j = 1 n m i j v j + 1 Γ ( q i ) j = 1 n m i j x j α Γ ( q i ) e α t α q i = | u i | + t q i Γ ( q i + 1 ) j = 1 n m i j v j + e α t α q i j = 1 n m i j x j α | u i | + Y α e α t j = 1 n m i j v j + e α t α q i j = 1 n m i j x j α
for
Y α = max i = 1 , , n , t R + t q i e α t Γ ( q i + 1 ) = max i = 1 , , n q i q i e q i α q i Γ ( q i + 1 ) .
Consequently, we arrive at
x i α | u i | + Y α j = 1 n m i j v j + 1 α q i j = 1 n m i j x j α ,
or setting w α = ( x 1 α , , x m α ) and | u | = ( | u 1 | , , | u n | ) , we obtain
w α ε | u | + Y α ε M v + 1 ρ ( M ) + ε M w α ε ,
for α ε = max i = 1 , , n ( ρ ( M ) + ε ) 1 / q i and ε > 0 is fixed. (10) is considered component-wise. From (10), we get
I 1 ρ ( M ) + ε M w α ε | u | + Y α ε M v
for the n × n identity matrix I. By the Neumann lemma, we have
I 1 ρ ( M ) + ε M 1 = i = 0 1 ρ ( M ) + ε M i ,
which is a positive matrix. Consequently, (11) implies
w α ε I 1 ρ ( M ) + ε M 1 | u | + Y α ε M v .
Letting ε 0 , we get
w α 0 I 1 ρ ( M ) M 1 | u | + Y α 0 M v
for α 0 = max i = 1 , , n ρ ( M ) 1 / q i . So we arrive at
| x i ( t ) | e α 0 t w i t R + , i = 1 , , n
for w = ( w 1 , , w n ) = I 1 ρ ( M ) M 1 | u | + Y α 0 M v , which implies
Λ ( u ) = lim sup t ln x ( t ) t lim sup t ln ( e α 0 t w ) t = lim sup t ln e α 0 t t + lim sup t w t = α 0 .
The proof is finished. ☐

4. Stability Result

This section is devoted to an autonomous version of (1)
D 0 + q i x i = j = 1 n a i j x j + g i ( x 1 , , x n ) , x i ( 0 ) = u i , i = 1 , , n
where a i , j R , g = ( g 1 , , g n ) : R n R n is locally Lipschitz with g ( 0 ) = 0 and lim x 0 g ( x ) x = 0 . We already know that (12) has a locally unique solution. First, we prove the following
Lemma 1.
Let λ > 0 and q , p ( 0 , 1 ] with q p . Then, it holds
S ( λ , q , p ) = sup t 0 ( t + 1 ) p 0 t ( t s ) q 1 E q , q ( ( t s ) q λ ) ( s + 1 ) p d s < .
Note that E q , q ( t ) 0 for any t R + by ([11] p. 85), where E q , q ( t ) is the two-parametric Mittag-Leffler function ([11] p. 56) .
Proof. 
By ([5] Proposition 2.4(i))or ([11] Formula (4.4.17)), there is a positive constant M ( q , λ ) such that
t q 1 E q , q ( t q λ ) M ( q , λ ) ( t + 1 ) q + 1 t 1 .
Here we note that one can use ([6] Lemma 3) to compute M ( q , λ ) . First, using ([11] Formula (4.4.4)), we derive
sup t [ 0 , 1 ] ( t + 1 ) p 0 t ( t s ) q 1 E q , q ( ( t s ) q λ ) ( s + 1 ) p d s sup t [ 0 , 1 ] ( t + 1 ) p 0 t ( t s ) q 1 E q , q ( ( t s ) q λ ) d s = sup t [ 0 , 1 ] ( t + 1 ) p 0 t z q 1 E q , q ( z q λ ) d z = sup t [ 0 , 1 ] ( t + 1 ) p t q E q , q + 1 ( λ t q ) = 2 p E q , q + 1 ( λ ) ,
since applying ([11] Formula (5.1.15)), we derive
d d t ( t + 1 ) p t q E q , q + 1 ( λ t q ) = p ( t + 1 ) p 1 t q E q , q + 1 ( λ t q ) + ( t + 1 ) p t q 1 E q , q ( λ t q ) 0
for t > 0 . Next, for t 1 , we derive
sup t 1 ( t + 1 ) p 0 t ( t s ) q 1 E q , q ( ( t s ) q λ ) ( s + 1 ) p d s sup t 1 ( t + 1 ) p 0 t 1 ( t s ) q 1 E q , q ( ( t s ) q λ ) ( s + 1 ) p d s + sup t 1 ( t + 1 ) p t 1 t ( t s ) q 1 E q , q ( ( t s ) q λ ) ( s + 1 ) p d s sup t 1 ( t + 1 ) p 0 t 1 M ( q , λ ) ( t s + 1 ) q + 1 ( s + 1 ) p d s + sup t 1 ( t + 1 ) p t p 0 1 z q 1 E q , q ( z q λ ) d z sup t 1 ( t + 1 ) p 0 t / 2 M ( q , λ ) ( t s + 1 ) q + 1 d s + sup t 1 2 p ( t + 1 ) p ( t + 2 ) p t / 2 t M ( q , λ ) ( t s + 1 ) q + 1 d s + 2 p E q , q + 1 ( λ ) sup t 1 2 q ( t + 1 ) p M ( q , λ ) q ( t + 2 ) q + sup t 1 2 p ( t + 1 ) p M ( q , λ ) q ( t + 2 ) p + 2 p E q , q + 1 ( λ ) ( 2 q + 2 p ) M ( q , λ ) q + 2 p E q , q + 1 ( λ )
Summarizing, we arrive at the estimate
S ( λ , q , p ) ( 2 q + 2 p ) M ( q , λ ) q + 2 p E q , q + 1 ( λ ) .
The proof is finished. ☐
Theorem 4.
Suppose a i i = λ i < 0 , i = 1 , , n . If
γ = max i = 1 , , n S ( λ i , q i , q n ) i j = 1 n | a i j | < 1 ,
then the zero solution of (12) is asymptotically stable.
Proof. 
The proof is motivated by the well-known Geršgoring type method [12,13]. Like above, we modify (12) outside of the unit ball B ( 1 ) such that the modified system is globally Lipschitz. Then, the solution of the modified (12) has a global unique solution on R + by Theorem 1. This solution of (12) has the fixed point form [8]
x i ( t ) = G i ( x , u ) , i = 1 , , n , G i ( x , u ) = E q i ( t q i λ i ) u i + 0 t ( t s ) q i 1 E q i , q i ( ( t s ) q i λ i ) ( i j = 1 n a i j x j ( s ) + g i ( x 1 ( s ) , , x n ( s ) ) ) d s ,
where E q ( t ) and E q , q ( t ) are the classical and two-parametric Mittag-Leffler functions, respectively ([11] p. 56). Fix T > 0 and set | | | x | | | = max t I | x ( t ) | ( t + 1 ) q n for x C ( I , R n ) . Then, (15) implies
| G i ( x , u ) ( t ) | ( t + 1 ) q n κ | u i | + ( t + 1 ) q n 0 t ( t s ) q i 1 E q i , q i ( ( t s ) q i λ i ) ( s + 1 ) q n d s × ( i j = 1 n | a i j | + h ( x ( s ) ) ) | | | x | | | ) κ u + S ( λ i , q i , q n ) ( i j = 1 n | a i j | + h ( | | | x | | | ) ) | | | x | | | )
for
κ = max i = 1 , , n sup t 0 E q i ( t q i λ i ) ( t + 1 ) q n <
(see ([11] Formula (3.4.30)) and h ( r ) = max x r | | g ( x ) | | x . Consequently, we obtain
| | | G ( x , u ) | | | κ u + ( γ + ω h ( | | | x | | | ) ) | | | x | | | )
for G ( x , u ) = ( G 1 ( x , u ) , , G n ( x , u ) ) and
ω = max i = 1 , , n i j = 1 n S ( λ i , q i , q n ) .
Taking r 1 ( 0 , 1 ) such that γ + ω h ( r 1 ) 1 + γ 2 , then (16) implies
| | | G ( x , u ) | | | κ u + 1 + γ 2 | | | x | | |
for any x B T ( r 1 ) = { x C ( I , R n ) | | | x | | | r 1 } . So, supposing
u 1 γ 2 κ r 1 ,
(17) gives
G : B T ( r 1 ) B T ( r 1 ) .
It is well-known that G : C ( I , R n ) C ( I , R n ) is continuous and compact [8]. Since B T ( r 1 ) is convex and bounded, by the Schauder fixed point theorem, G has a fixed point x B T ( r 1 ) . However, this a solution of a modified (12), which has a unique solution on R + . So, this unique solution satisfies x B T ( r 1 ) for any t > 0 ; i.e.,
x ( t ) r 1 ( t + 1 ) q n t 0 .
Certainly (19) gives x ( t ) 1 , so x ( t ) is also a unique solution of the original (12). Moreover, (17) leads to
x ( t ) 2 κ u ( 1 γ ) ( t + 1 ) q n t 0 ,
which determines the asymptotic stability of (12) at 0. The proof is finished. ☐
Remark 2.
Condition (14) is a Geršgoring type assumption [12,13]. Even for q 1 = = q n , condition (14) is new, since the general stability assumption (see [5] Theorem 2.2) is difficult to check. The above approach can be extended to FDEs on infinite lattices like in [14].

5. Examples

In this section, we give examples to demonstrate the validity of our theoretical results.
Example 1.
Consider
c D 0 + 0.6 x 1 = sin ( x 1 + 2 x 2 ) , t [ 0 , 1 ] , c D 0 + 0.5 x 2 = arctan ( 2 x 1 + x 2 ) , x 1 ( 0 ) = u 1 , x 2 ( 0 ) = u 2 .
Set q 1 = 0.6 > 0.5 = q 2 , T = 1 , f 1 ( t , x 1 , x 2 ) = sin ( x 1 + 2 x 2 ) and f 2 ( t , x 1 , x 2 ) = arctan ( 2 x 1 + x 2 ) . Clearly, f = ( f 1 , f 2 ) satisfying the uniformly Lipschitz condition with a Lipschitz constant L = 2 . By Theorem 1 or 2, (20) has a unique solution x C ( [ 0 , 1 ] , R 2 ) .
Next, | f 1 ( t , x 1 , x 2 ) |     m 11 | x 1 | + m 12 | x 2 | and | f 2 ( t , x 1 , x 2 ) |     m 21 | x 1 | + m 22 | x 2 | where m 11 = 1 = m 22 and m 12 = m 21 = 2 Set M = ( m i j ) R 2 × 2 , so
M = 1 2 2 1 ,
is a nonnegative matrix and ρ ( M ) = 3 . By Theorem 3, Λ ( u ) = 3 5 3 .
Example 2.
Consider
c D 0 + 0.6 x 1 = x 1 0.5 x 2 + x 1 + x 2 1 + ( x 1 + x 2 ) 2 sin ( x 1 + x 2 ) , t 0 , c D 0 + 0.5 x 2 = 0.5 x 1 x 2 + x 1 + x 2 1 + ( x 1 + x 2 ) 2 arctan ( x 1 + x 2 ) .
Set q 1 = 0.6 > 0.5 = q 2 , λ 1 = λ 2 = 1 , a 11 = a 22 = 1 and a 12 = 0.5 , a 21 = 0.5 , g 1 ( x 1 , x 2 ) = x 1 + x 2 1 + ( x 1 + x 2 ) 2 sin ( x 1 + x 2 ) , g 2 ( x 1 , x 2 ) = x 1 + x 2 1 + ( x 1 + x 2 ) 2 arctan ( x 1 + x 2 ) and u 1 = u 2 = 0 . Clearly, g = ( g 1 , g 2 ) satisfies the uniformly Lipschitz condition with a Lipschitz constant L = 2 and g ( ( 0 , 0 ) ) = ( 0 , 0 ) and lim ( x 1 , x 2 ) ( 0 , 0 ) g ( x 1 , x 2 ) ( x 1 , x 2 ) = 0 . We numerically derive
S ( 1 , 0.6 , 05 ) 1.0051 , S ( 1 , 0.5 , 0 , 5 ) 1
so γ 0.503 < 1 . By Theorem 4, the zero solution of (21) is asymptotically stable.

6. Conclusions

The existence and uniqueness of solutions of (1) are shown along with their stability and asymptotic properties. Theorems 1 and 2 extended the partial case of q i = q , i = 1 , 2 , , n which are known in the related literature. Theorems 3 and 4 are original results. Moreover, Theorems 1, 2, and 3 can be directly extended to infinite dimensional cases. The next step should be to derive an exact solution for scalar linear systems (i.e., (12) with g = 0 ) and to find an explicit variation of constant formula for nonhomogeneous systems. Of course, a more general stability criterion would also be interesting to find by generalizing ([5] Theorem 2.2), which is just for q 1 = = q n . Then, a derivation of a Gronwall-type inequality associated to MOFDEs would be challenging as well. This should extend Theorem 4 to infinite dimensional cases.

Acknowledgments

The authors thanks the referees for their careful reading of the manuscript and insightful comments, which help to improve the quality of the paper. This work is supported by National Natural Science Foundation of China (11661016) and the Slovak Research and Development Agency under the contract No. APVV-14-0378 and by the Slovak Grant Agency VEGA No. 2/0153/16 and No. 1/0078/17.

Author Contributions

M.F. and J.R.W. contributed equally in this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bonilla, B.; Rivero, M.; Trujillo, J.J. On systems of linear fractional differential equations with constant coefficients. Appl. Math. Comput. 2007, 187, 68–78. [Google Scholar] [CrossRef]
  2. Chalishajar, D.N.; Karthikeyan, K.; Trujillo, J.J. Existence of mild solutions for fractional impulsive semilinear integro-differential equations in Banach spaces. Commun. Appl. Nonlinear Anal. 2012, 4, 45–56. [Google Scholar]
  3. Chalishajar, D.N.; Karthikeyan, K. Existence and uniqueness results for boundary value problems of higher order fractional integro-differential equations involving Gronwall’s inequality in banach spaces. Acta Math. Sci. 2013, 33, 758–772. [Google Scholar] [CrossRef]
  4. Chalishajar, D.N.; Karthikeyan, K. Boundary value problems for impulsive fractional evolution integro-differential equations with Gronwall’s inequality in Banach spaces. Discontin. Nonlinearity Complex. 2014, 3, 33–48. [Google Scholar] [CrossRef]
  5. Cong, N.D.; Doan, T.S.; Siegmund, S.; Tuan, H.T. Linearized asymptotic stability for fractional differential equations. Electron. J. Qual. Theory Differ. Equ. 2016, 39, 1–13. [Google Scholar] [CrossRef]
  6. Cong, N.D.; Doan, T.S.; Siegmund, S.; Tuan, H.T. On stable manifolds for planar fractional differential equations. Appl. Math. Comput. 2014, 226, 157–168. [Google Scholar] [CrossRef]
  7. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  8. Zhou, Y. Fractional Evolution Equations and Inclusions: Analysis and Control; Academic Press: Amsterdam, The Netherlands, 2016. [Google Scholar]
  9. Škovránek, T.; Podlubny, I.; Petráš, I. Modeling of the national economies in state-space: A fractional calculus approach. Econ. Model. 2012, 29, 1322–1327. [Google Scholar] [CrossRef]
  10. Barreira, L.; Pesin, Y. Lectures on Lyapunov Exponents and Smooth Ergodic Theory; American Mathematical Society: Providence, RI, USA, 1998. [Google Scholar]
  11. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  12. Gil, M.I. Operator Functions and Localization of Spectra; Springer: Berlin, Germany, 2003. [Google Scholar]
  13. Varga, R.S. Geršgorin and His Circles; Springer: Berlin, Germany, 2011. [Google Scholar]
  14. Fečkan, M.; Kelemen, S. Multivalued integral manifolds in Banach spaces. Commun. Math. Anal. 2011, 10, 97–117. [Google Scholar]

Share and Cite

MDPI and ACS Style

Fečkan, M.; Wang, J. Mixed Order Fractional Differential Equations. Mathematics 2017, 5, 61. https://doi.org/10.3390/math5040061

AMA Style

Fečkan M, Wang J. Mixed Order Fractional Differential Equations. Mathematics. 2017; 5(4):61. https://doi.org/10.3390/math5040061

Chicago/Turabian Style

Fečkan, Michal, and JinRong Wang. 2017. "Mixed Order Fractional Differential Equations" Mathematics 5, no. 4: 61. https://doi.org/10.3390/math5040061

APA Style

Fečkan, M., & Wang, J. (2017). Mixed Order Fractional Differential Equations. Mathematics, 5(4), 61. https://doi.org/10.3390/math5040061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop