Next Article in Journal / Special Issue
Analytical Study of a ϕ− Fractional Order Quadratic Functional Integral Equation
Previous Article in Journal
On the Semi-Local Convergence of a Fifth-Order Convergent Method for Solving Equations
Previous Article in Special Issue
Analytical and Qualitative Study of Some Families of FODEs via Differential Transform Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and Uniqueness of Solutions to a Nabla Fractional Difference Equation with Dual Nonlocal Boundary Conditions

by
Nandhihalli Srinivas Gopal
and
Jagan Mohan Jonnalagadda
*
Department of Mathematics, Birla Institute of Technology and Science Pilani, Hyderabad 500078, India
*
Author to whom correspondence should be addressed.
Submission received: 30 December 2021 / Revised: 17 January 2022 / Accepted: 17 January 2022 / Published: 21 January 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
In this paper, we look at the two-point boundary value problem for a finite nabla fractional difference equation with dual non-local boundary conditions. First, we derive the associated Green’s function and some of its properties. Using the Guo–Krasnoselkii fixed point theorem on a suitable cone and under appropriate conditions on the non-linear part of the difference equation, we establish sufficient requirements for at least one and at least two positive solutions of the boundary value problem. Next, we discuss the existence and uniqueness of solutions to the considered problem. For this purpose, we use Brouwer and Banach fixed point theorem, respectively. Finally, we provide a few examples to illustrate the applicability of established results.

1. Introduction

Over the last few decades, the theory of fractional calculus has been extensively developed due to its properties, generalizing most results of differential calculus and its non-local nature of fractional derivatives. The contributions of several mathematicians over the span of three centuries have resulted in a robust theory of fractional differential equations for real functions. Its roots can be traced back to the Leibniz letter dated “30 September 1695”. Today fractional calculus has been successfully used for mathematical modelling in the fields of medical sciences, computational biology, economics, physics and several areas of engineering in the past three decades. For further applications and historical literature, we refer here to a few classical texts on fractional calculus by Miller and Ross [1], Samko et al. [2], Podlubny [3] and Kilbas et al. [4].
On the other side of the coin, nabla fractional calculus is a branch of mathematics, which deals with arbitrary order differences and sums in the backward sense. The theory of nabla fractional calculus is relatively young, with the most prominent works done in the past decade. The notion of nabla fractional difference and sum can be traced back to the work of Gray and Zhang [5], and Miller and Ross [6]. In this line, Atici and Eloe [7] developed nabla fractional Riemann–Liouville difference operator, initiated the study of nabla fractional initial value problem and established exponential law, product rule, and nabla Laplace transform. Following their works, the contributions of several mathematicians have made the theory of discrete fractional calculus a fruitful field of research in science and engineering, we refer here a few applications of discrete fractional equations [8,9,10]. We also refer here to a recent monograph by Goodrich and Peterson [11] and the references therein, which is an excellent source for all those who wish to work in this field.
The study of boundary value problems (BVPs) has a long past and can be followed back to the work of Euler and Taylor on vibrating strings. On the discrete fractional side, there is a sudden growth in interest for the development of nabla fractional BVPs. Many authors have studied nabla fractional BVPs recently. To name a few, Ahrendt [12], Goar [13], and Ikram [14] worked with self-adjoint Caputo nabla BVPs. Brackins [15] studied a particular class of self-adjoint Riemann–Liouville nabla BVPs and derived the Green’s function associated with it along with a few of its properties. Gholami et al. [16] obtained the Green’s function for a non-homogeneous Riemann–Liouville nabla BVP with Dirichlet boundary conditions. Jonnalagadda [17,18,19,20,21,22] analyzed some qualitative properties of two-point non-linear Riemann–Liouville nabla BVPs associated with a variety of boundary conditions. Goodrich [23,24] analyzed discrete fractional BVP’s with non-local conditions in the delta case. To the best of our knowledge, there are no results for the solutions of nabla fractional BVP’s with non-local conditions.
We consider the following boundary value problem with dual non-local conditions
( ρ ( o ) ϑ v ) ( z ) = q ( z , v ( z ) ) , z N o + 2 p , v ( o ) = g 1 ( v ) , v ( p ) = g 2 ( v ) ,
where o , p R with p o N 3 , 1 < ϑ < 2 , q : N o + 2 p × R R and the functionals g 1 , g 2 : R R . The present paper is organized as follows: Section 2 contains a few preliminaries on discrete fractional calculus. In Section 3, we construct the Green’s function corresponding to (1) and state a few of its properties. In Section 4, we study the existence of at least one and at least two positive solutions using Guo–Kranoselskii fixed point theorem on cones. In Section 5, we obtain sufficient conditions on the existence and uniqueness of solutions for the proposed class of boundary value problems using Brouwer and contraction mapping theorems, respectively. Finally, we conclude this article with a few examples.

2. Preliminaries

