Next Article in Journal
An Efficient and Robust Improved A* Algorithm for Path Planning
Previous Article in Journal
An Improvised SIMPLS Estimator Based on MRCD-PCA Weighting Function and Its Application to Real Data
Previous Article in Special Issue
Stability Analysis and Existence of Solutions for a Modified SIRD Model of COVID-19 with Fractional Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Separate Fractional (p,q)-Integrodifference Equations via Nonlocal Fractional (p,q)-Integral Boundary Conditions

by
Thongchai Dumrongpokaphan
1,
Sotiris K. Ntouyas
2,3 and
Thanin Sitthiwirattham
4,*
1
Research Group in Mathematics and Applied Mathematics, Department of Mathematics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
2
Department of Mathematics, University of Ioannina, 45110 Ioannina, Greece
3
Nonlinear Analysis and Applied Mathematics (NAAM)-Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
*
Author to whom correspondence should be addressed.
Submission received: 26 October 2021 / Revised: 10 November 2021 / Accepted: 17 November 2021 / Published: 19 November 2021

Abstract

:
In this paper, we study a boundary value problem involving ( p , q ) -integrodifference equations, supplemented with nonlocal fractional ( p , q ) -integral boundary conditions with respect to asymmetric operators. First, we convert the given nonlinear problem into a fixed-point problem, by considering a linear variant of the problem at hand. Once the fixed-point operator is available, existence and uniqueness results are established using the classical Banach’s and Schaefer’s fixed-point theorems. The application of the main results is demonstrated by presenting numerical examples. Moreover, we study some properties of ( p , q ) -integral that are used in our study.

1. Introduction

Quantum calculus or q-calculus is the modern name for the study of calculus without limits. q-calculus was first introduced by Jackson [1,2] in 1910. Quantum calculus has many applications in mathematics and physics, for example, in orthogonal polynomials, combinatorics, number theory, simple hypergeometric functions, dynamics and theory of relativity, to name a few.
An extension of quantum calculus, the ( p , q ) -calculus or post-quantum calculus, was introduced by Chakrabarti and Jagannathan in [3]. ( p , q ) -calculus is a generalization of q calculus including two independent quantum parameters p and q, reduced to q-calculus for the case p = 1 and to the classical q calculus when q 1 . Furthermore, ( p , q ) -calculus has many applications, such as physical sciences, combinatorics, hypergeometric functions, number theory, mechanics, Bézier curves and surfaces, etc. (for instances, see [4,5,6,7,8]). Many researchers have recently begun working on ( p , q ) -calculus and some results can be found in [9,10,11,12,13,14] and references cited therein. Recently, in [15], the authors introduced the fractional ( p , q ) -difference operators and studied their properties.
Recently, ( p , q ) -calculus was applied to establish several new types of inequalities (see [16,17] and references cited therein). In the literature, there exist few papers studying boundary value problems for ( p , q ) -difference equations, because the ( p , q ) -fractional operator has been introduced recently. In [18], the following ( p , q ) boundary value problem for second order ( p , q ) -difference equations with separated boundary conditions was studied:
D p , q 2 x ( t ) = f ( t , x ( p 2 t ) ) , t [ 0 , T / p 2 ] , α 1 x ( 0 ) + α 2 D p , q x ( 0 ) = α 3 , β 1 x ( T ) + β 2 D p , q x ( T / p ) = β 3 ,
where 0 < q < p 1 are two quantum numbers, D p , q 2 is the second order ( p , q ) -difference operator and f C ( [ 0 , T / p 2 ] × R , R ) , T > 0 , α i , β i , i = 1 , 2 , 3 are given real constants. A variety of new existence and uniqueness results were established using Banach’s, Schaefer’s and Krasnoselskii’s fixed-point theorems, as well as Leray–Schauder’s nonlinear alternative.
Some more results on ( p , q ) boundary value problems can be found in [19,20,21]. For existence results for boundary value problems fractional ( p , q ) -difference Schrödinger equations, we refer to [22].
Recently, in [15], the authors introduced the fractional ( p , q ) -integrodifference operators and studied their properties. Boundary value problems for fractional ( p , q ) -integrodifference equations with Robin boundary conditions were studied in [23], where the authors established existence and uniqueness results for the following problem:
D p , q α u ( t ) = F t , u ( t ) , Ψ p , q γ u ( t ) , D p , q ν u ( t ) , t I p , q T , λ 1 u ( η ) + λ 2 D p , q β u ( η ) = ϕ 1 ( u ) , η I p , q T 0 , T p , μ 1 u T p + μ 2 D p , q β u T p = ϕ 2 ( u ) ,
where I p , q T : = q k p k + 1 T : k N 0 { 0 } , 0 < q < p 1 , α ( 1 , 2 ] ; β , γ , ν ( 0 , 1 ] , λ 1 , λ 2 , μ 1 , μ 2 R + , F C I p , q T × R × R × R , R is given function, ϕ 1 , ϕ 2 : C I p , q T , R R are given functionals and
Ψ p , q γ u ( t ) : = I p , q γ φ u ( t ) = 1 p γ 2 Γ p , q ( γ ) 0 t ( t q s ) p , q γ 1 ̲ φ ( t , s ) u s p γ 1 d p , q s ,
are ( p , q ) -integral operators defined for φ C I p , q T × I p , q T , [ 0 , ) .
Motivated by the the aforementioned papers, our goal in this paper is to enrich the literature on boundary value problems for fractional ( p , q ) -integrodifference equations. More precisely, we introduce and study a nonlocal boundary value problem for ( p , q ) -integrodifference equations subject to fractional ( p , q ) -integral boundary conditions of the form
D p , q α u ( t ) = λ F t , u ( t ) , Ψ p , q γ u ( t ) + μ H t , u ( t ) , Υ p , q ν u ( t ) , t I p , q T , I p , q β g 1 ( η ) u ( η ) = ϕ 1 ( u ) , η I p , q T 0 , T p , I p , q β g 2 T p u T p = ϕ 2 ( u ) ,
where I p , q T : = q p k T p : k N 0 { 0 } ; 0 < q < p 1 ; α ( 1 , 2 ] , β , γ , ν ( 0 , 1 ] , λ , μ R + , F , H C I p , q T × R × R , R and g 1 , g 2 C I p , q T , R + are given functions, ϕ 1 , ϕ 2 : C I p , q T , R R are given functionals and
Ψ p , q γ u ( t ) : = I p , q γ φ u ( t ) = 0 t ( t q s ) p , q γ 1 ̲ p γ 2 Γ p , q ( γ ) φ t , s p γ 1 u s p γ 1 d p , q s ,
Υ p , q ν u ( t ) : = D p , q ν ψ u ( t ) = 0 t ( t q s ) p , q ν 1 ̲ p ν 2 Γ p , q ( ν ) ψ t , s p ν 1 u s p ν 1 d p , q s .
are operators defined for φ , ψ C I p , q T × I p , q T , [ 0 , ) .
We notice that the the boundary value problem (4) is of general type, concerning both ( p , q ) -fractional integral and ( p , q ) -fractional derivative operators. In addition, it contains nonlocal boundary conditions; it is well known that the study of nonlocal boundary value problems is of significance, since they have applications in physics and other areas of applied mathematics. We emphasize that the novelty of our paper lies in both the equation and the boundary conditions, contributing significantly to the existing literature on the topic. Our existence and uniqueness results rely on the standard tools of functional analysis. The methods used in our analysis are standard; however, their exposition in the framework of the boundary value problem (4) is new.
The remaining part of this manuscript is organized as follows. Section 2 contains some basic notions and known results of ( p , q ) -calculus. Furthermore, an auxiliary result is proved which plays a key role in transforming the given problem into a fixed-point problem. In Section 3, we prove the existence of a unique solution for the boundary value problem (4) via the Banach contraction mapping principle. An existence result is proved in Section 4, by using Schaefer’s fixed-point theorem. Finally, examples illustrating the applicability of the main results are presented in Section 5. The papers ends with a section that illustrates the conclusions.

