Search for Articles:

Contents

Solvability of discrete second-order two-point boundary value problems

Jagan Mohan Jonnalagadda1
1Department of Mathematics, Birla Institute of Technology and Science Pilani, Hyderabad, Telangana, India – 500078
Copyright © Jagan Mohan Jonnalagadda. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This article considers a second-order difference equation with constant coefficients in its standard form and two different classes of two-point homogeneous boundary conditions. First, we construct the corresponding Green functions and derive some important properties for further analysis. Next, we propose adequate conditions for the existence of solutions to the considered boundary value problems. Finally, we offer two examples to show the applicability of the main results.

Keywords: second-order difference equation, constant coefficients, standard form, boundary value problem, Green’s function, fixed point, existence

Introduction

Merdivenci [1] considered the following second-order boundary value problem for difference equations, and established criteria for the existence of solutions by an application of a fixed point theorem: \[ \left(\Delta^2 y\right)(\zeta – 1) = f(\zeta, y(\zeta)), \quad \zeta \in \mathbb{N}_{1}^{T}, \label{DE 1} \ \tag{1}\] \[A y(0) – B \left(\Delta y\right)(0) = 0, \quad C y(T) + D \left(\Delta y\right)(T) = 0, \label{BC 1} \tag{2}\] where \(A\), \(B\), \(C \geq 0\), \(A^{2} + B^{2} > 0\), \(D > 0\) with \(ACT + AD + BC > 0\); \(T \in \mathbb{N}_{3}\); \(y: \mathbb{N}^{T + 1}_{0} \rightarrow \mathbb{R}\), and \(f : \mathbb{N}^{T}_{1} \times \mathbb{R} \rightarrow \mathbb{R}\) is a continuous function with respect to its second argument. Also, Atici [2] proved the existence of solutions for the nonlinear discrete Sturm–Liouville equation \[\label{DE 2} \left(\Delta^2 y\right)(\zeta – 1) – \alpha y(\zeta) = f(\zeta, y(\zeta)), \quad \zeta \in \mathbb{N}_{1}^{T}, \tag{3}\] with \(\alpha \geq 0\) subject to two-point homogeneous separated boundary conditions (2) using some results of differentiable operators. Following these works, Aykut et al. [3] investigated the existence and uniqueness of solutions for the boundary value problem (3) – (2) using fixed point theory under the following assumptions: \(\alpha \geq 0\); \(A\), \(B\), \(C\), \(D \geq 0\), \(A + B > 0\), \(C + D > 0\), \(A + C > 0\) (if \(\alpha = 0\)). Further, in [46], Atici et al. proved the existence of solutions of the second-order nonlinear difference Eq. (3) for \(\alpha > 0\) with periodic boundary conditions \[\label{BC 2} y(0) = y(T), \quad \left(\Delta y\right)(0) = \left(\Delta y\right)(T), \tag{4}\] using suitable fixed point theorems.

In 2015, Lyons and Neugebauer [7] studied the second-order nonlinear difference Eq. (1) satisfying the anti-periodic boundary conditions \[\label{BC 3} y(0) + y(T) = 0, \quad \left(\Delta y\right)(0) + \left(\Delta y\right)(T) = 0. \tag{5}\]

It is obvious that the second-order difference Eq. (1) associated with (4) or any of the following pairs of two-point boundary conditions is at resonance. \[ \left(\Delta y\right)(0) = 0, \quad \left(\Delta y\right)(T) = 0, \label{BC 4} \ \tag{6}\] \[y(0) = 0, \quad \left(\Delta y\right)(0) = \left(\Delta y\right)(T). \label{BC 5} \tag{7}\]

At the same time, the boundary value problems (3) – (4) and (3) – (6) are nonresonant for \(\alpha > 0\). Recently, the authors of [810] observed that the second-order difference Eq. (1) associated with the following pairs of two-point boundary conditions is also at resonance. \[ y(0) = y(T), \quad\;\; \left(\Delta y\right)(0) + \left(\Delta y\right)(T) = 0, \label{BC 6} \ \tag{8}\] \[y(0) + y(T) = 0, \quad\qquad \left(\Delta y\right)(0) = \left(\Delta y\right)(T). \label{BC 7} \tag{9}\]

Inspired by these studies, in this article, we examine the second order difference equation \[\label{DE 3} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = f(\zeta, y(\zeta)), \quad \zeta \in \mathbb{N}_{1}^{T}, \tag{10}\] associated with the boundary conditions either (2) or \[\label{BC 8} A y(0) + B y(T) = 0, \quad C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0, \tag{11}\] where \(\alpha\), \(\beta \geq 0\); \(A\), \(B\), \(C\), \(D \in \mathbb{R}\) with \(A^{2} + B^{2}\), \(C^{2} + D^{2} > 0\); \(T \in \mathbb{N}_{3}\); \(y: \mathbb{N}^{T + 1}_{0} \rightarrow \mathbb{R}\), and \(f : \mathbb{N}^{T}_{1} \times \mathbb{R} \rightarrow \mathbb{R}\) is a continuous function with respect to its second argument.

Clearly, the boundary conditions (4), (5), (7), (8), and (9) are particular cases of (11). Moreover, a suitable choice of \(\alpha\) and \(\beta\) yield the boundary value problems (10) – (7), (10) – (8), and (10) – (9) non-resonant.

In particular, the following Tables 1 and 2 illustrate that the boundary value problems considered in [110] as special cases of the boundary value problems (10) – (2) and (10) – (11):

Table 1 Special cases of the boundary value problem (10) – (2) studied in the literature
\(\alpha\) \(\beta\) \(A\) \(B\) \(C\) \(D\) \(\text{Reference}\)
\(0\) \(0\) \(\geq 0\) \(\geq 0\) \(\geq 0\) \(> 0\) [1]
\(\geq 0\) \(0\) \(\geq 0\) \(\geq 0\) \(\geq 0\) \(> 0\) [2]
\(\geq 0\) \(0\) \(\geq 0\) \(\geq 0\) \(\geq 0\) \(\geq 0\) [3]
Table 2 Special cases of the boundary value problem (10) – (11) studied in the literature
\(\alpha\) \(\beta\) \(A\) \(B\) \(C\) \(D\) \(\text{Reference}\)
\(> 0\) \(0\) \(1\) \(- 1\) \(1\) \(- 1\) [46]
\(0\) \(0\) \(1\) \(1\) \(1\) \(1\) [7]
\(0\) \(0\) \(1\) \(- 1\) \(1\) \(1\) [810]
\(0\) \(0\) \(1\) \(1\) \(1\) \(- 1\) [810]

Consequently, the results established in this article not only encompass the works of [110] as special cases, but also extend to a broad range of second-order difference equations through appropriate choices of the parameters \(\alpha\) and \(\beta\). Moreover, they accommodate a wide variety of boundary conditions via suitable selections of the constants \(A\), \(B\), \(C\), and \(D\).

We organize the present article as follows: §2 contains preliminaries on discrete calculus [11] and fixed point theory [12]. In §3, we obtain the expressions for the Green functions associated with the boundary value problem (10) – (2) and (10) – (11). We also derive a few important properties of the Green functions, which we will use to deduce the main results. In §4, we establish sufficient conditions for the existence of solutions to the boundary value problems (10) – (2) and (10) – (11). §5 presents two examples to show the applicability of the main results.

2. Preliminaries

In this section, we present the following fundamentals of discrete calculus [11] and fixed point theory [12], which we will use throughout the article. Denote by \(\mathbb{N}_{a} = \{a, a + 1, a + 2, \ldots\}\) and \(\mathbb{N}^{b}_{a} = \{a, a + 1, a + 2, \ldots, b\}\) for any real numbers \(a\) and \(b\) such that \(b – a \in \mathbb{N}_{1}\).

Definition 1. [11] The forward jump operator \(\sigma : \mathbb{N}_{a} \rightarrow \mathbb{N}_{a + 1}\) is defined by \[\sigma(\zeta) = \zeta + 1, \quad \zeta \in \mathbb{N}_{a}.\]

Definition 2. [11] Let \(y: \mathbb{N}^{b}_{a} \rightarrow \mathbb{R}\) and \(M \in \mathbb{N}_1\). The first-order forward difference of \(y\) is defined by \[\left(\Delta y\right)(\zeta) = y(\sigma(\zeta)) – y(\zeta), \quad \zeta \in \mathbb{N}^{b – 1}_{a},\] and the \(M^{th}\)-order forward difference of \(y\) is defined recursively by \[\left(\Delta^{M} y\right)(\zeta) = \left(\Delta \left(\Delta^{M – 1} y\right)\right)(\zeta), \quad \zeta \in \mathbb{N}^{b – M}_{a}.\]

