Bifurcation analysis of Leslie-Gower predator-prey system with harvesting and fear effect

1.   Introduction

Recently, fractional calculus methods became of great interest, because it is a powerful tool for calculating the derivation of multiples systems. These methods study real world phenomena in many areas of natural sciences including biomedical, radiography, biology, chemistry, and physics ^{[1,2,3,4,5,6,7]}. Abundant publications focus on the Caputo fractional derivative (CFD) and the Caputo-Hadamard derivative. Additionally, other generalization of the previous derivatives, such as $\Psi$-Caputo, study the existence of solutions to some FDEs (see ^{[8,9,10,11,12,13,14]}).

In general, an $m$-point fractional boundary problem involves a fractional differential equation with fractional boundary conditions that are specified at m different points on the boundary of a domain. The fractional derivative is defined using the Riemann-Liouville fractional derivative or the Caputo fractional derivative. Solving these types of problems can be challenging due to the non-local nature of fractional derivatives. However, there are various numerical and analytical methods available for solving such problems, including the spectral method, the finite difference method, the finite element method, and the homotopy analysis method. The applications of $m$-point fractional boundary problems can be found in various fields, including physics, engineering, finance, and biology. These problems are useful in modeling and analyzing phenomena that exhibit non-local behavior or involve memory effects (see ^{[15,16,17,18]}).

Pantograph equations are a set of differential equations that describe the motion of a pantograph, which is a mechanism used for copying and scaling drawings or diagrams. The equations are based on the assumption that the pantograph arms are rigid and do not deform during operation, we can simply say that see ^[19]. One important application of the pantograph equations is in the field of drafting and technical drawing. Before the advent of computer-aided design (CAD) software, pantographs were commonly used to produce scaled copies of drawings and diagrams. By adjusting the lengths of the arms and the position of the stylus, a pantograph can produce copies that are larger or smaller than the original ^[20], electrodynamics ^[21] and electrical pantograph of locomotive ^[22].

Many authors studied a huge number of positive solutions for nonlinear fractional BVP using fixed point theorems (FPTs) such as SFPT, Leggett-Williams and Guo-Krasnosel'skii (see ^[23,24]). Some studies addressed the sign-changing of solution of BVPs ^{[25,26,27,28,29]}.

In this work, we use Schauder's fixed point theorem (SFPT) to solve the semipostone multipoint $\Psi$-Caputo fractional pantograph problem

where $\lambda\in\left(0, \frac{\mathcal{\Im}\mathfrak{-}r}{\mathcal{\Im}}\right), \mathfrak{\mathcal{D}}_{r}^{\nu; \psi}$ is $\Psi$-Caputo fractional derivative ($\Psi$-CFD) of order $\nu$, $1 < \nu\leq2, \ \zeta_{i}\in\mathbb{R}^{+} \left(1\leq i\leq m-2\right)$ such that $0 < \Sigma_{i = 1}^{m-2}\zeta_{i} < 1, \ \mathfrak{\eta}_{i}\in(r, \mathcal{\Im})$, and $\mathcal{F}:[r, \mathcal{\Im}]\times\mathbb{R\times R}\rightarrow \mathbb{R}$.

The most important aspect of this research is to prove the existence of a positive solution of the above $m$-point FBVP. Note that in ^[30], the author considered a two-point BVP using Liouville-Caputo derivative.

The article is organized as follows. In the next section, we provide some basic definitions and arguments pertinent to fractional calculus (FC). Section $3$ is devoted to proving the the main result and an illustrative example is given in Section $4$.

2.   Preliminaries

In the sequel, $\Psi$ denotes an increasing map $\Psi:[r_{1}, r_{2}]\rightarrow \mathbb{R}$ via $\Psi^{\prime }(\varsigma)\neq0$, $\forall$ $\varsigma$, and $[\alpha]$ indicates the integer part of the real number $\alpha$.

Definition 2.1. ^[4,5] Suppose the continuous function $\varkappa:(0, \infty)\rightarrow \mathbb{R}$. We define (RLFD) the Riemann-Liouville fractional derivative of order $\alpha > 0, n = [\alpha]+1$ by

where $n-1 < \alpha < n$.

