Next Article in Journal
General Properties on Differential Sets of a Graph
Next Article in Special Issue
Hybrid Method for Simulation of a Fractional COVID-19 Model with Real Case Application
Previous Article in Journal
Multi-Layer Fuzzy Sustainable Decision Approach for Outsourcing Manufacturer Selection in Apparel and Textile Supply Chain
Previous Article in Special Issue
Global Stability Condition for the Disease-Free Equilibrium Point of Fractional Epidemiological Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Periodic Fractional (p, q)-Integral Boundary Value Problems for Sequential Fractional (p, q)-Integrodifference Equations

by
Jarunee Soontharanon
1 and
Thanin Sitthiwirattham
2,*
1
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
2
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(4), 264; https://doi.org/10.3390/axioms10040264
Submission received: 8 September 2021 / Revised: 13 October 2021 / Accepted: 18 October 2021 / Published: 19 October 2021

Abstract

:
We study the existence results of a fractional (p, q)-integrodifference equation with periodic fractional (p, q)-integral boundary condition by using Banach and Schauder’s fixed point theorems. Some properties of (p, q)-integral are also presented in this paper as a tool for our calculations.

1. Introduction

The studies of quantum calculus with integer order were presented in the last three decades, and many researchers extensively studied calculus without a limit that deals with a set of nondifferentiable functions, the so-called quantum calculus. Many types of quantum difference operators are employed in several applications of mathematical areas, such as the calculus of variations, particle physics, quantum mechanics, and theory of relativity. The q-calculus, one type of quantum initiated by Jackson [1,2,3,4,5], was employed in several fields of applied sciences and engineering such as physical problems, dynamical system, control theory, electrical networks, economics, and so on [6,7,8,9,10,11,12,13,14].
For fractional quantum calculus, Agarwal [15] and Al-Salam [16] proposed fractional q-calculus, and Díaz and Osler [17] proposed fractional difference calculus. In 2017, Brikshavana and Sitthiwirattham [18] introduced fractional Hahn difference calculus. In 2019, Patanarapeelert and Sitthiwirattham [19] studied fractional symmetric Hahn difference calculus.
Later, the motivation of quantum calculus based on two parameters (p, q)-integer was presented. The (p, q)-calculus (postquantum calculus) was introduced by Chakrabarti and Jagannathan [20]. This calculus was used in many fields such as special functions, approximation theory, physical sciences, Lie group, hypergeometric series, Bézier curves, and surfaces. For some recent papers about (p, q)-differenceequations, we refer to [21,22,23,24,25,26,27,28,29,30,31,32,33] and the references therein. For example, the fundamental theorems of (p, q)-calculus and some (p, q)-Taylor formulas were studied in [21]. In [32], the (p, q)-Melin transform and its applications were studied. The Picard and Gauss–Weierstrass singular integral in (p, q)-calculus were introduced in [33]. For the boundary value problem for ( p , q ) -difference equations were studied in [34,35,36]. For example, the nonlocal boundary value problems for first-order (p, q)-difference equations were studied in [34]. The second-order (p, q)-difference equations with separated boundary conditions were studied in [35]. In [36], the authors studied the first-order and second-order (p, q)-difference equations with impulse.
Recently, Soontharanon and Sitthiwirattham [37] introduced the fractional (p, q)-difference operators and its properties. Now, this calculus was used in the inequalities [38,39] and the boundary value problems [40,41,42]. However, the study of the boundary value problems for fractional (p, q)-difference equation in the beginning, there are a few literature on this knowledge. In [40], the existence results of a fractional (p, q)-integrodifference equation with Robin boundary condition were studied in 2020. In 2021 [41], the authors investigated the boundary value problem of a class of fractional (p, q)-difference Schr o ¨ dinger equations. In the same year, the existence results of solution and positive solution for the boundary value problem of a class of fractional (p, q)-difference equations involving the Riemann–Liouville fractional derivative [42] were studied.
Motivated by the above papers, we seek to enrich the contributions in this new research area. In this paper, we introduce and study the boundary value problem involving function F, which depends on fractional (p, q)-integral and fractional (p, q)-difference, and the boundary condition is nonlocal. Our problem is sequential fractional (p, q)-integrodifference equation with periodic fractional (p, q)-integral boundary conditions of the form
D p , q α D p , q β u ( t ) = F t , u ( t ) , Ψ p , q γ u ( t ) , D p , q ν u ( t ) , t I p , q T , u 0 = u T p I p , q θ g ( η ) u ( η ) = φ u , η I p , q T 0 , T ,
where I p , q T : = q p k T p : k N 0 { 0 } ; 0 < q < p 1 ; α , β , γ , ν , θ ( 0 , 1 ] ; F C I p , q T × R × R × R , R , g I p , q T , R + are given functions; φ : C I p , q T , R R is given functional; and for ϕ C I p , q T × I p , q T , [ 0 , ) , we define an operator of the (p, q)-integral of the product of functions ϕ and u as
Ψ p , q γ u ( t ) : = I p , q γ ϕ u ( t ) = 1 p γ 2 Γ p , q ( γ ) 0 t ( t q s ) p , q γ 1 ̲ ϕ t , s p γ 1 u s p γ 1 d p , q s .
We aim to show the existence results to the problem (1). Firstly, we convert the given nonlinear problem (1) into a fixed point problem related to (1), by considering a linear variant of the problem at hand. Once the fixed point operator is available, we make use the classical Banach’s and Schauder’s fixed point theorems to establish existence results.
The paper is organized as follows: Section 2 contains some preliminary concepts related to our problem. We present the existence and uniqueness result in Section 2, and the existence of at least one solution in Section 4. To illustrate our results, we provide some examples in Section 5. Finally, Section 6 discusses our conclusions.

2. Preliminaries