The set of all real numbers and positive integers are denoted by R and Z + , respectively. We use the following notations, definitions, and known results of nabla fractional calculus [11]. We assume empty sums and products are 0 and 1, respectively.
Definition 1.
For o R , the sets N o and N o p , where p o Z + , are defined by
N o = { o , o + 1 , o + 2 , . . . } , N o p = { o , o + 1 , o + 2 , . . . , p } .
Definition 2.
We define the backward jump operator, ρ : N o + 1 N o , by
ρ ( z ) = z 1 , z N o + 1 .
Let v : N o R and M N 1 . The first order backward (nabla) difference of v is defined by ( v ) ( z ) = v ( z ) v ( z 1 ) , for z N o + 1 , and the Mth -order nabla difference of v is defined recursively by ( M v ) ( z ) = ( ( M 1 v ) ) ( z ) , for z N o + M .
Definition 3
(See [11]). Let z R \ { , 2 , 1 , 0 } and r R such that ( z + r ) R \ { , 2 , 1 , 0 } , the generalized rising function is defined by
z r ¯ = Γ ( z + r ) Γ ( z ) .
Here Γ ( · ) denotes the Euler gamma function. Also, if z { , 2 , 1 , 0 } and r R such that ( z + r ) R \ { , 2 , 1 , 0 } , then we use the convention that z r ¯ = 0 .
Definition 4
(See [11]). Let z, o R and ξ R \ { , 2 , 1 } . The ξ th-order nabla fractional Taylor monomial is given by
H ξ ( z , o ) = ( z o ) ξ ¯ Γ ( ξ + 1 ) ,
given that the right-hand side exists.
We observe the following properties of the nabla fractional Taylor monomials.
Lemma 1
(See [14,18]). Let ξ > 1 and t N o . Then the following hold:
i.
If z N ρ ( t ) , then H ξ ( z , ρ ( t ) ) 0 and if z N t , then H ξ ( z , ρ ( t ) ) > 0 .
ii.
If z N t and 1 < ξ < 0 , then H ξ ( z , ρ ( t ) ) is an increasing function of t.
iii.
If z N t + 1 and 1 < ξ < 0 , then H ξ ( z , ρ ( t ) ) is a decreasing function of z.
iv.
If z N ρ ( t ) and ξ > 0 , then H ξ ( z , ρ ( t ) ) is a decreasing function of t.
v.
If z N ρ ( t ) and ξ 0 , then H ξ ( z , ρ ( t ) ) is a non-decreasing function of z.
vi.
If z N t and ξ > 0 , then H ξ ( z , ρ ( t ) ) is an increasing function of z.
vii.
If 0 < v ξ , then H v ( z , o ) H ξ ( z , o ) , for each fixed z N o .
Lemma 2.
Let o , p be two real numbers such that 0 < o p and 1 < ϑ < 2 . Then ( o t ) ϑ 1 ¯ ( p t ) ϑ 1 ¯ is a decreasing function of t for t N 0 o 1 .
Proof. 
It is enough to show that ( o t ) ϑ 1 ¯ ( p t ) ϑ 1 ¯ < 0 .
Consider
( o t ) ϑ 1 ¯ ( p t ) ϑ 1 ¯ = ( p t ) ϑ 1 ¯ ( ϑ 1 ) ( o ρ ( t ) ) ϑ 2 ¯ + ( o t ) ϑ 1 ¯ ( ϑ 1 ) ( p ρ ( t ) ) ϑ 2 ¯ ( p t ) ϑ 1 ¯ ( p ρ ( t ) ) ϑ 1 ¯ = ( ϑ 1 ) ( p t ) ( p ρ ( t ) ) ϑ 2 ¯ ( o ρ ( t ) ) ϑ 2 ¯ + ( o t ) ( o ρ ( t ) ) ϑ 2 ¯ ( p ρ ( t ) ) ϑ 2 ¯ ( p t ) ϑ 1 ¯ ( p ρ ( t ) ) ϑ 1 ¯ = ( ϑ 1 ) ( p ρ ( t ) ) ϑ 2 ¯ ( o ρ ( t ) ) ϑ 2 ¯ ( p + t + o t ) ( p t ) ϑ 1 ¯ ( p ρ ( t ) ) ϑ 1 ¯ = ( ϑ 1 ) ( p ρ ( t ) ) ϑ 2 ¯ ( o ρ ( t ) ) ϑ 2 ¯ ( o p ) ( p t ) ϑ 1 ¯ ( p ρ ( t ) ) ϑ 1 ¯ .
Since p > o , it follows from Lemma 1 that ( o t ) ϑ 1 ¯ ( p t ) ϑ 1 ¯ < 0 . The proof is complete. □
Definition 5
(See [11]). Let v : N o + 1 R and ν > 0 . The ν th-order nabla sum of v is given by
( o ν v ) ( z ) = t = o + 1 z H ν 1 ( z , ρ ( t ) ) v ( t ) , z N o + 1 .
Definition 6
(See [11]). Let v : N o + 1 R , ν > 0 and choose A N 1 such that A 1 < ν A . The ν th-order Riemann–Liouville nabla difference of v is given by
( o ν v ) ( z ) = ( A ( o ( A ν ) v ) ) ( z ) , z N o + A .
Theorem 1
(See [11]). Assume μ > 0 and A 1 < μ A . Then a general solution of o μ v ( z ) = 0 is given by
v ( z ) = c 1 H μ 1 ( z , o ) + c 2 H μ 2 ( z , o ) + + c A H μ A ( z , o ) , for z N o .

3. Green’s Function