Definition 2.2. ^[4,5] The $\Psi$-Riemann-Liouville fractional integral ($\Psi$-RLFI) of order $\alpha > 0$ of a continuous function $\varkappa:\left[ r, \Im\right] \rightarrow \mathbb{R}$ is defined by

Definition 2.3. ^[4,5] The CFD of order $\alpha > 0$ of a function $\varkappa:\left[ 0, +\infty\right) \rightarrow \mathbb{R}$ is defined by

Definition 2.4. ^[4,5] We define the $\Psi$-CFD of order $\alpha > 0$ of a continuous function $\varkappa:\left[r, \Im\right] \rightarrow \mathbb{R}$ by

where $\partial_{\Psi}^{n} = \left(\frac{1}{\Psi^{\prime}(\varsigma)}\frac {d}{d\varsigma}\right) ^{n}, \, n\in\mathbb{N}$.

Lemma 2.1. ^[4,5] Suppose $q, \ell > 0$, and $\varkappa \mathit{\mbox{in}} \mathcal{C} ([r, \mathcal{\Im}], \mathbb{R})$. Then $\forall\varsigma\in\lbrack r, \mathcal{\Im}]$ and by assuming $F_{r}(\varsigma) = \Psi(\varsigma)-\Psi(r)$, we have

1) $\mathcal{I}_{r}^{q; \Psi}\mathcal{I}_{r}^{\ell; \Psi}\varkappa(\varsigma) = \mathcal{I}_{r}^{q+\ell; \Psi}\varkappa(\varsigma),$

2) $\mathfrak{\mathcal{D}}_{r}^{q; \Psi}\mathcal{I}_{r}^{q; \Psi}\varkappa (\varsigma) = \varkappa(\varsigma),$

3) $\mathcal{I}_{r}^{q; \Psi}(F_{r}(\varsigma))^{\ell-1} = \frac{\Gamma(\ell)}{\Gamma(\ell+q)}(F_{r}(\varsigma))^{\ell+q-1},$

4) $\mathfrak{\mathcal{D}}_{r}^{q; \Psi}(F_{r}(\varsigma))^{\ell-1} = \frac{\Gamma(\ell)}{\Gamma(\ell-q)}(F_{r}(\varsigma))^{\ell-q-1},$

5) $\mathfrak{\mathcal{D}}_{r}^{q; \Psi}(F_{r}(\varsigma))^{k} = 0, \ k = 0, \ldots, n-1, \ n\in\mathbb{N}, \ q \mathit{\mbox{in}} (n-1, n]$.

Lemma 2.2. ^[4,5] Let $n-1 < \alpha_{1}\leq n, \alpha _{2} > 0, \ r > 0, \ \varkappa\in L(r, \mathcal{\Im}), \ \mathfrak{\mathcal{D}}_{r}^{\alpha_{1};\Psi}\varkappa\in L(r, \mathcal{\Im})$. Then the differential equation

has the unique solution

and

with $\mathcal{W}_{\ell}\in \mathbb{R}, \ \ell \in \{0, 1, \ldots, n-1\}$.

Furthermore,

and

Here we will deal with the FDE solution of (1.1) and (1.2), by considering the solution of

bounded by the condition (1.2). We set

Lemma 2.3. Let $\nu\in(1, 2]$ and $\varsigma\in\lbrack r, \mathcal{\Im} ]$. Then, the FBVP (2.1) and (1.2) have a solution $\varkappa$ of the form

where

$i = 1, 2, ..., m-2$.

Proof. According to the Lemma 2.2 the solution of $\mathfrak{\mathcal{D}}_{r}^{\nu; \psi}\varkappa (\varsigma) = -h(\varsigma)$ is given by

where $c_{0}, c_{1}\in\mathbb{R}$. Since $\varkappa(r) = \vartheta_{1}$ and $\varkappa(\mathcal{\Im}) = \sum _{i = 1}^{m-2}\zeta_{i}\varkappa(\mathfrak{\eta}_{i})+\vartheta_{2}$, we get $c_{0} = \vartheta_{1}$ and

By substituting $c_{0}, c_{1}$ into Eq (2.3) we find,