Consider \[\mathcal{B}_{1} = \left\{y : \mathbb{N}^{T + 1}_{0} \rightarrow \mathbb{R} ~ | ~ A y(0) – B \left(\Delta y\right)(0) = 0, ~ C y(T) + D \left(\Delta y\right)(T) = 0\right\} \subseteq \mathbb{R}^{T + 2},\] and \[\mathcal{B}_{2} = \left\{y : \mathbb{N}^{T + 1}_{0} \rightarrow \mathbb{R} ~ | ~ A y(0) + B y(T) = 0, ~ C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0\right\} \subseteq \mathbb{R}^{T + 2}.\]

Clearly, \(\mathcal{B}_{1}\) and \(\mathcal{B}_{2}\) are Banach spaces equipped with the maximum norm defined by \[\left \|y\right\| = \max \limits_{\zeta \in \mathbb{N}^{T + 1}_{0}} \left|y(\zeta)\right|,\] for any \(y \in \mathcal{B}_{1}\) (or \(\mathcal{B}_{2}\)).

We apply the Leray–Schauder nonlinear alternative to establish sufficient conditions for the existence of solutions to the boundary value problem (10) – (11). We state this theorem as follows for convenience:

Theorem 1. [12] (Leray–Schauder Nonlinear Alternative) Let \(\mathcal{B} = \left(\mathcal{B}, \|\cdot\|\right)\) be a Banach space, \(C\) a closed, convex subset of \(\mathcal{B}\), \(U\) an open subset of \(C\) and \(0 \in U\). Suppose that \(T : \bar{U} \rightarrow C\) is a completely continuous map. Then, either

1. \(T\) has a fixed point in \(\bar{U}\); or

2. there exist a \(y \in \partial U\) and \(\lambda \in (0, 1)\) such that \(y = \lambda T y\).

3. Green functions & their properties

In this section, we obtain the expressions for the Green functions associated with the boundary value problems (10) – (2) and (10) – (11). We also derive a few important properties of the Green functions, which we will use to deduce the main results.

Denote by \[\begin{aligned} & \lambda = \frac{(2 + \alpha + \beta) + \sqrt{(\alpha + \beta)^{2} + 4 \beta}}{2}, \\ & \Lambda = \lambda – \frac{1 + \alpha}{\lambda} = \frac{\lambda^{2} – \alpha – 1}{\lambda}, \\ & g(x) = \lambda^{x} – \left(\frac{1 + \alpha}{\lambda}\right)^{x}, \quad x \in \mathbb{N}_{0}, \\ & h(x) = \lambda^{- x} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x}, \quad x \in \mathbb{N}_{0}, \\ & \omega(x) = g(x + 1) – (1 + \alpha) g(x), \quad x \in \mathbb{N}_{0}, \\ & w(x) = (1 + \alpha)^{-1} h(x) – h(x + 1), \quad x \in \mathbb{N}_{0}. \end{aligned}\]

Lemma 1. We have that

(1) \(\lambda \geq 1 + \alpha + \beta \geq 1 + \alpha\);

(2) \(\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha) = 0\);

(3) \(\lambda \geq 1\);

(4) \(\Lambda \geq 0\);

(5) \(g : \mathbb{N}_{0} \rightarrow \mathbb{R}\) is a nonnegative nondecreasing function;

(6) \(h : \mathbb{N}_{0} \rightarrow \mathbb{R}\) is a nonpositive nonincreasing function;

(7) \(\omega : \mathbb{N}_{0} \rightarrow \mathbb{R}\) is a nonnegative nondecreasing function;

(8) \(w : \mathbb{N}_{0} \rightarrow \mathbb{R}\) is a nonnegative nondecreasing function.

Proof. First, we prove (1). Consider \[\lambda = \frac{(2 + \alpha + \beta) + \sqrt{(\alpha + \beta)^{2} + 4 \beta}}{2} \geq \frac{(2 + \alpha + \beta) + (\alpha + \beta)}{2} = 1 + \alpha + \beta \geq 1 + \alpha.\]

The proof of (1) is complete. Next, we prove (2). The roots of the quadratic equation \(x^{2} – (2 + \alpha + \beta) x + (1 + \alpha) = 0\) are \[x = \frac{(2 + \alpha + \beta) \pm \sqrt{(2 + \alpha + \beta)^{2} – 4(1 + \alpha)}}{2} = \frac{(2 + \alpha + \beta) \pm \sqrt{(\alpha + \beta)^{2} + 4 \beta}}{2}.\]

Since \[\lambda = \frac{(2 + \alpha + \beta) + \sqrt{(\alpha + \beta)^{2} + 4 \beta}}{2},\] it satisfies the quadratic equation \(x^{2} – (2 + \alpha + \beta) x + (1 + \alpha) = 0\). The proof of (2) is complete. Now, we prove (3). From (1), we have \(\lambda \geq 1 + \alpha \geq 1\). The proof of (3) is complete. Next, we prove (4). From (1) and (3), we have \(\frac{1 + \alpha}{\lambda} \leq 1\) and \(\lambda \geq 1\) implying that \[\Lambda = \lambda – \frac{1 + \alpha}{\lambda} \geq 0.\]

The proof of (4) is complete. Now, we prove (5). Clearly, \(g(0) = 0\), \(g(1) = \Lambda \geq 0\), and \[g(2) = \lambda^{2} – \left(\frac{1 + \alpha}{\lambda}\right)^{2} = \Lambda \left(\lambda + \frac{1 + \alpha}{\lambda}\right) \geq 0.\]

For \(x \in \mathbb{N}_{0}\), consider \[\begin{aligned} \left(\Delta g\right)(x) & = g(x + 1) – g(x) \\ & = \left[\lambda^{x + 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{x + 1}\right] – \left[\lambda^{x} – \left(\frac{1 + \alpha}{\lambda}\right)^{x}\right] \\ & = \left[\lambda^{x + 1} – \lambda^{x}\right] – \left[\left(\frac{1 + \alpha}{\lambda}\right)^{x + 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{x}\right] \\ & = (\lambda – 1) \lambda^{x} – \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{x} \\ & = (\lambda – 1) \lambda^{x} + \left(1 – \frac{1 + \alpha}{\lambda}\right) \left(\frac{1 + \alpha}{\lambda}\right)^{x}. \end{aligned}\]

Clearly, \(\lambda – 1 \geq 0\) and \(\displaystyle 1 – \frac{1 + \alpha}{\lambda} \geq 0\). Since \(\lambda^{x} > 0\) and \(\left(\frac{1 + \alpha}{\lambda}\right)^{x} > 0\) for all \(x \in \mathbb{N}_{0}\), we obtain that \[\left(\Delta g\right)(x) \geq 0, \quad x \in \mathbb{N}_{0},\] implying that \(g\) is a nonnegative nondecreasing function. Next, we prove (6). Clearly, \(h(0) = 0\), and \[h(1) = \frac{1}{\lambda} – \frac{\lambda}{1 + \alpha} = \frac{1 + \alpha – \lambda}{\lambda (1 + \alpha)} \leq 0.\]

For \(x \in \mathbb{N}_{0}\), consider \[\begin{aligned} \left(\Delta h\right)(x) & = h(x + 1) – h(x) \\ & = \left[\lambda^{- x – 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1}\right] – \left[\lambda^{- x} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x}\right] \\ & = \left[\lambda^{- x – 1} – \lambda^{- x}\right] – \left[\left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x}\right] \\ & = (1 – \lambda) \lambda^{- x – 1} – \left(1 – \frac{1 + \alpha}{\lambda}\right) \left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1} \\ & = – (\lambda – 1) \lambda^{- x – 1} – \left(1 – \frac{1 + \alpha}{\lambda}\right) \left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1}. \end{aligned}\]

Clearly, \(\lambda – 1 \geq 0\) and \(\displaystyle 1 – \frac{1 + \alpha}{\lambda} \geq 0\). Since \(\lambda^{- x – 1} > 0\) and \(\left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1} > 0\) for all \(x \in \mathbb{N}_{0}\), we obtain that \[\left(\Delta h\right)(x) \leq 0, \quad x \in \mathbb{N}_{0},\] implying that \(h\) is a nonpositive nonincreasing function. Now, we prove (7). Clearly, \[\omega(0) = g(1) – (1 + \alpha) g(0) = \Lambda \geq 0.\]