In this section, we will build the Green’s function formula for our boundary value problem (1) and derive a few properties from it that will be used throughout the rest of the paper.
Theorem 2
(See [15]). The following nabla fractional boundary value problem
( ρ ( o ) ϑ v ) ( z ) = h ( z ) , z N o + 2 p , v ( o ) = v ( p ) = 0 ,
where o , p R with p o N 3 , 1 < ϑ < 2 and h : N o + 2 p R . Then (2) has the unique solution
v ( z ) = t = o + 2 p G ( z , t ) h ( t ) , z N o p ,
where Green’s function G ( z , t ) is given by
G ( z , t ) = G 1 ( z , t ) = H ϑ 1 ( z , o ) H ϑ 1 ( p , o ) H ϑ 1 ( p , ρ ( t ) ) , z N o t 1 , G 2 ( z , t ) = H ϑ 1 ( z , o ) H ϑ 1 ( p , o ) H ϑ 1 ( p , ρ ( t ) ) H ϑ 1 ( z , ρ ( t ) ) , z N t p .
Lemma 3.
The equivalent form of the following homogeneous nabla fractional boundary value problem with non-local conditions
( ρ ( o ) ϑ w ) ( z ) = 0 , for z N o + 2 p , w ( o ) = g 1 ( w ) , w ( p ) = g 2 ( w ) ,
is given by
w ( z ) = g 1 ( w ) ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) + g 2 ( w ) ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ , z N o p .
Proof. 
From Theorem 1, the general solution of the equation ( ρ ( o ) ϑ w ) ( z ) = 0 , is given by
w ( z ) = c 1 ( z o + 1 ) ϑ 1 ¯ + c 2 ( z o + 1 ) ϑ 2 ¯ , z N o p ,
where c 1 and c 2 are arbitrary constants. Using w ( o ) = g 1 ( w ) and w ( p ) = g 2 ( w ) , respectively in (7), we have
g 1 ( w ) Γ ( ϑ 1 ) = c 1 ( ϑ 1 ) + c 2 , g 2 ( w ) = c 1 ( p o + 1 ) ϑ 1 ¯ + c 2 ( p o + 1 ) ϑ 2 ¯ .
Now, solving the above system of equations for c 1 and c 2 , we have
c 1 = g 1 ( w ) ( p o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) ( p o ) ϑ 1 ¯ + g 2 ( w ) ( p o ) ϑ 1 ¯ , c 2 = g 1 ( w ) Γ ( ϑ 1 ) ( ϑ 1 ) g 1 ( w ) ( p o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) ( p o ) ϑ 1 ¯ + g 2 ( w ) ( p o ) ϑ 1 ¯ .
Substituting c 1 and c 2 in (7), we have
w ( z ) = g 2 ( w ) ( p o ) ϑ 1 ¯ g 1 ( w ) ( p o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) ( p o ) ϑ 1 ¯ ( z o + 1 ) ϑ 1 ¯ + g 1 ( w ) Γ ( ϑ 1 ) ( ϑ 1 ) g 1 ( w ) ( p o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) ( p o ) ϑ 1 ¯ + g 2 ( w ) ( p o ) ϑ 1 ¯ ( z o + 1 ) ϑ 2 ¯ = g 1 ( w ) ( p o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) ( p o ) ϑ 1 ¯ ( ϑ 1 ) ( z o + 1 ) ϑ 2 ¯ ( z o + 1 ) ϑ 1 ¯ + g 2 ( w ) ( p o ) ϑ 1 ¯ ( z o + 1 ) ϑ 1 ¯ ( ϑ 1 ) ( z o + 1 ) ϑ 2 ¯ + g 1 ( w ) Γ ( ϑ 1 ) ( z o + 1 ) ϑ 2 ¯ = g 1 ( w ) Γ ( ϑ 1 ) ( z o + 1 ) ϑ 2 ¯ ( p o + 1 ) ϑ 2 ¯ ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ + g 2 ( w ) ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ = g 1 ( w ) ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) + g 2 ( w ) ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ .
The proof is complete. □
Lemma 4.
w satisfies the following property:
max z N o p w ( z ) g 1 ( w ) + g 2 ( w ) ,
where w is given by (7).
Proof. 
Consider
w ( z ) = ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) g 1 ( w ) + ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ g 2 ( w ) .
Clearly, function ( p z p o ) is decreasing with respect to z for z N o p . It follows from Lemma 1 that function ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) is decreasing with respect to z and ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ is an increasing function of z for z N o p . Thus, we have
max z N o p ( p z p o ) = 1 ,
max z N o p ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) = ( o o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) = 1 ,
max z N o p ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ = 1 ,
implying that
w ( z ) g 1 ( w ) + g 2 ( w ) , for z N o p .
The proof is complete. □
Theorem 3
(See [17]). Let 1 < ϑ < 2 and q : N o p × R R . The equivalent form of (1) is given by
v ( z ) = A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) + t = o + 2 p G ( z , t ) q ( t , v ( t ) ) , z N o p ,
where A ( z ) = ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) , B ( z ) = ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ and Green’s function G ( z , t ) is given by (4), respectively.
Theorem 4
(See [15,17]). The Green’s function G ( z , t ) defined in (4) satisfies the following properties:
i.
G ( o , t ) = G ( p , t ) = 0 , for all s N o + 1 p .
ii.
G ( z , o + 1 ) = 0 , for all z N o p .
iii.
G ( z , t ) > 0 , for all ( z , t ) N o + 1 p 1 × N o + 2 p .
iv.
max z N o + 1 p 1 G ( z , t ) = G ( t 1 , t ) for all t N o + 2 p .
v.
t = o + 1 p G ( z , t ) Λ , for all ( z , t ) N o p × N o + 1 p , where
Λ = p o 1 ϑ Γ ( ϑ + 1 ) ( ϑ 1 ) ( p o ) + 1 ϑ ϑ 1 ¯ .

4. Existence of Positive Solution

