Next Article in Journal
Some New Inequalities and Extremal Solutions of a Caputo–Fabrizio Fractional Bagley–Torvik Differential Equation
Next Article in Special Issue
Well-Posedness and Global Attractors for Viscous Fractional Cahn–Hilliard Equations with Memory
Previous Article in Journal
Formal Verification of Fractional-Order PID Control Systems Using Higher-Order Logic
Previous Article in Special Issue
Multiplicity of Solutions for Fractional-Order Differential Equations via the κ(x)-Laplacian Operator and the Genus Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Sturm–Liouville Eigenvalue Problems, II

by
Mohammad Dehghan
* and
Angelo B. Mingarelli
School of Mathematics and Statistics, Carleton University, Ottawa, ON K1S 5B6, Canada
*
Author to whom correspondence should be addressed.
Submission received: 29 July 2022 / Revised: 21 August 2022 / Accepted: 25 August 2022 / Published: 30 August 2022
(This article belongs to the Special Issue Advances in Fractional Differential Operators and Their Applications)

Abstract

:
We continue the study of a non-self-adjoint fractional three-term Sturm–Liouville boundary value problem (with a potential term) formed by the composition of a left Caputo and left Riemann–Liouville fractional integral under Dirichlet type boundary conditions. We study the existence and asymptotic behavior of the real eigenvalues and show that for certain values of the fractional differentiation parameter α , 0 < α < 1 , there is a finite set of real eigenvalues and that, for α near 1 / 2 , there may be none at all. As α 1 we show that their number becomes infinite and that the problem then approaches a standard Dirichlet Sturm–Liouville problem with the composition of the operators becoming the operator of second order differentiation.

1. Introduction

This is a continuation of [1] where the results therein are extended to three-term Fractional Sturm–Liouville operators (with a potential term) formed by the composition of a left Caputo and left-Riemann–Liouville fractional integral. Similar kinds of spectral problems have been considered in [2,3,4,5,6,7,8,9,10,11]. Specifically, the boundary value problem considered here is of the form,
c D 0 + α D 0 + α y ( t ) + q ( t ) y ( t ) = λ y ( t ) , 1 / 2 < α < 1 , 0 t 1 ,
with boundary conditions
I 0 + 1 α y ( t ) | t = 0 = c 1 , and I 0 + 1 α y ( t ) | t = 1 = c 2 ,
where c 1 , c 2 are real constants and the real valued unspecified potential function, q L [ 0 , 1 ] . We note that these are not self-adjoint problems and so there may be a non-real spectrum, in general. A well-known property of the Riemann–Liouville integral gives that if the solutions are continuous on [ 0 , 1 ] then the boundary conditions (2) reduce to the usual fixed-end boundary conditions, y ( 0 ) = y ( 1 ) = 0 , as α 1 .
For the analogue of the Dirichlet problem described above we study the existence and asymptotic behaviour of the real eigenvalues and show that for each α , 0 < α < 1 , there is a finite set of real eigenvalues and that, for α near 1 / 2 , there may be none at all. As α 1 we show that their number becomes infinite and that the problem then approaches a standard Dirichlet Sturm–Liouville problem with the composition of the operators becoming the operator of second order differentiation acting on a suitable function space.
Our approach is different from most in this area. Specifically, we start with the existence and uniqueness of solutions of the Equation (1) along with the initial conditions (2), then we formulate the boundary value problem as an integral equation, after which we show that the solution of this integral equation as a function of λ is an entire function of λ of order of at most 1 / 2 α . Since α is between 1 / 2 and 1, this entire function is of fractional order and therefore must have an infinite number of zeros, some of which may be complex. These zeros are the eigenvalues of our problem and therefore we get their existence. Using asymptotic methods and hard analysis, we proved our bounds on the real eigenvalues of which there must be a finite number for each α . Finally, we show that as α tends to 1, the number of real eigenvalues becomes infinite and the original problem approaches the standard Sturm–Liouville eigenvalue problem.

2. Preliminaries

We recall some definitions from fractional calculus and refer the reader to our previous paper [1] for further details.
Definition 1.
The left and the right Riemann–Liouville fractional integrals I a + α and I b α of order α R + are defined by
I a + α f ( t ) : = 1 Γ ( α ) a t f ( s ) ( t s ) 1 α d s , t ( a , b ] ,
and
I b α f ( t ) : = 1 Γ ( α ) t b f ( s ) ( s t ) 1 α d s , t [ a , b ) ,
respectively. Here Γ ( α ) denotes Euler’s Gamma function. The following property is easily verified.
Property 1.
For a constant C, we have I a + α C = ( t a ) α Γ ( α + 1 ) · C .
The proof is by direct calculation.
Definition 2.
The left and the right Caputo fractional derivatives c D a + α and c D b α are defined by
c D a + α f ( t ) : = I a + n α D n f ( t ) = 1 Γ ( n α ) a t f ( n ) ( s ) ( t s ) α n + 1 d s , t > a ,
and
c D b α f ( t ) : = ( 1 ) n I b n α D n f ( t ) = ( 1 ) n Γ ( n α ) t b f ( n ) ( s ) ( s t ) α n + 1 d s , t < b ,
respectively, where f is sufficiently differentiable and n 1 α < n .
Definition 3.
Similarly, the left and the right Riemann–Liouville fractional derivatives D a + α and D b α are defined by
D a + α f ( t ) : = D n I a + n α f ( t ) = 1 Γ ( n α ) d n d t n a t f ( s ) ( t s ) α n + 1 d s , t > a ,
and
D b α f ( t ) : = ( 1 ) n D n I b n α f ( t ) = ( 1 ) n Γ ( n α ) d n d t n t b f ( s ) ( s t ) α n + 1 d s , t < b ,
respectively, where f is sufficiently differentiable and n 1 α < n .
Property 2.
For ( ν ) > 1 , 0 < α < 1 , and t > 0 , we have
D 0 + α ( t ν ) = Γ ( 1 + ν ) Γ ( 1 + ν α ) t ν α
Property 3.
For ( ν ) > 0 , 0 < α < 1 , and t > 0 , we have
c D 0 + α ( t ν ) = Γ ( 1 + ν ) Γ ( 1 + ν α ) t ν α
Property 4.
If y ( t ) L 1 ( a , b ) and I a + 1 α y , I b 1 α y A C [ a , b ] , then
I a + α D a + α y ( t ) = y ( t ) ( t a ) α 1 Γ ( α ) I a + 1 α y ( a ) , I b α D b α y ( t ) = y ( t ) ( b t ) α 1 Γ ( α ) I b 1 α y ( b ) .
Property 5.
If y ( t ) A C [ a , b ] and 0 < α 1 , then
I a + α c D a + α y ( t ) = y ( t ) y ( a ) , I b α c D b α y ( t ) = y ( t ) y ( b ) .
Property 6.
For 0 < α < 1 we have
D a + α f ( t ) = f ( a ) Γ ( 1 α ) ( t a ) α + c D a + α f ( t ) .