In this section, we provide some basic definitions, notations, and lemmas as follows. For 0 < q < p 1 , we define
[ k ] q : = 1 q k 1 q , k N 1 , k = 0 , [ k ] p , q : = p k q k p q = p k 1 [ k ] q p , k N 1 , k = 0 , [ k ] p , q ! : = [ k ] p , q [ k 1 ] p , q · · · [ 1 ] p , q = i = 1 k p i q i p q , k N 1 , k = 0 .
The (p, q)-forward jump and the (p, q)-backward jump operators are defined as
σ p , q k ( t ) : = q p k t and ρ p , q k ( t ) : = p q k t , for k N , respectively.
The q-analogue of the power function ( a b ) q n ̲ with n N 0 : = { 0 , 1 , 2 , } is given by
( a b ) q 0 ̲ : = 1 , ( a b ) q n ̲ : = i = 0 n 1 ( a b q i ) , a , b R .
The (p, q)-analogue of the power function ( a b ) p , q n ̲ with n N 0 is given by
( a b ) p , q 0 ̲ : = 1 , ( a b ) p , q n ̲ : = k = 0 n 1 ( a p k b q k ) , a , b R .
Generally, for α R , we define
( a b ) q α ̲ = a α i = 0 1 b a q i 1 b a q α + i , a 0 .
( a b ) p , q α ̲ = p α 2 ( a b ) q p α ̲ = a α i = 0 1 p α 1 b a q p i 1 b a q p i + α , a 0 .
In particular, a q α ̲ = a α , a p , q α ̲ = a p α and ( 0 ) q α ̲ = ( 0 ) p , q α ̲ = 0 for α > 0 .
The (p, q)-gamma and (p, q)-beta functions are defined by
Γ p , q ( x ) : = ( p q ) p , q x 1 ̲ ( p q ) x 1 = 1 q p p , q x 1 ̲ 1 q p x 1 , x R \ { 0 , 1 , 2 , } [ x 1 ] p , q ! , x N , B p , q ( x , y ) : = 0 1 t x 1 ( 1 q t ) p , q y 1 ̲ d p , q t = p 1 2 ( y 1 ) ( 2 x + y 2 ) Γ p , q ( x ) Γ q p , q ( y ) Γ p , q ( x + y ) ,
respectively.
Definition 1.
For 0 < q < p 1 and f : [ 0 , T ] R , we define the ( p , q ) -difference of f as
D p , q f ( t ) : = f ( p t ) f ( q t ) ( p q ) ( t ) , for t 0 f ( 0 ) , for t = 0
provided that f is differentiable at 0 and f is called ( p , q ) -differentiable on I p , q T if D p , q f ( t ) exists for all t I p , q T .
Observe that the function g ( t ) = D p , q f ( t ) is defined on [ 0 , T / p ] .
Definition 2.
Let I be any closed interval of R containing a , b and 0. Assuming that f : I R is a given function, we define ( p , q ) -integral of f from a to b by
a b f ( t ) d p , q t : = 0 b f ( t ) d p , q t 0 a f ( t ) d p , q t ,
where
I p , q f ( x ) = 0 x f ( t ) d p , q t = ( p q ) x k = 0 q k p k + 1 f q k p k + 1 x , x I ,
provided that the series converges at x = a and x = b and f is called ( p , q ) -integrable on [ a , b ] if it is ( p , q ) -integrable on [ a , b ] for all a , b I .
An operator I p , q N is defined as
I p , q 0 f ( x ) = f ( x ) a n d I p , q N f ( x ) = I p , q I p , q N 1 f ( x ) , N N .
The relations between (p, q)-difference and (p, q)-integral operators are given by
D p , q I p , q f ( x ) = f ( x ) a n d I p , q D p , q f ( x ) = f ( x ) f ( 0 ) .
Fractional (p, q)-integral and fractional (p, q)-difference of Riemann–Liouville type are defined as follows.
Definition 3.
For α > 0 , 0 < q < p 1 and f defined on I p , q T , the fractional ( p , q ) -integral is defined by
I p , q α f ( t ) : = 1 p α 2 Γ p , q ( α ) 0 t ( t q s ) p , q α 1 ̲ f s p α 1 d p , q s = ( p q ) t p α 2 Γ p , q ( α ) k = 0 q k p k + 1 t q p k + 1 t p , q α 1 ̲ f q k p k + α t ,
and ( I p , q 0 f ) ( t ) = f ( t ) .
Definition 4.
For α > 0 , 0 < q < p 1 and f defined on I p , q T , the fractional ( p , q ) -difference operator of Riemann–Liouville type of order α is defined by
D p , q α f ( t ) : = D p , q N I p , q N α f ( t ) = 1 p α 2 Γ p , q ( α ) 0 t ( t q s ) p , q α 1 ̲ f s p α 1 d p , q s ,
and D p , q 0 f ( t ) = f ( t ) , where N 1 < α < N , N N .
Lemma 1
([37]). Let α ( N 1 , N ) , N N , 0 < q < p 1 and f : I p , q T R . Then,
I p , q α D p , q α f ( t ) = f ( t ) + C 1 t α 1 + C 2 t α 2 + + C N t α N
for some C i R , i = 1 , 2 , , N .
Lemma 2
([37]). Let 0 < q < p 1 and f : I p , q T R be continuous at 0. Then,
0 x 0 s f ( τ ) d p , q τ d p , q s = 0 x p p q τ x f ( τ ) d p , q s d p , q τ .
Lemma 3
([37]). Let α , β > 0 , 0 < q < p 1 . Then,
( a ) 0 t ( t q s ) p , q α 1 ̲ s β d p , q s = t α + β B p , q ( β + 1 , α ) , ( b ) 0 t 0 x ( t q x ) p , q α 1 ̲ ( x q s ) p , q β 1 ̲ d p , q s d p , q x = B p , q ( β + 1 , α ) [ β ] p , q t α + β .
Lemma 4
([40]). Let α , β > 0 , 0 < q < p 1 and n Z . Then,
( a ) 0 t ( t q s ) p , q α 1 ̲ d p , q s = p α 2 Γ p , q ( α ) Γ p , q ( α + 1 ) t α , ( b ) 0 t 0 x p β 1 ( t q x ) p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x = p α 2 + β 2 Γ p , q ( α ) Γ p , q ( α + 1 ) t α + β , ( c ) 0 t ( t q s ) p , q β 1 ̲ s p β 1 α n d p , q s = p β 2 Γ p , q ( α n + 1 ) Γ p , q ( β ) Γ p , q ( α β n + 1 ) t α β n .
Lemma 5.
Let α , β , θ > 0 , 0 < q < p 1 and n Z . Then,
( a ) 0 t 0 x p θ 1 ( t q x ) p , q θ 1 ̲ x p θ 1 q s p , q β 1 ̲ s p α 1 p , q α n ̲ d p , q s d p , q x = p β 2 + θ 2 Γ p , q ( α n + 1 ) Γ p , q ( β ) Γ p , q ( θ ) Γ p , q ( α + β + θ n + 1 ) t α + β + θ n , ( b ) 0 t 0 y p θ 1 0 x p β 1 ( t q y ) p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x d p , q y = p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) Γ p , q ( α + β + θ + 1 ) t α + β + θ .
Proof. 
By Lemmas 2, 3 and 4 and definition of the (p, q)-beta function, we have
( a ) 0 t 0 x p θ 1 ( t q x ) p , q θ 1 ̲ x p θ 1 q s p , q β 1 ̲ s p α 1 p , q α n ̲ d p , q s d p , q x = 0 t ( t q x ) p , q θ 1 ̲ 0 x p θ 1 x p θ 1 q s p , q β 1 ̲ s p α 1 p , q α n ̲ d p , q s d p , q x = p β 2 p ( θ 1 ) ( α + β n ) · Γ p , q ( α n + 1 ) Γ p , q ( β ) Γ p , q ( α + β n + 1 ) 0 t ( t q x ) p , q θ 1 ̲ x α + β n d p , q x = p β 2 + θ 2 Γ p , q ( α n + 1 ) Γ p , q ( β ) Γ p , q ( θ ) Γ p , q ( α + β + θ n + 1 ) t α + β + θ n , ( b ) 0 t 0 y p θ 1 0 x p β 1 ( t q y ) p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x d p , q y = 0 t 0 y p θ 1 ( t q y ) p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ 0 x p β 1 x p β 1 q s p , q α 1 ̲ d p , q s d p , q x d p , q y = p α 2 Γ p , q ( α ) Γ p , q ( α + 1 ) 0 t 0 y p θ 1 ( t q y ) p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ x α p β 1 d p , q x d p , q y = p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) Γ p , q ( α + β + θ + 1 ) t α + β + θ .
The proof is complete. □
The following lemma, dealing with a linear variant of problem (1), plays an important role in the forthcoming analysis.
Lemma 6.
Let Ω 0 , α , β , θ ( 0 , 1 ] , 0 < q < p 1 , h C I p , q T , R and g C I p , q T , R + be given functions, φ : C I p , q T , R R be given functional. Then, the problem
D p , q α D p , q β u ( t ) = h ( t ) , t I p , q T ,
u ( 0 ) = u T p
I p , q θ g ( η ) u ( η ) = φ u , η I p , q T 0 , T p
has the unique solution:
u ( t ) = 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 t 0 x p β 1 t q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ h s p α 1 d p , q s d p , q x t β 1 Ω B η P [ h ] + A T φ ( u ) Q [ h ] + t α + β 1 Ω Γ p , q ( α + β ) A η P [ h ] + T p β 1 φ ( u ) Q [ h ]
where the functionals P [ h ] and Q [ h ] are defined by
P [ h ] : = 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ × h s p α 1 d p , q s d p , q x
Q [ h ] : = 1 p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) 0 η 0 y p θ 1 0 x p β 1 ( η q y ) p , q θ 1 ̲ × y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ g y p θ 1 h s p α 1 d p , q s d p , q x d p , q y
and the constants A T , A η , B η and Ω are defined by
A T : = 1 p β 2 Γ p , q ( β ) 0 T p T p q s p , q β 1 ̲ s p θ 1 α 1 ̲ d p , q s = T p α + β 1 Γ p , q ( α + β )
A η : = 1 p θ 2 Γ p , q ( θ ) 0 η ( η q s ) p , q θ 1 ̲ g s p θ 1 s p θ 1 β 1 ̲ d p , q s
B η : = 1 p β 2 + θ 2 Γ p , q ( β ) Γ p , q ( θ ) 0 η 0 x p θ 1 ( η q x ) p , q θ 1 ̲ x p θ 1 q s p , q β 1 ̲ g x p θ 1 × s p α 1 α 1 ̲ d p , q s d p , q x
Ω : = T p β 1 B η A T A η .
Proof. 
Taking fractional (p, q)-integral of order α for (2) and using Lemma 1, we then have
D p , q β u ( t ) = C 1 t α 1 + I p , q α h ( t ) = C 1 t α 1 + 1 p α 2 Γ p , q ( α ) 0 t ( t q s ) p , q α 1 ̲ h s p α 1 d p , q s .
Next, taking fractional (p, q)-difference of order β for (12), we have
u ( t ) = C 0 t β 1 + C 1 t α + β 1 Γ p , q ( α + β ) + 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) × 0 t 0 x p β 1 ( t q x ) p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ h s p α 1 d p , q s d p , q x .
Substituting t = 0 , T p into (13) and employing the condition (3), we have
C 0 T p β 1 + C 1 A T = P [ h ] .
By taking fractional (p, q)-integral of order θ for (13), we have
I p , q θ u ( t ) = C 0 t β + θ 1 Γ p , q ( β + θ ) + C 1 p β 2 + θ 2 Γ p , q ( β ) Γ p , q ( θ ) × 0 t 0 x p θ 1 ( t q x ) p , q θ 1 ̲ x p θ 1 q s p , q β 1 ̲ s p α 1 α 1 ̲ d p , q s d p , q x + 1 p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) 0 t 0 y p θ 1 0 x p β 1 ( t q y ) p , q θ 1 ̲ × y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ h s p α 1 d p , q s d p , q x d p , q y .
From the condition (4) we have
C 0 A η + C 1 B η = φ ( u ) Q [ h ]
Solving the system of linear Equations (14) and (16),we obtain
C 0 = B η P [ h ] A T φ ( u ) Q [ h ] Ω a n d C 1 = T p β 1 φ ( u ) Q [ h ] + A η P [ h ] Ω ,
where P [ h ] , Q [ h ] , A T , A η , B η and Ω are defined by (6)–(11), respectively.
After substituting C 0 , C 1 into (13), we obtain (15). We can prove the converse by direct computation. The proof is complete. □