In this section, we show the existence of at least one positive solution for the boundary value problem (1), using Guo–Krasnoselskii fixed point theorem on a suitable cone.
Definition 7.
Let B be a Banach space over R . A closed nonempty subset C of B is said to be a cone provided,
(i)
e x + f y C , for all x , y C and all e , f 0 , and
(ii)
x C and x C implies x = 0 .
Definition 8.
An operator is called completely continuous if it is continuous and maps bounded sets into pre-compact sets.
Lemma 5
(See [25]). [Guo–Krasnoselskii fixed point theorem] Let B be a Banach space and C B be a cone. Assume that Ω 1 and Ω 2 are open sets contained in B such that 0 Ω 1 and Ω ¯ 1 Ω 2 . Assume further that T : C ( Ω ¯ 2 \ Ω 1 ) C is a completely continuous operator. If either
i.
T v v for v C Ω 1 and T v v for v C Ω 2 ; or
ii.
T v v for v C Ω 1 and T v v for v C Ω 2 ;
holds, then T has at least one fixed point in C ( Ω ¯ 2 \ Ω 1 ) .
We establish the following lemmas, which will be used later in the proof of our main result.
Lemma 6.
There exits a number λ ( 0 , 1 ) , such that
min z N c d G ( z , t ) λ max z N o p G ( z , t ) = λ G ( t 1 , t ) ,
where, c , d N o + 1 p 1 , c = o + p o + 1 4 and d = o + 3 p o + 1 4 .
Proof. 
By using the properties of Green’s function and Taylor monomial from Definition 4, Lemma 1 and Theorem 4, respectively.
Consider, for t N o + 2 p ,
G ( z , t ) G ( t 1 , t ) = ( z o ) ϑ 1 ¯ ( t o 1 ) ϑ 1 ¯ , for t > z , ( z o ) ϑ 1 ¯ ( t o 1 ) ϑ 1 ¯ ( z t + 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ , for t z .
Now, for t > z and c z d . Function G 1 ( z , t ) is increasing with respect to z. Then,
min z N c d G 1 ( z , t ) = G 1 ( c , t ) = ( c o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ .
For z t and c z d , function G 2 ( z , t ) is decreasing with respect to z. Then,
min z N c d G 2 ( z , t ) = G 2 ( d , t ) , = ( d o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ ( d t + 1 ) ϑ 1 ¯ .
Thus,
min z N c d G ( z , t ) = G 2 ( d , t ) , for t N o + 2 c , min G 2 ( d , t ) , G 2 ( c , t ) , for t N c + 1 d 1 , G 1 ( c , t ) , for t N d p , = G 2 ( d , t ) , for t N o + 2 r , G 1 ( c , t ) , for t N r p ,
where c < r < d . Consider
min z N c d G ( z , t ) G ( t 1 , t ) = ( d o ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ ( d t + 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ , for t N o + 2 r , ( c o ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ , for t N r p .
Thus,
min z N c d G ( z , t ) λ ( t ) max z N o p G ( z , t ) ,
where
λ ( t ) = min ( c o ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ , ( d o ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ ( d t + 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ .
Let for t N r p , denote
λ 1 ( t ) = ( c o ) ϑ 1 ¯ ( t o + 1 ) ϑ 1 ¯ ( c o ) ϑ 1 ¯ ( p o + 1 ) ϑ 1 ¯ .
Similarly, for t N o + 2 r , we take
λ 2 ( t ) = 1 ( t o + 1 ) ϑ 1 ¯ ( d o ) ϑ 1 ¯ ( d t + 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ .
By Lemmas 1 and 2, we see that function ( d t + 1 ) ϑ 1 ¯ ( p t + 1 ) ϑ 1 ¯ is decreasing for t N o + 2 r . Then
λ 2 ( t ) 1 ( t o + 1 ) ϑ 1 ¯ ( d o ) ϑ 1 ¯ ( d o 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p o 1 ) ϑ 1 ¯ 1 ( d o ) ϑ 1 ¯ ( d o ) ϑ 1 ¯ ( d o 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p o 1 ) ϑ 1 ¯ .
Thus,
min z N c d G ( z , t ) λ max z N o p G ( z , t ) ,
where
λ = min ( c o ) ϑ 1 ¯ ( p o + 1 ) ϑ 1 ¯ , 1 ( d o 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p o 1 ) ϑ 1 ¯ ( d o ) ϑ 1 ¯ .
Since G 1 ( c , t ) > 0 and G 2 ( d , t ) > 0 , we have λ ( t ) > 0 for all t N o + 2 p , implying λ > 0 . It would be suffice to prove that one of the terms
( c o ) ϑ 1 ¯ ( p o + 1 ) ϑ 1 ¯ , 1 ( d o 1 ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ ( p o 1 ) ϑ 1 ¯ ( d o ) ϑ 1 ¯ < 1 .
It follows from Lemma 1 that
( c o ) ϑ 1 ¯ ( p o + 1 ) ϑ 1 ¯ < 1 .
Therefore, we conclude that λ ( 0 , 1 ) . The proof is complete. □
Lemma 7.
There exists a number λ 0 ( 0 , 1 ) such that
min z N c d w ( z ) λ 0 max z N o p A ( z ) g 1 ( w ) + max z N o p B ( z ) g 2 ( w ) .
where w is given by (6).
Proof. 
Clearly, function ( p z p o ) is decreasing with respect to z for z N o p . It follows from Lemma 1 that ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) is a decreasing function of z and function ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ is an increasing with respect to z for z N o p . Then, there exists M 1 , M 2 > 0 such that
min z N c d ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) = M 1 , min z N c d ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ = M 2 < 1 .
Take λ 0 = min ( M 1 , M 2 ) . Clearly λ 0 ( 0 , 1 ) . Thus, for all z N c d , we have
w ( z ) min z N c d ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) g 1 ( w ) + min z N c d ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ g 2 ( w ) = M 1 g 1 ( w ) + M 2 g 2 ( w ) λ 0 g 1 ( w ) + λ 0 g 2 ( w ) = λ 0 max z N o p ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) g 1 ( w ) + λ 0 max z N o p ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ g 2 ( w ) = λ 0 max z N o p A ( z ) g 1 ( w ) + max z N o p B ( z ) g 2 ( w ) .
The proof is complete. □
Lemma 8.
If g 1 , g 2 are non-negative, then there exists a constant λ ¯ ( 0 , 1 ) , such that
min z N c d t = o + 2 p G ( z , t ) q ( z , v ( t ) ) + min z N c d ( A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) )
λ ¯ max z N o p s = o + 2 p G ( z , t ) q ( z , v ( t ) ) + A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) .
Proof. 
It follows from Lemmas 6 and 7 have
min z N c d t = o + 2 p G ( z , t ) f ( t , v ( t ) ) + min z N c d ( A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) ) t = o + 2 p min z N c d G ( z , t ) q ( t , v ( t ) ) + λ 0 max z N o p A ( z ) g 1 ( v ) + max z N o p B ( z ) g 2 ( v ) λ t = o + 2 p max z N o p G ( z , t ) q ( t , v ( t ) ) + λ 0 max z N o p A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) λ ¯ max z N o p t = o + 2 p G ( z , t ) f ( t , v ( t ) ) + A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) .
where λ ¯ = min ( λ , λ 0 ) ( 0 , 1 ) . The proof is complete. □
We observe by Theorem 3 that the equivalent form of (1) is given by
v ( z ) = A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) + t = o + 2 p G ( z , t ) q ( t , v ( t ) ) , z N o p ,
where A ( z ) = ( p z p o ) ( z o + 1 ) ϑ 2 ¯ Γ ( ϑ 1 ) , B ( z ) = ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ and Green’s function G ( z , t ) is given by (4), respectively.
Note that any solution v : N o p R of (1) can be viewed as a real ( p o + 1 ) -tuple vector. Consequently, v R p o + 1 . Define the operator T : R p o + 1 R p o + 1 by
( T v ) ( z ) = A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) + t = o + 2 p G ( z , t ) q ( t , v ( t ) ) , z N o p .
Clearly, v is a fixed point of T if and only if v is a solution of (1). We use the fact that B = R p o + 1 is a Banach space equipped with the maximum norm
v = max z N o p | v ( z ) | ,
for any v B . We define the cone C by
C = v B : v ( z ) 0 and min z N c d v ( z ) λ ¯ v ( z ) .
Note that summation operator T is defined on a discrete finite set. Hence, operator T is trivially completely continuous.
Take
η = 1 t = o + 2 p G ( t 1 , t ) .
We state here the following hypothesis which will be used later
(H1)
q ( z , v ) 0 , ( z , v ) N o p × [ 0 , ) and g 1 ( v ) , g 2 ( v ) 0 , v C [ N o + 2 p R ] ,
(H2)
There exists a number r > 0 such that q ( z , v ) r η 2 , whenever 0 v r ,
(H3)
There exists a number r > 0 such that q ( z , v ) r η λ ¯ , whenever λ ¯ r v r ,
(H4)
Assume that lim v 0 + min z N o p q ( z , v ) v = , lim v min z N o p q ( z , v ) v = ,
(H5)
Assume that lim v 0 + min z N o p q ( z , v ) v = 0 , lim v min z N o p q ( z , v ) v = 0 ,
(G1)
The functionals g 1 ( v ) and g 2 ( v ) are linear. In particular, we assume that
g 1 ( v ) = i = o p m i v ( i ) and g 2 ( v ) = i = o p n i v ( i ) ,
where m i , n i R ,
(G2)
Assume i = o p m i + i = o p n i 1 2 .
Lemma 9.
Assume ( H 1 ) holds. Then, T : C C .
Proof. 
Let operator T be as defined in (18), then by using Lemma 8, we have
min z N c d ( T v ) ( z ) min z N c d t = o + 2 p G ( z , t ) q ( t , v ( t ) ) + min z N c d ( A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) ) λ ¯ max z N o p t = o + 2 p G ( z , t ) q ( z , v ( t ) ) + A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) . λ ¯ T v .
It is obvious that ( T v ) ( z ) 0 , whenever v C thus T : C C . □
Theorem 5.
Assume q ( z , v ) satisfy conditions { ( H 1 ) , ( H 2 ) , ( H 3 ) } . Also, g 1 ( v ) , g 2 ( v ) satisfy conditions { ( H 1 ) , ( G 1 ) , ( G 2 ) } . Then, the boundary value problem (1) has at least one positive solution.
Proof. 
We know that, T : C C is completely continuous. For 0 < r 1 < r 2 , define the set Ω 1 = { v C : v < r 1 } . Clearly, Ω 1 B is an open set with 0 Ω 1 . Since v = r 1 for v Ω 1 , ( H 2 ) , ( G 1 ) , ( G 2 ) holds for all v Ω 1 . Then we have for v C Ω 1 , by using Lemma 4.
T v max z N p o t = o + 2 p G ( z , t ) q ( t , v ( t ) ) + max z N o p A ( z ) g 1 ( v ) + max z N o p B ( z ) g 2 ( v ) t = o + 2 p max z N p o G ( z , t ) q ( t , v ( t ) ) + g 1 ( v ) + g 2 ( v ) η r 1 2 t = o + 2 p G ( t 1 , t ) + g 1 ( v ) + g 2 ( v ) 1 2 + i = o p m i + i = o p n i r 1 r 1 .
Therefore T v v , for v C Ω 1 . Similarly, set Ω 2 = { v C : v < r 2 } . Clearly, Ω 2 β is an open set and Ω ¯ 1 Ω 2 . Since v = r 2 for v Ω 2 , ( H 3 ) , ( G 1 ) , ( G 2 ) holds. We have
T v min z N c d | T v ( z ) | min z N c d t = o + 2 p G ( z , t ) q ( t , v ( t ) ) + min z N c d A ( z ) g 1 ( v ) + min z N c d B ( z ) g 2 ( v ) λ ¯ t = o + 2 p max z N p o G ( z , t ) q ( t , v ( t ) ) λ ¯ t = o + 2 p G ( t 1 , t ) q ( t , v ( t ) ) λ ¯ r 2 η λ ¯ t = o + 2 p G ( t 1 , t ) = r 2 .
Thus, we have T v v , for v Ω 2 C . By part (i) of Theorem 5, we conclude that operator T has a fixed point v 0 in C ( Ω ¯ 2 \ Ω 1 ) , satisfying r 1 < v 0 < r 2 . The proof is complete. □
Theorem 6.
Assume q ( z , v ) satisfy conditions { ( H 1 ) , ( H 2 ) , ( H 4 ) } . Also, g 1 ( v ) , g 2 ( v ) satisfy conditions { ( H 1 ) , ( G 1 ) , ( G 2 ) } . Then, the boundary value problem (1) has at least two positive solutions v 1 and v 2 with 0 < v 1 < r 2 < v 2 .
Proof. 
Let us choose a number N > 0 such that
N λ ¯ η > 1 ,
by condition (H4) there exists a number r * > 0 such that r * < r 1 < r 2 and q ( z , v ) N v for v [ 0 , r * ] and z N o p . Define the set Ω r * = { v C : v < r * } . Then we have
T v min z N c d | T v ( z ) | min z N c d t = o + 2 p G ( z , t ) q ( t , v ( t ) ) + min z N c d A ( z ) g 1 ( v ) + min z N c d B ( z ) g 2 ( v ) t = o + 2 p min z N c d G ( z , t ) q ( t , v ( t ) ) λ ¯ t = o + 2 p max z N p o G ( z , t ) q ( t , v ( t ) ) λ ¯ N v t = o + 2 p G ( t 1 , t ) > r * .
Thus, we have T v > v , for v Ω r * C . Next for the same N, we can find a number R 1 > 0 such that q ( z , v ) N v for v R 1 and z N o p . Choose R so that R = max r 2 , R 1 λ ¯ . Define the set Ω R = { v C : v < R } . We can show that T v v , for v Ω R C .
Finally define the set
Ω 2 = { v C : v < r 2 } .
Since v = r 2 condition (H2) holds for all v Ω 2 . Then, we have
T v max z N p o t = o + 2 p G ( z , t ) q ( t , v ( t ) ) + max z N o p A ( z ) g 1 ( v ) + max z N o p B ( z ) g 2 ( v ) t = o + 2 p max z N p o G ( z , t ) q ( t , v ( t ) ) + g 1 ( v ) + g 2 ( v ) r 2 η 2 t = o + 2 p G ( t 1 , t ) + g 1 ( v ) + g 2 ( v ) 1 2 + i = o p m i + i = o p n i r 2 r 2 .
Implying T v v , for v Ω r 2 C . Hence, we conclude that T has at least two fixed points say v 1 Ω 2 \ Ω ^ r * and v 2 Ω R \ Ω ^ 2 , where Ω ^ denoted the interior of the set Ω . In particular (1) has at least two positive solutions, say v 1 and v 2 satisfying 0 < v 1 < r 2 < v 2 . The proof is complete. □
Theorem 7.
Assume q ( z , v ) satisfy conditions { ( H 1 ) , ( H 3 ) , ( H 5 ) } . Also, g 1 ( v ) , g 2 ( v ) satisfy conditions { ( H 1 ) , ( G 1 ) , ( G 2 ) } . Then, the boundary value problem (1) has at least two positive solutions v 1 and v 2 with 0 < v 1 < r 2 < v 2 .

