Journal of Applied Mathematics and Physics
Vol.07 No.05(2019), Article ID:92572,10 pages
10.4236/jamp.2019.75072

Existence of Solutions for Boundary Value Problems of Vibration Equation with Fractional Derivative

Nan Yao, Zeyu Luo

College of Science, University of Shanghai for Science and Technology, Shanghai, China

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: April 8, 2019; Accepted: May 20, 2019; Published: May 23, 2019

ABSTRACT

In this paper, we investigate the solvability of boundary value problems for a class of vibration differential equation describing the fractional order damped system with signal stimulus. By presenting kernel function through the Laplace transform, and using the eigenvalue and the improved Leray-Schauder degree, the existence of solutions for boundary value problems is established.

Keywords:

Caputo Derivative, Boundary Value Problems for Vibration Equation, Laplace Transform, Eigenvalue, Improved Leray-Schauder Degree

1. Introduction

Fractional calculus has been applied more and more widely in many fields of science and engineering, many scholars have done a lot of research on it [1] - [8] . When describing fractional Brownian motion in Anomalous diffusion, if time fractional differential operator is introduced, fractional Langevin equation [9]

x ( t ) + γ d α 1 d t α 1 x ( t ) = Γ (t)

is obtained. The fractional Langevin equation describes both subdiffusion for 0 < α < 1 and superdiffusion for 1 < α < 2 .

Fractional differential equations are also used to describe damped vibrations in viscoelastic media. In [10] , Podlubny studied the initial value problem for the inhomogeneous Bagley-Torvik equation

{ A y ( t ) + B 0 D t 3 2 y ( t ) + C y ( t ) = f ( t ) , t > 0 , y ( 0 ) = 0 , y ( 0 ) = 0

and the numerical solutions are presented. The numerical solution is in agreement with the analytical solution, obtained with the help of the fractional Green’s function for a three-term fractional differential equation with constant coefficients.

Motivated by the above works, we study the following boundary value problems for a class of vibration differential equation describing the fractional order damped system with signal stimulus

{ x ( t ) b c D 0 + α x ( t ) = f ( t , x ( t ) ) , t ( 0 , 1 ) , x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) , (1.1)

where 0 < b < 1 , 1 < α < 2 , c 0 , c D 0 + α is the Caputo fractional derivative of order α , f : [ 0 , 1 ] × ( 0 , ) is continuous. Moreover, x ( t ) indicates the ratio of inertia force to mass, c D 0 + α x ( t ) is external damping term, namely, dissipative term, f represents the external force.

By using the Laplace transform, the kernel function is obtained. And then, by using the eigenvalue and the improved Leray-Schauder degree, the existence of the solutions to boundary value problem (1.1) is proved, see Theorem 1. So we can investigate the state of the oscillator motion under this system.

2. Preliminaries

In this part, we recall some definitions and lemmas which are critical to the existence result. The definitions of fractional integral and fractional derivative can be found in [11] [12] .

Definition 1. [10] Assume that function g ( t ) is defined in ( 0 , ) , then the Laplace transform of g ( t ) is defined as

g ˜ ( p ) = L [ g ( t ) ] ( p ) = 0 e p t g ( t ) d t , p > 0 ,

as long as the generalized integral is convergent.

Definition 2. [10] The original g ( t ) can be restored from the Laplace transform g ˜ ( p ) with the help of the inverse Laplace transform

g ( t ) = L 1 [ g ˜ ( p ) ] = 1 2 π i σ i σ + i g ˜ ( p ) e p t d p , σ = Re ( p ) > σ 0 .

where σ 0 lies in the right half plane of the absolute convergence of the Laplace integral.

Definition 3. [12] Let δ , β > 0 . The function

E δ , β ( z ) = k = 0 z k Γ ( δ k + β ) ,

whenever the series converges is called the two-parameter Mittag-Leffler function with parameters δ and β .

Lemma 1. [11] Let l 1 < δ l , l . The Laplace transform formula for c D 0 + δ g ( t ) is

L [ c D 0 + δ g ( t ) ] ( p ) = p δ L [ g ( t ) ] ( p ) j = 0 l 1 p δ 1 j g ( j ) ( 0 ) , p > 0.

Let

E δ , β ( k ) ( z ) = d k d z k E δ , β ( z ) , k = 0 , 1 , 2 , .

Lemma 2. [12] Let δ , β > 0 . The power series E δ , β ( z ) is convergent for all z . In other words, E δ , β ( z ) is an entire function.

Lemma 3. Let z , k = 0 , 1 , 2 , . Then

E δ , β ( k ) ( z ) = j = 0 z j Γ ( k + j + 1 ) Γ ( j + 1 ) Γ ( δ ( k + j ) + β ) .

Proof. By Definition 3 and Lemma 2, we can get

E δ , β ( k ) ( z ) = j = 0 d k d z k ( z j Γ ( δ j + β ) ) = j = k d k d z k ( Γ ( j + 1 ) z j k Γ ( j k + 1 ) Γ ( δ j + β ) ) = j = 0 z j Γ ( k + j + 1 ) Γ ( j + 1 ) Γ ( δ ( k + j ) + β ) .

Thus, the lemma can be obtained.

Lemma 4. [11] Let δ , β > 0 , E δ , β ( z ) be a two-parameter Mittag-Leffler function. Then

L [ t δ k + β 1 E δ , β ( k ) ( ± a t δ ) ] ( p ) = k ! p δ β ( p δ a ) k + 1 , p > | a | 1 δ .

Denote

g 1 ( t ) = t E 2 α , 2 ( b t 2 α ) , g 2 ( t ) = 1 b ( E 2 α , 1 ( b t 2 α ) 1 ) , t [ 0 , ) .

Lemma 5. The functions g 1 and g 2 , defined above, have the following properties.

1) g 1 ( t ) and g 2 ( t ) are represented by absolutely and uniformly convergent series and g 1 ( t ) , g 2 ( t ) 0 on [ 0 , 1 ] ;