3. Existence and Uniqueness Result

In this section, we prove the existence and uniqueness result for problem (1) by using Banach fixed point theorem as follows.
Lemma 7
([43] Banach fixed point theorem). Let a nonempty closed subset C of a Banach space X, then there is a unique fixed point for any contraction mapping P of C into itself.
Let C = C I p , q T , R be a Banach space of all function u with the norm defined by
u C = m a x u , D p , q ν u ,
where u = m a x t I p , q T | u ( t ) | and D p , q ν u C = m a x t I p , q T D p , q ν u ( t ) .
By Lemma 6, replacing h ( t ) by F t , u ( t ) , Ψ p , q γ u ( t ) , D p , q ν u ( t ) , we define an operator A : C C by
( A u ) ( t ) : = φ ( u ) Ω ( T p ) β 1 Γ p , q ( α + β ) t α + β 1 A T t β 1 + Q * [ F u ] Ω A T t β 1 ( T p ) β 1 Γ p , q ( α + β ) t α + β 1 + P * [ F u ] Ω A η t α + β 1 Γ p , q ( α + β ) B η t β 1 + 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 t 0 x p β 1 t q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ × F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 , D p , q ν u s p α 1 d p , q s d p , q x
where the functionals P * [ F u ] and Q * [ F u ] are defined by
P * [ F u ] : = 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ × F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 , D p , q ν u s p α 1 d p , q s d p , q x
Q * [ F u ] : = 1 p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) × 0 η 0 y p θ 1 0 x p β 1 η q y p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ g y p θ 1 × F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 , D p , q ν u s p α 1 d p , q s d p , q x d p , q y
and the constants A T , A η , B η and Ω are defined by (8)–(11), respectively.
We see that the problem (1) has solution if and only if the operator A has fixed point.
Theorem 1.
Assume that F : I p , q T × R × R × R R is continuous, ϕ : I p , q T × I p , q T [ 0 , ) is continuous with ϕ 0 = m a x ϕ ( t , s ) : ( t , s ) I p , q T × I p , q T , and φ : C I p , q T , R R is given functional. Suppose that the following conditions hold:
( H 1 )
There exist positive constants L 1 , L 2 , L 3 such that for each t I p , q T and u i , v i R , i = 1 , 2 , 3 ,
| F t , u 1 , u 2 , u 3 F t , v 1 , v 2 , v 3 | L 1 | u 1 v 1 | + L 2 | u 2 v 2 | + L 3 | u 3 v 3 | .
( H 2 )
There exists a positive constant ω such that for each u , v C ,
| φ ( u ) φ ( v ) | ω u v C .
( H 3 )
For each t I p , q T , 0 < g < g ( t ) < G .
( H 4 )
X : = ω O T + L + L 3 Θ 1 ,
where
L : = L 1 + L 2 ϕ 0 ( T p ) γ Γ p , q ( γ + 1 ) ,
Θ : = O T G η α + β + θ Γ p , q ( α + β + θ 1 ) + ( T p ) α + β Γ p , q ( α + β + 1 ) O η + 1 ,
O T : = T p α + β 1 Γ p , q ( α + β ) + A T T p β 1 m i n | Ω | ,
O η : = T p α Γ p , q ( α + β ) m a x   A η + m a x   B η T p β 1 m i n | Ω | .
Then, problem (1) has a unique solution in I p , q T .
Proof. 
For each t I p , q T and u , v C ,
| Ψ p , q γ u ( t ) Ψ p , q γ v ( t ) | ϕ 0 p γ 2 Γ p , q ( γ ) 0 t t q s p , q γ 1 ̲ u s p γ 1 v s p γ 1 d p , q s ϕ 0 p γ 2 Γ p , q ( γ ) u v 0 T p T p q s p , q γ 1 ̲ d p , q s = ϕ 0 T p γ Γ p , q ( γ + 1 ) u v .
Denote that
F | u v | ( t ) : = | F t , u ( t ) , Ψ p , q γ u ( t ) , D p , q ν u ( t ) F t , v ( t ) , Ψ p , q γ v ( t ) , D p , q ν v ( t ) | .
By using Lemma 5(a), we obtain
| P * [ F u ] P * [ F v ] | 1 p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ × F | u v | s p α 1 d p , q s d p , q x L 1 | u v | + L 2 | Ψ p , q γ u Ψ p , q γ v | + L 3 | D p , q ν u D p , q ν v | p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ × x p β 1 q s p , q α 1 ̲ s p α 1 d p , q s d p , q x L 1 + L 2 ϕ 0 T p γ Γ p , q ( γ + 1 ) | u v | + L 3 | D p , q ν u D p , q ν v | T p α + β Γ p , q ( α + β + 1 ) L + L 3 T p α + β Γ p , q ( α + β + 1 ) u v C ,
and by using Lemma 5(b), we have
| Q * [ F u ] Q * [ F v ] | G p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) 0 η 0 y p θ 1 0 x p β 1 η q y p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ × x p β 1 q s p , q α 1 ̲ F | u v | s p α 1 d p , q s d p , q x d p , q y G L 1 | u v | + L 2 | Ψ p , q γ u Ψ p , q γ v | + L 3 | D p , q ν u D p , q ν v | p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) × 0 η 0 y p θ 1 0 x p β 1 η q y p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x d p , q y L 1 + L 2 ϕ 0 T p γ Γ p , q ( γ + 1 ) | u v | + L 3 | D p , q ν u D p , q ν v | G η α + β + θ Γ p , q ( α + β + θ + 1 ) G L + L 3 η α + β + θ Γ p , q ( α + β + θ + 1 ) u v C .
Then,
| ( A u ) ( t ) ( A v ) ( t ) | ω u v C T p β 1 | Ω | T p α + β 1 Γ p , q ( α + β ) + A T + ( L + L 3 ) u v C G η α + β + θ | Ω | Γ p , q ( α + β + θ 1 ) . T p β 1 T p α + β 1 Γ p , q ( α + β ) + A T + ( L + L 3 ) u v C T p α + β | Ω | Γ p , q ( α + β + 1 ) . T p β 1 T p α Γ p , q ( α + β ) A η + B η + ( L + L 3 ) u v C p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x { O T ω + ( L + L 3 ) G η α + β + θ Γ p , q ( α + β + θ 1 ) + O η ( L + L 3 ) T p α + β Γ p , q ( α + β + 1 ) + ( L + L 3 ) T p α + β Γ p , q ( α + β + 1 ) } u v C = X u v C .
Taking fractional (p, q)-difference of order ν for (17), we get
( D p , q ν A u ) ( t ) = φ ( u ) Ω ( T p ) β 1 Γ p , q ( α + β ) Γ p , q ( α + β ) Γ p , q ( α + β ν ) t α + β ν 1 A T Γ p , q ( β ) Γ p , q ( β ν ) t β ν 1 + Q * [ F u ] Ω A T Γ p , q ( β ) Γ p , q ( β ν ) t β ν 1 ( T p ) β 1 Γ p , q ( α + β ) Γ p , q ( α + β ) Γ p , q ( α + β ν ) t α + β ν 1 + P * [ F u ] Ω A η Γ p , q ( α + β ) Γ p , q ( α + β ν ) Γ p , q ( α + β ) t α + β ν 1 B η Γ p , q ( β ) Γ p , q ( β ν ) t β ν 1 + 1 p α 2 + β 2 + ν 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( ν ) × 0 t 0 y p ν 1 0 x p β 1 t q y p , q ν 1 ̲ y p ν 1 q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ × F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 , D p , q ν u s p α 1 d p , q s d p , q x d p , q y .
Thus,
| ( D p , q ν A u ) ( t ) ( D p , q ν A v ) ( t ) | ω u v C T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) O T ( L + L 3 ) u v C G η α + β + θ Γ p , q ( α + β + θ 1 ) . T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) O T + ( L + L 3 ) u v C T p α + β Γ p , q ( α + β + 1 ) . T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) O η + ( L + L 3 ) u v C ( T p ) α + β ν Γ p , q ( α + β ν + 1 ) { O T ω + ( L + L 3 ) G η α + β + θ Γ p , q ( α + β + θ 1 ) T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) + O η ( L + L 3 ) T p α + β Γ p , q ( α + β + 1 ) T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) + ( L + L 3 ) T p α + β ν Γ p , q ( α + β ν + 1 ) } u v C < X u v C .
From (26) and (28), we have
A u A v C X u v C .
By ( H 4 ) , we can conclude that A is a contraction. Thus, by using Banach fixed point theorem in lemma 7, A has a fixed point, which is a unique solution of problem (1) on I p , q T . □