2. Preliminaries

In this section, some concepts regarding our study are recalled. Let 0 < q < p 1 . The following notations are used:
[ 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 .
By
σ p , q k ( t ) : = q p k t and ρ p , q k ( t ) : = p q k t , for k N ,
we denote, respectively, the ( p , q ) -forward jump and the ( p , q ) -backward jump operators. For the power function ( a b ) q n ̲ with n N 0 : = { 0 , 1 , 2 , } , the q-analogue is given by
( a b ) q 0 ̲ : = 1 , ( a b ) q n ̲ : = i = 0 n 1 ( a b q i ) , a , b R .
and
( a b ) p , q 0 ̲ : = 1 , ( a b ) p , q n ̲ : = k = 0 n 1 ( a p k b q k ) , a , b R ,
is the ( p , q ) -analogue of the power function ( a b ) p , q n ̲ , with n N 0 .
For α R , we define a general form
( 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 .
Note that 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 ) , f o r t 0 f ( 0 ) , f o r t = 0
if f is differentiable at 0. If D p , q f ( t ) exists for all t I p , q T , then f is called ( p , q ) -differentiable on I p , q T .
Definition 2.
Let us assume that I is a closed interval of R containing a , b and 0 and f : I R is a given function. The ( p , q ) -integral of the function f from a to b is defined by
a b f ( s ) d p , q s : = 0 b f ( s ) d p , q s 0 a f ( s ) d p , q s ,
where
I p , q f ( t ) = 0 t f ( s ) d p , q s = ( p q ) t k = 0 q k p k + 1 f q k p k + 1 t , t I ,
provided that the series converges at t = a and t = b . If f is ( p , q ) -integrable on [ a , b ] for all a , b I , then f is called ( p , q ) -integrable on [ a , b ] .
Next, we define an operator I p , q N as
I p , q 0 f ( x ) = f ( x ) and 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 ) and I p , q D p , q f ( x ) = f ( x ) f ( 0 ) .
Next, we introduce the Riemann–Liouville type of fractional ( p , q ) -integral and fractional ( p , q ) -difference.
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 ) , where the notation α 2 is a combination.
Definition 4.
For α > 0 , 0 < q < p 1 and f defined on I p , q T . The Riemann–Liouville type fractional ( p , q ) -difference operator of 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 .
Next, we introduce lemmas that are used in the main results.
Lemma 1
([15]). 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
([15]). Let 0 < q < p 1 and f : I p , q T R a continuous at 0 function. Then, we have
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
([15]). 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
([23]). 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 ( 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 .
We employ the above lemmas to obtain the new results as follows.
Lemma 5.
Let α , β > 0 , 0 < q < p 1 and n Z . Then,
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 ( β ) Γ p , q ( α + β + 1 ) t α + β .
Proof. 
From Lemma 3 ( a ) and the definition of the ( p , q ) -beta function, we have
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 = B p , q ( 1 , α ) p α ( β 1 ) 0 t ( t q x ) p , q β 1 ̲ x α d p , q x = B p , q ( 1 , α ) p α ( β 1 ) B p , q ( α + 1 , β ) t α + β = p α 2 + β 2 Γ p , q ( α ) Γ p , q ( β ) Γ p , q ( α + β + 1 ) t α + β .
The proof is complete. □
The following lemma, concerning a linear variant of problem (4), plays a significant role in the forthcoming analysis.
Lemma 6.
Let Λ 0 , α ( 1 , 2 ] , β ( 0 , 1 ] , 0 < q < p 1 , h C I p , q T , R and g 1 , g 2 C I p , q T , R + be given functions and ϕ 1 , ϕ 2 : C I p , q T , R R be given functionals. Then, the boundary value problem
D p , q α u ( t ) = h ( t ) , t I p , q T ,
I p , q β g 1 ( η ) u ( η ) = ϕ 1 ( u ) , η I p , q T 0 , T p ,
I p , q β g 2 T p u T p = ϕ 2 ( u ) ,
has the unique solution
u ( t ) = 1 p α 2 Γ p , q ( α ) 0 t t q s p , q α 1 ̲ h s p α 1 d p , q s t α 1 Λ B T O η [ ϕ 1 , h ] B η O T [ ϕ 2 , h ] + t α 2 Λ A T O η [ ϕ 1 , h ] A η O T [ ϕ 2 , h ]
where the functionals O η [ ϕ 1 , h ] , O T [ ϕ 2 , h ] are defined by
O η [ ϕ 1 , h ] : = ϕ 1 ( u ) 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 1 x p β 1 h s p α 1 d p , q s d p , q x ,
O T [ ϕ 2 , h ] : = ϕ 2 ( 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 ̲ g 2 x p β 1 h s p α 1 d p , q s d p , q x ,
and the constants A η , A T , B η , B T and Λ are defined by
A η : = 1 p β 2 Γ p , q ( β ) 0 η η q s p , q β 1 ̲ g 1 s p β 1 s p β 1 α 1 d p , q s
A T : = 1 p β 2 Γ p , q ( β ) 0 T p T p q s p , q β 1 ̲ g 2 s p β 1 s p β 1 α 1 d p , q s
B η : = 1 p β 2 Γ p , q ( β ) 0 η η q s p , q β 1 ̲ g 1 s p β 1 s p β 1 α 2 d p , q s
B T : = 1 p β 2 Γ p , q ( β ) 0 T p T p q s p , q β 1 ̲ g 2 s p β 1 s p β 1 α 2 d p , q s
Λ : = A T B η A η B T .
Proof. 
Taking fractional ( p , q ) -integral of order α for (7) and using Lemma 1, we have
u ( t ) = C 1 t α 1 + C 2 t α 2 + 1 p α 2 Γ p , q ( α ) 0 t ( t q s ) p , q α 1 ̲ h s p α 1 d p , q s .
Then, we take fractional ( p , q ) -integral of order β for (18); we have
I p , q β u ( t ) = C 1 t α + β 1 Γ p , q ( α + β ) + C 2 t α + β 2 Γ p , q ( α + β 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 ̲ h s p α 1 d p , q s d p , q x .
Substituting t = η into ( 19 ) and employing the condition ( 8 ) , we have
C 1 A η + C 2 B η = O η [ ϕ 1 , h ] .
Taking t = T p into ( 19 ) and employing the condition ( 9 ) , we have
C 1 A T + C 2 B T = O T [ ϕ 2 , h ] .
Solving (20) and (21), we find that
C 1 = B η O T B T O η Λ and C 2 = A T O η A η O T Λ ,
where O η [ ϕ 1 , h ] , O T [ ϕ 2 , h ] , A η , A T , B η , B T and Λ are defined by (11)–(17), respectively. Substituting the constants C 1 , C 2 into ( 18 ) , we obtain ( 10 ) . We can prove the converse by direct computation. This completes the proof. □

3. Existence and Uniqueness Result

In this section, an existence and uniqueness result for the problem ( 4 ) is proved, via the Banach contraction mapping principle. By C = C I p , q T , R , we denote the Banach space furnished with the norm
u C = u + D p , q ν u ,
where u = max t I p , q T | u ( t ) | and D p , q ν u = max 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 ) + μ H t , u ( t ) , Υ p , q ν u ( t ) , we define an operator A : C C by
( A u ) ( t ) : = 1 p α 2 Γ p , q ( α ) 0 t t q s p , q α 1 ̲ ( λ F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 + μ H s p α 1 , u s p α 1 , Υ p , q ν u s p α 1 ) d p , q s t α 1 Λ B T O η * [ ϕ 1 , F u + H u ] B η O T * [ ϕ 2 , F u + H u ] + t α 2 Λ A T O η * [ ϕ 1 , F u + H u ] A η O T * [ ϕ 2 , F u + H u ]
where the functionals O η * [ ϕ 1 , F u + H u ] , O T * [ ϕ 2 , F u + H u ] are defined by
O η * [ ϕ 1 , F u + H u ] : = ϕ 1 ( u ) 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 1 x p β 1 ( λ F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 + μ H s p α 1 , u s p α 1 , Υ p , q ν u s p α 1 ) d p , q s d p , q x ,
O T * [ ϕ 2 , F u + H u ] : = ϕ 2 ( 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 ̲ × g 2 x p β 1 ( λ F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 + μ H s p α 1 , u s p α 1 , Υ p , q ν u s p α 1 ) d p , q s d p , q x ,
and the constants A η , A T , B η , B T , Λ are defined by (13)–(17), respectively.
Notice that a fixed point of the operator A is a solution of the problem (4).
Theorem 1.
Let us assume that F , H : I p , q T × R × R R are continuous, φ , ψ : I p , q T × I p , q T [ 0 , ) are continuous with φ 0 = max φ ( t , s ) : ( t , s ) I p , q T × I p , q T and ψ 0 = max ψ ( t , s ) : ( t , s ) I p , q T × I p , q T . In addition, we suppose that:
( H 1 )
There exist constants M i > 0 such that, for each t I p , q T and u i , v i R , i = 1 , 2
| F t , u 1 , u 2 F t , v 1 , v 2 | M 1 | u 1 v 1 | + M 2 | u 2 v 2 | .
( H 2 )
There exist constants N i > 0 such that, for each t I p , q T and u i , v i R , i = 1 , 2
| H t , u 1 , u 2 H t , v 1 , v 2 | N 1 | u 1 v 1 | + N 2 | u 2 v 2 | .
( H 3 )
There exist constants ω 1 , ω 2 > 0 such that, for each u , v C
| ϕ 1 ( u ) ϕ 1 ( v ) | ω 1 u v C and | ϕ 2 ( u ) ϕ 2 ( v ) | ω 2 u v C .
( H 4 )
For each t I p , q T , g 1 g 1 ( t ) G 1 and g 2 g 2 ( t ) G 2 .
( H 5 )
X : = L Φ + G 1 η α + β Θ T * + G 2 T p α + β Θ η * Γ p , q ( α + β + 1 ) + ω 1 Θ T * + ω 2 Θ η * < 1 ,
where
L : = λ M 1 + M 2 φ 0 T p γ Γ p , q ( γ + 1 ) + μ N 1 + N 2 ψ 0 T p ν Γ p , q ( 1 ν ) ,
Φ : = T p α Γ p , q ( α + 1 ) + T p α + γ Γ p , q ( α + γ + 1 ) + T p α ν Γ p , q ( α ν + 1 ) ,
Θ T * : = Θ T + Θ ¯ T ,
Θ η * : = Θ η + Θ ¯ η ,
Θ T : = 1 min | Λ | max B T T p α 1 + max A T T p α 2 ,
Θ ¯ T : = 1 min | Λ | max B T T p α ν 1 Γ p , q ( α ν ) + max A T T p α ν 2 Γ p , q ( α ν 1 ) ,
Θ η : = 1 min | Λ | max B η T p α 1 + max A η T p α 2 ,
Θ ¯ η : = 1 min | Λ | max B η T p α ν 1 Γ p , q ( α ν ) + max A η T p α ν 2 Γ p , q ( α ν 1 ) .
Then, the boundary value problem ( 4 ) has a unique solution on I p , q T .
Proof. 
For each t I p , q T and u , v C , we have
| ( Ψ 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 .
Similarly, we have | ( Υ p , q ν u ) ( t ) ( Υ p , q ν v ) ( t ) | ψ 0 T p ν Γ p , q ( 1 ν ) u v .
We set
F | u v | ( t ) : = | F t , u ( t ) , ( Ψ p , q γ u ) ( t ) F t , v ( t ) , ( Ψ p , q γ v ) ( t ) | , H | u v | ( t ) : = | H t , u ( t ) , ( Υ p , q ν u ) ( t ) H t , v ( t ) , ( Υ p , q ν v ) ( t ) | .
Then, we obtain
| O η * [ ϕ 1 , F u + H u ] O η * [ ϕ 1 , F v + H v ] | | ϕ 1 ( u ) ϕ 1 ( v ) | + 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 1 x p β 1 λ F | u v | s p α 1 + μ H | u v | s p α 1 d p , q s d p , q x ω 1 u v C + λ M 1 + M 2 φ 0 T p γ Γ p , q ( γ + 1 ) + μ N 1 + N 2 ψ 0 T p ν Γ p , q ( 1 ν ) × G 1 η α + β Γ p , q ( α + β + 1 ) u v C ω 1 + L G 1 η α + β Γ p , q ( α + β + 1 ) u v C .
Similarly, we obtain
| O T * [ ϕ 1 , F u + H u ] O T * [ ϕ 1 , F v + H v ] | ω 2 + L G 2 T p α + β Γ p , q ( α + β + 1 ) u v C ,
and
| ( A u ) ( t ) ( A v ) ( t ) | 1 p α 2 Γ p , q ( α ) 0 T p T p q s p , q α 1 ̲ λ F | u v | s p α 1 + μ H | u v | s p α 1 d p , q s + T p α 1 | Λ | { | B T | | O η * [ ϕ 1 , F u + H u ] O η * [ ϕ 1 , F v + H v ] | + | B η | | O T * [ ϕ 2 , F u + H u ] O T * [ ϕ 2 , F v + H v ] | } + T p α 2 | Λ | { | A T | | O η * [ ϕ 1 , F u + H u ] O η * [ ϕ 1 , F v + H v ] | + | A η | | O T * [ ϕ 2 , F u + H u ] O T * [ ϕ 2 , F v + H v ] | } L T p α Γ p , q ( α + 1 ) + G 1 η α + β Θ T + G 2 T p α + β Θ η Γ p , q ( α + β + 1 ) + ω 1 Θ T + ω 2 Θ η u v C .
Next, we consider ( D p , q ν A u ) . We have
( D p , q ν A 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 ̲ × ( λ F s p α 1 , u s p α 1 , Ψ p , q γ u s p α 1 + μ H [ s p α 1 , u s p α 1 , Υ p , q ν u s p α 1 ] ) d p , q s d p , q x B T O η * [ ϕ 1 , F u + H u ] B η O T * [ ϕ 2 , F u + H u ] Λ p ν 2 Γ p , q ( ν ) × 0 t t q s p , q ν 1 ̲ s p ν 1 α 1 d p , q s + A T O η * [ ϕ 1 , F u + H u ] A η O T * [ ϕ 2 , F u + H u ] Λ p ν 2 Γ p , q ( ν ) × 0 t t q s p , q ν 1 ̲ s p ν 1 α 2 d p , q s .
Similarly as above, we have
| ( D p , q ν F u ) ( t ) ( D p , q ν F v ) ( t ) | L T p α ν Γ p , q ( α ν + 1 ) + G 1 η α + β Θ ¯ T + G 2 T p α + β Θ ¯ η Γ p , q ( α + β + 1 ) + ω 1 Θ ¯ T + ω 2 Θ ¯ η u v C .
From (33) and (35), we obtain
F u F v C X u v C .
Thus, by ( H 5 ) the operator A is a contraction. By the Banach contraction mapping principle we deduce that A has a fixed point which is the unique solution of problem (4) on I p , q T . The proof is finished. □

4. An Existence Result

In this section, we present an existence result for the boundary value problem (4) by using the Schaefer’s fixed-point theorem [24].
Theorem 2.
Let us assume that F , H : I p , q T × R × R R are continuous functions and φ 1 , φ 2 : C I p , q T , R R are given functionals. Let us suppose that the following conditions hold:
( H 6 )
There exist positive constants F , H such that, for each t I p , q T and u i R , i = 1 , 2 ,
| F t , u 1 , u 2 | F a n d | H t , u 1 , u 2 | H .
( H 7 )
There exist positive constants O 1 , O 2 such, that for each u C ,
| φ 1 ( u ) | O 1 a n d | φ 2 ( u ) | O 2 .
Then, the boundary problem ( 4 ) has at least one solution on I p , q T .
Proof. 
We need to show that the operator A is compact by applying the well-known Arzelá–Ascoli theorem. So, we show that the operator A ( B R ) is a uniformly bounded set, where B R = { u C : u C R , R > 0 } and an equicontinuous set.
(i) For each t I p , q T and u B R , we have
| O η * [ ϕ 1 , F u + H u ] | O 1 + [ λ F + μ H ] 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 1 x p β 1 d p , q s d p , q x O 1 + G 1 η α + β [ λ F + μ H ] Γ p , q ( α + β + 1 ) .
Similarly, we have
| O T * [ ϕ 2 , F u + H u ] | O 2 + G 2 T p α + β [ λ F + μ H ] Γ p , q ( α + β + 1 ) .
From (36) and (37) and for each t I p , q T , we find that
| ( A u ) ( t ) | [ λ F + μ H ] p α 2 Γ p , q ( α ) 0 T p T p q s p , q α 1 ̲ d p , q s + T p α 1 | Λ | | B T | | O η * [ ϕ 1 , F u + H u ] | + | B η | | O T * [ ϕ 2 , F u + H u ] | + T p α 2 | Λ | | A T | | O η * [ ϕ 1 , F u + H u ] | + | A η | | O T * [ ϕ 2 , F u + H u ] | [ λ F + μ H ] T p α Γ p , q ( α + 1 ) + G 1 η α + β Θ T + G 2 T p α + β Θ η Γ p , q ( α + β + 1 ) + O 1 Θ T + O 2 Θ η .
In addition, we obtain
| D p , q ν A u ( t ) | [ λ F + μ H ] T p α ν Γ p , q ( α ν + 1 ) + G 1 η α + β Θ ¯ T + G 2 T p α + β Θ ¯ η Γ p , q ( α + β + 1 ) + O 1 Θ ¯ T + O 2 Θ ¯ η .
Form (38) and (39) we obtain
A u C [ λ F + μ H ] Φ + G 1 η α + β Θ T * + G 2 T p α + β Θ η * Γ p , q ( α + β + 1 ) + O 1 Θ T * + O 2 Θ η * < ,
which implies that A ( B R ) is uniformly bounded.
(ii) We show that A ( B R ) is equicontinuous. For any t 1 , t 2 I p , q T with t 1 < t 2 , we have
| ( A u ) ( t 2 ) ( A u ) ( t 1 ) | [ λ F + μ H ] Γ p , q ( α + 1 ) | t 2 α t 1 α | + | t 2 α 1 t 1 α 1 | | Λ | | B T | O η * [ ϕ 1 , F u + H u ] + | B η | O T * [ ϕ 2 , F u + H u ] + | t 2 α 2 t 1 α 2 | | Λ | | A T | O η * [ ϕ 1 , F u + H u ] + | A η | O T * [ ϕ 2 , F u + H u ] ,
and
| ( D p , q ν A u ) ( t 2 ) ( D p , q ν A u ) ( t 1 ) | [ λ F + μ H ] Γ p , q ( α ν + 1 ) | t 2 α ν t 1 α ν | + Γ p , q ( α ) | t 2 α ν 1 t 1 α ν 1 | | Λ | Γ p , q ( α ν ) | B T | O η * [ ϕ 1 , F u + H u ] + | B η | O T * [ ϕ 2 , F u + H u ] + Γ p , q ( α 1 ) | t 2 α ν 2 t 1 α ν 2 | | Λ | Γ p , q ( α ν 1 ) | A T | O η * [ ϕ 1 , F u + H u ] + | A η | O T * [ ϕ 2 , F u + H u ] .
The right-hand side of (40) and (41) tends to zero as t 1 t 2 , independently of u , which implies that A ( B R ) is an equicontinuous set. By using the Arzelá–Ascoli theorem, the set A ( B R ) is compact.
(iii) Finally, we show that W = { u C : u = ξ A u , 0 < ξ < 1 } is a bounded set. Let u W . Then, as in (i), we have
| u ( t ) | ξ A x C [ λ F + μ H ] Φ + G 1 η α + β Θ T * + G 2 T p α + β Θ η * Γ p , q ( α + β + 1 ) + O 1 Θ T * + O 2 Θ η * ,
which yields
A u C [ λ F + μ H ] Φ + G 1 η α + β Θ T * + G 2 T p α + β Θ η * Γ p , q ( α + β + 1 ) + O 1 Θ T * + O 2 Θ η * .
Therefore, W is bounded.
Hence, by Schaefer’s fixed-point theorem, we deduce that the operator A has a fixed point, which is a solution of boundary value problem (4). The proof is finished. □

5. Examples

Example 1.
Let us consider the fractional ( p , q ) -integrodifference equation
D 2 3 , 1 2 4 3 u ( t ) = e [ cos 2 ( 2 π t ) + π ] 200 + e sin 2 ( 2 π t ) | u ( t ) | + e ( 10 t + π ) | Ψ 2 3 , 1 2 1 3 u ( t ) | 1 + | u ( t ) | + e [ sin 2 ( 2 π t ) + 10 ] ( t + 20 ) 2 | u ( t ) | + e ( 5 t + π 2 ) | Υ 2 3 , 1 2 1 4 u ( t ) | 1 + | u ( t ) | , t I 2 3 , 1 2 10
subject to fractional ( p , q ) -integral boundary condition
I 2 3 , 1 2 1 2 e + cos 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 ) , I 2 3 , 1 2 1 2 π + sin 15 2 u ( 15 ) = i = 0 D i | u ( t i ) | 1 + | u ( t i ) | , t i = σ 2 3 , 1 2 i ( 10 ) ,
where φ ( t , s ) = e | t s | ( 2 t + π ) 3 , ψ ( t , s ) = e | t s | ( t + e ) 3 and C i , D i are given constants, with 1 1000 i = 0 C i e 1000 and 1 500 i = 0 D i π 500 .
Here, p = 2 3 , q = 1 2 , α = 4 3 , β = 1 2 , γ = 1 3 , ν = 1 4 ,   T = 10 , η = σ 2 3 , 1 2 4 ( 10 ) = 1215 256 , λ = e π , μ = e 10 , ϕ 1 ( u ) = i = 0 C i | u ( t i ) | 1 + | u ( t i ) | , ϕ 2 = i = 0 D i | u ( t i ) | 1 + | u ( t i ) | , g 1 ( t ) = ( e + cos t ) 2 , g 2 ( t ) = ( π + sin t ) 2 ,
F t , u ( t ) , Ψ p , q γ u ( t ) = e cos 2 ( 2 π t ) 200 + e sin 2 ( 2 π t ) | u ( t ) | + e ( 10 t + π ) | Ψ 2 3 , 1 2 1 3 u ( t ) | 1 + | u ( t ) |
and
H t , u ( t ) , Υ p , q ν u ( t ) = e sin 2 ( 2 π t ) ( t + 20 ) 2 | u ( t ) | + e ( 5 t + π 2 ) | Υ 2 3 , 1 2 1 4 u ( t ) | 1 + | u ( t ) | .
For all t I 2 3 , 1 2 10 and u , v R , we have
F t , u , Ψ p , q γ u F t , v , Ψ p , q γ v 1 201 | u v | + 1 201 e π Ψ p , q γ u Ψ p , q γ v , H t , u , Υ p , q ν u F t , v , Υ p , q ν v 1 400 | u v | + 1 400 e π 2 Υ p , q ν u Υ p , q ν v .
Thus, ( H 1 ) and ( H 2 ) hold with M 1 = 0.004975 , M 2 = 0.0002145 , N 1 = 0.0025 and N 2 = 0.0005197 .
For all u , v C ,
| ϕ 1 ( u ) ϕ 1 ( v ) | e 1000 u v C , | ϕ 2 ( u ) ϕ 2 ( v ) | π 500 u v C .
So, ( H 3 ) holds with ω 1 = 0.00272 and ω 2 = 0.00628 .
Moreover, ( H 4 ) holds with g 1 = 2.9525 , G 1 = 13.8256 , g 2 = 4.5864 and G 2 = 17.1528 .
After calculating, we find that
| A η | 20.9307 , | A T | 32.7815 , | B η | 39.839 , | B T | 62.3957 , | Λ | 410.0706 , φ 0 = 0.0323   a n d   ψ 0 = 0.04979 .
We can show that
L 0.0002159 , Θ T 0.38884 , Θ η 0.24799 , Θ ¯ T 0.19883 , Θ ¯ η 0.12504 , Θ T * 1.5769 a n d Θ η * 1.0066 .
So, ( H 5 ) holds with
X 0.64592 < 1 .
Hence, by Theorem 1, the boundary value problem (43) and (44) has a unique solution on I 2 3 , 1 2 10 .
Example 2.
Let us consider the fractional ( p , q ) -integrodifference equation
D 2 3 , 1 2 4 3 u ( t ) = 1 10 t + 1 3 e ( t + π ) u ( t ) + Ψ 2 3 , 1 2 1 3 u ( t ) + 1 15 t + 2 5 e ( t + 10 ) u ( t ) + Υ 2 3 , 1 2 1 4 u ( t ) , t I 2 3 , 1 2 10
with fractional ( p , q ) -integral boundary condition
I 2 3 , 1 2 1 2 e + cos 1215 256 2 u 1215 256 = i = 0 C i e | u ( t i ) | , t i = σ 2 3 , 1 2 i ( 10 ) , I 2 3 , 1 2 1 2 π + sin 15 2 u ( 15 ) = i = 0 D i e | u ( t i ) | , t i = σ 2 3 , 1 2 i ( 10 ) ,
where C i and D i are given constants with 1 1000 i = 0 C i π 1000 and 1 500 i = 0 D i e 500 .
Here, p = 2 3 , q = 1 2 , α = 4 3 , β = 1 2 , γ = 1 3 , ν = 1 4 ,   T = 10 , η = σ 2 3 , 1 2 4 ( 10 ) = 1215 256 , λ = e π , μ = e 10 .
It is clear that F t , u , Ψ p , q γ u 23 15 = F , H t , u , Υ p , q ν u 77 75 = H for t I 2 3 , 1 2 10 and φ 1 ( u ) π 1000 = O 1 , φ 2 ( u ) e 500 = O 2 for u C .
Hence, ( H 5 ) and ( H 6 ) hold. Therefore, the boundary value problem (45) and (46) has at least one solution on I 2 3 , 1 2 10 by Theorem 2.

6. Conclusions

A nonlocal fractional ( p , q ) -integral boundary value problem for separate fractional ( p , q ) -integrodifference Equation ( 4 ) is studied. Our problem contains two fractional ( p , q ) -difference operators and two fractional ( p , q ) -integral operators. The existence of a unique solution is established via the Banach contraction mapping principle, while the existence result is proved using the Schaefer’s fixed-point theorem. In addition, some properties of the ( p , q ) -integral are also studied. It is imperative to mention that our results are new in the given configuration and enrich the literature on boundary value problems involving ( p , q ) -integrodifference equations. In the future, we plan to extend this work by considering new boundary value problems.

Author Contributions

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

Funding

This research study was funded by the Science, Research and Innovation Promotion Fund under Basic Research Plan—Saun Dusit University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research study was supported by Chiang Mai University.

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. 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]
  4. 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]
  5. Mursaleen, M.; Nasiruzzaman, M.; Khan, A.; Ansari, K.J. Some approximation results on Bleimann-Butzer-Hahn operators defined by (p,q)-integers. Filomat 2016, 30, 639–648. [Google Scholar] [CrossRef] [Green Version]
  6. 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]
  7. 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]
  8. Khan, K.; Lobiyal, D.K. Bézier curves based on Lupaş (p,q)-analogue of Bernstein functions in CAGD. J. Comput. Appl. Math. 2017, 317, 458–477. [Google Scholar] [CrossRef]
  9. Jagannathan, R.; Rao, K.S. Two-parameter quantum algebras, twin-basic number, and associated generalized hypergeometric series. arXiv 2006, arXiv:math/0602613. [Google Scholar]
  10. Sahai, V.; Yadav, S. Representations of two parameter quantum algebras and (p,q)-special functions. J. Math. Anal. Appl. 2007, 335, 268–279. [Google Scholar] [CrossRef]
  11. Sadjang, P.N. On the fundamental theorem of (p,q)-calculus and some (p,q)-Taylor formulas. arXiv 2018, arXiv:1309.3934. [Google Scholar]
  12. Duran, U. Post Quantum Calculus. Master’s Thesis, University of Gaziantep, Şehitkamil/Gaziantep, Turkey, 2016. [Google Scholar]
  13. 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]
  14. Cheng, W.T.; Zhang, W.H.; Cai, Q.B. (p,q)-gamma operators which preserve x2. J. Inequal. Appl. 2019, 2019, 108. [Google Scholar] [CrossRef] [Green Version]
  15. Soontharanon, J.; Sitthiwirattham, T. On fractional (p,q)-Calculus. Adv. Differ. Equ. 2020, 2020, 35. [Google Scholar] [CrossRef] [Green Version]
  16. Pheak, N.; Nonlaopon, K.; Tariboon, J.; Ntouyas, S.K. Fractional (p,q)-calculus on finite intervals and some Integral inequalities. Symmetry 2021, 3, 504. [Google Scholar]
  17. 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, 333. [Google Scholar]
  18. Promsakon, C.; Kamsrisuk, N.; Ntouyas, S.K.; Tariboon, J. On the second-order (p,q)-difference equations with separated boundary conditions. Adv. Math. Phys. 2018, 2018, 9089865. [Google Scholar] [CrossRef] [Green Version]
  19. 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]
  20. Nuntigrangjana, T.; Putjuso, S.; Ntouyas, S.K.; Tariboon, J. Impulsive quantum (p,q)-difference equations. Adv. Differ. Equ. 2020, 2020, 98. [Google Scholar] [CrossRef] [Green Version]
  21. Qin, Z.; Sun, S. Positive solutions for fractional (p,q)-difference boundary value problems. J. Appl. Math. Comput. 2021, 1–8. [Google Scholar] [CrossRef]
  22. Qin, Z.; Sun, S. On a nonlinear fractional (p,q)-difference Schrödinger equation. J. Appl. Math. Comput. 2021. [Google Scholar] [CrossRef]
  23. Soontharanon, J.; Sitthiwirattham, T. Existence results of nonlocal Robin boundary value problems for fractional (p,q)-integrodifference equations. Adv. Differ. Equ. 2020, 2020, 342. [Google Scholar] [CrossRef]
  24. Smart, D.R. Fixed Point Theorems; Cambridge University Press: Cambridge, UK, 1980. [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

Dumrongpokaphan, T.; Ntouyas, S.K.; Sitthiwirattham, T. Separate Fractional (p,q)-Integrodifference Equations via Nonlocal Fractional (p,q)-Integral Boundary Conditions. Symmetry 2021, 13, 2212. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13112212

AMA Style

Dumrongpokaphan T, Ntouyas SK, Sitthiwirattham T. Separate Fractional (p,q)-Integrodifference Equations via Nonlocal Fractional (p,q)-Integral Boundary Conditions. Symmetry. 2021; 13(11):2212. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13112212

Chicago/Turabian Style

Dumrongpokaphan, Thongchai, Sotiris K. Ntouyas, and Thanin Sitthiwirattham. 2021. "Separate Fractional (p,q)-Integrodifference Equations via Nonlocal Fractional (p,q)-Integral Boundary Conditions" Symmetry 13, no. 11: 2212. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13112212

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