The Mittag-Leffler Function

The function E δ ( z ) defined by
E δ ( z ) : = k = 0 z δ Γ ( δ k + 1 ) , ( z C , ( δ ) > 0 ) ,
was introduced by Mittag-Leffler [12]. In particular, when δ = 1 and δ = 2 , we have
E 1 ( z ) = e z , E 2 ( z ) = cosh ( z ) .
The generalized Mittag-Leffler function E δ , θ ( z ) is defined by
E δ , θ ( z ) = k = 0   z k Γ ( δ k + θ ) ,
where z , θ C and Re   ( δ ) > 0 . When θ = 1 , E δ , θ ( z ) coincides with the Mittag-Leffler function (9):
E δ , 1 ( z ) = E δ ( z ) .
Two other particular cases of (11) are as follows:
E 1 , 2 ( z ) = e z 1 z , E 2 , 2 ( z ) = sinh ( z ) z .
Property 7.
For any δ with ( δ ) > 0 and for any z 0 we have
E δ , δ ( z ) = 1 z E δ , 0 ( z )
Further properties of this special function may be found in [13].
Property 8
(See [14], p. 43). If 0 < δ < 2 and μ ( δ π 2 , min ( π , δ π ) ) , then the function E δ , θ ( z ) has the following exponential expansion as | z |
E δ , θ ( z ) = 1 δ z 1 θ δ exp ( z 1 δ ) k = 1 N 1 Γ ( θ δ k ) 1 z k + O ( 1 z N + 1 ) , | arg ( z ) | μ , k = 1 N 1 Γ ( θ δ k ) 1 z k + O ( 1 z N + 1 ) , μ | arg ( z ) | π .

3. Existence and Uniqueness of the Solution of SLPs