4. Existence of at Least One Solution

In this section, we prove the existence of at least one solution to (1). The following lemmas reviewing the Schauder’s fixed point theorem are also provided.
Lemma 8
([43] Arzelá-Ascoli theorem). A collection of functions in C [ a , b ] with the sup norm, is relatively compact if and only if it is uniformly bounded and equicontinuous on [ a , b ] .
Lemma 9
([43]). If a set is closed and relatively compact, then it is compact.
Lemma 10
([44] Schauder’s fixed point theorem). Let ( D , d ) be a complete metric space, U be a closed convex subset of D, and T : D D be the map such that the set T u : u U is relatively compact in D. Then, the operator T has at least one fixed point u * U : T u * = u * .
Theorem 2.
Assume that F : I p , q T × R × R × R R is continuous, and φ : C I p , q T , R R is given functional. Suppose that the following conditions hold:
( H 5 )
There exists a positive constant M such that for each t I p , q T and u i R , i = 1 , 2 , 3 ,
| F t , u 1 , u 2 , u 3 | M .
( H 6 )
There exists a positive constant N such that for each u C ,
| φ ( u ) | N .
Then, problem ( 1 ) has at least one solution on I p , q T .
Proof. 
To prove this theorem, we proceed as follows.
Step I. Verify A maps bounded sets into bounded sets in B R = { u C : u C R } . Let us prove that for any R > 0 , there exists a positive constant L such that for each x B R , we have A u C L . By using Lemma 5, for each t I p , q T and u B R , we have
| P * [ F u ] | M p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x M T p α + β Γ p , q ( α + β + 1 ) ,
| Q * [ F u ] | G M p α 2 + β 2 + θ 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( θ ) 0 η 0 y p θ 1 0 x p β 1 η q y p , q θ 1 ̲ y p θ 1 q x p , q β 1 ̲ × x p β 1 q s p , q α 1 ̲ d p , q s d p , q x d p , q y G M η α + β + θ Γ p , q ( α + β + θ + 1 ) .
From (29) and (30), we have
| ( A u ) ( t ) | N O T + G M η α + β + θ Γ p , q ( α + β + θ 1 ) O T + M T p α + β Γ p , q ( α + β + 1 ) O η + M p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) 0 T p 0 x p β 1 T p q x p , q β 1 ̲ x p β 1 q s p , q α 1 ̲ d p , q s d p , q x N O T + M O T G η α + β + θ Γ p , q ( α + β + θ 1 ) + T p α + β Γ p , q ( α + β + 1 ) ( O η + 1 ) N O T + M Θ : = L .
We find that
| D p , q ν A u ( t ) | N O T T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) + M [ O T G η α + β + θ Γ p , q ( α + β + θ 1 ) T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) + T p α + β Γ p , q ( α + β + 1 ) ( O η + 1 ) T p ν Γ p , q ( α + β ) Γ p , q ( α + β ν ) ] < L .
Thus, ( A u ) C L , which implies that A is uniformly bounded.
Step II. Since F is continuous, we can conclude that the operator A is continuous on B R .
Step III. For any t 1 , t 2 I p , q T with t 1 < t 2 , we find that
| ( A u ) ( t 1 ) ( A u ) ( t 2 ) | | t 2 β 1 t 1 β 1 | | Ω | A T ( N + Q * [ F u ] ) + B η P * [ F u ] + | t 2 α + β 1 t 1 α + β 1 | | Ω | Γ p , q ( α + β ) T p β 1 ( N + Q * [ F u ] ) + A η P * [ F u ] + M Γ p , q ( α + β ) | t 2 α + β t 1 α + β | ,
and
| ( D p , q ν A u ) ( t 2 ) ( D p , q ν A u ) ( t 1 ) | | t 2 α + β t 1 α + β | | Ω | Γ p , q ( β ν ) A T ( N + Q * [ F u ] ) + B η P * [ F u ] + | t 2 α + β ν 1 t 1 α + β ν 1 | | Ω | Γ p , q ( α + β ν ) × T p β 1 ( N + Q * [ F u ] ) + A η P * [ F u ] + M Γ p , q ( α + β ν + 1 ) | t 2 α + β ν t 1 α + β ν | .
We see that the right-hand side of (34) and (35) tends to be zero when | t 2 t 1 | 0 . Thus, A is relatively compact on B R . This implies that A ( B R ) is an equicontinuous set. By Arzelá-Ascoli theorem in Lemma 8, Lemma 9, and the above steps, we see that A : C C is completely continuous. Hence, we can conclude from Schauder fixed point theorem in Lemma 10 that problem (1) has at least one solution. □