2) g 1 ( t ) = E 2 α , 1 ( b t 2 α ) > 1 ;

3) g 1 ( t ) b g 2 ( t ) = 1 .

Proof. 1) By Lemma 2, we can show g 1 ( t ) is an entire function. Thus, g 1 ( t ) is represented by absolutely and uniformly convergent series on [ 0 , 1 ] .

Similarly, g 2 ( t ) is also represented by absolutely and uniformly convergent series on [ 0 , 1 ] . And we can easily have g 1 ( t ) , g 2 ( t ) 0 .

2) In view of g 1 ( t ) is represented by absolutely and uniformly convergent series on [ 0 , 1 ] ,

g 1 ( t ) = n = 0 b n t ( 2 α ) n Γ ( ( 2 α ) n + 1 ) = E 2 α , 1 ( b t 2 α ) > 1 .

3) From the definition of g 1 ( t ) , we have

g 1 ( t ) b g 2 ( t ) = E 2 α , 1 ( b t 2 α ) b 1 b ( E 2 α , 1 ( b t 2 α ) 1 ) = 1 .

The proof is complete.

Lemma 6. For y ( t ) is continuous on [ 0 , 1 ] and c g 1 ( 1 ) 0 , the unique solution of

{ x ( t ) b c D 0 + α x ( t ) = y ( t ) , t ( 0 , 1 ) , x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) (2.1)

is

x ( t ) = 0 1 G ( t , s ) y ( s ) d s , t [ 0 , 1 ] , (2.2)

where