In this section we convert (1) and (2) to an integral equation and prove that it has a solution that satisfies the relevant equations and initial conditions. First, we proceed formally. Separating terms in (1), we get
c D 0 + α D 0 + α y ( t ) = ( q ( t ) λ ) y ( t ) , 1 / 2 < α < 1 , 0 t 1 .
Taking the left Riemann–Liouville fractional integrals I a + α on both sides of the above equation and using Property 5, we have
D 0 + α y ( t ) D 0 + α y ( t ) | t = 0 = I 0 + α ( ( q ( t ) λ ) y ( t ) ) .
Taking the left Riemann–Liouville fractional integrals I a + α from both sides of the above equation once again and using Property 4, we get
y ( t ) t α 1 Γ ( α ) I 0 + 1 α y ( t ) | t = 0 I 0 + α ( D 0 + α y ( t ) | t = 0 ) = I 0 + α ( I 0 + α ( ( q ( t ) λ ) y ( t ) ) )
Using Property 1, we can write
y ( t ) = c 1 t α 1 Γ ( α ) + c 2 t α Γ ( α + 1 ) + I 0 + α ( I 0 + α ( ( q ( t ) λ ) y ( t ) ) )
in which
c 1 = I 0 + 1 α y ( t ) | t = 0 , c 2 = D 0 + α y ( t ) | t = 0 .
We obtain, through the double fractional integral in the above equation, the following:
y ( t ) = c 1 t α 1 Γ ( α ) + c 2 t α Γ ( α + 1 ) + 1 Γ 2 ( α ) 0 t ( t s ) α 1 0 s ( q ( r ) λ ) y ( r ) ( s r ) 1 α d r d s .
By changing the order of integrals in the above equation we get
y ( t ) = c 1 t α 1 Γ ( α ) + c 2 t α Γ ( α + 1 ) + 1 Γ 2 ( α ) 0 t ( q ( r ) λ ) y ( r ) r t ( t s ) α 1 ( s r ) α 1 d s d r
Solving the inner integral gives us
y ( t , λ ) = c 1 t α 1 Γ ( α ) + c 2 t α Γ ( α + 1 ) + 1 Γ ( 2 α ) 0 t ( q ( s ) λ ) y ( s , λ ) ( t s ) 2 α 1 d s .
We will now show that (15) has a solution that exists in a neighbourhood of t = 0 and is unique there. Working backwards will then provide us with a unique solution to (1) and (2). Although this result already appears in [15], we give a shorter proof part of which will be required later.
To this end, let t > 0 . Define
y n ( t , λ ) = y 0 ( t , λ ) + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( q ( s ) λ ) y n 1 ( s , λ ) d s ,
where
y 0 ( t , λ ) = c 1 t α 1 Γ ( α ) + c 2 t α Γ ( α + 1 ) .
Let λ C , | λ | < Λ , where Λ > 0 is arbitrary but fixed. Then,
| y 1 ( t , λ ) y 0 ( t , λ ) | 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ | | y 0 ( s , λ ) | d s | | q | | + Λ Γ ( 2 α ) 0 t ( t s ) 2 α 1 | y 0 ( s , λ ) | d s ,
in which | | q | | = sup t [ 0 , 1 ] | q ( t ) | . Substituting (17) in (18) and using the fact that,
a t ( t s ) α 1 ( s a ) β 1 d s = ( t a ) α + β 1 Γ ( α ) Γ ( β ) Γ ( α + β ) ,
we have
| y 1 ( t , λ ) y 0 ( t , λ ) | ( | | q | | + Λ ) c 1 Γ ( 3 α ) t 3 α 1 + c 2 t 3 α Γ ( 3 α + 1 ) ,
Now, for n = 2 in (16) we get
| y 2 ( t , λ ) y 1 ( t , λ ) | 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ | | y 1 ( s , λ ) y 0 ( s , λ ) | d s 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ | ( | | q | | + Λ ) ( c 1 Γ ( 3 α ) s 3 α 1 + c 2 s 3 α Γ ( 3 α + 1 ) ) d s ( | | q | | + Λ ) 2 c 1 Γ ( 5 α ) t 5 α 1 + c 2 Γ ( 5 α + 1 ) t 5 α .
Continuing in this way we get that the series
y ( t , λ ) : = y 0 ( t , λ ) + n = 1 ( y n ( t , λ ) y n 1 ( t , λ ) )
where
n = 1 | y n ( t , λ ) y n 1 ( t , λ ) | c 1 t 1 n = 1 ( | | q | | + Λ ) n Γ ( 2 n α + α ) t 2 n α + α + c 2 n = 1 ( | | q | | + Λ ) n Γ ( 2 n α + α + 1 ) t 2 n α + α
converges uniformly on compact subsets of ( 0 , 1 ] . Denote the sum of the infinite series in (20) by y ( t , λ ) . So, by virtue of (17) and (21), (20) gives us,
| y ( t , λ ) | | y 0 ( t , λ ) | + n = 1 | y n ( t , λ ) y n 1 ( t , λ ) | c 1 Γ ( α ) t α 1 + c 2 Γ ( α + 1 ) t α + c 1 t 1 n = 1 ( | | q | | + Λ ) n Γ ( 2 n α + α ) t 2 n α + α + c 2 n = 1 ( | | q | | + Λ ) n Γ ( 2 n α + α + 1 ) t 2 n α + α = c 1 t 1 n = 0 ( | | q | | + Λ ) n Γ ( 2 n α + α ) t 2 n α + α + c 2 n = 0 ( | | q | | + Λ ) n Γ ( 2 n α + α + 1 ) t 2 n α + α = c 1 t α 1 E 2 α , α ( ( | | q | | + Λ ) t 2 α ) + c 2 t α E 2 α , α + 1 ( ( | | q | | + Λ ) t 2 α ) .
Note that for a solution y ( t , λ ) of (15) to be C ( [ 0 , 1 ] ) , it is necessary and sufficient that c 1 = 0 , i.e., I 0 + 1 α y ( t ) | t = 0 = 0 . This then proves the global existence of a solution of (15) on [ δ , 1 ] , δ > 0 , since q L [ 0 , 1 ] for given c 1 and c 2 , as defined in (2).
From the proof comes the following a priori estimate when c 1 = 0 , that is,
| y ( t , λ ) | c 2 t α 1 Γ ( α + 1 ) + | E 2 α , α + 1 ( ( | | q | | + Λ ) t 2 α ) | c 2 1 Γ ( α + 1 ) + | E 2 α , α + 1 ( ( | | q | | + Λ ) t 2 α ) |
valid for each t [ 0 , 1 ] and all | λ | < Λ .
The previous bound can be made into an absolute constant by taking the sup over all t and | λ | < Λ . Of course, the bound goes to infinity as | λ | over non-real values, as it must. Thus,
| y ( t , λ ) | c 2 1 Γ ( α + 1 ) + sup | λ | < Λ , t [ 0 , 1 ] | E 2 α , α + 1 ( ( | | q | | + Λ ) t 2 α ) | = c 2 1 Γ ( α + 1 ) + | E 2 α , α + 1 ( ( | | q | | + Λ ) ) | : = c 3 .
for all | λ | < Λ , t [ 0 , 1 ] . Uniqueness follows easily by means of Gronwall’s inequality, as usual. Let ε > 0 . Assume that (15) has two solutions y ( t , λ ) , z ( t , λ ) . Since q L [ 0 , 1 ] and | λ | < Λ we can derive that,
| y ( t , λ ) z ( t , λ ) | ε e 1 Γ ( 2 α ) ( | | q | | + Λ ) t 2 α 2 α .
and since t [ 0 , 1 ] , we get
| y ( t , λ ) z ( t , λ ) |   O ( ε )
where the O-term can be made independent of both t , λ . Letting ε 0 yields uniqueness for t [ 0 , 1 ] and | λ | < Λ .

4. Another Integral Equation