For \(x \in \mathbb{N}_{0}\), consider \[\begin{aligned} \left(\Delta \omega\right)(x) & = \omega(x + 1) – \omega(x) \\ & = \left[g(x + 2) – (1 + \alpha) g(x + 1)\right] – \left[g(x + 1) – (1 + \alpha) g(x)\right] \\ & = g(x + 2) – (2 + \alpha) g(x + 1) + (1 + \alpha) g(x) \\ & = \left[\lambda^{x + 2} – \left(\frac{1 + \alpha}{\lambda}\right)^{x + 2}\right] – (2 + \alpha) \left[\lambda^{x + 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{x + 1}\right] \\ & \quad + (1 + \alpha) \left[\lambda^{x} – \left(\frac{1 + \alpha}{\lambda}\right)^{x}\right] \\ & = \left[\lambda^{x + 2} – (2 + \alpha) \lambda^{x + 1} + (1 + \alpha)\lambda^{x}\right] \\ & \quad – \left[\left(\frac{1 + \alpha}{\lambda}\right)^{x + 2} – (2 + \alpha)\left(\frac{1 + \alpha}{\lambda}\right)^{x + 1} + (1 + \alpha)\left(\frac{1 + \alpha}{\lambda}\right)^{x}\right] \\ & = \lambda^{x} \left[\lambda^{2} – (2 + \alpha) \lambda + (1 + \alpha)\right] \\ & \quad – \frac{1}{\lambda} \left(\frac{1 + \alpha}{\lambda}\right)^{x + 1} \left[\lambda^{2} – (2 + \alpha) \lambda + (1 + \alpha)\right] \\ & = \frac{1}{\lambda} \left[\lambda^{2} – (2 + \alpha) \lambda + (1 + \alpha)\right] g(x + 1) \\ & \geq \frac{1}{\lambda} \left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha)\right] g(x + 1) \\ & = 0, \end{aligned}\] implying that \(\omega\) is a nonnegative nondecreasing function. Finally, we prove (8). Clearly, \[w(0) = (1 + \alpha)^{-1} h(0) – h(1) = \frac{\lambda – \alpha – 1}{\lambda (1 + \alpha)} \geq 0.\]

For \(x \in \mathbb{N}_{0}\), consider \[\begin{aligned} \left(\Delta w\right)(x) & = w(x + 1) – w(x) \\ & = \left[(1 + \alpha)^{-1} h(x + 1) – h(x + 2)\right] – \left[(1 + \alpha)^{-1} h(x) – h(x + 1)\right] \\ & = – h(x + 2) + \left(1 + \frac{1}{1 + \alpha}\right) h(x + 1) – \frac{h(x)}{(1 + \alpha)} \\ & = – \left[\lambda^{- x – 2} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x – 2}\right] + \left(1 + \frac{1}{1 + \alpha}\right) \left[\lambda^{- x – 1} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1}\right] \\ & \quad – \frac{1}{(1 + \alpha)}\left[\lambda^{- x} – \left(\frac{1 + \alpha}{\lambda}\right)^{- x}\right] \\ & = – \left[\lambda^{- x – 2} – \left(1 + \frac{1}{1 + \alpha}\right) \lambda^{- x – 1} + \frac{1}{(1 + \alpha)}\lambda^{- x}\right] \\ & \quad + \left[\left(\frac{1 + \alpha}{\lambda}\right)^{- x – 2} – \left(1 + \frac{1}{1 + \alpha}\right)\left(\frac{1 + \alpha}{\lambda}\right)^{- x – 1} + \frac{1}{(1 + \alpha)}\left(\frac{1 + \alpha}{\lambda}\right)^{- x}\right] \\ & = – \frac{1}{\lambda (1 + \alpha)} \left[\lambda^{2} – (2 + \alpha) \lambda + (1 + \alpha)\right] h(x + 1) \\ & \geq – \frac{1}{\lambda (1 + \alpha)} \left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha)\right] h(x + 1) = 0, \end{aligned}\] implying that \(w\) is a nonnegative nondecreasing function. The proof is complete. ◻

Remark 1. To the end of this article, we assume that \(\Lambda \neq 0\).

Lemma 2. A general solution of the homogeneous second-order difference equation \[\label{BVP L 0} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = 0, \quad \zeta \in \mathbb{N}_{1}, \tag{12}\] is given by \[\label{GS L 0} y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}, \quad \zeta \in \mathbb{N}_{0}, \tag{13}\] where \(C_{1}\) and \(C_{2}\) are arbitrary constants.

Proof. In order to prove that (13) is a general solution of (12), it is enough to show that \(\displaystyle y_{1}(\zeta) = \lambda^{\zeta}\) and \(\displaystyle y_{2}(\zeta) = \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}\) are two linearly independent solutions of (12) on \(\mathbb{N}_{0}\). For \(\zeta \in \mathbb{N}_{1}\), consider \[\begin{aligned} \left(\Delta^2 y_{1}\right)(\zeta – 1) – \alpha \left(\Delta y_{1}\right)(\zeta – 1) – \beta y_{1}(\zeta) & = y_{1}(\zeta + 1) – (2 + \alpha + \beta) y_{1}(\zeta) + (1 + \alpha) y_{1}(\zeta – 1) \\ & = \lambda^{\zeta + 1} – (2 + \alpha + \beta) \lambda^{\zeta} + (1 + \alpha) \lambda^{\zeta – 1} \\ & = \lambda^{\zeta – 1} \left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha)\right] \\ & = 0, \end{aligned}\] and \[\begin{aligned} \left(\Delta^2 y_{2}\right)(\zeta – 1) – \alpha \left(\Delta y_{2}\right)(\zeta – 1) – \beta y_{2}(\zeta) & = y_{2}(\zeta + 1) – (2 + \alpha + \beta) y_{2}(\zeta) + (1 + \alpha) y_{2}(\zeta – 1) \\ & = \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta + 1} – (2 + \alpha + \beta) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + (1 + \alpha) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – 1} \\ & = \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – 1} \left[\left(\frac{1 + \alpha}{\lambda}\right)^{2} – (2 + \alpha + \beta) \left(\frac{1 + \alpha}{\lambda}\right) + (1 + \alpha)\right] \\ & = \frac{1}{\lambda} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} \left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha)\right] \\ & = 0, \end{aligned}\] implying that \(\displaystyle y_{1}(\zeta) = \lambda^{\zeta}\) and \(\displaystyle y_{2}(\zeta) = \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}\) are solutions of (12) on \(\mathbb{N}_{0}\). Further, the Wronskian of \(y_{1}\) and \(y_{2}\) is given by \[\begin{aligned} W\left(y_{1}, y_{2}\right)(\zeta) & = \begin{vmatrix} y_{1}(\zeta) & y_{2}(\zeta) \\ \left(\Delta y_{1}\right)(\zeta) & \left(\Delta y_{2}\right)(\zeta) \end{vmatrix} \\ & = \begin{vmatrix} \lambda^{\zeta} & \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} \\ (\lambda – 1) \lambda^{\zeta} & \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} \end{vmatrix} \\ & = \left(\frac{1 + \alpha}{\lambda} – 1\right) \lambda^{\zeta} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} – (\lambda – 1) \lambda^{\zeta} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} \\ & = – \Lambda (1 + \alpha)^{\zeta} \neq 0, \quad \zeta \in \mathbb{N}_{1}, \end{aligned}\] implying that \(\displaystyle y_{1}(\zeta) = \lambda^{\zeta}\) and \(\displaystyle y_{2}(\zeta) = \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}\) are linearly independent on \(\mathbb{N}_{0}\). The proof is complete. ◻

Lemma 3. Assume \(k\) is a real-valued function defined on a discrete set \(\mathbb{N}_{1}\). A general solution of the nonhomogeneous second-order difference equation \[\label{BVP L 00} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = k(\zeta), \quad \zeta \in \mathbb{N}_{1}, \tag{14}\] is given by \[ y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta – 1}_{s = 1}\left[\lambda^{\zeta – s} – \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s}\right]k(s) \label{GS L 1} \ \tag{15}\] \[ = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s), \label{GS L 1A} \tag{16}\] for \(\zeta \in \mathbb{N}_{0}\). Here \(C_{1}\), \(C_{2}\) are arbitrary constants and \[g(\zeta – s) = \lambda^{\zeta – s} – \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s}, \quad s \in \mathbb{N}^{\zeta – 1}_{1}.\]

Proof. In view of Lemma 2, it is enough to show that \[v(\zeta) = \frac{1}{\Lambda}\sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s), \quad \zeta \in \mathbb{N}_{0},\] is a particular solution of (14). For this purpose, we claim that \[\label{BVP L 11} \left(\Delta^2 v\right)(\zeta – 1) – \alpha \left(\Delta v\right)(\zeta – 1) – \beta v(\zeta) = k(\zeta), \quad \zeta \in \mathbb{N}_{1}. \tag{17}\]