G ( t , s ) = { g 1 ( t ) g 1 ( 1 s ) c g 1 ( 1 ) + g 1 ( t s ) , 0 s t 1 , g 1 ( t ) g 1 ( 1 s ) c g 1 ( 1 ) , 0 t < s 1. (2.3)

Proof. By Lemma 1, we have

L [ x ( t ) ] ( p ) = p 2 L [ x ( t ) ] ( p ) x ( 0 ) p x ( 0 ) , L [ c D 0 + α x ( t ) ] ( p ) = p α L [ x ( t ) ] ( p ) x ( 0 ) p α 1 .

Apply Laplace transform to both sides of x ( t ) b c D 0 + α x ( t ) = y ( t ) , we can easily obtain

p 2 L [ x ( t ) ] ( p ) x ( 0 ) p x ( 0 ) b ( p α L [ x ( t ) ] ( p ) x ( 0 ) p α 1 ) = L [ y ( t ) ] ( p ) ,

namely,

L [ x ( t ) ] ( p ) = x ( 0 ) p + x ( 0 ) b x ( 0 ) p α 1 + L [ y ( t ) ] ( p ) p 2 b p α . (2.4)

If | b p α 2 | < 1 , we have

1 p 2 b p α = p 2 1 b p α 2 = p 2 n = 0 ( b p α 2 ) n = n = 0 b n p n ( α 2 ) 2 ,

By virtue of Lemma 4, we can show

1 p 2 b p α = n = 0 b n p n ( α 2 ) 2 = L [ g 1 ( t ) ] ( p ) , p > 0.

Similarly, if b p α 2 < 1 and p > 0 , we have

p p 2 b p α = n = 0 b n p n ( α 2 ) 1 = L [ g 1 ( t ) ] ( p ) , p α 1 p 2 b p α = n = 0 b n p n ( α 2 ) + α 3 = L [ g 2 ( t ) ] ( p ) .

So (2.4) is equivalent to

L [ x ( t ) ] ( p ) = x ( 0 ) L [ g 1 ( t ) ] ( p ) + x ( 0 ) L [ g 1 ( t ) ] ( p ) b x ( 0 ) L [ g 2 ( t ) ] ( p ) + L [ y ( t ) ] ( p ) L [ g 1 ( t ) ] ( p ) . (2.5)

Furthermore,

L [ y ( t ) ] ( p ) L [ g 1 ( t ) ] ( p ) = L [ 0 t g 1 ( t s ) y ( s ) d s ] ( p ) ,

then we can get the inverse Laplace transform for (2.5) is

x ( t ) = x ( 0 ) g 1 ( t ) + x ( 0 ) g 1 ( t ) b x ( 0 ) g 2 ( t ) + 0 t g 1 ( t s ) y ( s ) d s = x ( 0 ) ( g 1 ( t ) b g 2 ( t ) ) + x ( 0 ) g 1 ( t ) + 0 t g 1 ( t s ) y ( s ) d s = x ( 0 ) + x ( 0 ) g 1 ( t ) + 0 t g 1 ( t s ) y ( s ) d s . (2.6)

Because x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) and c g 1 ( 1 ) 0 , we can show

x ( 0 ) = 0 , x ( 0 ) = 0 1 g 1 ( 1 s ) y ( s ) d s c g 1 ( 1 ) . (2.7)

Substituting (2.7) into (2.6), we get

x ( t ) = 0 1 G ( t , s ) y ( s ) d s ,

where G ( t , s ) is defined by (2.3).

On the other hand, by using the above proof, if x = x ( t ) satisfies (2.2), we obtain that x satisfies x ( t ) b c D 0 + α x ( t ) = y ( t ) and x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) . The proof is complete.

Lemma 7. The eigenfunction of

{ x ( t ) b c D 0 + α x ( t ) = λ x ( t ) , t ( 0 , 1 ) , x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) (2.8)

is x ( t ) = K n = 0 λ n n ! t 2 n + 1 E 2 α , 2 + α n ( n ) ( b t 2 α ) and its corresponding eigenvalue λ is the solution of equation

c = n = 0 λ n n ! E 2 α , 1 + α n ( n ) ( b ) ,

where K is a constant and λ > 0 .

Proof. Let x = x ( t ) is the solution of boundary value problem (2.8). Apply Laplace transform to both sides of x ( t ) b c D 0 + α x ( t ) = λ x ( t ) , we can easily obtain

L [ x ( t ) ] ( p ) = x ( 0 ) p + x ( 0 ) b x ( 0 ) p α 1 p 2 λ b p α . (2.9)

If | λ p 2 b p α | < 1 , we have

1 p 2 λ b p α = p α p 2 α b 1 1 λ p α p 2 α b = n = 0 λ n p α ( n + 1 ) ( p 2 α b ) n + 1 .

By virtue of Lemma 4, we can show, if p > λ 1 2 α ,

1 p 2 λ b p α = n = 0 λ n p α ( n + 1 ) ( p 2 α b ) n + 1 = L [ n = 0 λ n n ! t 2 n + 1 E 2 α , 2 + α n ( n ) ( b t 2 α ) ] ( p ) .

Similarly, if | λ p 2 b p α | < 1 and p > λ 1 2 α ,

p p 2 λ b p α = n = 0 λ n p 1 α ( n + 1 ) ( p 2 α b ) n + 1 = L [ n = 0 ( λ t 2 ) n n ! E 2 α , 1 + α n ( n ) ( b t 2 α ) ] ( p ) , p α 1 p 2 λ b p α = n = 0 λ n p α n 1 ( p 2 α b ) n + 1 = L [ n = 0 λ n n ! t 2 n + 2 α E 2 α , 3 α + α n ( n ) ( b t 2 α ) ] ( p ) .