In the previous section we showed that (15) has a solution that, for each λ C , exists on [ 0 , 1 ] , is unique, and is continuous there if and only if c 1 = 0 . On the other hand, if c 1 0 then the solution is merely continuous on all compact subsets of ( 0 , 1 ] . In this section we find another expression for the integral equation which is equivalent to both (15) and the problem (1) with boundary conditions (2).
Lemma 1.
For 0 < α < 1 and 0 < t < 1 , we have
c D 0 + α D 0 + α t α 1 E 2 α , α ( λ t 2 α ) = λ t α 1 E 2 α , α ( λ t 2 α )
Proof. 
Using properties of the Mittag-Leffler function we can write
D 0 + α t α 1 E 2 α , α ( λ t 2 α ) = D 0 + α k = 0 ( λ ) k t 2 α k + α 1 Γ ( 2 α k + α ) = k = 0 ( λ ) k D 0 + α t 2 α k + α 1 Γ ( 2 α k + α ) = k = 0 ( λ ) k t 2 α k 1 Γ ( 2 α k ) = t 1 E 2 α , 0 ( λ t 2 α ) = λ t 2 α 1 E 2 α , 2 α ( λ t 2 α ) ,
in which the third and the last equalities come from Property 2 and Property 7, respectively. Now, taking the left Caputo fractional derivative of both sides of (23) we get
c D 0 + α D 0 + α ( t α 1 E 2 α , α ( λ t 2 α ) ) = c D 0 + α ( λ t 2 α 1 E 2 α , 2 α ( λ t 2 α ) ) = λ c D 0 + α k = 0 ( λ ) k t 2 α k + 2 α 1 Γ ( 2 α k + 2 α ) = λ k = 0 ( λ ) k c D 0 + α t 2 α k + 2 α 1 Γ ( 2 α k + 2 α ) = λ t α 1 k = 0 ( λ ) k t 2 α k Γ ( 2 α k + α ) = λ t α 1 E 2 α , α ( λ t 2 α )
as required. □
Lemma 2.
For 0 < α < 1 and 0 < t < 1 , we have
c D 0 + α D 0 + α t α E 2 α , α + 1 ( λ t 2 α ) = λ t α E 2 α , α + 1 ( λ t 2 α )
Proof. 
Once again, using the properties of the Mittag-Leffler function we can write
D 0 + α t α E 2 α , α + 1 ( λ t 2 α ) = D 0 + α k = 0 ( λ ) k t 2 α k + α Γ ( 2 α k + α + 1 ) = k = 0 ( λ ) k D 0 + α t 2 α k + α Γ ( 2 α k + α + 1 ) = k = 0 ( λ ) k t 2 α k Γ ( 2 α k + 1 ) = E 2 α , 1 ( λ t 2 α ) .
in which the third equality comes from Property 3. Now, taking the left Caputo fractional derivative of both sides of (24) we get
c D 0 + α D 0 + α t α E 2 α , α + 1 ( λ t 2 α ) = c D 0 + α E 2 α , 1 ( λ t 2 α ) = c D 0 + α k = 0 ( λ ) k t 2 α k Γ ( 2 α k + 1 ) = k = 0 ( λ ) k c D 0 + α t 2 α k Γ ( 2 α k + 1 ) = k = 1 ( λ ) k t α ( 2 k 1 ) Γ ( 1 + α ( 2 k 1 ) ) = k = 0 ( λ ) k + 1 t α ( 2 k + 1 ) Γ ( 1 + α ( 2 k + 1 ) ) = λ t α k = 0 ( λ ) k t 2 α k Γ ( 2 k α + α + 1 ) = λ t α E 2 α , α + 1 ( λ t 2 α )
as desired. □
Lemma 3.
For 0 < α < 1 and 0 < t < 1 , we have
c D 0 + α D 0 + α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s = q ( t ) y ( t ) + λ 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s
Proof. 
Let c 4 = 1 / Γ ( 1 α ) . Observe that,
I 0 + 1 α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s = c 4 0 t 0 r ( r s ) 2 α 1 E 2 α , 2 α ( λ ( r s ) 2 α ) q ( s ) y ( s ) d s ( t r ) α d r = c 4 0 t q ( s ) y ( s ) s t ( r s ) 2 α 1 ( t r ) α E 2 α , 2 α ( λ ( r s ) 2 α ) d r d s = c 4 0 t q ( s ) y ( s ) k = 0 ( λ ) k Γ ( 2 α k + 2 α ) s t ( r s ) 2 α 1 + 2 α k ( t r ) α d r d s = 0 t q ( s ) y ( s ) k = 0 ( λ ) k ( t s ) 2 α k + α Γ ( 2 α k + α + 1 ) d s = 0 t q ( s ) y ( s ) ( t s ) α E 2 α , α + 1 ( λ ( t s ) 2 α ) d s .
Next, differentiating both sides of (25) with respect to t and noting that D 0 + α = D ( I 0 + 1 α ) we find,
D 0 + α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s = 0 t ( t s ) α 1 E 2 α , α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s .
as d d t ( t α E 2 α , α + 1 ( λ t 2 α ) ) = t α 1 E 2 α , α ( λ t 2 α ) . Next, we are going to take the left Caputo fractional derivative of both sides of (26). However, since the right hand side of (26) as a function of t is zero at t = 0 , we can use Property 6 and replace the Caputo fractional derivative c D 0 + α by the Riemann–Liouville one D 0 + α . In order to do so, first we need to apply I 0 + 1 α followed by the classical derivative of the right hand side of (26) as follows,
I 0 + 1 α 0 t ( t s ) α 1 E 2 α , α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s = c 4 0 t 0 r ( r s ) α 1 E 2 α , α ( λ ( r s ) 2 α ) q ( s ) y ( s ) d s ( t r ) α d r = c 4 0 t q ( s ) y ( s ) s t ( r s ) α 1 ( t r ) α E 2 α , α ( λ ( r s ) 2 α ) d r d s = c 4 0 t q ( s ) y ( s ) k = 0 ( λ ) k Γ ( 2 α k + α ) s t ( r s ) α 1 + 2 α k ( t r ) α d r d s = 0 t q ( s ) y ( s ) k = 0 ( λ ) k ( t s ) 2 α k Γ ( 2 α k + 1 ) d s = 0 t q ( s ) y ( s ) E 2 α , 1 ( λ ( t s ) 2 α ) d s .
Taking the derivative of the previous equation and using the fact stated in the previous paragraph, we get
c D 0 + α 0 t ( t s ) α 1 E 2 α , α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s = q ( t ) y ( t ) + 0 t q ( s ) y ( s ) ( t s ) 1 E 2 α , 0 ( λ ( t s ) 2 α ) d s = q ( t ) y ( t ) λ 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s ) y ( s ) d s .
where we used Property 7 to arrive at the second equality above. Combining (26) and (27) completes the proof. □
Theorem 1.
For 1 / 2 < α < 1 , the integral equation
y ( t , λ ) = c 1 t α 1   E 2 α , α ( λ t 2 α ) + c 2 t α E 2 α , α + 1 ( λ t 2 α ) + 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s
satisfies (1) with initial conditions I 0 + 1 α y ( t ) | t = 0 = c 1 and D 0 + α y ( t ) | t = 0 = c 2 in which c 1 and c 2 are given constants, and that this solution is unique.
Proof. 
We apply c D 0 + α D 0 + α on both sides of (28) to find,
c D 0 + α D 0 + α y ( t , λ ) = c D 0 + α D 0 + α c 1 t α 1   E 2 α , α ( λ t 2 α ) + c 2 t α E 2 α , α + 1 ( λ t 2 α ) + c D 0 + α D 0 + α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s = λ c 1 t α 1   E 2 α , α ( λ t 2 α ) + λ c 2 t α E 2 α , α + 1 ( λ t 2 α ) q ( t ) y ( t ) + λ 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s = q ( t ) y ( t ) + λ c 1 t α 1   E 2 α , α ( λ t 2 α ) + c 2 t α E 2 α , α + 1 ( λ t 2 α ) + λ   0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s = q ( t ) y ( t ) + λ y ( t , λ ) ,
in which second equality come from Lemma 1, Lemma 2, and Lemma 3. We verify the initial conditions. Taking I 0 + 1 α of both sides (28), we get,
I 0 + 1 α y ( t , λ ) = I 0 + 1 α c 1 t α 1   E 2 α , α ( λ t 2 α ) + c 2 t α E 2 α , α + 1 ( λ t 2 α ) + I 0 + 1 α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s = c 1 E 2 α , 1 ( λ t 2 α ) + c 2 t E 2 α , 2 ( λ t 2 α ) + 0 t ( t s ) α E 2 α , 2 α + 1 ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s ,
where the third term of the second equality comes from (25). Since E 2 α , 1 ( λ t 2 α ) | t = 0 = 1 and the other two terms of the above equality vanish when t = 0 , we have verified the first initial condition. Again Taking D 0 + α on both sides (28), we can find,
D 0 + α y ( t , λ ) = D 0 + α c 1 t α 1   E 2 α , α ( λ t 2 α ) + c 2 t α E 2 α , α + 1 ( λ t 2 α ) + D 0 + α 0 t ( t s ) 2 α 1 E 2 α , 2 α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s = c 1 λ t 2 α 1 E 2 α , 2 α ( λ t 2 α ) + c 2 E 2 α , 1 ( λ t 2 α ) + 0 t ( t s ) α 1 E 2 α , α ( λ ( t s ) 2 α ) q ( s )   y ( s , λ )   d s ,
where the second equality above comes from (23), (24) and (26). The second initial condition can readily be obtained by substituting t = 0 in (31). □