5. Existence of Solutions

In this section, we present the existence and uniqueness results of the boundary value problem (1) using Brouwer fixed point theorem and contraction mapping theorem, respectively and also, we construct an example to illustrate the same results.
Theorem 8
(See [25]). [Brouwer fixed point theorem]. Let C 0 be a nonempty compact convex subset of R n and T be a continuous mapping of C 0 into itself. Then, T has a fixed point in C 0 .
Theorem 9.
Assume q ( z , v ) , g 1 ( v ) and g 2 ( v ) is continuous with respect to v , for each z N o p . Assume there exist a positive constant L, such that
max L v L { g 1 ( v ) , g 2 ( v ) } L 2 ( Λ + 1 )
and
max ( z , v ) N o p × [ L , L ] { q ( z , v ) } L Λ + 1 ,
where Λ is given by (9). Then, the boundary value problem (1) has a solution.
Proof. 
Consider
C L = { v N o p R , v L } .
Clearly, C L is a non-empty compact convex subset of R p o + 1 . Let T be a operator as defined in (18). It is clear that T is a continuous operator. Therefore the main objective is to show that T : C L C L , then Theorem 8 can be invoked.
Let
Ω 0 = L Λ + 1 ,
By using (22), (23) and Lemma 4
T v max z N o p | A ( z ) g 1 ( v ) + B ( z ) g 2 ( v ) | + max z N o p t = o + 2 p G ( z , t ) q ( t , v ( t ) ) g 1 ( v ) + g 2 ( v ) + ( L Λ + 1 ) t = o + 2 p max z N o p G ( z , t ) ( L 2 ( Λ + 1 ) ) + ( L 2 ( Λ + 1 ) ) + ( L Λ + 1 ) Λ = Ω 0 ( 1 + Λ ) = L .
Thus, T v L and T : C L C L . It follows at once by Brouwer fixed point theorem, that there exist a fixed point of T, say v 0 C , such that | v 0 ( z ) | L . The proof is complete. □
Theorem 10
(See [25]). [Contraction Mapping Theorem]. Let S be a closed subset of R n . Assume T : S S is a contraction mapping, i.e., there exists a number ξ , 0 ξ 1 , such that
T v T u ξ v u ,
for all v , u S . Then, T has a unique fixed point v 0 S .
Theorem 11.
Assume that q ( z , v ) , g 1 ( v ) and g 2 ( v ) are Lipschitz with respect to v , i.e., there exists a , b , c > 0 , such that | q ( z , v 1 ) q ( z , v 2 ) | a v 1 v 2 , | g 1 ( v 1 ) g 1 ( v 2 ) | b v 1 v 2 and | g 2 ( v 1 ) g 2 ( v 2 ) | c v 1 v 2 , whenever v 1 , v 2 C [ N o p R ] . Then the boundary value problem (1) has a unique solution provided
a Λ + b + c < 1 ,
holds.
Proof. 
Consider
T v 1 T v 2 = max z N o p | T v 1 T v 2 | max z N o p t = o + 2 p G ( z , t ) q ( t , v 1 ) max z N o p t = o + 2 p G ( z , t ) f ( t , v 2 ) + max z N o p ( p z p o ) H ϑ 2 ( z , ρ ( o ) ) g 1 ( v 1 ) max z N o p ( p z p o ) H ϑ 2 ( z , ρ ( o ) ) g 1 ( v 2 ) + max z N o p ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ g 2 ( v 1 ) max z N o p ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ g 2 ( v 2 ) t = o + 2 p max z N o p G ( z , t ) | q ( t , v 1 ) q ( t , v 2 ) | + max z N o p ( p z p o ) H ϑ 2 ( z , ρ ( o ) ) | g 1 ( v 1 ) g 1 ( v 2 ) | + max z N o p ( z o ) ϑ 1 ¯ ( p o ) ϑ 1 ¯ | g 2 ( v 1 ) g 2 ( v 2 ) | Λ a v 1 v 2 + b v 1 v 2 + c v 1 v 2 ( Λ a + b + c ) v 1 v 2 .
Thus, using (24) T is a contraction on R p o + 1 . Hence, by Theorem 10, the result follows. The proof is complete. □
Example 1.
Suppose, ϑ = 1.1 and o = 0 , p = 10 , q ( z , v ) = s i n ( v ) 15 + z , g 1 ( v ) = t = o p v ( t ) 20 and g 2 ( v ) = t = o p v ( t ) 10 . Then, (1) becomes
( ρ ( 0 ) 1.1 v ) ( z ) = s i n ( v ) 15 + z , z N 2 10 , v ( 0 ) = t = o p v ( t ) 20 , v ( 10 ) = t = o p v ( t ) 10 .
Clearly, q ( z , v ) , g 1 ( v ) and g 2 ( v ) are Lipschitz with respect to v with Lipschitz constant a , b and c , respectively. Here a = 1 15 , b = 1 20 and c = 1 10 . Then
Λ = ( ( p o ) ( ϑ 1 ) + 1 ϑ ) ϑ 1 ¯ ( ( p o 1 ) ϑ Γ ( ϑ + 1 ) ) = 8.05 .
and
( a Λ + b + c ) = 0.6866 < 1 .
Thus, by Theorem 11 the boundary value problem (25) has a unique solution.
Example 2.
Suppose, ϑ = 1.5 and o = 0 , p = 6 , q ( z , v ) = 1 20 v + v 2 , g 1 ( v ) = 1 11 v ( 1 ) 1 23 v ( 4 ) and g 2 ( v ) = 1 9 v ( 5 ) 1 10 v ( 2 ) . Then, (1) becomes
( ρ ( 0 ) 1.5 v ) ( z ) = 1 20 v + v 2 , z N 2 6 , v ( 0 ) = 1 11 v ( 1 ) 1 23 v ( 4 ) , v ( 6 ) = 1 9 v ( 5 ) 1 10 v ( 2 ) .
Here for a given values of o , p we get λ ¯ = 0.1191 and η = 0.77 . We see that q ( z , v ) satisfies condition (H1) and (H4), also by taking r 2 = 2 we have that
q ( z , v ) = 1 20 v + v 2 1 20 r 2 + r 2 2 = 0.270 < η r 2 2 = 0.77 ,
and g 1 , g 2 satisfy conditions (G1) and (G2). Thus, all conditions of Theorem 6 are satisfied, hence (26) has at least two positive solutions v 1 and v 2 such that 0 < v 1 < 2 < v 2 .