5. Examples

In this section, to illustrate our results, we consider some examples.
Example 1.
Consider the following fractional ( p , q ) ) -integrodifference equation as
D 2 3 , 1 2 3 4 D 2 3 , 1 2 1 2 u ( t ) = 1 100 e 2 + t 3 ( 1 + | u ( t ) | ) [ e 3 t u 2 + 2 | u | + e ( π + s i n 2 π t ) Ψ 2 3 , 1 2 1 3 u ( t ) + e ( 2 π + c o s 2 π t ) D 2 3 , 1 2 1 4 u ( t ) ] , t I 2 3 , 1 2 10 = 10 1 2 k 2 3 k + 1 : k N 0 { 0 }
with periodic fractional ( p , q ) -integral boundary condition
u ( 0 ) = u ( 15 ) I 2 3 , 1 2 2 3 2 e + s i n 1215 256 2 u 1215 256 = i = 0 C i | u ( t i ) | 1 + | u ( t i ) | , t i = σ 2 3 , 1 2 i ( 10 ) ,
where C i is given constants with 1 1000 i = 0 C i e 1000 and ϕ ( t , s ) = e | t s | ( t + e ) 3 .
Letting α = 3 4 , β = 1 2 , γ = 1 3 , ν = 1 4 , θ = 2 3 , p = 2 3 , q = 1 2 , T = 10 , η = σ 2 3 , 1 2 4 ( 10 ) = 1215 256 , g ( t ) = ( 20 e + s i n   t ) 2 and F t , u ( t ) , Ψ p , q γ u ( t ) , D p , q ν u ( t ) = 1 100 e 2 + t 3 ( 1 + | u ( t ) | ) ×
e 3 t u 2 + 2 | u | + e ( π + s i n 2 π t ) Ψ 2 3 , 1 2 1 3 u ( t ) + e ( 2 π + c o s 2 π t ) D 2 3 , 1 2 1 4 u ( t ) .
Using above values, we find that
ϕ 0 = 0.0498 , | A T | = 2.06344 , | A η | 264.588 , | B η | 196.777 and | Ω | 283.525 .
For all t I 2 3 , 1 2 10 and u , v R , we find that
F t , u , Ψ p , q γ u , D p , q ν u F t , v , Ψ p , q γ v , D p , q ν v 1 100 e 2 | u v | + 1 1000 e 2 + π Ψ p , q γ u Ψ p , q γ v + 1 100 e 2 + 2 π D p , q ν u D p , q ν v .
Thus, ( H 1 ) holds with L 1 = 0.001353 , L 2 = 5.848 × 10 5 and L 3 = 2.5273 × 10 6 . So L = 0.00136 .
For all u , v C ,
| φ ( u ) φ ( v ) | e 1000 u v C .
Thus, ( H 2 ) holds with ω = 0.002718 .
In addition, ( H 3 ) holds with g = 19.6831 , G = 41.42935 .
Since
O T = 0.001885 O η = 2.10484 and Θ = 89.5277 ,
therefore, ( H 4 ) holds with
X = 0.121989 < 1 .
Hence, by Theorem 1 this problem has a unique solution.
Example 2.
Consider the following fractional ( p , q ) -integrodifference equation as
D 2 3 , 1 2 3 4 D 2 3 , 1 2 1 2 u ( t ) = 1 10 t + 1 3 e ( t + 1 10 ) u ( t ) + Ψ 2 3 , 1 2 1 3 u ( t ) + D 2 3 , 1 2 1 4 u ( t ) , t I 2 3 , 1 2 10
with periodic fractional ( p , q ) -integral boundary condition
u ( 0 ) = u ( 15 ) I 2 3 , 1 2 2 3 2 e + s i n 1215 256 2 u 1215 256 = i = 0 C i e | u ( t i ) | , t i = σ 2 3 , 1 2 i ( 10 ) ,
where D i is given constants with 1 500 i = 0 D i e 500 .
Letting α = 3 4 , β = 1 2 , γ = 1 3 , ν = 1 4 , θ = 2 3 , p = 2 3 , q = 1 2 , T = 10 , η = 1215 256 . It is clear that F t , u , Ψ p , q γ u , D p , q ν u 23 15 = M for t I 2 3 , 1 2 10 , and φ ( u ) e 500 = N for u C . Thus, we can conclude from Theorem 2 that our problem has at least one solution.