5. Analyticity of Solutions with Respect to the Parameter λ

In this section we show that the solutions (15) or (28) are, generally speaking, entire functions of the parameter λ for each t under consideration and λ C . First, we show continuity with respect to said parameter. Consider the case where c 1 = 0 , i.e., y C [ 0 , 1 ] .
Lemma 4.
Let y C [ 0 , 1 ] , λ C . Then, for each fixed t [ 0 , 1 ] , y ( t , λ ) is continuous with respect to λ.
Proof. 
Let Λ > 0 be arbitrary but fixed, and let | λ | , | λ 0 | < Λ . Using (28),
y ( t , λ ) y ( t , λ 0 ) = 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( q ( s ) λ ) y ( s , λ ) ( q ( s ) λ 0 ) y ( s , λ 0 ) d s = 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( λ 0 λ ) y ( s , λ ) + ( q ( s ) λ 0 ) ( y ( s , λ ) y ( s , λ 0 ) ) d s .
So,
y ( t , λ ) y ( t , λ 0 ) = ( λ λ 0 ) 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 y ( s , λ ) d s + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( q ( s ) λ 0 ) ( y ( s , λ ) y ( s , λ 0 ) ) d s .
Now, let ϵ > 0 and | λ λ 0 | < δ where δ > 0 is to be chosen later. Then,
| y ( t , λ ) y ( t , λ 0 ) | δ 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | y ( s , λ ) | d s + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ 0 | | y ( s , λ ) y ( s , λ 0 ) | d s .
Using (22) and Gronwall’s inequality, we get
| y ( t , λ ) y ( t , λ 0 ) | δ c 3 t 2 α 2 α   Γ ( 2 α ) + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ 0 | | y ( s , λ ) y ( s , λ 0 ) | d s δ c 3 Γ ( 2 α + 1 )   e 1 Γ ( 2 α ) 0 1 ( t s ) 2 α 1 | q ( s ) λ 0 |   d s   δ c 3 Γ ( 2 α + 1 )     e 1 Γ ( 2 α ) 0 1 ( 1 s ) 2 α 1 | q ( s ) λ 0 |   d s : = C δ
where
C = c 3 Γ ( 2 α + 1 )     e 1 Γ ( 2 α ) 0 1 ( 1 s ) 2 α 1 | q ( s ) λ 0 |   d s
is a function of α and λ 0 only as q L ( 0 , 1 ) . Thus, for any t [ 0 , 1 ] , the continuity of y ( t , λ ) follows by choosing δ < ε C . It also follows from this that,
sup t [ 0 , 1 ] | y ( t , λ ) y ( t , λ 0 ) | < ε , | λ λ 0 | < δ .
 □
Next, we consider the differentiability of y ( t , λ ) with respect to λ .
Lemma 5.
Let y C [ 0 , 1 ] , λ C . Then, for each fixed t [ 0 , 1 ] , y ( t , λ ) is differentiable with respect to λ .
Proof. 
As before let | λ | < Λ , t [ 0 , 1 ] . Equation (32) can be rewritten as
y ( t , λ ) y ( t , λ 0 ) λ λ 0 = 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 y ( s , λ )   d s + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( q ( s ) λ 0 ) y ( s , λ ) y ( s , λ 0 ) λ λ 0   d s .
As y ( t , λ 0 ) is given, we define h ( t , λ 0 ) to be the unique solution of the Volterra integral equation of the second kind,
h ( t , λ 0 ) = 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 y ( s , λ 0 )   d s + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 ( q ( s ) λ 0 ) h ( s , λ 0 )   d s .
So,
| y ( t , λ ) y ( t , λ 0 ) λ λ 0 h ( t , λ 0 ) | 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | y ( s , λ ) y ( s , λ 0 ) | d s + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ 0 | | y ( s , λ ) y ( s , λ 0 ) λ λ 0 h ( s , λ 0 ) | d s .
Let ε > 0 and choose δ > 0 as in (33). Using Gronwall’s inequality and (33) we get, for t [ 0 , 1 ] ,
| y ( t , λ ) y ( t , λ 0 ) λ λ 0 h ( t , λ 0 ) | ε 2 α Γ ( 2 α ) + 1 Γ ( 2 α ) 0 t ( t s ) 2 α 1 | q ( s ) λ 0 | | y ( s , λ ) y ( s , λ 0 ) λ λ 0 h ( s , λ 0 ) | d s , ε Γ ( 2 α + 1 )   e 1 Γ ( 2 α ) 0 1 ( t s ) 2 α 1 | q ( s ) λ 0 |   d s = O ( ε ) .
for λ near λ 0 since, for t [ 0 , 1 ] , 0 1 ( t s ) 2 α 1 | q ( s ) λ 0 |   d s = O ( 1 ) . Thus,
y ( t , λ ) λ | λ = λ 0 : = lim λ λ 0 y ( t , λ ) y ( t , λ 0 ) λ λ 0 = h ( t , λ 0 ) ,
exists at λ 0 . Since λ 0 is arbitrary y λ ( t , λ ) exists for all λ with | λ | < Λ , real or complex and the result follows. □
Theorem 2.
For each t [ 0 , 1 ] , y ( t , λ ) is an entire function of λ.
Proof. 
This follows from Lemma 5 since λ C and | λ | < Λ where Λ > 0 is arbitrary. □