6. Conclusions

In this paper, we constructed the Green’s function corresponding to the boundary value problem (1), and stated a few of its properties. We studied the existence of at least one and at least two positive solutions for the boundary value problem (1) using Guo–Kranoselskii fixed point theorem on cones, and obtained sufficient conditions on the existence and uniqueness of solutions for the proposed class of boundary value problems using Brouwer and contraction mapping theorems, respectively. Finally, we demonstrated the applicability of established results with a few examples. To the best of our knowledge, there is no prior progress made with non-local conditions in the nabla case.

Author Contributions

All authors equally contributed this paper and approved the final version. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to express their full thanks to the respected editorial office and also author N. S. Gopal acknowledges support received from CSIR-HRDG New Delhi, Government of India.

Conflicts of Interest

We declare that no competing interest exist regarding this manuscript.

References

  1. Miller, K.S.; Ross, B. Fractional difference calculus. In Proceedings of the International Symposium on Univalent Functions, Fractional Calculus, and Their Applications, Koriyama, Japan, 1–5 May 1988; pp. 139–152. [Google Scholar]
  2. Kilbas, A.A.; Marichev, O.I.; Samko, S.G. Fractional Integrals and Derivatives. Theory and Applications; Translated from the 1987 Russian Origina. Revised by the authors; Gordon and Breach Science Publishers: Yverdon, Switzerland, 1993; Volume xxxvi, 976p. [Google Scholar]
  3. Podlubny, I. Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, 198. Mathematics in Science and Engineering; Academic Press, Inc.: San Diego, CA, USA, 1999. [Google Scholar]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations. 204. North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Gray, H.L.; Zhang, N.F. On a new definition of the fractional difference. Math. Comput. 1988, 50, 513–529. [Google Scholar] [CrossRef]
  6. Miller, K.S.; Ross, B. Ellis Horwood Series: Mathematics and Its Applications; Srivastava, H.M., Owa, S., Eds.; Ellis Horwood Ltd.: Chichester, UK; Halsted Press [John Wiley and Sons, Inc.]: New York, NY, USA, 1989; 404p. [Google Scholar]
  7. Atici, F.M.; Eloe, P.W. Discrete fractional calculus with the nabla operator. Electron. J. Qual. Theory Differ. Equ. 2009, I, 12. [Google Scholar]
  8. Wu, F. Nabla Fractional Calculus and Its Application in Analyzing Tumor Growth of Cancer. Master’s Thesis, Western Kentucky University, Bowling Green, KY, USA, 2012. [Google Scholar]
  9. Huang, L.; Bae, Y. Chaotic Dynamics of the Fractional-Love Model with an External Environment. Entropy 2018, 20, 53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Huang, L.; Bae, Y.-C. Non-linear Behavior in Fractional-Order Romeo and Juliet’s Love Model Influenced by External Force with Fuzzy Function. Int. J. Fuzzy Syst. 2018, 21, 2018. [Google Scholar]
  11. Goodrich, C.; Peterson, A.C. Discrete Fractional Calculus; Springer: Cham, Switzerland, 2015. [Google Scholar]
  12. Ahrendt, K.; De Wolf, L.; Mazurowski, L.; Mitchell, K.; Rolling, T.; Veconi, D. Initial and Boundary Value Problems for the Caputo Fractional Self-Adjoint Difference Equations. Enlight. Pure Appl. Math. 2016, 2, 1. [Google Scholar]
  13. Goar, J.S. A Caputo Boundary Value Problem in Nabla Fractional Calculus. Ph.D. Thesis, The University of Nebraska, Lincoln, NE, USA, 2016; 112p. [Google Scholar]
  14. Ikram, A. Lyapunov inequalities for nabla Caputo boundary value problems. J. Differ. Equ. Appl. 2019, 25, 757–775. [Google Scholar] [CrossRef] [Green Version]
  15. Brackins, A. Boundary Value Problems of Nabla Fractional Difference Equations. Ph.D. Thesis, The University of Nebraska, Lincoln, NE, USA, 2014; 92p. [Google Scholar]
  16. Gholami, Y.; Ghanbari, K. Coupled systems of fractional -difference boundary value problems. Differ. Equ. Appl. 2016, 8, 459–470. [Google Scholar] [CrossRef] [Green Version]
  17. Jonnalagadda, J.M. On two-point Riemann-Liouville type nabla fractional boundary value problems. Adv. Dyn. Syst. Appl. 2018, 13, 141–166. [Google Scholar]
  18. Jonnalagadda, J.M. On a nabla fractional boundary value problem with general boundary conditions. AIMS Math. 2020, 5, 204–215. [Google Scholar] [CrossRef]
  19. Jonnalagadda, J.M. Existence results for solutions of nabla fractional boundary value problems with general boundary conditions. Adv. Theory Non-Linear Anal. Appl. 2020, 4, 29–42. [Google Scholar] [CrossRef] [Green Version]
  20. Jonnalagadda, J.M. An ordering on Green’s function and a Lyapunov-type inequality for a family of nabla fractional boundary value problems. Fract. Differ. Calc. 2019, 9, 109–124. [Google Scholar] [CrossRef] [Green Version]
  21. Jonnalagadda, J.M. Discrete fractional Lyapunov-type inequalities in nabla sense. Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 2020, 27, 397–419. [Google Scholar]
  22. Jonnalagadda, J.M. Lyapunov-type inequalities for discrete Riemann-Liouville fractional boundary value problems. Int. J. Differ. Equ. 2018, 13, 85–103. [Google Scholar]
  23. Goodrich, C. Existence and uniqueness of solutions to a fractional difference equation with non-local conditions. Comput. Math. Appl. 2011, 61, 191–202. [Google Scholar] [CrossRef]
  24. Goodrich, C.S. On positive solutions to non-local fractional and integer-order difference equations. Appl. Anal. Discrete Math. 2011, 5, 122–132. [Google Scholar] [CrossRef]
  25. Agarwal, R.P.; Meehan, M.; O’regan, D. Fixed Point Theory and Applications. 141. Cambridge Tracts in Mathematics; Cambridge University Press: Cambridge, UK, 2001. [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

Gopal, N.S.; Jonnalagadda, J.M. Existence and Uniqueness of Solutions to a Nabla Fractional Difference Equation with Dual Nonlocal Boundary Conditions. Foundations 2022, 2, 151-166. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010009

AMA Style

Gopal NS, Jonnalagadda JM. Existence and Uniqueness of Solutions to a Nabla Fractional Difference Equation with Dual Nonlocal Boundary Conditions. Foundations. 2022; 2(1):151-166. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010009

Chicago/Turabian Style

Gopal, Nandhihalli Srinivas, and Jagan Mohan Jonnalagadda. 2022. "Existence and Uniqueness of Solutions to a Nabla Fractional Difference Equation with Dual Nonlocal Boundary Conditions" Foundations 2, no. 1: 151-166. https://0-doi-org.brum.beds.ac.uk/10.3390/foundations2010009

Article Metrics

Back to TopTop