To see this, for \(\zeta \in \mathbb{N}_{1}\), consider \[\begin{aligned} & \left(\Delta^2 v\right)(\zeta – 1) – \alpha \left(\Delta v\right)(\zeta – 1) – \beta v(\zeta) \\ & = v(\zeta + 1) – (2 + \alpha + \beta) v(\zeta) + (1 + \alpha) v(\zeta – 1) \\ & = \frac{1}{\Lambda} \sum^{\zeta}_{s = 1}g(\zeta – s + 1)k(s) – \left(\frac{2 + \alpha + \beta}{\Lambda}\right)\sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s) + \left(\frac{1 + \alpha}{\Lambda}\right) \sum^{\zeta – 2}_{s = 1}g(\zeta – s – 1)k(s) \\ & = \frac{1}{\Lambda} \left[g(1) k(\zeta) + g(2) k(\zeta – 1)\right] – \left(\frac{2 + \alpha + \beta}{\Lambda}\right) g(1) k(\zeta – 1) \\ & \quad + \frac{1}{\Lambda} \sum^{\zeta – 2}_{s = 1}\left[\lambda^{\zeta – s + 1} – (2 + \alpha + \beta) \lambda^{\zeta – s} + (1 + \alpha) \lambda^{\zeta – s – 1} \right]k(s) \\ & \quad – \frac{1}{\Lambda} \sum^{\zeta – 2}_{s = 1}\Bigg[\left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s + 1} – (2 + \alpha + \beta) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s} + (1 + \alpha) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s – 1}\Bigg]k(s) \\ & = k(\zeta) + \left(\frac{\lambda^{2} + \alpha + 1}{\lambda}\right) k(\zeta – 1) – (2 + \alpha + \beta) k(\zeta – 1) \\ & \quad + \frac{1}{\Lambda} \lambda^{\zeta – s – 1} \sum^{\zeta – 2}_{s = 1}\left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha) \right]k(s) \\ & \quad – \frac{1}{\lambda \Lambda} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s} \sum^{\zeta – 2}_{s = 1}\left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha) \right]k(s) \\ & = k(\zeta) + \frac{1}{\lambda}\left[\lambda^{2} – (2 + \alpha + \beta) \lambda + (1 + \alpha) \right]k(\zeta – 1) \\ & = k(\zeta), \end{aligned}\] implying that (17) holds. The proof is complete. ◻

Remark 2. Consider (15). Then, for \(\zeta \in \mathbb{N}_{0}\), \[\begin{aligned} \nonumber \left(\Delta y\right)(\zeta) & = y(\zeta + 1) – y(\zeta) \\ \nonumber & = C_{1} \lambda^{\zeta + 1} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta + 1} + \frac{1}{\Lambda} \sum^{\zeta}_{s = 1}g(\zeta – s + 1)k(s) – \left[C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda} \sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s)\right] \\ \nonumber & = C_{1} (\lambda – 1) \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda} g(1) k(\zeta) + \frac{1}{\Lambda} \sum^{\zeta – 1}_{s = 1}\left[g(\zeta – s + 1) – g(\zeta – s)\right]k(s) \\ \nonumber & = C_{1} (\lambda – 1) \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda} \sum^{\zeta}_{s = 1}\left[g(\zeta – s + 1) – g(\zeta – s)\right]k(s) \\ \nonumber & = C_{1} (\lambda – 1) \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta}_{s = 1}\left[(\lambda – 1)\lambda^{\zeta – s} – \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta – s}\right]k(s). \end{aligned}\]

Using Lemmas 2 and 3, we obtain the expressions for the Green functions associated with the linear boundary value problems \[\label{BVP L 1} \begin{cases} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = k(\zeta), \quad \zeta \in \mathbb{N}_{1}^{T}, \\ A y(0) – B \left(\Delta y\right)(0) = 0, \quad C y(T) + D \left(\Delta y\right)(T) = 0, \end{cases} \tag{18}\] and \[\label{BVP L 2} \begin{cases} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = k(\zeta), \quad \zeta \in \mathbb{N}_{1}^{T}, \\ A y(0) + B y(T) = 0, \quad C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0. \end{cases} \tag{19}\]

For convenience, we use the following notations: \[\begin{aligned} & E_{1} = A + B \lambda^{T}, \\ & E_{2} = A + B \left(\frac{1 + \alpha}{\lambda}\right)^{T}, \\ & E_{3} = C + D \lambda^{T}, \\ & E_{4} = C + D \left(\frac{1 + \alpha}{\lambda}\right)^{T}, \\ & F_{1} = A – B (\lambda – 1), \\ & F_{2} = A – B \left(\frac{1 + \alpha}{\lambda} – 1\right), \\ & F_{3} = C \lambda^{T} + D (\lambda – 1) \lambda^{T}, \\ & F_{4} = C \left(\frac{1 + \alpha}{\lambda}\right)^{T} + D \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{T}, \\ – &\xi = \Lambda \left[AC + BD (1 + \lambda)^{T}\right] + AD \left[g(T + 1) – g(T)\right] + BC \omega(T – 1), \\ -& \chi = AC g(T) + AD \left[g(T + 1) – g(T)\right] + BC \omega(T – 1) + BD (\lambda – 1)\left(1 – \frac{1 + \alpha}{\lambda}\right) g(T). \end{aligned}\]

Remark 3. To the end of this article, we assume that \(\xi \neq 0\) and \(\chi \neq 0\).

Lemma 4. Assume \(k\) is a real-valued function defined on a discrete finite set \(\mathbb{N}^{T}_{1}\). Then, the linear boundary value problem (19) has a unique solution given in the form \[\label{H Sol} y(\zeta) = \sum^{T}_{s = 1}\mathcal{H}(\zeta, s) k(s), \quad \zeta \in \mathbb{N}^{T + 1}_{0}, \tag{20}\] where \[\label{Green 2} \mathcal{H}(\zeta, s) = \frac{1}{\Lambda} \begin{cases} \mathcal{H}_{1}(\zeta, s), \quad s \in \mathbb{N}^{\zeta}_{1}, \\ \mathcal{H}_{2}(\zeta, s), \quad s \in \mathbb{N}^{T}_{\zeta}, \end{cases} \tag{21}\] \[\begin{aligned} \xi \mathcal{H}_{2}(\zeta, s) =& AD g(\zeta) \left[g(T – s + 1) – g(T – s)\right] \\ & + BC g(T – s) \omega(\zeta – 1) + BD \Lambda (1 + \alpha)^{T} h(s – \zeta), \quad (\zeta, s) \in \mathbb{N}^{T + 1}_{0} \times \mathbb{N}^{\zeta}_{1}, \end{aligned}\] and \[\begin{aligned} \xi \mathcal{H}_{1}(\zeta, s) =& AD g(\zeta) \left[g(T – s + 1) – g(T – s)\right] + BC g(T – s) \omega(\zeta – 1) + BD \Lambda (1 + \alpha)^{T} g(\zeta – s) \\& + \xi g(\zeta – s), \quad (\zeta, s) \in \mathbb{N}^{T + 1}_{0} \times \mathbb{N}^{T}_{\zeta}. \end{aligned}\]

Proof. From Lemma 3, a general solution of the second-order difference equation in (19) is given by \[\label{G 1} y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s), \quad \zeta \in \mathbb{N}^{T + 1}_{0}, \tag{22}\] where \(C_{1}\) and \(C_{2}\) are arbitrary constants. It follows from Remark 2 that \[\begin{aligned} \label{G 2} \left(\Delta y\right)(\zeta) = C_{1} (\lambda – 1) \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta}_{s = 1}\left[g(\zeta – s + 1) – g(\zeta – s)\right]k(s), \quad \zeta \in \mathbb{N}^{T}_{0}. \end{aligned} \tag{23}\]

Using \(A y(0) + B y(T) = 0\) in (22), we obtain \[\label{G 3} C_{1} E_{1} + C_{2} E_{2} = – \frac{B}{\Lambda}\sum^{T – 1}_{s = 1}g(T – s)k(s). \tag{24}\]

Using \(C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0\) in (23), we obtain \[\label{G 4} (\lambda – 1) C_{1} E_{3} + \left(\frac{1 + \alpha}{\lambda} – 1\right) C_{2} E_{4} = – \frac{D}{\Lambda}\sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s)\right]k(s). \tag{25}\]

From (24) and (25), we have \[\label{G 5} C_{1} = – \frac{1}{\Lambda \xi} \sum^{T}_{s = 1}\left[\left(\frac{1 + \alpha}{\lambda} – 1\right) B E_{4} g(T – s) – D E_{2} \left[g(T – s + 1) – g(T – s)\right]\right]k(s), \tag{26}\] and \[\label{G 6} C_{2} = – \frac{1}{\Lambda \xi} \sum^{T}_{s = 1}\left[- B (\lambda – 1)E_{3} g(T – s) + D E_{1} \left[g(T – s + 1) – g(T – s)\right] \right]k(s). \tag{27}\]

Substituting the equalities (26) and (27) in (22), we obtain (20). The proof is complete. ◻

Remark 4. Imposing the boundary conditions \(A y(0) + B y(T) = 0\) and \(C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0\) on the general solution \(y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}\) of the associated homogeneous second-order difference equation \[\left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = 0, \quad \zeta \in \mathbb{N}_{1}^{T},\] yields a linear algebraic system for the constants \(C_{1}\) and \(C_{2}\): \[\begin{pmatrix} E_{1} & E_{2} \\ (\lambda – 1) E_{3} & \left(\frac{1 + \alpha}{\lambda} – 1\right) E_{4} \end{pmatrix} \begin{pmatrix} C_{1} \\ C_{2} \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}.\]