6. A Dirichlet-Type Problem

Let y C [ 0 , 1 ] , λ C be fixed. In this case we note that the first of the boundary conditions (2) is equivalent to the usual fixed end (Dirichlet) boundary conditions, that is,
y C [ 0 , 1 ] I 0 + 1 α y ( t , λ ) | t = 0 = 0 y ( 0 , λ ) = 0 .
The continuity assumption implies that there is a number M such that | y ( t , λ ) | M , for all t [ 0 , 1 ] . Thus,
| I 0 + 1 α y ( t , λ ) | M Γ ( α ) 0 t ( t s ) α   d s = M   t 1 α ( 1 α ) Γ ( α ) ,
and so I 0 + 1 α y ( t , λ ) | t = 0 = 0 . On the other hand (15) now implies that c 1 = 0 , i.e., y ( 0 , λ ) = 0 , so that y C [ 0 , 1 ] . However, the condition y ( 1 , λ ) = 0 is independent of the statement that I 0 + 1 α y ( t , λ ) | t = 1 = 0 .
Since, for any z 0 , the Mittag-Leffler functions satisfy
E δ , δ ( z ) = 1 z E δ , 0 ( z ) ,
we get
t 2 α 1 E 2 α , 2 α ( λ t 2 α ) = 1 λ t   E 2 α , 0 ( λ t 2 α ) .
Hence, using (28) and (35) we get
y ( t , λ ) = c 1   t α 1   E 2 α , α ( λ t 2 α ) + c 2   t α E 2 α , α + 1 ( λ t 2 α ) 0 t E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) q ( s ) y ( s , λ )   d s .
Remark 1.
When α 1 , the integral Equation (36) becomes
y ( t , λ ) = y ( 0 , λ ) cos ( λ t ) + y ( 0 , λ ) sin ( λ t ) λ + 0 t sin ( λ ( t s ) ) λ q ( s ) y ( s , λ ) d s ,
which is exactly the integral equation equivalent of the classical Sturm–Liouville equation y + q ( t ) y = λ y for λ > 0 .
Remark 2.
Observe that, for each α,
lim s t E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) = 0 , i f α ( 1 / 2 , 1 ] , 1 , i f α = 1 / 2 .
and so, for each 1 / 2 < α < 1 , the kernel appearing in (36) is uniformly bounded on [ 0 , 1 ] . This agrees with the equivalent result for the classical case (37).

7. Existence and Asymptotic Distribution of the Eigenvalues