6. Conclusions

A fractional (p, q)-integrodifference equation with periodic fractional (p, q)-integral boundary condition (1) is studied. Our problem contains three fractional (p, q)-difference operators, and two fractional (p, q)-integral operators. We establish the conditions for the existence and uniqueness of solution for problem (1) by using the Banach fixed point theorem, and this result is shown in Theorem 1. We also established the conditions of at least one solution by using the Schauder’s fixed point theorem, and this result is shown in Theorem 2. The choice to use of Theorems 1 or 2 depends on the conditions of the assumptions. The main results are illustrated by a numerical example. Some properties of fractional (p, q)-integral needed in our study are also discussed. The results of the paper are new and enrich the subject of boundary value problems for fractional (p, q)-difference equations. In the future work, we may extend this work by considering new boundary value problems.

Author Contributions

Conceptualization, J.S. and T.S.; methodology, J.S. and T.S.; formal analysis, J.S. and T.S.; funding acquisition, J.S and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok. Contract no.KMUTNB-62-KNOW-22.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their gratitude to anonymous referees for very helpful suggestions and comments which led to improvements of our original manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jackson, F.H. On q-difference equations. Am. J. Math. 1910, 32, 305–314. [Google Scholar] [CrossRef]
  2. Jackson, F.H. On q-difference integrals. Q. J. Pure Appl. Math. 1910, 41, 193–203. [Google Scholar]
  3. Ernst, T. A Comprehensive Treatment of q-Calculus; Springer: Basel, Switzerland, 2012. [Google Scholar]
  4. Kac, V.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002. [Google Scholar]
  5. Ernst, T. A New Notation for q-Calculus and a New q-Taylor Formula; U.U.D.M. Report 1999:25; Department of Mathematics, Uppsala University: Uppsala, Sweden, 1999; ISSN 1101-3591. [Google Scholar]
  6. Finkelstein, R.J. Symmetry group of the hydrogen atom. J. Math. Phys. 1967, 8, 443–449. [Google Scholar] [CrossRef]
  7. Finkelstein, R.J. q-gauge theory. Int. J. Mod. Phys. A 1996, 11, 733–746. [Google Scholar] [CrossRef] [Green Version]
  8. Ilinski, K.N.; Kalinin, G.V.; Stepanenko, A.S. q-functional field theory for particles with exotic statistics. Phys. Lett. A 1997, 232, 399–408. [Google Scholar] [CrossRef] [Green Version]
  9. Finkelstein, R.J. q-field theory. Lett. Math. Phys. 1995, 34, 169–176. [Google Scholar] [CrossRef]
  10. Finkelstein, R.J. The q-Coulomb problem. J. Math. Phys. 1996, 37, 2628–2636. [Google Scholar] [CrossRef]
  11. Cheng, H. Canonical Quantization of Yang-Mills Theories. Perspectives in Mathematical Physics; International Press: Boston, MA, USA, 1996. [Google Scholar]
  12. Finkelstein, R.J. q-gravity. Lett. Math. Phys. 1996, 38, 53–62. [Google Scholar] [CrossRef]
  13. Negadi, T.; Kibler, M. A q-deformed Aufbau Prinzip. J. Phys. A Math. Gen. 1992, 25, 157–160. [Google Scholar] [CrossRef]
  14. Siegel, W. Introduction to String Field Theory, Volume 8 of Advanced Series in Mathematical Physics; World Scientific Publishing: Teaneck, NJ, USA, 1988. [Google Scholar]
  15. Agarwal, R.P. Certain fractional q-integrals and q-derivatives. Proc. Camb. Phil. Soc. 1969, 66, 365–370. [Google Scholar] [CrossRef]
  16. Al-Salam, W.A. Some fractional q-integrals and q-derivatives. Proc. Edin. Math. Soc. 1966, 15, 135–140. [Google Scholar] [CrossRef] [Green Version]
  17. Díaz, J.B.; Osler, T.J. Differences of fractional order. Math. Comp. 1974, 28, 185–202. [Google Scholar] [CrossRef]
  18. Brikshavana, T.; Sitthiwirattham, T. On Fractional Hahn Calculus. Adv. Differ. Equ. 2017, 354, 1–15. [Google Scholar] [CrossRef] [Green Version]
  19. Patanarapeelert, N.; Sitthiwirattham, T. On Fractional Symmetric Hahn Calculus. Mathematics 2019, 7, 873. [Google Scholar] [CrossRef] [Green Version]
  20. Chakrabarti, R.; Jagannathan, R. A (p, q)-oscillator realization of two-parameter quantum algebras. J. Math. Phys. Math. Gen. 1991, 24, 5683–5701. [Google Scholar] [CrossRef]
  21. Sadjang, P.N. On the fundamental theorem of (p, q)-calculus and some (p, q)-Taylor formulas. Results Math. 2018, 73, 1–21. [Google Scholar]
  22. Mursaleen, M.; Ansari, K.J.; Khan, A. On (p, q)-analogue of Bernstein operators. Appl. Math. Comput. 2015, 266, 874–882. [Google Scholar] [CrossRef]
  23. Mursaleen, M.; Ansari, K.J.; Khan, A. Some approximation results by (p, q)-analogue of Bernstein-Stancu operators. Appl. Math. Comput. 2015, 264, 392–402, Corrigendum in Appl. Math. Comput 2015, 269, 744 –746. [Google Scholar] [CrossRef]
  24. Mursaleen, M.; Khan, F.; Khan, A. Approximation by (p, q)-Lorentz polynomials on a compact disk. Complex Anal. Oper. Theory 2016, 10, 1725–1740. [Google Scholar] [CrossRef] [Green Version]
  25. Rahman, S.; Mursaleen, M.; Alkhaldi, A.H. Convergence of iterates of q-Bernstein and (p, q)-Bernstein operators and the Kelisky-Rivlin type theorem. Filomat 2018, 32, 4351–4364. [Google Scholar] [CrossRef]
  26. Jebreen, H.B.; Mursaleen, M.; Ahasan, M. On the convergence of Lupas (p, q)-Bernstein operators via contraction principle. J. Ineq. Appl. 2019, 2019. [Google Scholar] [CrossRef] [Green Version]
  27. Khan, A.; Sharma, V. Statistical approximation by (p, q)-analogue of Bernstein-Stancu operators. Azerbaijan J. Math. 2018, 8, 100–121. [Google Scholar]
  28. Khan, K.; Lobiyal, D.K. Bezier curves based on Lupas (p, q)-analogue of Bernstein functions in CAGD. J. Comput. Appl. Math. 2017, 317, 458–477. [Google Scholar] [CrossRef]
  29. Milovanovic, G.V.; Gupta, V.; Malik, N. (p, q)-Beta functions and applications in approximation. Bol. Soc. Mat. Mex. 2018, 24, 219–237. [Google Scholar] [CrossRef] [Green Version]
  30. Cheng, W.T.; Zhang, W.H.; Cai, Q.B. (p, q)-gamma operators which preserve x2. J. Inequal. Appl. 2019, 2019, 1–14. [Google Scholar] [CrossRef] [Green Version]
  31. Nasiruzzaman, M.; Mukheimer, A.; Mursaleen, M. Some Opial-type integral inequalities via (p, q)-calculus. J. Ineq. Appl. 2019, 2019, 1–15. [Google Scholar] [CrossRef] [Green Version]
  32. Jain, P.; Basu, C.; Panwar, V. On the (p, q)-Melin Transform and Its Applications. Acta Math. Sci. Ser. B 2021, 41, 1719–1732. [Google Scholar] [CrossRef]
  33. Aral, A.; Deniz, E.; Erbay, H. The Picard and Gauss-Weierstrass singular integral in (p, q)-calculus. Bull. Malays. Math. Sci. Soc. 2020, 43, 1569–1583. [Google Scholar] [CrossRef]
  34. Kamsrisuk, N.; Promsakon, C.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for (p, q)-difference equations. Differ. Equ. Appl. 2018, 10, 183–195. [Google Scholar] [CrossRef] [Green Version]
  35. Promsakon, C.; Kamsrisuk, N.; Ntouyas, S.K.; Tariboon, J. On the second-order quantum (p, q)-difference equations with separated boundary conditions. Adv. Math. Phys. 2018, 2018, 9089865. [Google Scholar] [CrossRef] [Green Version]
  36. Nuntigrangjana, T.; Putjuso, S.N.; Ntouyas, S.K.; Tariboon, J. Impulsive quantum (p, q)-difference equations. Adv. Differ. Equ. 2020, 2020, 1–20. [Google Scholar] [CrossRef] [Green Version]
  37. Soontharanon, J.; Sitthiwirattham, T. On Fractional (p, q)-Calculus. Adv. Differ. Equ. 2020, 2020, 1–18. [Google Scholar] [CrossRef] [Green Version]
  38. Pheak, N.; Nonlaopon, K.; Tariboon, J.; Ntouyas, S.K. Fractional (p, q)-Calculus on Finite Intervals and Some Integral Inequalities. Symmetry 2021, 13, 504. [Google Scholar]
  39. Pheak, N.; Nonlaopon, K.; Tariboon, J.; Ntouyas, S.K. Praveen Agarwal, Some trapezoid and midpoint type inequalities via fractional (p, q)-calculus. Adv. Differ. Equ. 2021, 2021, 1–22. [Google Scholar]
  40. Soontharanon, J.; Sitthiwirattham, T. Existence Results of Nonlocal Robin Boundary Value Problems for Fractional (p, q)-Integrodifference Equations. Adv. Differ. Equ. 2020, 2020, 1–17. [Google Scholar] [CrossRef]
  41. Qin, Z.; Sun, S. On a nonlinear fractional (p, q)-difference Schrödinger equation. J. Appl. Math. Comput. 2021. [Google Scholar] [CrossRef]
  42. Qin, Z.; Sun, S. Positive solutions for fractional (p, q)-difference boundary value problems. J. Appl. Math. Comput. 2021. [Google Scholar] [CrossRef]
  43. Griffel, D.H. Applied Functional Analysis; Ellis Horwood Publishers: Chichester, UK, 1981. [Google Scholar]
  44. Guo, D.; Lakshmikantham, V. Nonlinear Problems in Abstract Cone; Academic Press: Orlando, FL, USA, 1988. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Soontharanon, J.; Sitthiwirattham, T. On Periodic Fractional (p, q)-Integral Boundary Value Problems for Sequential Fractional (p, q)-Integrodifference Equations. Axioms 2021, 10, 264. https://doi.org/10.3390/axioms10040264

AMA Style

Soontharanon J, Sitthiwirattham T. On Periodic Fractional (p, q)-Integral Boundary Value Problems for Sequential Fractional (p, q)-Integrodifference Equations. Axioms. 2021; 10(4):264. https://doi.org/10.3390/axioms10040264

Chicago/Turabian Style

Soontharanon, Jarunee, and Thanin Sitthiwirattham. 2021. "On Periodic Fractional (p, q)-Integral Boundary Value Problems for Sequential Fractional (p, q)-Integrodifference Equations" Axioms 10, no. 4: 264. https://doi.org/10.3390/axioms10040264

APA Style

Soontharanon, J., & Sitthiwirattham, T. (2021). On Periodic Fractional (p, q)-Integral Boundary Value Problems for Sequential Fractional (p, q)-Integrodifference Equations. Axioms, 10(4), 264. https://doi.org/10.3390/axioms10040264

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