The determinant of the co-efficient matrix of this system is given by \(- \xi\). The condition \(\xi \neq 0\) guarantees that the homogeneous boundary value problem \[\begin{cases} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = 0, \quad \zeta \in \mathbb{N}_{1}^{T}, \\ A y(0) + B y(T) = 0, \quad C \left(\Delta y\right)(0) + D \left(\Delta y\right)(T) = 0. \end{cases}\] has only the trivial solution. Consequently, the linear boundary value problem (19) has a unique solution given in the form (20).

Lemma 5. Assume \(k\) is a real-valued function defined on a discrete finite set \(\mathbb{N}^{T}_{1}\). Then, the linear boundary value problem (18) has a unique solution given in the form \[\label{G Sol} y(\zeta) = \sum^{T}_{s = 1}\mathcal{G}(\zeta, s) k(s), \quad \zeta \in \mathbb{N}^{T + 1}_{0}, \tag{28}\] where \[\label{Green 1} \mathcal{G}(\zeta, s) = \frac{1}{\Lambda} \begin{cases} \mathcal{G}_{1}(\zeta, s), \quad s \in \mathbb{N}^{\zeta}_{1}, \\ \mathcal{G}_{2}(\zeta, s), \quad s \in \mathbb{N}^{T}_{\zeta}, \end{cases} \tag{29}\] \[\begin{aligned} \chi \mathcal{G}_{2}(\zeta, s) = \left[C g(T – s) + D \left(g(T – s + 1) – g(T – s)\right)\right] \times \left[A g(\zeta) + B \omega(\zeta – 1)\right], \quad (\zeta, s) \in \mathbb{N}^{T + 1}_{0} \times \mathbb{N}^{\zeta}_{1}, \end{aligned} \tag{30}\] and \[\mathcal{G}_{1}(\zeta, s) = \mathcal{G}_{2}(\zeta, s) + g(\zeta – s), \quad (\zeta, s) \in \mathbb{N}^{T + 1}_{0} \times \mathbb{N}^{T}_{\zeta}. \tag{31}\]

Proof. From Lemma 3, a general solution of the second-order difference equation in (18) is given by \[\label{H 1} y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta – 1}_{s = 1}g(\zeta – s)k(s), \quad \zeta \in \mathbb{N}^{T + 1}_{0}, \tag{32}\] where \(C_{1}\) and \(C_{2}\) are arbitrary constants. It follows from Remark 2 that \[\begin{aligned} \label{H 2} \left(\Delta y\right)(\zeta) = C_{1} (\lambda – 1) \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda} – 1\right) \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta} + \frac{1}{\Lambda}\sum^{\zeta}_{s = 1}\left[g(\zeta – s + 1) – g(\zeta – s)\right]k(s), \quad \zeta \in \mathbb{N}^{T}_{0}. \end{aligned} \tag{33}\]

Using \(A y(0) – B \left(\Delta y\right)(0) = 0\) in (32) – (33), we obtain \[\label{H 3} C_{1} F_{1} + C_{2} F_{2} = 0. \tag{34}\]

Using \(C y(T) + D \left(\Delta y\right)(T) = 0\) in (32) – (33), we obtain \[\label{H 4} C_{1} F_{3} + C_{2} F_{4} = – \frac{1}{\Lambda}\sum^{T}_{s = 1}\left[C g(T – s) + D \left(g(T – s + 1) – g(T – s)\right)\right]k(s). \tag{35}\]

From (34) and (35), we have \[\label{H 5} C_{1} = \frac{F_{2}}{\Lambda \chi}\sum^{T}_{s = 1}\left[C g(T – s) + D \left(g(T – s + 1) – g(T – s)\right)\right]k(s), \tag{36}\] and \[\label{H 6} C_{2} = – \frac{F_{1}}{\Lambda \chi}\sum^{T}_{s = 1}\left[C g(T – s) + D \left(g(T – s + 1) – g(T – s)\right)\right]k(s). \tag{37}\]

Substituting the equalities (36) and (37) in (32), we obtain (28). The proof is complete. ◻

Remark 5. Imposing the boundary conditions \(A y(0) – B \left(\Delta y\right)(0) = 0\) and \(C y(T) + D \left(\Delta y\right)(T) = 0\) on the general solution \(y(\zeta) = C_{1} \lambda^{\zeta} + C_{2} \left(\frac{1 + \alpha}{\lambda}\right)^{\zeta}\) of the associated homogeneous second-order difference equation \[\left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = 0, \quad \zeta \in \mathbb{N}_{1}^{T},\] yields a linear algebraic system for the constants \(C_{1}\) and \(C_{2}\): \[\begin{pmatrix} F_{1} & F_{2} \\ F_{3} & F_{4} \end{pmatrix} \begin{pmatrix} C_{1} \\ C_{2} \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}.\]

The determinant of the co-efficient matrix of this system is given by \(- \chi\). The condition \(\chi \neq 0\) guarantees that the homogeneous boundary value problem \[\begin{cases} \left(\Delta^2 y\right)(\zeta – 1) – \alpha \left(\Delta y\right)(\zeta – 1) – \beta y(\zeta) = 0, \quad \zeta \in \mathbb{N}_{1}^{T}, \\ A y(0) – B \left(\Delta y\right)(0) = 0, \quad C y(T) + D \left(\Delta y\right)(T) = 0, \end{cases}\] has only the trivial solution. Consequently, the linear boundary value problem (18) has a unique solution given in the form (28).

Lemma 6. The Green’s function \(\mathcal{H}(\zeta, s)\) given by (21) satisfies the following property: \[\sum^{T}_{s = 1}\left|\mathcal{H}(\zeta, s)\right| \leq \Upsilon, \quad \zeta \in \mathbb{N}^{T + 1}_{0},\] where \[\begin{aligned} \Upsilon & = \frac{1}{\left|\xi\right|} \max_{\zeta \in \mathbb{N}_{0}^{T + 1}} \Bigg(\Bigg[|A| |D| g(\zeta) \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s)\right] \\ & \quad + |B| |C| \omega(\zeta – 1) \sum^{T}_{s = 1} g(T – s) + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{T}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg). \end{aligned}\]

Proof. For \(\zeta \in \mathbb{N}^{T + 1}_{0}\), consider \[\begin{aligned} \sum^{T}_{s = 1}\left|\mathcal{H}(\zeta, s)\right| & = \sum^{T}_{s = \zeta}\left|\mathcal{H}_{2}(\zeta, s)\right| + \sum^{\zeta}_{s = 1}\left|\mathcal{H}_{1}(\zeta, s)\right| \\ & = \frac{1}{\left|\xi\right|} \sum^{T}_{s = \zeta}\left|\xi \mathcal{H}_{2}(\zeta, s)\right| + \frac{1}{\left|\xi\right|} \sum^{\zeta}_{s = 1}\left|\xi \mathcal{H}_{1}(\zeta, s)\right| \\ & \leq \frac{1}{\left|\xi\right|} \sum^{T}_{s = 1}\Bigg|AD g(\zeta) \left[g(T – s + 1) – g(T – s)\right] + BC g(T – s) \omega(\zeta – 1)\Bigg| \\ & \quad + \frac{1}{\left|\xi\right|} \sum^{\zeta}_{s = 1} \Bigg|BD \Lambda (1 + \alpha)^{T} g(\zeta – s)\Bigg| + \frac{1}{\left|\xi\right|} \sum^{T}_{s = \zeta} \Bigg|BD \Lambda (1 + \alpha)^{T} h(s – \zeta)\Bigg| + \sum^{\zeta}_{s = 1}g(\zeta – s) \\ & \leq \frac{1}{\left|\xi\right|} \Bigg[|A| |D| g(\zeta) \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s)\right] \\ & \quad + |B| |C| \omega(\zeta – 1) \sum^{T}_{s = 1} g(T – s) + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{T}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s) \\ & \leq \Upsilon. \end{aligned}\]

The proof is complete. ◻

Lemma 7. The Green’s function \(\mathcal{G}(\zeta, s)\) given by (29) satisfies the following property: \[\sum^{T}_{s = 1}\left|\mathcal{G}(\zeta, s)\right| \leq \Theta, \quad \zeta \in \mathbb{N}^{T + 1}_{0},\] where \[\begin{aligned} \Theta & = \max_{\zeta \in \mathbb{N}_{0}^{T + 1}} \Bigg(\frac{\left[\left|A\right| g(\zeta) + \left|B\right| \omega(\zeta – 1)\right]}{\left|\chi\right|} \times \Bigg[|C| \sum^{T}_{s = 1} g(T – s) + |D| \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s) \right]\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg). \end{aligned}\]