Without loss of generality we may assume that c 2 = 1 in (36) and y ( t , λ ) is the corresponding solution. In the sequel we always assume that 1 / 2 < α < 1 .
Lemma 6.
For each t [ 0 , 1 ] , 1 / 2 < α < 1 , and | arg ( λ ) | μ where μ ( α π , π ) , we have t α E 2 α , α + 1 ( λ t 2 α ) 0 as | λ | .
Proof. 
By (14) we can write
t α E 2 α , α + 1 ( λ t 2 α ) = t α 1 2 α ( λ t 2 α ) 1 ( α + 1 ) 2 α exp ( λ t 2 α ) 1 2 α + O 1 λ = i 2 α λ exp ( λ ) 1 2 α t + O 1 λ = i 2 α λ exp | λ | 1 2 α cos ( arg ( λ ) 2 α ) + i sin ( arg ( λ ) 2 α ) t + O 1 λ .
Therefore,
| t α E 2 α , α + 1 ( λ t 2 α ) | = 1 2 α λ exp | λ | 1 2 α cos ( arg ( λ ) 2 α ) t .
Regarding the assumption on arg ( λ ) , we have cos ( arg ( λ ) 2 α ) < 0 and it completes the proof. □
Lemma 7.
For each t [ 0 , 1 ] , s [ 0 , t ] , 1 / 2 < α < 1 , and | arg ( λ ) | μ where μ ( α π , π ) , we have E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) 0 as | λ | .
Proof. 
By (14) we can write
E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) = 1 2 α ( λ ( t s ) 2 α ) 1 2 α exp ( λ ( t s ) 2 α ) 1 2 α + O 1 λ λ ( t s ) = 1 2 α ( λ ) 1 2 α λ exp ( λ ) 1 2 α ( t s ) + O 1 λ 2 .
Then,
E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) = 1 2 α | λ | ( 2 α 1 ) / 2 α   exp ( t s ) | λ | 1 / 2 α cos arg ( λ ) 2 α + O 1 | λ | 2 .
Arguing as in the previous lemma we reach the desired conclusion. □
Lemma 8.
For each t [ 0 , 1 ] , s [ 0 , t ] , 1 / 2 < α < 1 , and | arg ( λ ) | μ where μ ( α π , π ) , we have E 2 α , 1 α ( λ ( t s ) 2 α ) λ ( t s ) α 0 as | λ |
Proof. 
By (14) we can write
E 2 α , 1 α ( λ ( t s ) 2 α ) λ ( t s ) α = 1 2 α ( λ ( t s ) 2 α ) 1 ( 1 α ) 2 α exp ( λ ( t s ) 2 α ) 1 2 α + O 1 λ λ ( t s ) α = i λ 2 α exp ( λ ) 1 2 α ( t s ) + O 1 λ 2 .
Then,
| E 2 α , 1 α ( λ ( t s ) 2 α ) λ ( t s ) α | = 1 2 α | λ |   exp ( t s ) | λ | 1 / 2 α cos arg ( λ ) 2 α .
The result follows since the exponential term is uniformly bounded. □
Lemma 9.
For each t [ 0 , 1 ] , and 1 / 2 < α < 1 , the solution y ( t , λ ) is an entire function of λ of order at most 1 / 2 α .
Proof. 
Let λ C . Define f by
y ( t , λ ) = exp t | λ | 1 / 2 α cos arg ( λ ) 2 α f ( t ) .
Then, using (36),
f ( t ) = t α E 2 α , α + 1 ( λ t 2 α ) exp t | λ | 1 / 2 α cos arg ( λ ) 2 α 0 t E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) exp ( t s ) | λ | 1 / 2 α cos arg ( λ ) 2 α   q ( s ) f ( s )   d s
Applying Lemma 6 there exists Λ R + such that for all | λ | > Λ we have
| f ( t ) | 1 + 1 2 α | λ | ( 2 α 1 ) / 2 α 0 t | q ( s ) |   | f ( s ) |   d s
which, on account of Gronwall’s inequality, gives us
| f ( t ) | exp 1 2 α | λ | ( 2 α 1 ) / 2 α 0 1 | q ( s ) |   d s
for all sufficiently large | λ | . Thus, f L [ 0 , 1 ] so that (39) yields, for some M,
| y ( t , λ ) | M exp | λ | 1 / 2 α
and the order claim is verified. □
Lemma 10.
For each t [ 0 , 1 ] , I 0 + 1 α y ( t , λ ) is an entire function of λ of order at most 2 α .
Proof. 
This is clear from the definition, the possible values of α , and since y ( t , λ ) is itself entire and of order at most 1 / 2 α , from Lemma 9. □
Lemma 11.
The boundary value problem (1) and (2) has infinitely many complex eigenvalues (real eigenvalues are not to be excluded here).
Proof. 
By Lemma 10, we know that I 0 + 1 α y ( t ,   λ ) is entire for each t [ 0 , 1 ] , and 1 / 2 < α < 1 as well. So, the eigenvalues of our problem are given by the zeros of I 0 + 1 α y ( 1 ,   λ ) , which must be countably infinite in number since the latter function is of fractional order 1 / 2 α (on account of the restriction on α ). This gives us the existence of infinitely many eigenvalues, generally in C . □
Next, we give the asymptotic distribution of these eigenvalues when α is either very close to 1 / 2 from the right or very close to 1 from the left. Recall (36) with c 2 = 1 , so that
y ( t , λ ) = t α E 2 α , α + 1 ( λ t 2 α ) 0 t E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) q ( s ) y ( s , λ )   d s .
An iterative method for solving for approximate solutions of (41) maybe found in [16]. Keeping in mind the boundary condition (2) at t = 1 , we calculate I 0 + 1 α y ( t ,   λ ) and then evaluate this at t = 1 in order to find the dispersion relation for the eigenvalues. However, our derivation is theoretical in nature. A straightforward though lengthy calculation using (41) and the definition of the Mittag-Leffler functions show that
I 0 + 1 α y ( t ,   λ ) = I 0 + 1 α { t α E 2 α , α + 1 ( λ t 2 α ) } + I 0 + 1 α 0 t E 2 α , 0 ( λ ( t s ) 2 α ) λ ( t s ) q ( s ) y ( s , λ )   d s = t E 2 α , 2 ( λ t 2 α ) + 1 λ 0 t E 2 α , 1 α ( λ ( t s ) 2 α ) ( t s ) α q ( s ) y ( s , λ )   d s
so that the eigenvalues of (1) and (2) are given by those λ C such that
E 2 α , 2 ( λ ) + 1 λ 0 1 E 2 α , 1 α ( λ ( 1 s ) 2 α ) ( 1 s ) α q ( s ) y ( s , λ )   d s = 0 .
Let us consider first the case where λ R . Lemma 8 implies that the right side of (38) tends to 0 as λ . Indeed this, combined with (39), implies that
| E 2 α , 1 α ( λ ( t s ) 2 α ) λ ( t s ) α   y ( s , λ ) | = O 1 | λ |
for all sufficiently large λ .
Thus, the real eigenvalues of the problem (1) and (2) become the zeros of a transcendental equation of the form,
E 2 α , 2 ( λ ) + O 1 λ = 0 .
We are concerned with the asymptotic behaviour of these real zeros. Recall the distribution of the real zeros of E 2 α , 2 ( λ ) in [1]. There we showed that, for each n = 0 , 1 , 2 , , N 1 , where N depends on α , the interval
I n ( α ) : = ( 2 n + 1 2 + 1 2 α ) π sin ( π 2 α ) 2 α , ( 2 n + 3 2 + 1 2 α ) π sin ( π 2 α ) 2 α ,
always contains at least two real zeros of E 2 α , 2 ( λ ) . For α 1 , these intervals approach the intervals
( 2 n + 1 ) 2   π 2 , ( 2 n + 2 ) 2   π 2 ,
whose end-points are each eigenvalues of the Dirichlet problem for the classical equation y = λ   y on [ 0 , 1 ] . Since each interval I n contains two zeros we can denote the first of these two zeros by λ 2 n ( α ) . Equation (44) now gives the a-priori estimate
( 2 n + 1 2 + 1 2 α ) π sin ( π 2 α ) 2 α λ 2 n ( α ) ( 2 n + 3 2 + 1 2 α ) π sin ( π 2 α ) 2 α .
For each α < 1 , and close to 1, and for large λ , the real zeros of the preceding equation approach those of E 2 α , 2 ( λ ) and spread out towards the end-points of intervals of the form (44). For α close to 1 / 2 there are no zeros, the first two zeros appearing only when α 0.7325 . For α larger than this critical value, the zeros appear in pairs and in intervals of the form (44).
Next, recall that for α < 1 there are only finitely many such real zeros, (see [1]) their number growing without bound as α 1 . It also follows from Lemma 11 that, for each α , the remaining infinitely many eigenvalues must be non-real. As α 1 these non-real eigenvalues tend to the real axis thereby forming more and more real eigenvalues until the spectrum is totally real when α = 1 and the problem then reduces to a (classical) regular Sturm–Liouville problem.
Finally, for α close to 1, (45) leads to the approximation,
λ 2 n ( α ) ( 2 n + 2 ) π sin ( π 2 α ) 2 α ,
from which this, in conjunction with (44) and α 1 , we can derive the classical eigenvalue asymptotics, λ n n 2 π 2 as n .