So (2.9) is equivalent to

L [ x ( t ) ] ( p ) = x ( 0 ) L [ n = 0 λ n n ! t 2 n + 1 E 2 α , 2 + α n ( n ) ( b t 2 α ) ] ( p ) + x ( 0 ) L [ n = 0 ( λ t 2 ) n n ! E 2 α , 1 + α n ( n ) ( b t 2 α ) ] ( p ) b x ( 0 ) L [ n = 0 λ n n ! t 2 n + 2 α E 2 α , 3 α + α n ( n ) ( b t 2 α ) ] ( p ) . (2.10)

Furthermore, we can get the inverse Laplace transform for (2.10) is

x ( t ) = x ( 0 ) n = 0 λ n n ! t 2 n + 1 E 2 α , 2 + α n ( n ) ( b t 2 α ) + x ( 0 ) n = 0 ( λ t 2 ) n n ! E 2 α , 1 + α n ( n ) ( b t 2 α ) b x ( 0 ) n = 0 λ n n ! t 2 n + 2 α E 2 α , 3 α + α n ( n ) ( b t 2 α ) .

Because x ( 0 ) = 0 , x ( 1 ) = c x ( 0 ) , we can show

x ( t ) = K n = 0 λ n n ! t 2 n + 1 E 2 α , 2 + α n ( n ) ( b t 2 α ) , c = n = 0 λ n n ! E 2 α , 1 + α n ( n ) ( b ) .

The proof is complete.

Lemma 8. The function G ( t , s ) defined by (2.3) is continuous on [ 0 , 1 ] × [ 0 , 1 ] .

Proof. By the definition of G ( t , s ) and Lemma 5, we get G ( t , s ) is continuous for ( t , s ) [ 0 , 1 ] × [ 0 , 1 ] .The proof is complete.

3. The Existence of the Solutions

Let

β = lim x sup t [ 0 , 1 ] f ( t , x ) x .

Throughout this paper, we always suppose that the following conditions are satisfied.

(H1) There exists constant r > 0 , τ ( 0 , 1 ) such that f ( , x ) > r τ g 1 ( 1 τ ) for any | x | < r .

(H2) There exists λ > 0 such that λ β , here λ satisfies

c = n = 0 λ n n ! E 2 α , 1 + α n ( n ) ( b ) .

Let E = C [ 0 , 1 ] , with the norm x = max t [ 0 , 1 ] | x ( t ) | . Obviously, ( E , ) is a Banach space.

Define the operators T , A : E E ,

T x ( t ) = 0 1 G ( t , s ) x ( s ) d s , A x ( t ) = 0 1 G ( t , s ) f ( s , x ( s ) ) d s .

By virtue of Lemma 6, the solution of boundary value problem (1.1) is equivalent to the fixed point of the operator A; Boundary value problem (2.8) is equivalent to the following integral equation

x ( t ) = λ 0 1 G ( t , s ) x ( s ) d s . (3.1)

Therefore, λ T x = x , we have 1 λ is the eigenvalue of operator T corresponding to the eigenfunction (3.1).

Lemma 9. A : E E is completely continuous.

Proof. Let G 0 = max t , s [ 0 , 1 ] | G ( t , s ) | . Obviously, A : E E .

Let { x n } E , x E such that x n x 0 as n . So there exists ρ > 0 such that x n ρ , x ρ .

Let D = { x E ; x ρ } , M = max t [ 0 , 1 ] , x [ ρ , ρ ] f ( t , x ) . Then

| A x n ( t ) A x ( t ) | = | 0 1 G ( t , s ) ( f ( s , x n ( s ) ) f ( s , x ( s ) ) ) d s | 0 1 | G ( t , s ) f ( s , x n ( s ) ) f ( s , x ( s ) ) | d s ,

| G ( t , s ) f ( s , x n ( s ) ) f ( s , x ( s ) ) | 2 M G 0

By virtue of Lebesgue’s dominated convergence theorem, we have

lim n 0 1 | G ( t , s ) f ( s , x n ( s ) ) f ( s , x ( s ) ) | d s = 0 1 | G ( t , s ) | lim n | f ( s , x n ( s ) ) f ( s , x ( s ) ) | d s = 0 ,

so A x n A x 0 as n + . Hence, the operator A is continuous.