where $\mathfrak{\varpi}(\varsigma, \tau)$ is given by (2.2). Hence the required result.

Lemma 2.4. If $0 < \sum_{i = 1}^{m-2}\zeta_{i} < 1$, then

$i)$ $\Delta > 0$,

$ii)$ $(\Psi\left(\mathcal{\Im}\right) -\Psi\left(\tau\right))^{\nu -1}-\sum_{j = i}^{m-2}\zeta_{j}(\Psi\left(\mathfrak{\eta}_{j}\right) -\Psi\left(\tau\right))^{\nu-1} > 0$.

Proof. $i)$ Since $\mathfrak{\eta}_{i} < \mathcal{\Im}$, we have

If $1-\Sigma_{i = 1}^{m-2}\zeta_{i} > 0$, then $\left(\Psi\left(\mathcal{\Im }\right) -\Psi\left(r\right) \right) -\Sigma_{i = 1}^{m-2}\zeta_{i} (\Psi\left(\mathfrak{\eta}_{i}\right) -\Psi\left(r\right)) > 0$. So we have $\Delta > 0$.

$ii)$ Since $0 < \nu-1\leq1$, we have $(\Psi\left(\mathfrak{\eta}_{i}\right) -\Psi\left(\tau\right))^{\nu-1} < (\Psi\left(\mathcal{\Im}\right) -\Psi\left(\tau\right))^{\nu-1}$. Then we obtain

and so

Remark 2.1. Note that $\int_{r}^{\mathcal{\Im}} \mathfrak{\varpi}(\varsigma, \tau)\Psi^{\prime}\left(\tau\right) d\tau$ is bounded $\forall \varsigma\in[r, \mathcal{\Im}]$. Indeed

Remark 2.2. Suppose $\Upsilon(\varsigma)\in L^{1}[r, \mathcal{\Im}]$, and $w(\varsigma)$ verify

then $w(\varsigma) = \int_{r}^{\mathcal{\Im}}\mathfrak{\varpi}(\varsigma, \tau)\Upsilon(\tau)\Psi^{\prime}\left(\tau\right) d\tau$.

Next we recall the Schauder fixed point theorem.

Theorem 2.1. ^[23] [SFPT] Consider the Banach space $\Omega$. Assume $\aleph$ bounded, convex, closed subset in $\Omega$. If $\digamma:\aleph\rightarrow \aleph$ is compact, then it has a fixed point in $\aleph$.

3.   Existence result

We start this section by listing two conditions which will be used in the sequel.

● $(\Sigma1)$ There exists a nonnegative function $\Upsilon\in L^{1} [r, \mathcal{\Im}]$ such that $\int_{r}^{\mathcal{\Im}}\Upsilon(\varsigma)d\varsigma > 0$ and $\mathcal{F}(\varsigma, \varkappa, v)\geq-\Upsilon (\varsigma)$ for all $(\varsigma, \varkappa, v)\in\lbrack r, \mathcal{\Im} ]\times\mathbb{R}\times\mathbb{R}$.

● $(\Sigma2)$ $\mathcal{G} (\varsigma, \varkappa, v)\neq0$, for $(\varsigma, \varkappa, v)\in\lbrack r, \mathcal{\Im}]\times\mathbb{R}\times\mathbb{R}$.

Let $\aleph = \mathcal{C}([r, \mathcal{\Im}], \mathbb{R})$ the Banach space of CFs (continuous functions) with the following norm

First of all, it seems that the FDE below is valid

Here the existence of solution satisfying the condition (1.2), such that $\mathcal{G}:[r, \mathcal{\Im}]\times\mathbb{R}\times\mathbb{R}\rightarrow \mathbb{R}$

and $\varkappa^{\ast}(\varsigma) = \max\{(\varkappa-w)(\varsigma), 0\},$ hence the problem (2.5) has $w$ as unique solution. The mapping $Q:\aleph\rightarrow \aleph$ accompanied with the (3.1) and (1.2) defined as

where the relation (2.2) define $\mathfrak{\varpi} (\varsigma, \tau)$. The existence of solution of the problems (3.1) and (1.2) give the existence of a fixed point for $Q$.