Proof. For \(\zeta \in \mathbb{N}^{T + 1}_{0}\), consider \[\begin{aligned} \sum^{T}_{s = 1}\left|\mathcal{G}(\zeta, s)\right| & = \sum^{T}_{s = \zeta}\left|\mathcal{G}_{2}(\zeta, s)\right| + \sum^{\zeta}_{s = 1}\left|\mathcal{G}_{1}(\zeta, s)\right| \\ & = \frac{1}{\left|\chi\right|} \sum^{T}_{s = \zeta}\left|\chi \mathcal{G}_{2}(\zeta, s)\right| + \frac{1}{\left|\chi\right|} \sum^{\zeta}_{s = 1}\left|\chi \mathcal{G}_{1}(\zeta, s)\right| \\ & \leq \frac{1}{\left|\chi\right|} \sum^{T}_{s = 1}\Bigg|\left[C g(T – s) + D \left(g(T – s + 1) – g(T – s)\right)\right] \left[A g(\zeta) + B \omega(\zeta – 1)\right]\Bigg| + \sum^{\zeta}_{s = 1}g(\zeta – s) \\ & \leq \frac{\left[\left|A\right| g(\zeta) + \left|B\right| \omega(\zeta – 1)\right]}{\left|\chi\right|} \times \Bigg[|C| \sum^{T}_{s = 1} g(T – s) + |D| \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s) \right]\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s) \\ & \leq \Theta. \end{aligned}\]

The proof is complete. ◻

4. Main Results

This section establishes sufficient conditions on the existence of solutions to the boundary value problems (10) – (2) and (10) – (11). Lemma 5 implies the equivalence between the solutions of (10) – (2) and the solutions of the summation equation \[y(\zeta) = \sum^{T}_{s = 1}\mathcal{G}(\zeta, s) f(s, y(s)), \quad \zeta \in \mathbb{N}^{T + 1}_{0}.\]

Similarly, Lemma 4 implies the equivalence between the solutions of (10) – (11) and the solutions of the summation equation \[y(\zeta) = \sum^{T}_{s = 1}\mathcal{H}(\zeta, s) f(s, y(s)), \quad \zeta \in \mathbb{N}^{T + 1}_{0}.\]

Define the operators \(\mathfrak{R}: \mathcal{B}_{1} \rightarrow \mathcal{B}_{1}\) and \(\mathfrak{S} : \mathcal{B}_{2} \rightarrow \mathcal{B}_{2}\) by \[\left(\mathfrak{R} y\right)(\zeta) = \sum^{T}_{s = 1}\mathcal{G}(\zeta, s)f(s, y(s)), \quad \zeta \in \mathbb{N}_{0}^{T + 1},\] and \[\left(\mathfrak{S} y\right)(\zeta) = \sum^{T}_{s = 1}\mathcal{H}(\zeta, s)f(s, y(s)), \quad \zeta \in \mathbb{N}_{0}^{T + 1}.\]

It follows from Lemma 5 that \(y\) is a fixed point of \(\mathfrak{R}\) if, and only if, \(y\) is a solution of (10) – (2). Also, from Lemma 4, \(y\) is a fixed point of \(\mathfrak{S}\) if, and only if, \(y\) is a solution of (10) – (11). Let \[L = \left\{y \in \mathcal{B}_{1} : \left\|y\right\| \leq n\right\},\] and \[K = \left\{y \in \mathcal{B}_{2} : \left\|y\right\| \leq r\right\}.\]

Clearly, \(L\) and \(K\) are nonempty bounded closed convex subsets of the finite dimensional normed space \(\mathcal{B}_{1}\) and \(\mathcal{B}_{2}\), respectively.

Now, we apply Theorem 1 to discuss the existence of solutions to (10) – (2) and (10) – (11).

Theorem 2. Assume the following conditions hold:

(C1) There exist \(p: \mathbb{N}^{T}_{1} \rightarrow [0, \infty)\) and a nondecreasing function \(q : [0, \infty) \rightarrow [0, \infty)\) such that \[\left|f(\zeta, y)\right| \leq p(\zeta) q\left(\|y\|\right), \quad (\zeta, y) \in \mathbb{N}^{T}_{1} \times \mathbb{R}.\]

(C2) There exists \(N > 0\) such that \[\frac{N}{\Upsilon \Omega q\left(N\right)} > 1,\] where \[\Omega = \max_{\zeta \in \mathbb{N}^{T}_{1}} p(\zeta).\]

Then, (10) – (11) has a solution in \(\mathcal{B}_{2}\).

Proof. Since \(\mathbb{N}^{T + 1}_{0}\) is a discrete set, \(\mathcal{B}_{2}\) is finite dimensional normed space and \(\mathfrak{S}\) is a continuous operator, it follows immediately that \(\mathfrak{S}\) is completely continuous.

Next, we suppose \(y \in \mathcal{B}_{2}\) and that for some \(0 < \nu < 1\), \(y = \nu \mathfrak{S} y\). Then, for \(\zeta \in \mathbb{N}^{T + 1}_{0}\), and again by (C1), \[\begin{aligned} \left|y(\zeta)\right| & = \left|\nu (\mathfrak{S}y)(\zeta)\right| \\ & \leq \sum_{s = 1}^{T}\left|\mathcal{H}(\zeta, s)\right|\left|f(s, y(s))\right| \\ & \leq \sum_{s = 1}^{T}\left|\mathcal{H}(\zeta, s)\right|p(s) q\left(|y(s)|\right) \\ & \leq q\left(\|y\|\right) \sum_{s = 1}^{T}\left|\mathcal{H}(\zeta, s)\right|p(s) \\ & \leq \Upsilon \Omega q \left(\|y\|\right), \end{aligned}\] implying that \[\frac{\|y\|}{\Upsilon \Omega q \left(\|y\|\right)} \leq 1.\]

It follows from (C2) that \(\|y\| \neq N\). If we set \[U = \Big\{y \in \mathcal{B}_{2} : \|y\| < N\Big\},\] then the operator \(\mathfrak{S} : \bar{U} \rightarrow \mathcal{B}_{2}\) is completely continuous. From the choice of \(U\), there is no \(y \in \partial U\) such that \(y = \nu \mathfrak{S} y\) for some \(0 < \nu < 1\). It follows from Theorem 1 that \(\mathfrak{S}\) has a fixed point \(y_0 \in \bar{U}\), which is a desired solution of (10) – (11). ◻

Theorem 3. Assume the following conditions hold:

(C1) There exist \(p: \mathbb{N}^{T}_{1} \rightarrow [0, \infty)\) and a nondecreasing function \(q : [0, \infty) \rightarrow [0, \infty)\) such that \[\left|f(\zeta, y)\right| \leq p(\zeta) q\left(\|y\|\right), \quad (\zeta, y) \in \mathbb{N}^{T}_{1} \times \mathbb{R}.\]

(C3) There exists \(M > 0\) such that \[\frac{M}{\Theta \Omega q\left(M\right)} > 1,\] where \[\Omega = \max_{\zeta \in \mathbb{N}^{T}_{1}} p(\zeta).\]

Then, (10) – (2) has a solution in \(\mathcal{B}_{1}\).

Proof. Since \(\mathbb{N}^{T + 1}_{0}\) is a discrete set, \(\mathcal{B}_{1}\) is finite dimensional normed space and \(\mathfrak{R}\) is a continuous operator, it follows immediately that \(\mathfrak{R}\) is completely continuous.

Next, we suppose \(y \in \mathcal{B}_{1}\) and that for some \(0 < \nu < 1\), \(y = \nu \mathfrak{R} y\). Then, for \(\zeta \in \mathbb{N}^{T + 1}_{0}\), and again by (C1), \[\begin{aligned} \left|y(\zeta)\right| & = \left|\nu (\mathfrak{R}y)(\zeta)\right| \\ & \leq \sum_{s = 1}^{T}\left|\mathcal{G}(\zeta, s)\right|\left|f(s, y(s))\right| \\ & \leq \sum_{s = 1}^{T}\left|\mathcal{G}(\zeta, s)\right|p(s) q\left(|y(s)|\right) \\ & \leq q\left(\|y\|\right) \sum_{s = 1}^{T}\left|\mathcal{G}(\zeta, s)\right|p(s) \\ & \leq \Theta \Omega q \left(\|y\|\right), \end{aligned}\] implying that \[\frac{\|y\|}{\Theta \Omega q \left(\|y\|\right)} \leq 1.\]

It follows from (C2) that \(\|y\| \neq M\). If we set \[V = \Big\{y \in \mathcal{B}_{1} : \|y\| < M\Big\},\] then the operator \(\mathfrak{R} : \bar{V} \rightarrow \mathcal{B}_{1}\) is completely continuous. From the choice of \(V\), there is no \(y \in \partial V\) such that \(y = \nu \mathfrak{R} y\) for some \(0 < \nu < 1\). It follows from Theorem 1 that \(\mathfrak{R}\) has a fixed point \(y_0 \in \bar{V}\), which is a desired solution of (10) – (2). ◻