For each x in the bounded area D,

| A x ( t ) | = | 0 1 G ( t , s ) f ( s , x ( s ) ) d s | G 0 M .

Consequently, the operator A is uniformly bounded.

By the continuity of G ( t , s ) on [ 0 , 1 ] × [ 0 , 1 ] , ε > 0 , δ > 0 for any t 1 , t 2 , s 1 , s 2 [ 0 , 1 ] , if | t 1 t 2 | < δ , | s 1 s 2 | < δ , then we have

| G ( t 1 , s 1 ) G ( t 2 , s 2 ) | < ε .

If | t 1 t 2 | < δ , x D , we obtain

| A x ( t 2 ) A x ( t 1 ) | = | 0 1 ( G ( t 2 , s ) G ( t 1 , s ) ) f ( s , x ( s ) ) d s | ε M .

Then, through the Arzela-Ascoli theorem, the operator A is compact on D.

To summarize, A : E E is completely continuous. The proof is complete.

Lemma 10. The operator A is Frechét differentiable at , and A ( ) = β T .

Proof. Since A ( β ) = lim x sup t [ 0 , 1 ] f ( t , x ) x , then for any ε > 0 and | x | > N , there exists N > 0 such that | f ( t , x ) x β | < ε 2 , for any t [ 0 , 1 ] . Namely, | f ( t , x ) β x | < ε 2 | x | . Let M 0 = max t [ 0 , 1 ] , | x | N | f ( t , x ) β x | . Then for any x , | f ( t , x ) β x | < ε 2 | x | + M 0 .

So we have

| ( A x β T x ) ( t ) | = | 0 1 G ( t , s ) ( f ( s , x ( s ) ) β x ( s ) ) d s | ( ε 2 x + M 0 ) G 0 .

Thus, ( A x β T x ) x ( ε 2 x + M 0 ) G 0 x = ε 2 + M 0 G 0 x .

For the above ε , there exists N 0 > 0 such that M 0 G 0 x < ε 2 , namely, ( A x β T x ) x ε for x > N 0 . Thus, we can show

lim x A x β T x x = 0.

The proof is complete.

Lemma 11. [13] Let Ω be a bounded open set in infinite dimensional real Banach space E, θ Ω and A : Ω ¯ E be completely continuous. Suppose that A x x , A x x , x Ω . Then deg ( I A , Ω , θ ) = 0 .

Lemma 12. [14] Let A be a completely continuous operator which is defined on a Banach space E. Assume that 1 is not an eigenvalue of the asymptotic derivative. The completely continuous vector field I A is then nonsingular on spheres S ρ = { x | x = ρ } of sufficiently large radius ρ and

deg ( I A , B ( θ , ρ ) , θ ) = ( 1 ) k ,

where k is the sum of the algebraic multiplicities of the real eigenvalues of A ( ) in ( 1 , + ) .

Generalizing the previous lemmas, we obtain the following result.

Theorem 1. If (H1) and (H2) hold, then boundary value problem (1.1) has at least one nontrivial solution.

Proof. Obviously B ( θ , r ) is bounded open set, and θ ( B ( θ , r ) ) . Via Lemma 9, we get A : B ( θ , r ) ¯ E is completely continuous.

Combining (H2) and Lemma 10, we obtain the eigenvalue of A ( ) is β λ 1 . Therefore, through Lemma 12, we get

deg ( I A , B ( θ , ρ ) , θ ) = ( 1 ) k .

By Lemma 7, we have c = n = 0 λ n n ! E 2 α , 1 + α n ( n ) ( b ) . Therefore,

c = n = 0 λ n n ! E 2 α , 1 + α n ( n ) ( b ) = n = 0 λ n n ! j = 0 b j Γ ( n + j + 1 ) Γ ( j + 1 ) Γ ( ( 2 α ) j + 2 n + 1 ) > n = 0 b n Γ ( ( 2 α ) n + 1 ) = g 1 ( 1 ) .

Through the definition of G ( t , s ) and Lemma 5, we get G ( t , s ) 0 for ( t , s ) [ 0 , 1 ] × [ 0 , 1 ] . Considering (H1), for any x = r ,

| ( A x ) ( 1 ) | = 0 1 G ( 1 , s ) f ( s , x ( s ) ) d s 0 τ g 1 ( 1 s ) f ( s , x ( s ) ) d s τ g 1 ( 1 τ ) r τ g 1 ( 1 τ ) = x .