Theorem 3.1. Suppose the conditions $(\Sigma1)$ and $(\Sigma2)$ hold. If there exists $\rho > 0$ such that

where $L\geq\max\{\left\vert \mathcal{G}(\varsigma, \varkappa, v)\right\vert :\varsigma\in\lbrack r, \mathcal{\Im}],$ $|\varkappa |, |v|\leq\rho\}$ and $M$ is defined in (2.4), then, the problems (3.1) and (3.2) have a solution $\varkappa (\varsigma)$.

Proof. Since $P: = \{\varkappa\in\aleph:\Vert\varkappa\Vert\leq\rho\}$ is a convex, closed and bounded subset of $B$ described in the Eq (3.3), the SFPT is applicable to $P$. Define $Q:P\rightarrow \aleph$ by (3.3). Clearly $Q$ is continuous mapping. We claim that range of $Q$ is subset of $P$. Suppose $\varkappa\in P$ and let $\varkappa^{\ast}(\varsigma)\leq\varkappa(\varsigma)\leq\rho$, $\forall \varsigma\in\lbrack r, \mathcal{\Im}]$. So

for all $\varsigma\in\lbrack r, \mathcal{\Im}]$. This indicates that $\Vert Q\varkappa\Vert\leq\rho$, which proves our claim. Thus, by using the Arzela-Ascoli theorem, $Q:\aleph\rightarrow \aleph$ is compact. As a result of SFPT, $Q$ has a fixed point $\varkappa$ in $P$. Hence, the problems (3.1) and (1.2) has $\varkappa$ as solution.

Lemma 3.1. $\varkappa^{\ast}(\varsigma)$ is a solution of the FBVP (1.1), (1.2) and $\varkappa(\varsigma) > w(\varsigma)$ for every $\varsigma\in\lbrack r, \mathcal{\Im}]$ iff the positive solution of FBVP (3.1) and (1.2) is $\varkappa = \varkappa^{\ast}+w$.

Proof. Let $\varkappa(\varsigma)$ be a solution of FBVP (3.1) and (1.2). Then

So,

Then we get the existence of the solution with the condition

For the converse, if $\varkappa^{\ast}$ is a solution of the FBVP (1.1) and (1.2), we get

which leads to

We easily see that

i.e., $\varkappa(r) = \vartheta_{1}$ and

So,

Thus $\varkappa(\varsigma)$ is solution of the problem FBVP (3.1) and (3.2).

4.   Example

We propose the given FBVP as follows

Let $\Psi\left(\varsigma\right) = \log\varsigma,$ where $\mathcal{F} (\varsigma, \varkappa(\varsigma), \varkappa(1+\frac{1}{2}\varsigma)) = \frac{\varsigma}{1+\varsigma}\arctan(\varkappa(\varsigma)+\varkappa (1+\frac{1}{2}\varsigma))$.

Taking $\Upsilon(\varsigma) = \varsigma$ we get $\int_{1}^{e}\varsigma d\varsigma = \frac{e^{2}-1}{2} > 0$, then the hypotheses $(\Sigma1)$ and $(\Sigma2)$ hold. Evaluate $\Delta\cong0.366$, $M\cong3.25$ we also get $\left\vert \mathcal{G}(\varsigma, \varkappa, v)\right\vert < \pi+e = L$ such that $|\varkappa|\leq\rho$, $\rho = 17$, we could just confirm that

By applying the Theorem 3.1 there exit a solution $\varkappa (\varsigma)$ of the problem (4.1) and (4.2).

5.   Conclusions

In this paper, we have provided the proof of BVP solutions to a nonlinear $\Psi$-Caputo fractional pantograph problem or for a semi-positone multi-point of (1.1) and(1.2). What's new here is that even using the generalized $\Psi$-Caputo fractional derivative, we were able to explicitly prove that there is one solution to this problem, and that in our findings, we utilize the SFPT. The results obtained in our work are significantly generalized and the exclusive result concern the semi-positone multi-point $\Psi$-Caputo fractional differential pantograph problem (1.1) and (1.2).

Acknowledgment

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through Small Groups (RGP.1/350/43).

Conflict of interest

The authors declare no conflict of interest.