5. Examples

In this section, we provide two examples to demonstrate the applicability of established results.

Example 1. Consider (10) – (11) with \(A = B = C = D = 1\), \(T = 5\), \(\alpha = 1\), \(\beta = 0\) and \(f(\zeta, s) = \zeta s^2\). Clearly, \[\left|f(\zeta, s)\right| \leq p(\zeta) q\left(|s|\right), \quad (\zeta, s) \in \mathbb{N}^{5}_{1} \times \mathbb{R},\] where \[p(\zeta) = \zeta, \quad \zeta \in \mathbb{N}^{5}_{1},\] and \[q\left(|s|\right) = |s|^2 = s^2, \quad s \in \mathbb{R}.\]

Also, \(p : \mathbb{N}^{5}_{1} \rightarrow [0, \infty)\) and \(q : [0, \infty) \rightarrow [0, \infty)\) is a nondecreasing function. Thus, the assumption (C1) of Theorem 2 holds. Further, we have \[\Omega = \max_{\zeta \in \mathbb{N}^{5}_{1}} p(\zeta) = 5.\]

Now, we calculate \(\Upsilon\). We have \(\lambda = 2\), \(\Lambda = 1\), \[\begin{aligned} & g(x) = 2^{x} – 1, \quad x \in \mathbb{N}_{0}, \\ & h(x) = 2^{- x} – 1, \quad x \in \mathbb{N}_{0}, \\ & \omega(x) = 1, \quad x \in \mathbb{N}_{0}, \\ & w(x) = \frac{1}{2}, \quad x \in \mathbb{N}_{0}. \end{aligned}\]

Also, \[\begin{aligned} – \xi & = \Lambda \left[AC + BD (1 + \lambda)^{T}\right] + AD \left[g(T + 1) – g(T)\right] + BC \omega(T – 1) \\ & = 1 + 3^{5} + 2^{6} – 2^{5} + 1 \\ & = 277, \end{aligned}\] implying that \(\xi = – 277\). The corresponding Green’s function is given by \[\mathcal{H}(\zeta, s) = \begin{cases} \mathcal{H}_{1}(\zeta, s), \quad s \in \mathbb{N}^{\zeta}_{1}, \\ \mathcal{H}_{2}(\zeta, s), \quad s \in \mathbb{N}^{5}_{\zeta}, \end{cases}\] where \[\mathcal{H}_{2}(\zeta, s) = – \frac{1}{277}\left(2^{\zeta – s + 6} – 33\right),\] and \[\mathcal{H}_{1}(\zeta, s) = \frac{213}{277} 2^{\zeta – s} – \frac{244}{277}.\]

Consequently, \[\begin{aligned} \Upsilon & = \frac{1}{\left|\xi\right|} \max_{\zeta \in \mathbb{N}_{0}^{T + 1}}\Bigg(\Bigg[|A| |D| g(\zeta) \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s)\right] \\ & \quad + |B| |C| \omega(\zeta – 1) \sum^{T}_{s = 1} g(T – s) + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{T}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg) \\ & = \frac{1}{277} \max_{\zeta \in \mathbb{N}_{0}^{6}}\Bigg(\Bigg[g(\zeta) \sum^{5}_{s = 1}\left[g(6 – s) – g(5 – s)\right] + \omega(\zeta – 1) \sum^{5}_{s = 1} g(5 – s) + (32) \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + (32) \sum^{5}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg). \end{aligned}\]

We have \[\begin{aligned} \sum^{\zeta}_{s = 1}g(\zeta – s) & = \sum^{\zeta}_{s = 1} 2^{\zeta – s} – \sum^{\zeta}_{s = 1} 1 \\ & = \left[\frac{2^{\zeta} – 1}{2 – 1}\right] – \zeta \\ & = 2^{\zeta} – \zeta – 1,\\ \sum^{T}_{s = \zeta}\left|h(s – \zeta)\right| & = \sum^{5}_{s = \zeta}\left[1 – 2^{\zeta – s}\right] \\ & = \sum^{5}_{s = \zeta} 1 – \sum^{5}_{s = \zeta} 2^{\zeta – s} \\ & = 6 – \zeta – \sum^{5 – \zeta}_{s = 0} \left(\frac{1}{2}\right)^{s} \\ & = 6 – \zeta – \left[\frac{1 – \left(\frac{1}{2}\right)^{6 – \zeta}}{1 – \left(\frac{1}{2}\right)}\right] \\ & = 4 – \zeta + 2^{\zeta – 5},\\ \sum^{5}_{s = 1} g(5 – s) &= 2^{5} – 5 – 1 \\ &= 26, \end{aligned}\] and \[\begin{aligned} \sum^{5}_{s = 1}\left[g(6 – s) – g(5 – s)\right] & = \sum^{5}_{s = 1}g(6 – s) – \sum^{5}_{s = 1}g(5 – s) \\ & = \sum^{5}_{s = 1}\left[2^{6 – s} – 1\right] – 26 \\ & = \sum^{5}_{s = 1} 2^{6 – s} – \sum^{5}_{s = 1} 1 – 26 \\ & = 2 \left[\frac{2^{5} – 1}{2 – 1}\right] – 5 – 26 \\ & = 62 – 31 = 31. \end{aligned}\]

Then, \[\begin{aligned} \Upsilon & = \frac{1}{277} \max_{\zeta \in \mathbb{N}_{0}^{6}}\Bigg(\Bigg[31 \left(2^{\zeta} – 1\right) + 26 + (32) \left(2^{\zeta} – \zeta – 1\right) \\ & \quad + (32) \left(4 – \zeta + 2^{\zeta – 5}\right)\Bigg] + \left(2^{\zeta} – \zeta – 1\right)\Bigg) = \frac{3860}{277}. \end{aligned}\]

Then, there exists \(0 < N < \frac{277}{19300}\) such that \[\frac{N}{\Upsilon \Omega q\left(N\right)} > 1,\]implying that the assumption (C2) of Theorem 2 holds. Therefore, by Theorem 2, the boundary value problem (10) – (11) has a solution in \(\mathcal{B}_{2}\).

Example 2. Consider (10) – (11) with \(A = B = C = D = 1\), \(T = 5\), \(\alpha = 0\), \(\beta = \frac{1}{2}\) and \(f(\zeta, s) = \zeta s^2\). Clearly, \[\left|f(\zeta, s)\right| \leq p(\zeta) q\left(|s|\right), \quad (\zeta, s) \in \mathbb{N}^{5}_{1} \times \mathbb{R},\] where \[p(\zeta) = \zeta, \quad \zeta \in \mathbb{N}^{5}_{1},\] and \[q\left(|s|\right) = |s|^2 = s^2, \quad s \in \mathbb{R}.\]

Also, \(p : \mathbb{N}^{5}_{1} \rightarrow [0, \infty)\) and \(q : [0, \infty) \rightarrow [0, \infty)\) is a nondecreasing function. Thus, the assumption (C1) of Theorem 2 holds. Further, we have \[\Omega = \max_{\zeta \in \mathbb{N}^{5}_{1}} p(\zeta) = 5.\]

Now, we calculate \(\Upsilon\). We have \(\lambda = 2\), \(\Lambda = \frac{3}{2}\), \[\begin{aligned} & g(x) = 2^{x} – 2^{- x}, \quad x \in \mathbb{N}_{0}, \\ & h(x) = 2^{- x} – 2^{x}, \quad x \in \mathbb{N}_{0}, \\ & \omega(x) = 2^{x} + 2^{- x – 1}, \quad x \in \mathbb{N}_{0}, \\ & w(x) = 2^{x} + 2^{- x – 1}, \quad x \in \mathbb{N}_{0}. \end{aligned}\]

Also, \[\begin{aligned} – \xi & = \Lambda \left[AC + BD (1 + \lambda)^{T}\right] + AD \left[g(T + 1) – g(T)\right] + BC \omega(T – 1) \\ & = \frac{3}{2} \left[1 + 3^{5}\right] + \left[2^{6} – 2^{- 6}\right] – \left[2^{5} – 2^{- 5}\right] + 2^{4} + 2^{- 5} \\ & = \frac{26499}{64}, \end{aligned}\] implying that \(\xi = – \frac{26499}{64}\). The corresponding Green’s function is given by \[\mathcal{H}(\zeta, s) = \frac{2}{3} \begin{cases} \mathcal{H}_{1}(\zeta, s), \quad s \in \mathbb{N}^{\zeta}_{1}, \\ \mathcal{H}_{2}(\zeta, s), \quad s \in \mathbb{N}^{5}_{\zeta}, \end{cases}\] where \[\begin{aligned} \xi \mathcal{H}_{2}(\zeta, s) = \left(2^{\zeta} – 2^{- \zeta}\right) \left(2^{\zeta – s} + 2^{s – \zeta – 1}\right) + \left(2^{5 – s} – 2^{s – 5}\right)\left(2^{\zeta – 1} + 2^{- \zeta}\right) + 3 \left(2^{\zeta – s – 1} – 2^{s – \zeta – 1}\right), \end{aligned}\] and \[\begin{aligned} \nonumber \xi \mathcal{H}_{1}(\zeta, s) =& \left(2^{\zeta} – 2^{- \zeta}\right) \left(2^{\zeta – s} + 2^{s – \zeta – 1}\right) \\& + \left(2^{5 – s} – 2^{s – 5}\right)\left(2^{\zeta – 1} + 2^{- \zeta}\right) + 3 \left(2^{\zeta – s – 1} – 2^{s – \zeta – 1}\right) + \xi \left(2^{\zeta – s} – 2^{s – \zeta}\right). \end{aligned}\]