8. Conclusions

We consider the fractional eigenvalue problem,
c D 0 + α D 0 + α y ( t ) + q ( t ) y ( t ) = λ y ( t ) , 0 t 1 ,
where α is a real parameter, 1 / 2 < α < 1 , λ is a generally unspecified complex parameter, with mixed Caputo and Riemann–Liouville derivatives and q an essentially bounded function, subject to the following boundary conditions involving the Riemann–Liouville integrals,
I 0 + 1 α y ( t ) | t = 0 = 0 , and I 0 + 1 α y ( t ) | t = 1 = 0 .
We show that this problem admits, for each α under consideration, and for eigenfunctions that are in C [ 0 , 1 ] , a finite number of real eigenvalues and an infinite number of non-real eigenvalues. The real eigenvalues, though finite in number for each α , are approximated by (44) and (45), which as α 1 gives the classical asymptotic relation λ n n 2 π 2 as n .
As α 1 we observe that the spectrum obtained approaches the Sturm–Liouville spectrum of the classical problem
y + q ( t ) y = λ y , y ( 0 ) = y ( 1 ) = 0 .
The same results hold if the eigenfunctions are merely C ( 0 , 1 ] (i.e., c 1 0 ) except that now the latter have an infinite discontinuity at t = 0 for each α . The proofs are identical and are therefore omitted.

Author Contributions

Conceptualization, M.D. and A.B.M.; methodology, M.D. and A.B.M.; software, M.D. and A.B.M.; validation, M.D. and A.B.M.; formal analysis, M.D. and A.B.M.; investigation, M.D. and A.B.M.; resources, M.D. and A.B.M.; data curation, Not applicable; writing—original draft preparation, M.D. and A.B.M.; writing—review and editing, M.D. and A.B.M.; visualization, Not applicable; supervision, M.D. and A.B.M.; project administration, M.D. and A.B.M.; funding acquisition, M.D. and A.B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by a grant from the Office of the Dean of Science, Carleton University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data were used to support this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dehghan, M.; Mingarelli, A.B. Fractional Sturm-Liouville eigenvalue problems, I. Rev. Real Acad. Cienc. Exactas Fis. Nat. Ser. A Mat. 2020, 114, 46. [Google Scholar]
  2. Aleroev, T.S.; Aleroeva, H.T. Problems of Sturm-Liouville type for differential equations with fractional derivatives. In Volume 2 Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 21–46. [Google Scholar] [CrossRef]
  3. Al-Mdallal, Q.M. An efficient method for solving fractional Sturm-Liouville problems. Chaos Solitons Fractals 2009, 40, 183–189. [Google Scholar]
  4. Blaszczyk, T.; Ciesielski, M. Numerical solution of fractional Sturm-Liouville equation in integral form. Fract. Calc. Appl. Anal. 2014, 17, 307–320. [Google Scholar]
  5. Heydarpour, Z.; Izadi, J.; George, R.; Ghaderi, M.; Rezapour, S. On a Partial Fractional Hybrid Version of Generalized Sturm-Liouville-Langevin Equation. Fractal Fract. 2022, 6, 269. [Google Scholar]
  6. Klimek, M.; Agrawal, O.P. On a regular fractional Sturm-Liouville problem with derivatives of order (0, 1). In Proceedings of the 13th International Carpathian Control Conference, High Tatras, Slovakia, 28–31 May 2012. [Google Scholar] [CrossRef]
  7. Klimek, M.; Agrawal, O.P. Fractional Sturm-Liouville problem. Comput. Math. Appl. 2013, 66, 795–812. [Google Scholar]
  8. Klimek, M.; Odzijewicz, T.; Malinowska, A.B. Variational methods for the fractional Sturm-Liouville problem. J. Math. Anal. Appl. 2014, 416, 402–426. [Google Scholar] [CrossRef]
  9. Klimek, M.; Ciesielski, M.; Blaszczyk, T. Exact and Numerical Solution of the Fractional Sturm-Liouville Problem with Neumann Boundary Conditions. Entropy 2022, 24, 143. [Google Scholar] [CrossRef] [PubMed]
  10. Min, D.; Chen, F. Variational methods to the p-Laplacian type nonlinear fractional order impulsive differential equations with Sturm-Liouville boundary-value problem. Fract. Calc. Appl. Anal. 2021, 24, 1069–1093. [Google Scholar] [CrossRef]
  11. Zayernouri, M.; Karniadakis, G.E. Fractional Sturm-Liouville eigen-problems: Theory and numerical approximation. J. Comput. Phys. 2013, 252, 495–517. [Google Scholar]
  12. Mittag-Leffler, G. Sur la nouvelle function Eα(x). C. R. Acad. Sci. Paris 1903, 137, 554–558. [Google Scholar]
  13. Erdélyi, A.; Bateman, H. Higher Transcendental Functions; McGraw-Hill Book Company, Inc.: New York, NY, USA, 1953; Volume 3. [Google Scholar]
  14. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Application of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  15. Maralani, E.M.; Saei, F.D.; Akbarfam, A.J.; Ghanbari, K. Eigenvalues of fractional Sturm-Liouville problems by successive method. Comput. Methods Differ. Equ. 2021, 9, 1163–1175. [Google Scholar]
  16. Ma, H.J.; Wang, Y. Full Information Control of Borel-Measurable Markov Jump. Mathematics 2022, 10, 37. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dehghan, M.; Mingarelli, A.B. Fractional Sturm–Liouville Eigenvalue Problems, II. Fractal Fract. 2022, 6, 487. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6090487

AMA Style

Dehghan M, Mingarelli AB. Fractional Sturm–Liouville Eigenvalue Problems, II. Fractal and Fractional. 2022; 6(9):487. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6090487

Chicago/Turabian Style

Dehghan, Mohammad, and Angelo B. Mingarelli. 2022. "Fractional Sturm–Liouville Eigenvalue Problems, II" Fractal and Fractional 6, no. 9: 487. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6090487

Article Metrics

Back to TopTop