i.e. A x x . By Lemma 11, we get

deg ( I A , B ( θ , r ) , θ ) = 0.

In conclusion, deg ( I A , B ( θ , ρ ) \ B ( θ , r ) , θ ) = ( 1 ) k 0 = ( 1 ) k 0 .

So we get at least one x B ( θ , ρ ) \ B ( θ , r ) is a fixed point of the operator A. That is to say, x is one nontrivial solution of nonlinear problem (1.1). The proof is complete.

4. Conclusion

Theorem 1 is the main result of this paper. By Theorem 1, boundary value problem (1.1) has at least one nontrivial solution under the conditions of (H1) and (H2). Because boundary value problem (1.1) has at least one nontrivial solution, we can investigate the state of the oscillator motion under this system in the later research.

Acknowledgements

We thank the Editor and the referee for their comments.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Yao, N. and Luo, Z.Y. (2019) Existence of Solutions for Boundary Value Problems of Vibration Equation with Fractional Derivative. Journal of Applied Mathematics and Physics, 7, 1067-1076. https://doi.org/10.4236/jamp.2019.75072

References

  1. 1. Bai, Z., Chen, Y., Lian, H., et al. (2014) On the Existence of Blow up Solutions for a Class of Fractional Differential Equations. Fractional Calculus and Applied Analysis, 17, 1175-1187. https://doi.org/10.2478/s13540-014-0220-2

  2. 2. Bai, Z., Dong, X. and Yin, C. (2016) Existence Results for Impulsive Nonlinear Fractional Differential Equation with Mixed Boundary Conditions. Boundary Value Problems, 2016, Article ID: 63. https://doi.org/10.1186/s13661-016-0573-z

  3. 3. Liu, Y. (2016) New Existence Results for Positive Solutions of Boundary Value Problems for Coupled Systems of Multi-Term Fractional Differential Equations. Hacettepe Journal of Mathematics and Statistics, 45, 391-416. https://doi.org/10.15672/HJMS.20164512499

  4. 4. Liu, X. and Jia, M. (2018) The Method of Lower and Upper Solutions for the General Boundary Value Problems of Fractional Differential Equations with P-Laplacian. Advances in Difference Equations, 2018, Article ID: 28. https://doi.org/10.1186/s13662-017-1446-1

  5. 5. Han, Z., Lu, H. and Zhang, C. (2015) Positive Solutions for Eigenvalue Problems of Fractional Differential Equation with Generalized P-Laplacian. Applied Mathematics and Computation, 257, 526-536. https://doi.org/10.1016/j.amc.2015.01.013

  6. 6. Liu, X., Jia, M. and Ge, W. (2017) The Method of Lower and Upper Solutions for Mixed Fractional Four-Point Boundary Value Problem with P-Laplacian Operator. Applied Mathematics Letters, 65, 56-62. https://doi.org/10.1016/j.aml.2016.10.001

  7. 7. Jia, M. and Liu, X. (2014) Multiplicity of Solutions for Integral Boundary Value Problems of Fractional Differential Equations with Upper and Lower Solutions. Applied Mathematics and Computation, 232, 313-323. https://doi.org/10.1016/j.amc.2014.01.073

  8. 8. Liu, X., Jia, M. and Ge, W. (2013) Multiple Solutions of a P-Laplacian Model Involving a Fractional Derivative. Advances in Difference Equations, 126, 1-12. https://doi.org/10.1186/1687-1847-2013-126

  9. 9. Lutz, E. (2001) Fractional Langevin Equation. Physical Review E, 64, Article ID: 051106. https://doi.org/10.1103/PhysRevE.64.051106

  10. 10. Podluny, I. (1999) Fractional Differential Equations. Academic Press, San Diego.

  11. 11. Kilbas, A.A., Srivastava, H.M. and Trujillo, J.J. (2006) Theory and Applications of Fractional Differential Equations. Elsevier B V, Amsterdam.

  12. 12. Diethelm, K. (2010) The Analysis of Fractional Differential Equations. Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-642-14574-2

  13. 13. Guo, D. and Lakshmikantham, V. (2014) Nonlinear Problems in Abstract Cones. Academic Press, Cambridge.

  14. 14. Krasnosel’skii, M.A. and Zabreiko, P.P. (1984) Geometrical Methods of Nonlinear Analysis. Springer-Verlag, Berlin, Heidelberg.