Consequently, \[\begin{aligned} \Upsilon & = \frac{1}{\left|\xi\right|} \max_{\zeta \in \mathbb{N}_{0}^{T + 1}}\Bigg(\Bigg[|A| |D| g(\zeta) \sum^{T}_{s = 1}\left[g(T – s + 1) – g(T – s)\right] \\ & \quad + |B| |C| \omega(\zeta – 1) \sum^{T}_{s = 1} g(T – s) + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + |B| |D| \Lambda (1 + \alpha)^{T} \sum^{T}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg) \\ & = \frac{64}{26499} \max_{\zeta \in \mathbb{N}_{0}^{6}}\Bigg(\Bigg[g(\zeta) \sum^{5}_{s = 1}\left[g(6 – s) – g(5 – s)\right] + \omega(\zeta – 1) \sum^{5}_{s = 1} g(5 – s) + (32) \sum^{\zeta}_{s = 1} g(\zeta – s) \\ & \quad + (32) \sum^{5}_{s = \zeta} \left|h(s – \zeta)\right|\Bigg] + \sum^{\zeta}_{s = 1}g(\zeta – s)\Bigg). \end{aligned}\]

We have \[\begin{aligned} \sum^{\zeta}_{s = 1}g(\zeta – s) & = \sum^{\zeta}_{s = 1} 2^{\zeta – s} – \sum^{\zeta}_{s = 1} \left(\frac{1}{2}\right)^{\zeta – s} \\ & = \left[\frac{2^{\zeta} – 1}{2 – 1}\right] – \left[\frac{1 – \left(\frac{1}{2}\right)^{\zeta}}{1 – \left(\frac{1}{2}\right)}\right] \\ & = 2^{\zeta} + 2^{1 – \zeta} – 3, \end{aligned}\] \[\begin{aligned} \sum^{T}_{s = \zeta}\left|h(s – \zeta)\right| & = \sum^{5}_{s = \zeta}\left[\left(\frac{1}{2}\right)^{\zeta – s} – 2^{\zeta – s}\right] \\ & = \sum^{5}_{s = \zeta} \left(\frac{1}{2}\right)^{\zeta – s} – \sum^{5}_{s = \zeta} 2^{\zeta – s} \\ & = \sum^{5 – \zeta}_{s = 0} 2^{s} – \sum^{5 – \zeta}_{s = 0} \left(\frac{1}{2}\right)^{s} \\ & = \left[\frac{2^{6 – \zeta} – 1}{2 – 1}\right] – \left[\frac{1 – \left(\frac{1}{2}\right)^{6 – \zeta}}{1 – \left(\frac{1}{2}\right)}\right] \\ & = 2^{6 – \zeta} + 2^{\zeta – 5} – 3, \\ \sum^{5}_{s = 1} g(5 – s) &= 2^{5} + 2^{- 4} – 3 = \frac{465}{16}, \end{aligned}\] and \[\begin{aligned} \sum^{5}_{s = 1}\left[g(6 – s) – g(5 – s)\right] & = \sum^{5}_{s = 1}g(6 – s) – \sum^{5}_{s = 1}g(5 – s) \\ & = \sum^{5}_{s = 1}\left[2^{6 – s} – \left(\frac{1}{2}\right)^{6 – s}\right] – \frac{465}{16} \\ & = \sum^{5}_{s = 1} 2^{6 – s} – \sum^{5}_{s = 1} \left(\frac{1}{2}\right)^{6 – s} – \frac{465}{16} \\ & = 2 \left[\frac{2^{5} – 1}{2 – 1}\right] – \frac{1}{2} \left[\frac{1 – \left(\frac{1}{2}\right)^{5}}{1 – \left(\frac{1}{2}\right)}\right] – \frac{465}{16} \\ & = \frac{1023}{32}. \end{aligned}\]

Then, \[\begin{aligned} \Upsilon & = \frac{64}{26499} \max_{\zeta \in \mathbb{N}_{0}^{6}}\Bigg(\Bigg[\frac{1023}{32} \left(2^{\zeta} – 2^{- \zeta}\right) + \frac{465}{16} \left(2^{\zeta – 1} + 2^{- \zeta}\right) + (32) \left(2^{\zeta} + 2^{1 – \zeta} – 3\right) \\ & \quad + (32) \left(2^{6 – \zeta} + 2^{\zeta – 5} – 3\right)\Bigg] + \left(2^{\zeta} + 2^{1 – \zeta} – 3\right)\Bigg) = \frac{319360}{26499}. \end{aligned}\]

Then, there exists \(0 < N < \frac{85}{5122}\) such that \[\frac{N}{\Upsilon \Omega q\left(N\right)} > 1,\]implying that the assumption (C2) of Theorem 2 holds. Therefore, by Theorem 2, the boundary value problem (10) – (11) has a solution in \(\mathcal{B}_{2}\).

6. Conclusion

In this article, we established sufficient conditions for the existence of solutions to the boundary value problems (10) – (2) and (10) – (11) using Leray–Schauder nonlinear alternative. The results established in this article not only subsume the works of [110] as particular cases, but also extend to a broad class of second-order difference equations through appropriate choices of the parameters \(\alpha\) and \(\beta\). Moreover, they accommodate a wide variety of boundary conditions via suitable selections of the constants \(A\), \(B\), \(C\), and \(D\). To the best of our knowledge, no existing studies in the literature address the existence of solutions for the discrete boundary value problems (10) – (2) and(10) – (11) by employing nonlinear analytical techniques such as fixed point theory, fixed point index theory, coincidence degree theory, critical point theory, or variational methods.

Conflicts of Interest

The author declares no conflict of interest.

Data Availability

No data is required for this research.

Funding Information

No funds, grants, or other support was received.

References

  1. Merdivenci, F. (1995). Two positive solutions of a boundary value problem for difference equations. Journal of Difference Equations and Applications, 1(3), 263-270.

  2. Atici, F. M. (2000). Existence of positive solutions of nonlinear discrete Sturm-Liouville problems. Mathematical and Computer Modelling, 32(5-6), 599-607.

  3. Aykut, N., & Guseinov, G. S. (2003). On positive solutions of boundary value problems for nonlinear second order difference equations. Turkish Journal of Mathematics, 27(4), 481-508.

  4. Atici, F. M., & Cabada, A. (2003). Existence and uniqueness results for discrete second-order periodic boundary value problems. Computers & Mathematics With Applications, 45(6-9), 1417-1427.

  5. Atici, F. M., & Guseinov, G. S. (1999). Positive periodic solutions for nonlinear difference equations with periodic coefficients. Journal of Mathematical Analysis and Applications, 232(1), 166-182.

  6. Atici, F. M., Cabada, A., & Otero-Espinar, V. (2003). Criteria for existence and nonexistence of positive solutions to a discrete periodic boundary value problem. Journal of Difference Equations and Applications, 9(9), 765-775.

  7. Lyons, J. W., & Neugebauer, J. T. (2015). A Difference Equation with Anti-Periodic Boundary Conditions. Dynamics of Continuous, Discrete and Impulsive Systems, Series A: Mathematical Analysis, 22(1), 47.

  8. Graef, J. R., Kong, L., & Wang, M. (2021). A variational framework for a second order discrete boundary value problem with mixed periodic boundary conditions. Results in Mathematics, 76(2), 98.

  9. Kong, L., & Wang, M. (2020). Multiple and particular solutions of a second order discrete boundary value problem with mixed periodic boundary conditions. Electronic Journal of Qualitative Theory of Differential Equations, 2020, 47.

  10. Kong, L., & Wang, M. (2020). Existence of solutions for a second order discrete boundary value problem with mixed periodic boundary conditions. Applied Mathematics Letters, 102, 106138.

  11. Bohner, M., & Peterson, A. (2001). Dynamic Equations on Time Scales: An Introduction With Applications. Springer Science & Business Media.

  12. Agarwal, R. P., Meehan, M., & O’regan, D. (2001). Fixed Point Theory and Applications (Vol. 141). Cambridge university press.