An Outlook on Hybrid Fractional Modeling of a Heat Controller with Multi-Valued Feedback Control

Abstract

In this study, we extend the investigations of fractional-order models of thermostats and guarantee the solvability of hybrid Caputo fractional models for heat controllers, satisfying some nonlocal hybrid multi-valued conditions with multi-valued feedback control, which involves the Chandrasekhar kernel, by using hybrid Dhage’s fixed point theorem. A part of this study is dedicated to transforming this problem into an equivalent integral representation and then proving some existence results to achieve our aims. Furthermore, the continuous dependence of the unique solution on the control variable and on the set of selections will be discussed. Moreover, we provide an illustration to support our results.

1. Introduction

Hundreds of years ago, humans tried to find a tool or an instrument that could control heat exchange or temperature. A device that could control heat transfer or temperature. The object was to have an easier life, allowing technology to take care of this virtual task. The result was the thermostat, which is found in furnaces, air conditioners, refrigerators, cars, etc. Many scholars have discussed mathematical models for thermostats, for example, refs. [1,2,3,4,5,6,7,8,9,10].

In 1997 [2], two new mathematical models characterizing the dynamic behavior of motor vehicle thermostats were introduced; these models including delay-differential equations have been solved. Another modern mathematical model of the energetic behavior of indoor regulators found in an engine’s cooling framework was presented, in addition to a calculation of numerical solutions [3].

Webb [4] introduced a mathematical treatment for thermostats in 2005. The first mathematical model for thermostat control was developed by Webb [4] in the form:

$\begin{array}{c}\hfill \left\{\begin{array}{c}{\nu }^{″}\left(\tau \right)+h\left(\tau \right)\mathcal{G}(\tau ,\nu \left(\tau \right))=0\hfill \\ {\nu }^{\prime }\left(0\right)=0,\mathcal{b}{\nu }^{\prime }\left(\tau \right)+\nu \left(\tau \right)=0,\tau \in [0,1],\mathcal{b}>0\hfill \end{array}\right.\end{array}$

A model for a thermostat with a second-order nonlocal boundary value problem was established in 2012 by Webb [5]. The sensors act linearly; one gives feedback to a controller at one endpoint from a portion of the interval, and the other provides feedback to a controller at the other extremity. Numerous positive solutions and a few nonexistent solutions were established by the proof of certain useful characteristics.

Second-order differential equations for inclusion and fractional hybrid versions of thermostat models have also been published [11]. On this issue, hybrid boundary value criteria have also been taken into consideration [11]. The system’s historical memory is protected by the Caputo–Fabrizio fractional-order derivative, which has been used to model and study the complications of childhood mumps-related hearing loss in [12].

To characterize the energetic behavior of a car indoor regulator, two modern models including delay-differential conditions with hysteresis were formulated in 1997 [2], and it was discovered that these two models were solvable. An entirely novel mathematical representation of the dynamic behavior of a thermostat installed in an engine cooling circuit was demonstrated, along with a calculation of numerical results [13].

Hybrid differential equations have received great attention [1,9,14,15]. Dhage and Lakshmikantham [14] initiated and presented a discussion of hybrid differential equations. A generalized version of the hybrid Dhage’s fixed point results was used by Baleanu et al. [15].

Shen et al. [6] established a fractional order model for a thermostat using the same boundary conditions as in [5].

A further extension of the second-order differential equation of a thermostat model to a fractional hybrid equation with nonlocal hybrid conditions has been considered in [11]:

(1)$\begin{array}{c}\hfill {-}^c{\mathcal{D}}^{\gamma }\left(\frac{z\left(\mathfrak{r}\right)}{h(\mathfrak{r},z(\mathfrak{r}\left)\right)}\right)\in \mathrm{\Psi }(\mathfrak{r},z(\mathfrak{r}\left)\right),\mathfrak{r}\in [0,1]\end{array}$

with the hybrid conditions

(2)$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{D}\left(\frac{z\left(\mathfrak{r}\right)}{h(\mathfrak{r},z(\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=0}=0,\hfill \\ {\lambda }_1{}^c{\mathcal{D}}^{\gamma -1}\left(\frac{z\left(\mathfrak{r}\right)}{h(\mathfrak{r},z(\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=1}+{\lambda }_2\left(\frac{z\left(\mathfrak{r}\right)}{h(\mathfrak{r},z(\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=\eta }=0,\hfill \end{array}\right.\end{array}$

Some existence results have been investigated. Moreover, two examples are illustrated in [11].

Recently, the authors of [9] have established a model for thermostats involving hybrid integro-differential inclusions:

$-\frac{d^2}{d{\mathfrak{r}}^2}\left(\frac{\nu \left(\mathfrak{r}\right)}{h(\mathfrak{r},\nu (\mathfrak{r}\left)\right)}\right)\in {\int }_0^1\frac{\mathfrak{r}}{\mathfrak{r}+\tau }\mathrm{\Phi }\left(\tau ,{\int }_0^1\frac{\tau }{\tau +\varrho }\psi (\varrho ,\nu \left(\varrho \right))d\varrho \right)d\tau ,\mathfrak{r}\in [0,1]$

satisfying the hybrid nonlocal conditions:

$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{D}\left(\frac{\nu \left(\mathfrak{r}\right)}{h(\mathfrak{r},\nu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=0}=0,\hfill \\ \lambda {}^c{\mathcal{D}}^{\gamma }\left(\frac{\nu \left(\mathfrak{r}\right)}{h(\mathfrak{r},\nu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=\alpha }+\left(\frac{\nu \left(\mathfrak{r}\right)}{h(\mathfrak{r},\nu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=\eta }=0,\gamma \in (0,1],\alpha \in (0,1],\eta \in (0,1]\hfill \end{array}\right.\end{array}$

and proved some existence and continuous dependency results.

Inspired by the investigated mathematical models in [9,11], we establish some existence results for hybrid fractional modeling of thermostats.

(3)$\begin{array}{c}\hfill {-}^c{\mathcal{D}}^{\gamma }\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right)\in {\int }_0^{\mathfrak{r}}\frac{\tau }{\mathfrak{r}+\tau }\mathrm{\Phi }\left(\tau ,\nu \left(\tau \right),{\int }_0^1\frac{\tau }{\tau +\varrho }\psi (\varrho ,\mu \left(\varrho \right))d\varrho \right)d\tau ,\mathfrak{r}\in I=[0,1]\end{array}$

with a nonlocal hybrid multi-valued boundary condition

(4)$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{D}\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=0}=0,\hfill \\ {\lambda }_1{}^c{\mathcal{D}}^{\gamma -1}\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=1}+{\lambda }_2\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=\eta }\in \mathrm{\Theta }(\mathfrak{r},\mu \left(\mathfrak{r}\right)),\hfill \end{array}\right.\end{array}$

and with a multi-valued control variable in the form of

(5)$\nu \left(\mathfrak{r}\right)\in \mathrm{\Omega }\left(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\right),$

where $\gamma >0$ is a real number with $n-1<\gamma <n,$ and ${\lambda }_i,i=1,2$ are positive real parameters, $\mathcal{D}=\frac{d}{dt}$, ${}^c{\mathcal{D}}^{\gamma }$ is the Caputo derivative of order $\gamma ,$ where $\mathrm{\Phi },\mathrm{\Omega }:I\times \mathbb{R}\times \mathbb{R}\to P\left(\mathbb{R}\right)$, $\mathrm{\Theta }:I\times \mathbb{R}\to P\left(\mathbb{R}\right)$ are multi-valued maps, $\psi :I\times \mathbb{R}\to \mathbb{R}$ is continuous, and $\mathcal{F}\in C(I\times \mathbb{R},\mathbb{R}\backslash \{0\left\}\right)$.

This study is the first attempt to discuss the solvability of the hybrid fractional model of thermostats (3) satisfying the nonlocal hybrid multi-valued condition (4) under multi-valued constraints (5) in $C(I,\mathbb{R})$. Furthermore, it will be established that the solution of this problem is unique and it depends continuously on the control variable (5) and on the set of selections $S_{\mathrm{\Phi }}.$ Finally, an example is presented to clarify our results.

To reach our goal, we need to investigate the single-valued problem that corresponds to the mentioned problem (3) and (4)

(6)$\begin{array}{cc}\hfill {-}^c{\mathcal{D}}^{\gamma }\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right)& ={\int }_0^{\mathfrak{r}}\frac{\tau }{\mathfrak{r}+\tau }\varphi \left(\tau ,\nu \left(\tau \right),{\int }_0^1\frac{\tau }{\tau +\varrho }\psi (\tau ,\mu \left(\varrho \right))d\varrho )\right)d\tau ,\mathfrak{r}\in I\hfill \end{array}$

with nonlocal hybrid condition

(7)$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{D}\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=0}=0,\hfill \\ {\lambda }_1{}^c{\mathcal{D}}^{\gamma -1}\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=1}+{\lambda }_2\left(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)}\right){|}_{\mathfrak{r}=\eta }=\theta (\mathfrak{r},\mu \left(\mathfrak{r}\right)),\hfill \end{array}\right.\end{array}$

and the control variable is provided as

(8)$\nu \left(\mathfrak{r}\right)=\varpi \left(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\right),$

with $\varphi \in S_{\mathrm{\Phi }}$, $\theta \in S_{\mathrm{\Theta }},$ and $\varpi \in S_{\mathrm{\Omega }}$.

Our problem (3) and (4) involves Chandrasekhar’s kernel; integral equations containing this kernel have been treated and discussed by many scholars in different classes and by various techniques due to their usage in numerous branches of research and engineering, including traffic theory, neutron transport theory, kinetic theory of gases, and radiative transfer theory (for examples, see [16,17]).

2. Single-Valued Problem

Consider the nonlocal problem (6) and (7) with feedback control (8), assuming the following:

• $\left({\mathcal{H}}_1\right)$

  $\varphi :I\times \mathbb{R}\times \mathbb{R}\to \mathbb{R}$ is continuous in $\nu ,\tau$ for every $\mathfrak{r}\in I$, and measurable for almost all $\mathfrak{r}\in I$ and $\forall \nu ,\tau \in \mathbb{R}$. There are two integrable functions $m,k_1:I\to I$ with

  $\left|\varphi (\mathfrak{r},\nu ,\tau )\right|\le m\left(\mathfrak{r}\right)+k_1\left(\mathfrak{r}\right)\left(\right|\nu |+\left|\tau \right|),\mathfrak{r}\in I$

  and

  ${\int }_0^1\frac{1}{\tau +\varrho }\left|m\left(\varrho \right)\right|d\varrho \le m,\mathrm{and}{\int }_0^1\frac{1}{\tau +\varrho }\left|k_1\left(\varrho \right)\right|d\varrho \le k.$

• $\left({\mathcal{H}}_2\right)$

  $\psi \in C(I\times \mathbb{R},\mathbb{R})$, and there exists a continuous function $k_2:[0,1]\times [0,1]\to \mathbb{R}$, and a non-decreasing continuous map $\chi :[0,\infty )\to (0,\infty ),$ with

  $\begin{array}{c}\hfill \left|\psi \right(\mathfrak{r},\tau \left)\right|\le k_2\left(\mathfrak{r}\right)\chi (\parallel \tau \parallel ),\end{array}$

  and

  ${\int }_0^1\frac{1}{\tau +\varrho }\left|k_2\left(\varrho \right)\right|d\varrho \le k^{*}.$

• $\left({\mathcal{H}}_3\right)$

  Let $\theta :I\times \mathbb{R}\to \mathbb{R}$ be a Lipschitzian function, with

  $|\theta (\mathfrak{r},{\mu }_1)-\theta (\mathfrak{r},{\mu }_2)|\le k_3\left(\mathfrak{r}\right)|{\mu }_1-{\mu }_2|.$

• $\left({\mathcal{H}}_4\right)$

  $\mathcal{F}\in C(I\times \mathbb{R},\mathbb{R}\backslash \{0\left\}\right)$ and there is a continuous function $\omega :I\to I$, where

  $|\mathcal{F}(\mathfrak{r},{\mu }_1)-\mathcal{F}(\mathfrak{r},{\mu }_2)|\le \omega \left(\mathfrak{r}\right)|{\mu }_1-{\mu }_2|,$

  $\forall {\mu }_1,{\mu }_2\in \mathbb{R}$ and $\mathfrak{r}\in I$.

• $\left({\mathcal{H}}_5\right)$

  $\omega \in C(I\times \mathbb{R}\times \mathbb{R},\mathbb{R})$, and there exists a measurable and bounded function $\delta :I\to \mathbb{R},$ which has norm $\parallel \delta \parallel ,$ with

  $\left|\varpi \right(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\left)\right|\le \delta \left(\mathfrak{r}\right),\mathfrak{r}\in I,$

  where $\delta ={max}_{\mathfrak{r}\in I}\{\parallel \delta \parallel \}$.

• $\left({\mathcal{H}}_6\right)$

  The real number r is the positive root of

  $\begin{array}{cc}\hfill & -K_3\parallel \omega \parallel r^2+\left(1-\left[m\parallel \omega \parallel +\parallel \omega \parallel \theta +k\left(k^{*}\chi (\parallel \vartheta \parallel )+\parallel \delta \parallel \right)]\mathrm{\Lambda }+k_{3}G\right]\right)r\\ & -G(\frac{1}{{\lambda }_2}(\mathrm{\Theta }+k_{3}1+m)+[m+k{k}^{*}\chi (\parallel \mu \parallel )]\mathrm{\Lambda })=0,\end{array}$

  where $G={sup}_{\mathfrak{r}\in I}\left|\mathcal{F}(\mathfrak{r},0)\right|$, and

  (9)$\mathrm{\Lambda }=\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\lambda }_2{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}.$

Remark 1.

Using assumptions $\left({\mathcal{H}}_3\right)$ and $\left({\mathcal{H}}_4\right)$, we have

$\left|\theta (\mathfrak{r},\mu )\right|\le k_3\left(\mathfrak{r}\right)\left|\mu \right|+\mathrm{\Theta },\mathrm{\Theta }=\underset{\mathfrak{r}\in I}{sup}\left|\theta (\mathfrak{r},0)\right|.$

and

$\begin{array}{cc}\hfill \left|\mathcal{F}\right(\mathfrak{r},\mu \left)\right|& \le \omega \left|\mu \left(\mathfrak{r}\right)\right|+G,withG=\underset{\mathfrak{r}\in I}{sup}\left|\mathcal{F}(\mathfrak{r},0)\right|.\hfill \end{array}$

Lemma 1.

$\mu \in C(I,\mathbb{R})$ is a solution of the hybrid differential equation

(10) $\begin{array}{c}\hfill {}^c{\mathcal{D}}^{\gamma }(\frac{\mu \left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)})+{\int }_0^{\mathfrak{r}}\frac{\mathfrak{r}}{\mathfrak{r}+\tau }\chi \left(\tau \right)d\tau =0\mathfrak{r}\in [0,1],\gamma \in (1,2],\end{array}$

with the condition (7) and feedback control (8) if $\mu \in C(I,\mathbb{R})$ is a solution of the following equation:

(11) $\begin{array}{cc}\hfill \mu \left(\mathfrak{r}\right)=& \mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))[-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau ]-\frac{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)\theta (\mathfrak{r},\mu (\mathfrak{r}\left)\right)}{{\lambda }_2}.\hfill \end{array}$

Proof. 

Assume that ${\mu }_0$ is a solution of (10). Then, we find constants ${\alpha }_0,{\alpha }_1\in \mathbb{R}$ that satisfy

(12)${\mu }_0=\mathcal{F}(\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))\left[-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +{\alpha }_0+{\alpha }_1\mathfrak{r}\right].$

Then,

$\mathcal{D}(\frac{{\mu }_0\left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))})=-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -2}}{\mathrm{\Gamma }(\gamma -1)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +{\alpha }_1$

and

${}^c{\mathcal{D}}^{\gamma -1}(\frac{{\mu }_0\left(\mathfrak{r}\right)}{\mathcal{F}(\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))})=-{\int }_0^{\mathfrak{r}}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +{\alpha }_1\frac{{\mathfrak{r}}^{2-\gamma }}{\mathrm{\Gamma }(3-\gamma )}.$

Therefore, ${\alpha }_1=0$ and

${\alpha }_0=\frac{\theta (\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))}{{\lambda }_2}+\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau .$

We get, by replacing the values ${\alpha }_0$ and ${\alpha }_1$ in (12),

$\begin{array}{cc}\hfill {\mu }_0\left(\mathfrak{r}\right)=& \mathcal{F}(\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))[-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau ]-\frac{\mathcal{F}(\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))\theta (\mathfrak{r},{\mu }_0\left(\mathfrak{r}\right))}{{\lambda }_2}\hfill \end{array}$

This indicates that for the fractional integral Equation (11), ${\mu }_0$ is the solution. Conversely, it is obvious that for the fractional hybrid problem (7) and (10), ${\mu }_0$ is a solution of (11).  □

Corollary 1.

Let $\mu \in C(I,\mathbb{R})$ be a solution of problem (6) and (7) with feedback control (8). Then, it satisfies

(13) $\begin{array}{cc}\hfill \mu \left(\mathfrak{r}\right)=& \mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))[-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau ]\hfill \\ \hfill & -\frac{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)\theta (\mathfrak{r},\mu (\mathfrak{r}\left)\right)}{{\lambda }_2}.\hfill \end{array}$

Proof. 

From Lemma 1, we have

$\begin{array}{cc}\hfill \mu \left(\mathfrak{r}\right)=& \mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))[-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho ds\hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\chi \left(\varrho \right)d\varrho d\tau ]-\frac{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)\theta (\mathfrak{r},\mu (\mathfrak{r}\left)\right)}{{\lambda }_2}.\hfill \end{array}$

Now, for

$\chi \left(\varrho \right)=\varphi (\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma ),\varrho ,\varsigma \in I$

we obtain the result. □

2.1. Existence of Solutions

Theorem 1.

Let (${\mathcal{H}}_1$)–(${\mathcal{H}}_6$) be verified. Therefore, a solution for (13) exists.

Proof. 

Consider the ball ${\mathfrak{V}}_{ϵ}\left(0\right)={\{\tau \in X:\parallel \tau \parallel }_X\le ϵ\}.$

Clearly, ${\mathfrak{V}}_{ϵ}\left(0\right)$ is a closed, convex, and bounded subset of the Banach space X. Regard the operators $\mathcal{A}:X\to X,\mathcal{B}:{\mathfrak{V}}_{ϵ}\left(0\right)\to X$ defined by:

(14)$\left(\mathcal{A}\mu \right)\left(\mathfrak{r}\right)=ϝ(\mathfrak{r},\mu (\mathfrak{r}\left)\right),\mathfrak{r}\in I$

(15)$\begin{array}{cc}\hfill \left(\mathcal{B}\mu \right)\left(\mathfrak{r}\right)=& -{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau ,\hfill \end{array}$

and

(16)$\begin{array}{cc}\hfill \left(\mathcal{C}\mu \right)\left(\mathfrak{r}\right)=& \frac{\mathcal{F}(\mathfrak{r},\mu (\mathfrak{r}\left)\right)\theta (\mathfrak{r},\mu (\mathfrak{r}\left)\right)}{{\lambda }_2}.\hfill \end{array}$

As a solution to problem (6) and (7) with feedback control (8) exists, it is evident that X satisfies the operator equation $\mathcal{A}\mu \mathcal{B}\mu +\mathcal{C}\mu =\mu$. By utilizing the presumptions of the theorem of three operators with Banach algebra, due to Dhage [18] and the problem (6) and (7) with (8), we show that such a solution exists.

In the beginning, we demonstrate that operators $\mathcal{A},\mathcal{C}$ are Lipschitzian with a constant $\parallel \omega \parallel$ on normed space X. For evidence of this, take ${\mu }_1,{\mu }_2\in X,$ then

$\begin{array}{ccc}\hfill |\left(\mathcal{A}{\mu }_1\right)\left(\mathfrak{r}\right)-\left(\mathcal{A}{\mu }_2\right)\left(\mathfrak{r}\right)|& =& |\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\hfill \\ & \le & \omega \left(\mathfrak{r}\right)|{\mu }_1\left(\mathfrak{r}\right)-{\mu }_2\left(\mathfrak{r}\right)|,\hfill \end{array}$

∀${\mu }_1,{\mu }_2\in {\mathfrak{V}}_{ϵ}\left(0\right)$, and $\mathfrak{r}\in I$, and then

$\parallel \mathcal{A}{\mu }_1-\mathcal{A}{\mu }_2{\parallel }_{\mathcal{X}}\le \parallel \omega \parallel \parallel {\mu }_1-{\mu }_2{\parallel }_{\mathcal{X}},$

The operator $\mathcal{A}$ is then Lipschitzian on ${\mathfrak{V}}_{ϵ}\left(0\right)$ with the constant $\parallel \omega \parallel$.

Similarly, we have ∀${\mu }_1,{\mu }_2\in {\mathfrak{V}}_{ϵ}\left(0\right)$

$\begin{array}{ccc}& & |\left(\mathcal{C}{\mu }_1\right)\left(\mathfrak{r}\right)-\left(\mathcal{C}{\mu }_2\right)\left(\mathfrak{r}\right)|\hfill \\ & =& \frac{1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))\theta (\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))\theta (\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\hfill \\ & \le & \frac{1}{{\lambda }_2}\left(|\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))\left|\right|\theta (\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\theta (\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|+|\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))\left|\right|\theta (\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\right)\hfill \\ & \le & \frac{1}{{\lambda }_2}\left([\parallel \omega \parallel |{\mu }_1\left(\mathfrak{r}\right)|+G]k_3|{\mu }_1\left(\mathfrak{r}\right)-{\mu }_2\left(\mathfrak{r}\right)|+\omega \left(\mathfrak{r}\right)|{\mu }_1\left(\mathfrak{r}\right)-{\mu }_2\left(\mathfrak{r}\right)\left|\right[\mathrm{\Theta }+k_3\left(\mathfrak{r}\right)\parallel u_2\parallel ]\right)\hfill \\ & \le & \frac{\left(\right[\parallel \omega \parallel \parallel {\mu }_1\parallel +G]k_3+\parallel \omega \parallel [\mathrm{\Theta }+k_3\parallel u_2\parallel \left]\right)}{{\lambda }_2}|{\mu }_1\left(\mathfrak{r}\right)-{\mu }_2\left(\mathfrak{r}\right)|.\hfill \end{array}$

When applying the supremum to $\mathfrak{r}\in I$, we obtain

$\parallel \mathcal{C}{\mu }_1-\mathcal{C}{\mu }_2{\parallel }_X\le \frac{\left(\right[\parallel \omega \parallel \parallel {\mu }_1\parallel +G]k_3+\parallel \omega \parallel [\mathrm{\Theta }+k_3\parallel u_2\parallel \left]\right)}{{\lambda }_2}{\parallel {\mu }_1-{\mu }_2\parallel }_X,$

Then, $\mathcal{C}$ is a Lipschitz mapping on X with the Lipschitz constant $\frac{\left(\right[\parallel \omega \parallel \parallel {\mu }_1\parallel +G]k_3+\parallel \omega \parallel [\mathrm{\Theta }+k_3\parallel u_2\parallel \left]\right)}{{\lambda }_2}$.

In second step, the aim is to show that $\mathcal{B}$ is continuous and compact on ${\mathfrak{V}}_{ϵ}\left(0\right)$ into X; thus, the continuity of $\mathcal{B}$ on ${\mathfrak{V}}_{ϵ}\left(0\right)$ is demonstrated. Taking $\left\{{\mu }_n\right\}$ as the series converges to $\mu \in {\mathfrak{V}}_{ϵ}\left(0\right)$ and considering that $\mathfrak{r}\in I,$ the continuous functions $\varphi (\mathfrak{r},\mu (\mathfrak{r}\left)\right)$ and $\psi (\mathfrak{r},\mu (\mathfrak{r}\left)\right)$ in $X$ give $\varphi (\mathfrak{r},{\mu }_n\left(\mathfrak{r}\right))\to \varphi (\mathfrak{r},\mu \left(\mathfrak{r}\right)),$ and $\psi (\mathfrak{r},{\mu }_n\left(\mathfrak{r}\right))\to \psi (\mathfrak{r},\mu \left(\mathfrak{r}\right))$ (from $\left({\mathcal{H}}_2\right)$) and $\left({\mathcal{H}}_3\right)$), by the Lebesgue Dominated Convergence Theorem. Then, we have

$\begin{array}{cc}\hfill \underset{n\to \infty }{lim}\left(\mathcal{B}{\mu }_n\right)\left(\mathfrak{r}\right)& =-\underset{n\to \infty }{lim}{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\underset{n\to \infty }{lim}\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\underset{n\to \infty }{lim}{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & =-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\underset{n\to \infty }{lim}\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\underset{n\to \infty }{lim}\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\underset{n\to \infty }{lim}\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & =-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\underset{n\to \infty }{lim}\psi (\varsigma ,{\mu }_n\left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)d\varrho d\tau =\left(\mathcal{B}\mu \right)\left(\mathfrak{r}\right),\hfill \end{array}$

for all $\mathfrak{r}\in I$. Thus, $\mathcal{B}{\mu }_n\to \mathcal{B}\mu$ as $n\to \infty$ uniformly on $\mathbb{R}$; then, the operator $\mathcal{B}$ is continuous on ${\mathfrak{V}}_{ϵ}\left(0\right)$.

Next, we demonstrate that the operator $\mathcal{B}$ is compact on ${\mathfrak{V}}_{ϵ}\left(0\right)$. It is sufficient to prove that $\mathcal{B}\left({\mathfrak{V}}_{ϵ}\left(0\right)\right)$ is a uniformly bounded and equicontinuous set in ${\mathfrak{V}}_{ϵ}\left(0\right)$. Using $\left({\mathcal{H}}_2\right)$ and $\left({\mathcal{H}}_3\right)$, then

$\begin{array}{cc}\hfill & |(\mathcal{B}\mu \left)\right(\mathfrak{r})|\hfill \\ \hfill \le & {\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[m\left(\varrho \right)+k_1\left(\varrho \right)(|\nu \left(\varrho \right)|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,\mu \left(\varsigma \right))|d\varsigma )]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\nu \left(\varrho \right)|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,\mu \left(\varsigma \right))|d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\nu \left(\varrho \right)|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,\mu \left(\varsigma \right))|d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill \le & {\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|\hfill \\ \hfill & +{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel \left)d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^1\frac{\tau }{\tau +\varrho }[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|\hfill \\ \hfill & +{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel \left)d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|\left(|\delta \left(\mathfrak{r}\right)|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel )d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill \le & \left[\frac{m}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}m+\frac{{\lambda }_{2}m{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}\right]+k\left[\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}\right]\left(k^{*}\chi (\parallel \mu \parallel )+\parallel \delta \parallel \right)\hfill \\ \hfill & \le [m+k\left(k^{*}\chi (\parallel \mu \parallel )+\parallel \delta \parallel \right)]\mathrm{\Lambda },\hfill \end{array}$

$\forall \mathfrak{r}\in I$ and $\mu \in {\mathfrak{V}}_{ϵ}\left(0\right)$. As a result, $\parallel \mathcal{B}\tau \parallel \le [m+k\left(k^{*}\chi (\parallel \tau \parallel )+\parallel \delta \parallel \right)]\mathrm{\Lambda }$, such that $\mathrm{\Lambda }$ is given in (9). As a consequence, it could be concluded that the set $\mathcal{B}\left({\mathfrak{V}}_{ϵ}\left(0\right)\right)$ in the normed space X is uniformly bounded. Hence, the equicontinuity of $\mathcal{B}$ is investigated. To achieve our goal, assuming ${\mathfrak{r}}_1,{\mathfrak{r}}_2\in I$ with ${\mathfrak{r}}_1<{\mathfrak{r}}_2$, then

$\begin{array}{ccc}& & |\left(\mathcal{B}\mu \right)\left({\mathfrak{r}}_2\right)-\left(\mathcal{B}\mu \right)\left({\mathfrak{r}}_1\right)|\hfill \\ & \le & {\int }_0^{{\mathfrak{r}}_1}\left(\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}-\frac{{({\mathfrak{r}}_1-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}\right)\hfill \\ & \times & {\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel )d\varsigma \right)\right]d\varrho d\tau \hfill \\ & +& {\int }_{{\mathfrak{r}}_1}^{{\mathfrak{r}}_2}\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[m(\varrho )+k_1\left(\varrho \right)(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|\hfill \\ & +& {\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel \left)d\varsigma \right)]d\varrho d\tau \hfill \\ & \le & {\int }_0^{{\mathfrak{r}}_1}\left(\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}-\frac{{({\mathfrak{r}}_1-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}\right){\int }_0^1\frac{\tau }{\tau +\varrho }[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|(|\delta \left(\varrho \right)|\hfill \\ & +& {\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel \left)d\varsigma \right)]d\varrho d\tau \hfill \\ & +& {\int }_{{\mathfrak{r}}_1}^{{\mathfrak{r}}_2}\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|\left(|\delta \left(\varrho \right)|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \mu \parallel )d\varsigma \right)\right]d\varrho d\tau \hfill \\ & \le & \left[m+k(\parallel \delta \parallel +k^{*}\chi (\parallel \mu \parallel ))\right]\left[{\int }_0^{t_1}\left(\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}-\frac{{({\mathfrak{r}}_1-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}\right)d\tau +{\int }_{{\mathfrak{r}}_1}^{{\mathfrak{r}}_2}\frac{{({\mathfrak{r}}_2-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}d\tau \right],\hfill \end{array}$

This does not depend on $\mu \in {\mathfrak{V}}_{ϵ}\left(0\right).$ Then, ∀$\epsilon >0$, and we find $\rho >0$, where

$|{\mathfrak{r}}_2-{\mathfrak{r}}_1|<\rho \stackrel{\hspace{1em}\hspace{1em}}{\Rightarrow }|\left(\mathcal{B}\mu \right)\left({\mathfrak{r}}_2\right)-\left(\mathcal{B}\mu \right)\left({\mathfrak{r}}_1\right)|<\epsilon ,$

$\forall {\mathfrak{r}}_1,{\mathfrak{r}}_2\in I,\mu \in {\mathfrak{V}}_{ϵ}\left(0\right).$ Then, $\forall \epsilon >0$. This demonstrates $\mathcal{B}$ in X is an equicontinuous set. The Arzela–Ascoli theorem states that $\mathcal{B}$ is compact because it is equicontinuous and uniformly bounded on set X. Consequently, the operator $\mathcal{B}$ on ${\mathfrak{V}}_{ϵ}\left(0\right)$ is completely continuous.

However, by utilizing (${\mathcal{H}}_3$), then

$\begin{array}{c}\hfill {\mathcal{M}}_0^{*}=\parallel \mathcal{B}\left({\mathfrak{V}}_{ϵ}\left(0\right)\right){\parallel }_X=[m+k(\parallel \delta \parallel +k^{*}\chi (\parallel \mu \parallel )\left)\right]\mathrm{\Lambda }.\end{array}$

Putting $L^{*}=\parallel \omega \parallel$, then we get $L^{*}{\mathcal{M}}_0^{*}<1$. Consequently, the Dhage hybrid fixed point theorem’s [19] presumptions hold, and if either condition (a) or (b) is valid, then Dhage’s hybrid fixed point theorem [19] is justified. In order for $\mu =\mathcal{A}\mu \mathcal{B}\vartheta +\mathcal{C}\mu$, let $\mu \in X$ and $\nu \in S$ be random elements. Then, there is

$\begin{array}{cc}\hfill & \left|\mu \right(\mathfrak{r}\left)\right|\hfill \\ & =\left|\mathcal{A}\mu \right(\mathfrak{r}\left)\right|\left|\mathcal{B}\nu \right(\mathfrak{r}\left)\right|+\left|\mathcal{C}\mu \right(\mathfrak{r}\left)\right|\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))|(-{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\vartheta \left(\varsigma \right))d\varsigma \right)|d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)|d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\vartheta \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau )\hfill \\ \hfill & +\frac{\left(\right[\parallel \omega \parallel \parallel \mu \parallel +G\left]\right[\mathrm{\Theta }+k_3\parallel \mu \parallel \left]\right)}{{\lambda }_2}\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))|\left({\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\right[m\left(\varrho \right)+k_1\left(\varrho \right)(\left|\varpi (\varrho ,\nu \left(\varrho \right),\mu \left(\varrho \right))\right|\hfill \\ \hfill & +{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)\left|\chi \right(\parallel \vartheta \parallel \left)d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|\left(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \vartheta \parallel )d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \vartheta \parallel )d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left[|m\left(\varrho \right)|+|k_1\left(\varrho \right)|\left(|\varpi (\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)|\chi (\parallel \vartheta \parallel )d\varsigma \right)\right]d\varrho d\tau )\hfill \\ \hfill & +\frac{\left(\right[\parallel \omega \parallel \parallel \mu \parallel +G\left]\right[\mathrm{\Theta }+k_3\parallel \mu \parallel \left]\right)}{{\lambda }_2}.\hfill \end{array}$

Taking supremum over $\mathfrak{r}\in I$, we have

$\begin{array}{ccc}\hfill \parallel \mu \parallel & \le & \left[\parallel \mu \parallel \parallel \omega \parallel +G\right][m+k\left(k^{*}\chi (\parallel \vartheta \parallel )+\parallel \delta \parallel \right)]\mathrm{\Lambda }+\frac{\left(\right[\parallel \mu \parallel \parallel \omega \parallel +G\left]\right[\mathrm{\Theta }+k_3\parallel \mu \parallel \left]\right)}{{\lambda }_2}\hfill \\ & \le & ϵ.\hfill \end{array}$

Therefore, all of the requirements of Dhage’s hybrid fixed point theorem [19] are held. Therefore, $\mu =\mathcal{A}\mu \mathcal{B}\vartheta +\mathcal{C}\mu$ has a solution in S. Hence, problem (6) and (7) with feedback control (8) is solvable in S on I.  □

2.2. Existence of the Unique Solution

With the aim of proving some uniqueness results of the problem (6) and (7) involving (8), assume that:

• ${\left({\mathcal{H}}_1\right)}^{*}$

  Let $\varphi :I\times \mathbb{R}\times \mathbb{R}\to \mathbb{R}$ be a Lipschitzian mapping, where

  $|\varphi (\mathfrak{r},{\nu }_1,{\mu }_1)-\varphi (\mathfrak{r},,{\nu }_2,{\mu }_2)|\le k_1\left(\mathfrak{r}\right)\left(\right|{\nu }_1-{\nu }_2|+|{\mu }_1-{\mu }_2\left|\right).$

  From this assumption, we see that the assumption $\left({\mathcal{H}}_1\right)$ is valid; then,

  $\left|\varphi (\mathfrak{r},\nu ,\mu )\right|\le k_1\left(\mathfrak{r}\right)\left(\right|\nu |+\left|\mu \right|)+m,m=\underset{\mathfrak{r}\in I}{sup}\left|\varphi (\mathfrak{r},0,0)\right|.$

• ${\left({\mathcal{H}}_2\right)}^{*}$

  $\psi (\mathfrak{r},\mu \left(\mathfrak{r}\right))=k_2\left(\mathfrak{r}\right)\mu \left(\mathfrak{r}\right)$.

Theorem 2.

Let the assumptions of Theorem 1 hold, and replace assumption $\left({\mathcal{H}}_1\right)$ with ${\left({\mathcal{H}}_1\right)}^{*}$ and $\left({\mathcal{H}}_2\right)$ with ${\left({\mathcal{H}}_2\right)}^{*}$ with

$\mathrm{\Lambda }\left[\parallel \omega \parallel (m+k_1(\parallel \delta \parallel +k_{2}r\left)\right)+(\parallel \omega \parallel r+G)k_1{k}_2\right]+\frac{1}{{\lambda }_2}\left[k_3\left[\parallel \omega \parallel r+G\right]+\parallel \omega \parallel \left[k_{3}r+\mathrm{\Theta }\right]\right]<1.$

Then, the hybrid problem (6) and (7) with feedback control (8) has a unique continuous solution.

Proof. 

Let ${\mu }_1,{\mu }_2$ be two solutions of Equation (11), so

$\begin{array}{cc}\hfill & |{\mu }_1\left(\mathfrak{r}\right)-{\mu }_2\left(\mathfrak{r}\right)|\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\hfill \\ \hfill & \times {\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\varpi (\varrho ,\nu \left(\varrho \right),{\mu }_1\left(\varrho \right))|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))|d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }{k}_1\left(\varrho \right){\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))-\psi (\varsigma ,{\mu }_2\left(\varsigma \right))|d\varsigma d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\hfill \\ \hfill & \times {\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\varpi (\varrho ,\nu \left(\varrho \right),{\mu }_1\left(\varrho \right))|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }{k}_1\left(\varsigma \right)|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))|d\varsigma )\right)\right]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }{k}_1\left(\varrho \right){\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))-\psi (\varsigma ,{\mu }_2\left(\varsigma \right))|d\varsigma d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|\hfill \\ \hfill & \times {\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[m\left(\varrho \right)+k_1\left(\varrho \right)\left(|\varpi (\varrho ,\nu \left(\varrho \right),{\mu }_1\left(\varrho \right))|+{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))\right)|d\varsigma )\right]d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }{k}_1\left(\varrho \right){\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))-\psi (\varsigma ,{\mu }_2\left(\varsigma \right))|d\varsigma d\varrho d\tau \hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[|\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))||\theta (\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))-\theta (\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))|+|\mathcal{F}(\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))||\theta (\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))|\right]\hfill \\ \hfill & \le \frac{\parallel \omega \parallel \parallel {\mu }_1-{\mu }_2\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel {\mu }_1\parallel \right)\right]+\frac{\left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}{k}_1{k}_2\parallel {\mu }_1-{\mu }_2\parallel \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\parallel \omega \parallel \parallel {\mu }_1-{\mu }_2\parallel \left[m+k_1\left(\parallel \delta \parallel +k_2\parallel {\mu }_1\parallel \right)\right]\hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }_2\left(\mathfrak{r}\right)|+G\right]k_1{k}_2\parallel {\mu }_1-{\mu }_2\parallel \hfill \\ \hfill & +\frac{{\eta }^{\gamma }\parallel \omega \parallel \parallel {\mu }_1-{\mu }_2\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel {\mu }_1\parallel \right)\right]+\frac{{\eta }^{\gamma }\left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}{k}_1{k}_2\parallel {\mu }_1-{\mu }_2\parallel \hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[k_3\parallel {\mu }_1-{\mu }_2\parallel \left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G\right]+\parallel \omega \parallel \parallel {\mu }_1-{\mu }_2\parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]\right].\hfill \end{array}$

Taking the supremum over $\mathfrak{r}\in I$, we have

$\begin{array}{cc}\hfill \parallel {\mu }_1-{\mu }_2\parallel & \le [\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}]\parallel \omega \parallel \parallel {\mu }_1-{\mu }_2\parallel \left[m+k_1\left(\parallel \delta \parallel +k_{2}r\right)\right]\hfill \\ \hfill & +[\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}]\left[\parallel \omega \parallel r+G\right]k_1{k}_2\parallel {\mu }_1-{\mu }_2\parallel \hfill \\ \hfill & +\parallel {\mu }_1-{\mu }_2\parallel \left[k_3\left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]\right]\hfill \\ \hfill & \le \parallel {\mu }_1-{\mu }_2\parallel (\mathrm{\Lambda }\left[\omega (m+k_1\left(\parallel \delta \parallel +k_{2}r\right))+(\parallel \omega \parallel r+G)k_1{k}_2\right]\hfill \\ \hfill & +\frac{1}{{\lambda }_2}[k_3\left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G+\parallel \omega \parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]\right]),\hfill \end{array}$

and

$\begin{array}{cc}\hfill & [1-(\mathrm{\Lambda }\left[\parallel \omega \parallel (m+k_1(\parallel \delta \parallel +k_{2}r\left)\right)+(\parallel \omega \parallel r+G)k_1{k}_2\right]\hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[k_3\left[\parallel \omega \parallel r+G\right]+\parallel \omega \parallel \left[k_{3}r+\mathrm{\Theta }\right]\right])]\parallel {\mu }_1-{\mu }_2\parallel \le 0,\hfill \end{array}$

which denotes

${\mu }_1\left(\mathfrak{r}\right)={\mu }_2\left(\mathfrak{r}\right).$

  □

2.3. Continuous Dependency on the Control Variable

Here, we study the continuous dependence on the control variable $\nu .$

Definition 1.

The solution of problem (6) and (7) with feedback control (8) depends continuously on $\nu$ if $\forall ϵ>0,\exists \mathrm{{\rm Y}}>0$, where

$\begin{array}{ccc}\hfill |\nu \left(\mathfrak{r}\right)-{\nu }^{*}\left(\mathfrak{r}\right)|=|\varpi (\mathfrak{r},\nu ,\mu )-\varpi (\mathfrak{r},{\nu }^{*},\mu )|& <& \mathrm{{\rm Y}},\mathfrak{r}\in [0,1].\hfill \end{array}$

Then, $\parallel \mu -{\mu }^{*}\parallel <ϵ.$

Now, we prove the result.

Theorem 3.

Suppose that Theorem 5 is verified. Thus, the solution of (13) depends continuously on the control variable $\nu$.

Proof. 

For the two solutions $\mu$ and ${\mu }^{*}$ of (13), corresponding to the control variables $\nu ,{\nu }^{*}$, we get

$\begin{array}{cc}\hfill & |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)|\hfill \\ \hfill & +|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)-\varphi \left(\varrho ,{\nu }^{*}\left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)|]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))\right|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left[\right|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)|\hfill \\ \hfill & +|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)-\varphi \left(\varrho ,{\nu }^{*}\left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)|]d\varrho d\tau \hfill \\ \hfill & +|ϝ(\mathfrak{r},\mu \left(\mathfrak{r}\right))-ϝ(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)|\hfill \\ \hfill & +|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)-\varphi (\varrho ,{\nu }^{*}\left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma )\left|\right]d\varrho d\tau \hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[\right|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))\theta (\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))\theta (\mathfrak{r},\mu \left(\mathfrak{r}\right))\left|\right]\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[m\left(\varrho \right)+k_1\left(\varrho \right)(\left|\varpi (\varrho ,\nu \left(\varrho \right),\mu \left(\varrho \right))\right|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\left|\psi (\varsigma ,\mu \left(\varsigma \right))\right|d\varsigma \left)\right]d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left(k_1\left(\varrho \right)\right(|\nu \left(\varrho \right)-{\nu }^{*}\left(\varrho \right)|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,\mu \left(\varsigma \right))-\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))|d\varsigma \left)d\varrho \right)d\tau +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|\hfill \\ \hfill & \times {\int }_0^{\varrho }\frac{{(\varrho -\tau )}^{1-\gamma }}{\mathrm{\Gamma }(2-\gamma )}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[m\left(\varrho \right)+k_1\left(\varrho \right)(\left|\varpi (\varrho ,\nu \left(\varrho \right),\mu \left(\varrho \right))\right|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }{k}_1\left(\varsigma \right)\left|\psi (\varsigma ,\mu \left(\varsigma \right))\right|d\varsigma \left)\right)]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left(k_1\left(\varrho \right)\right(|\nu \left(\varrho \right)-\nu {\left(\varrho \right)}^{*}|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,{\mu }_1\left(\varsigma \right))-\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))|d\varsigma \left)d\varrho \right)ds\hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[m\left(\varrho \right)+k_1\left(\varrho \right)(\left|\varpi (\varrho ,\nu \left(\varrho \right),\mu \left(\varrho \right))\right|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\left|\psi (\varsigma ,\mu \left(\varsigma \right))\right)\left|d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^s\frac{s}{s+\varrho }\left(k_1\left(\varrho \right)\right(|\nu \left(\varrho \right)-\nu {\left(\varrho \right)}^{*}|\hfill \\ \hfill & +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|\psi (\varsigma ,\mu \left(\varsigma \right))-\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))|d\varsigma \left)d\varrho \right)d\tau \hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))\left|\right|\theta (\mathfrak{r},{\mu }_2\left(\mathfrak{r}\right))-\theta (\mathfrak{r},{\mu }_1\left(\mathfrak{r}\right))|+|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))\left|\right|\theta (\mathfrak{r},\mu \left(\mathfrak{r}\right))|\right]\hfill \\ \hfill & \le \parallel \omega \parallel |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }[\left|m\left(\varrho \right)\right|\hfill \\ \hfill & +|k_1\left(\varrho \right)|\left(\left|\delta \left(\varrho \right)\right|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)\left|\right|\mu \left(\varsigma \right)|d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left(k_1\left(\varrho \right)\left(\gamma +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }{k}_2\left(\varsigma \right)|\mu \left(\varsigma \right)-{\mu }^{*}\left(\varsigma \right)|d\varsigma \right)d\varrho \right)d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\parallel \omega \parallel |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|\hfill \\ \hfill & \times {\int }_0^1{\int }_0^1\frac{\tau }{\tau +\varrho }\left[\left|m\left(\varrho \right)\right|+|k_1\left(\varrho \right)|\left(\left|\delta \left(\varrho \right)\right|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)\left|\right|\mu \left(\varsigma \right)|d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left(|k_1\left(\varrho \right)|\left(\gamma +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)\left|\right|\mu \left(\varsigma \right)-{\mu }^{*}\left(\varsigma \right)|d\varsigma \right)d\varrho \right)d\tau \hfill \\ \hfill & +\parallel \omega \parallel |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|\hfill \\ \hfill & \times {\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left[\left|m\left(\varrho \right)\right|+|k_1\left(\varrho \right)|\left(\left|\delta \left(\varrho \right)\right|+{\int }_0^1\frac{\varrho }{\varrho +\varsigma }|k_2\left(\varsigma \right)\left|\right|\mu \left(\varsigma \right)|d\varsigma \right)\right]d\varrho d\tau \hfill \\ \hfill & +\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\frac{\tau }{\tau +\varrho }\left(k_1\left(\varrho \right)\left(\gamma +{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }{k}_2\left(\varsigma \right)|\mu \left(\varsigma \right)-{\mu }^{*}\left(\varsigma \right)|d\varsigma \right)d\varrho \right)d\tau \hfill \\ \hfill & +\frac{k_3|{\mu }^{*}\left(\mathfrak{r}\right)-\mu \left(\mathfrak{r}\right)|}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]+\parallel \omega \parallel |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|\left[k_3\left|\mu \left(\mathfrak{r}\right)\right|+\mathrm{\Theta }\right]\hfill \\ \hfill & \le \frac{\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]+\frac{\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}(k_1\left(\mathrm{{\rm Y}}+k_2\parallel \mu -{\mu }^{*}\parallel \right)\hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]\hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }^{*}\left(t\right)|+G\right]\left(k_1\left(\mathrm{{\rm Y}}+k_2\parallel \mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\frac{{\eta }^{\gamma }\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]+\frac{{\eta }^{\gamma }\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}\left(k_1\left(\mathrm{{\rm Y}}+k_2|\mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\frac{1}{{\lambda }_1}\left[k_3\parallel \mu -{\mu }^{*}\parallel \left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[k_3\parallel \mu \parallel +\mathrm{\Theta }\right]\right].\hfill \end{array}$

Taking the supremum over $\mathfrak{r}\in I,$ we have

$\begin{array}{cc}\hfill \parallel \mu -{\mu }^{*}\parallel & \le [\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}]\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[m+k_1\left(\parallel \delta \parallel +k_{2}r\right)\right]\hfill \\ \hfill & +[\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2\mathrm{\Gamma }(3-\gamma )}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}]\left[\parallel \omega \parallel r+G\right]\left(k_1\left(\mathrm{{\rm Y}}+k_2\parallel \mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\parallel \mu -{\mu }^{*}\parallel \left[k_3\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel \mu \parallel +\mathrm{\Theta }\right]\right]\hfill \\ \hfill & \le \parallel \mu -{\mu }^{*}\parallel (\mathrm{\Lambda }\left[\parallel \omega \parallel (m+k_1\left(\parallel \delta \parallel +k_{2}r\right))+(\parallel \omega \parallel r+G)k_1{k}_2\right]\hfill \\ \hfill & +\frac{[k_3\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]}{{\lambda }_2})+\left[\parallel \omega \parallel r+G\right]\mathrm{{\rm Y}}\mathrm{\Lambda }{k}_1\hfill \end{array}$

and

$\begin{array}{cc}\hfill \parallel \mu -{\mu }^{*}\parallel & \le \frac{\parallel \omega \parallel r+G]\mathrm{{\rm Y}}\mathrm{\Lambda }{k}_1}{1-\left(\mathrm{\Lambda }\left[\parallel \omega \parallel (m+k_1\left(\parallel \delta \parallel +k_{2}r\right))+(\parallel \omega \parallel r+G)k_1{k}_2\right]+\frac{k_3\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]}{{\lambda }_2}\right)}\hfill \\ \hfill & =ϵ.\hfill \end{array}$

The previous inequality leads to the following result:

$\parallel \mu -{\mu }^{*}\parallel \le ϵ.$

This demonstrates the solution’s continuous dependence on the control variable function $\nu$. □

3. Set-Valued Problem

The study of inclusion problems has drawn much interest based on their extensive applications and actual problems [15,20,21]. Regarding the differential inclusion problems and some results of existence, see [22,23,24,25].

• $\left({\mathcal{H}}_1^{**}\right)$

  Let $\mathrm{\Phi }:I\times \mathbb{R}\times \mathbb{R}\to 2^{\mathbb{R}}$ be non-empty and convex and let subset $\forall (\mathfrak{r},\nu ,\mu )\in I\times \mathbb{R}\times \mathbb{R},$ where

  • (i)

    $\mathrm{\Phi }(\mathfrak{r},.,.)$ is upper semicontinuous in $(\nu ,\mu )\in \mathbb{R}\times \mathbb{R}$, $\forall \mathfrak{r}\in I$.

  • ($ii$)

    $\mathrm{\Phi }(.,\nu ,\mu )$ is measurable in $\mathfrak{r}\in I,$ $\forall (\nu ,\mu )\in \mathbb{R}\times \mathbb{R}$.

  • ($iii$)

    There exist $m,k_1:I\to I$, where $m,k_1\in L^1\left(I\right)$ with

    $\left|\mathrm{\Phi }(\mathfrak{r},\nu ,\tau )\right|=sup\left\{\right|\varphi |:\varphi \in \mathrm{\Phi }(\mathfrak{r},\nu ,\tau )\}\le m\left(\mathfrak{r}\right)+k_1\left(\mathfrak{r}\right)\left|\right(\nu |+\left|\tau \right|),\mathfrak{r}\in I$

    and

    ${\int }_0^1\frac{1}{\tau +\varrho }\left|m\left(\varrho \right)\right|d\varrho \le m,and{\int }_0^1\frac{1}{\tau +\varrho }\left|k_1\left(\varrho \right)\right|d\varrho \le k.$

  $\left({\mathcal{H}}_3^{*}\right)$ Let $\mathrm{\Theta }:I\times \mathbb{R}\to 2^{\mathbb{R}}$ be a Lipschitzian set-valued function with a nonempty compact convex subset of $2^{\mathbb{R}}$, where

  $\parallel \mathrm{\Theta }(\mathfrak{r},{\mu }_1)-\mathrm{\Theta }(\mathfrak{r},{\mu }_2)\parallel \le k_3\left(\mathfrak{r}\right)|{\mu }_1-{\mu }_2|.$

Remark 2.

Obviously, we can deduce that, as shown in Remark 1 [9], there exists a Carathéodory mapping $\varphi \in \mathrm{\Phi }$ [26] that is measurable in $\mathfrak{r}\in I,\forall \nu ,\mu \in \mathbb{R}$ and continuous in ν, $\mu \in \mathbb{R}$, $\forall \mathfrak{r}\in I,$

$\left|\varphi (\mathfrak{r},\nu .\mu )\right|\le m\left(\mathfrak{r}\right)+k_1\left(\mathfrak{r}\right)\left(\right|\nu |+\left|\mu \right|),\mathfrak{r}\in I.$

In addition, the set of Lipschitzian selections $S_{\mathrm{\Theta }}$ is nonempty [26] and $\theta \in S_{\mathrm{\Theta }}$ meets

$|\theta (\mathfrak{r},{\mu }_1)-\theta (\mathfrak{r},{\mu }_2)|\le k_3\left(\mathfrak{r}\right)|{\mu }_1-{\mu }_2|,$

Then,

$\left|\theta (\mathfrak{r},\mu )\right|\le k_1\left(\mathfrak{r}\right)\left|\mu \right|+\mathrm{\Theta },\mathrm{\Theta }=\underset{\mathfrak{r}\in I}{sup}\left|\theta (\mathfrak{r},0)\right|$

and satisfies the nonlocal problem (6) and (7) with feedback control (8).

Let $\mathrm{\Omega }(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)):I\times {\mathbb{R}}^{+}\times {\mathbb{R}}^{+}\to 2^{{\mathbb{R}}^{+}}$ meet the mentioned conditions:

• (a) 

  The set $\mathrm{\Omega }(\mathfrak{r},\nu (\mathfrak{r}),\mu (\mathfrak{r}\left)\right)$ is a non-empty, closed and convex subset for all

  $(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right))\in I\times {\mathbb{R}}^{+}\times {\mathbb{R}}^{+}$.

• (b) 

  $\mathrm{\Omega }(\mathfrak{r},\nu (\mathfrak{r}),\mu (\mathfrak{r}\left)\right)$ is upper semicontinuous in $\nu ,\mu \in {\mathbb{R}}^{+}$ for each $\mathfrak{r}\in I$.

• (c) 

  $\mathrm{\Omega }(.,\nu (\mathfrak{r}),\mu (\mathfrak{r}\left)\right)$ is measurable in $\mathfrak{r}\in I$ for each $\nu ,\mu \in {\mathbb{R}}^{+}$.

• (d) 

  There exists two measurable and bounded functions $\delta :I\to \mathbb{R},$ with norm $\parallel \delta \parallel ,$ where

  $\left|\mathrm{\Omega }\right(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\left)\right|=sup\left\{\right|\eta |:\eta \in \mathrm{\Omega }(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\left)\right\}\le \delta \left(\mathfrak{r}\right),\mathfrak{r}\in I,$

  with $\delta ={max}_{\mathfrak{r}\in I}\{\parallel \delta \parallel \}$.

Remark 3.

We can infer from assumption (i) that the set of selections $S_{\mathrm{\Omega }}$$(i=1,2,\dots ,k)$ of the set-valued function Ω is nonempty and there exists a Carathéodory function $\varpi \in \mathrm{\Omega }$ (see [26]) such that

$\left|\varpi \right(\mathfrak{r},\nu \left(\mathfrak{r}\right),\mu \left(\mathfrak{r}\right)\left)\right|\le \delta \left(\mathfrak{r}\right),$

fulfilling the implicit equation

(17)$\nu \left(\tau \right)=\varpi (\mathfrak{r},\nu (\mathfrak{r}),\mu (\mathfrak{r}\left)\right),\mathfrak{r}\in I.$

Thus, any solution of problem (6) and (7) with multi-valued feedback control (8) is a solution of problem (6) and (7) with feedback control (17).

3.1. Existence Results

Now, based on the main findings in Section 2, we present in the following the results obtained for the nonlocal hybrid modeling of a heat controller (3) via the multi-valued condition (4) with feedback control (5).

Theorem 4.

Let the assumptions $\left({\mathcal{H}}_1^{**}\right)$, $\left({\mathcal{H}}_3^{*}\right)$ and $\left({\mathcal{H}}_2^{*}\right)$, $\left({\mathcal{H}}_4\right)$, $\left({\mathcal{H}}_5\right)$ hold. Then, the problem (3) and (4) has one solution, $\mu \in C(I,\mathbb{R})$.

In the aim of demonstrating uniqueness result of (3)–(5), we replace the assumption $\left({\mathcal{H}}_1\right)$ by

• $\left({\mathcal{H}}_1^{***}\right)$. Let $\mathrm{\Phi }:I\times \mathbb{R}\times \mathbb{R}\to 2^{\mathbb{R}}$ be a Lipschitzian multi-valued mapping with a nonempty convex compact subset of $2^R$, with

  $\parallel \mathrm{\Phi }(\mathfrak{r},{\nu }_1,{\mu }_1)-\mathrm{\Phi }(\mathfrak{r},{\nu }_2,{\mu }_2)\parallel \le k_1\left(\mathfrak{r}\right)\left(\right|{\nu }_1-{\nu }_2|+|{\mu }_1-{\mu }_2\left|\right).$

From these assumptions, we can observe that $\left({\mathcal{H}}_1^{**}\right)$ is held. In addition, the Lipschitzian selection set $S_{\mathrm{\Phi }}$ is nonempty ([26]) and $\varphi \in S_{\mathrm{\Phi }}$ meets

$|\varphi (\mathfrak{r},,{\nu }_1,{\mu }_1)-\varphi (\mathfrak{r},,{\nu }_2,{\mu }_2)|\le k_1\left(\mathfrak{r}\right)\left(\right|{\nu }_1-{\nu }_2|+|{\mu }_1-{\mu }_2\left|\right),$

then, we have

$\left|\varphi (\mathfrak{r},\nu ,\mu )\right|\le k_1\left(\mathfrak{r}\right)\left(\right|\nu |+\left|\mu \right|)+m,m=\underset{\mathfrak{r}\in I}{sup}\left|\varphi (\mathfrak{r},0,0)\right|$

Theorem 5.

Let Theorem 4 be verified and replace assumption $\left({\mathcal{H}}_1^{**}\right)$ with assumption $\left({\mathcal{H}}_1^{***}\right)$. Then, inclusion problem (3)–(5) has a unique solution, $x\in C(I,\mathbb{R})$.

3.2. Continuous Dependency on the Sets of Selections

Definition 2.

The solution of problem (7) and (6), with multi-valued feedback control (17) depends continuously on the set $S_{\mathrm{\Phi }}$. If $\forall ϵ>0,\exists \mathrm{{\rm Y}}>0$, such that

$\begin{array}{ccc}\hfill |\varphi (\mathfrak{r},\nu ,\mu )-{\varphi }^{*}(\mathfrak{r},\nu ,\mu )|& <& \mathrm{{\rm Y}},\varphi ,{\varphi }^{*}\in S_{\mathrm{\Phi }},\mathfrak{r}\in [0,1],\hfill \end{array}$

then, $\parallel \mu -{\mu }^{*}\parallel <ϵ.$

Theorem 6.

Let Theorem 5 be verified. Then, the solution of (13) depends continuously on the set $S_{\mathrm{\Phi }}$ of all Lipschitzian selections of ϕ.

Proof. 

Let the functions $\mu \left(\mathfrak{r}\right)$ and ${\mu }^{*}\left(\mathfrak{r}\right)$ of (13), correspond to $\varphi ,{\varphi }^{*}\in S_{\mathrm{\Phi }}$, respectively; then,

$\begin{array}{cc}\hfill & |\mu \left(\mathfrak{r}\right)-{\mu }^{*}\left(\mathfrak{r}\right)|\hfill \\ \hfill & \le |\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varrho \right)\right|d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -{\varphi }^{*}\left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))\right|{\int }_0^1{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\tau }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -{\varphi }^{*}\left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)]d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }\left|\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\right|d\varrho d\tau \hfill \\ \hfill & +|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\frac{\tau }{\tau +\varrho }[\varphi \left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,\mu \left(\varsigma \right))d\varsigma \right)\hfill \\ \hfill & -{\varphi }^{*}\left(\varrho ,\nu \left(\varrho \right),{\int }_0^{\varrho }\frac{\varrho }{\varrho +\varsigma }\psi (\varsigma ,{\mu }^{*}\left(\varsigma \right))d\varsigma \right)]d\varrho d\tau \hfill \end{array}$

$\begin{array}{cc}\hfill & +\frac{1}{{\lambda }_2}|\mathcal{F}(\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))\theta (\mathfrak{r},{\mu }^{*}\left(\mathfrak{r}\right))-\mathcal{F}(\mathfrak{r},\mu \left(\mathfrak{r}\right))\theta (\mathfrak{r},\mu \left(\mathfrak{r}\right))|\hfill \\ \hfill & \le \parallel \omega \parallel \left|\mu \left(\mathfrak{r}\right)\right)-{\mu }^{*}\left(\mathfrak{r}\right)|{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\left[\left|m\left(\varrho \right)\right|+|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]d\varrho d\tau \hfill \\ \hfill & +\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^{\mathfrak{r}}\frac{{(\mathfrak{r}-\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^{\tau }\left(\gamma +|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)d\varrho \right)d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\parallel \omega \parallel \left|\mu \left(\mathfrak{r}\right)\right)-{\mu }^{*}\left(\mathfrak{r}\right)|{\int }_0^1{\int }_0^{\mathfrak{r}}\left[\left|m\left(\varrho \right)\right|+|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel {\mu }_1\parallel \right)\right]d\varrho d\tau \hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^1{\int }_0^{\tau }\left(\mathrm{{\rm Y}}+|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)d\varrho \right)d\tau \hfill \\ \hfill & +\parallel \omega \parallel \left|\mu \left(\mathfrak{r}\right)\right)-{\mu }^{*}\left(\mathfrak{r}\right)|{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\left[\left|m\left(\varrho \right)\right|+|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]d\varrho d\tau \hfill \\ \hfill & +\left[\parallel \omega \parallel |{\mu }^{*}\left(\mathfrak{r}\right)|+G\right]{\int }_0^{\eta }\frac{{(\eta -\tau )}^{\gamma -1}}{\mathrm{\Gamma }\left(\gamma \right)}{\int }_0^1\left(\mathrm{{\rm Y}}+|k_1\left(\varrho \right)|\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)d\varrho \right)d\tau \hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[k_3\parallel \mu -{\mu }^{*}\parallel \left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[k_3\parallel \mu \parallel +\mathrm{\Theta }\right]\right]\hfill \\ \hfill & \le \frac{\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]+\frac{\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}(\mathrm{{\rm Y}}+k_1\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)\hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]\hfill \\ \hfill & +\frac{{\lambda }_1}{{\lambda }_2}\left[\parallel \omega \parallel |{\mu }^{*}\left(t\right)|+G\right]\left(\mathrm{{\rm Y}}+k_1\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\frac{{\eta }^{\gamma }\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel }{\mathrm{\Gamma }(\gamma +1)}\left[m+k_1\left(\parallel \delta \parallel +k_2\parallel \mu \parallel \right)\right]+\frac{{\eta }^{\gamma }\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]}{\mathrm{\Gamma }(\gamma +1)}\left(\mathrm{{\rm Y}}+k_1\left(\parallel \delta \parallel +k_2|\mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\frac{1}{{\lambda }_2}\left[k_3\parallel \mu -{\mu }^{*}\parallel \left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[k_3\parallel \mu \parallel +\mathrm{\Theta }\right]\right].\hfill \end{array}$

For $\mathfrak{r}\in [0,1]$, we obtain

$\begin{array}{cc}\hfill \parallel \mu -{\mu }^{*}\parallel & \le \left[\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}\right]\parallel \omega \parallel \parallel \mu -{\mu }^{*}\parallel \left[m+k_1\left(\parallel \delta \parallel +k_{2}r\right)\right]\hfill \\ \hfill & +\left[\frac{1}{\mathrm{\Gamma }(\gamma +1)}+\frac{{\lambda }_1}{{\lambda }_2\mathrm{\Gamma }(3-\gamma )}+\frac{{\eta }^{\gamma }}{\mathrm{\Gamma }(\gamma +1)}\right]\left[\parallel \omega \parallel r+G\right]\left(\mathrm{{\rm Y}}+k_1\left(\parallel \delta \parallel +k_2\parallel \mu -{\mu }^{*}\parallel \right)\right)\hfill \\ \hfill & +\frac{\parallel \mu -{\mu }^{*}\parallel }{{\lambda }_2}\left[k_3\left[\parallel \omega \parallel \parallel {\mu }^{*}\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel \mu \parallel +\mathrm{\Theta }\right]\right]\hfill \\ \hfill & \le \parallel \mu -{\mu }^{*}\parallel (\mathrm{\Lambda }\left[\omega (m+k_1\left(\parallel \delta \parallel +k_{2}r\right))+(\parallel \omega \parallel r+G)\left(\mathrm{{\rm Y}}+k_1\left(\parallel \delta \parallel +k_2\right)\right)\right]\hfill \\ \hfill & +\frac{\left[k_3\left[\parallel \omega \parallel \parallel {\mu }_2\parallel +G\right]+\parallel \omega \parallel \left[k_3\parallel {\mu }_1\parallel +\mathrm{\Theta }\right]\right]}{{\lambda }_2})+\parallel \omega \parallel r+G][\mathrm{{\rm Y}}+k_1\parallel \delta \parallel ]\mathrm{\Lambda }\hfill \end{array}$

and

$\begin{array}{cc}\hfill \parallel \mu -{\mu }^{*}\parallel & \le \frac{\left[\parallel \omega \parallel r+G\right][\mathrm{{\rm Y}}+k_1\parallel \delta \parallel ]\mathrm{\Lambda }}{1-\left(\mathrm{\Lambda }[\parallel \omega \parallel (m+k_1(\parallel \delta \parallel +k_{2}r\left)\right)+(\parallel \omega \parallel r+G)k_1{k}_2]+\frac{[k_3[\parallel \omega \parallel r+G]+\parallel \omega \parallel [k_{3}r+\mathrm{\Theta }]]}{{\lambda }_2}\right)}\hfill \\ \hfill & =ϵ\hfill \end{array}$

The previous inequality leads to the following result:

$\parallel \mu -{\mu }^{*}\parallel \le ϵ.$

This demonstrates the solution of the problem (7) and (6) with multi-valued feedback control (17) is continuously dependent on the set $S_{\mathrm{\Phi }}$. □

3.3. Example

In accordance with the indicated hybrid BVP (3) and (4), we take into account the fractional order hybrid inclusion problem

(18)$\begin{array}{cc}\hfill {-}^c{\mathcal{D}}^{\gamma }\left(\frac{\mu \left(\mathfrak{r}\right)}{\frac{e^{-l{n}^2(\mathfrak{r}+1)}\left|\mu \left(\mathfrak{r}\right)\right|}{1+\left|\mu \right(\mathfrak{r}\left)\right|}+8}\right)& \in \left[{\int }_0^{\mathfrak{r}}\frac{\tau }{\mathfrak{r}+\tau }[\tau +\frac{{\int }_0^{\tau }\nu \left(\tau \right)\sin \frac{\left|\mu \right(\varrho \left)\right|}{1+\left|\mu \right(\varrho \left)\right|}d\varrho }{2(1+{\int }_0^{\tau }\nu \left(\tau \right)\sin \frac{\left|\mu \right(\varrho \left)\right|}{1+\left|\mu \right(\varrho \left)\right|}d\varrho )}]d\tau ,0\right],\mathfrak{r}\in I\hfill \end{array}$

via the nonlocal conditions

(19)$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{D}\left(\frac{\mu \left(\mathfrak{r}\right)}{\frac{e^{-l{n}^2(\mathfrak{r}+1)}\left|\mu \right|}{100(1+|\mu \left|\right)}+8}\right){|}_{\mathfrak{r}=0}=0,\\ {\lambda }_1{}^c{\mathcal{D}}^{\gamma -1}\left(\frac{\mu \left(\mathfrak{r}\right)}{\frac{e^{-l{n}^2(\mathfrak{r}+1)}\left|\mu \right|}{100(1+|\mu \left|\right)}+8}\right){|}_{\mathfrak{r}=1}+{\lambda }_2\left(\frac{\mu \left(\mathfrak{r}\right)}{\frac{e^{-l{n}^2(\mathfrak{r}+1)}\left|\mu \right|}{100(1+|\mu \left|\right)}+8}\right){|}_{\mathfrak{r}=\eta }\in \left[\frac{\mathfrak{r}}{2}+\frac{\mu \left(\mathfrak{r}\right)}{4+\mathfrak{r}},0\right],\end{array}\right.\end{array}$

with multi-valued feedback

(20)$\mu \left(\mathfrak{r}\right)\in [0.1\mu \left(\mathfrak{r}\right)+\frac{1}{200}\cos \left(\mathfrak{r}\right)+e^{\frac{-3}{2}\mathfrak{r}}\nu \left(\mathfrak{r}\right),0].$

Let $\gamma =\frac{7}{4},\gamma -1=\frac{3}{4},\eta =0.89$, and ${\lambda }_1={\lambda }_2=\frac{9}{5}$.

Define $g(\mathfrak{r},\mu \left(\mathfrak{r}\right))=\frac{e^{-l{n}^2(\mathfrak{r}+1)}\left|u\right|}{100(1+|u\left|\right)}+8$ and the multi-valued map $\mathrm{\Phi }$ by

$\mathrm{\Phi }(\mathfrak{r},\nu ,\mu )=\left[\mathfrak{r}+\frac{{\int }_0^{\mathfrak{r}}\nu \left(\mathfrak{r}\right)\sin \frac{\left|\mu \right(\mathfrak{r}\left)\right|}{1+\left|\mu \right(\mathfrak{r}\left)\right|}d\tau }{2(1+{\int }_0^{\mathfrak{r}}\nu \left(\mathfrak{r}\right)\sin \frac{\left|\mu \right(\tau \left)\right|}{1+\left|\mu \right(\tau \left)\right|}d\tau )},0\right]\mathrm{and}\theta (\mathfrak{r},x\left(\mathfrak{r}\right))=\frac{\mathfrak{r}}{2}+\frac{\mu \left(\mathfrak{r}\right)}{4+\mathfrak{r}}.$

If $w\left(\mathfrak{r}\right)=\frac{e^{-l{n}^2(\mathfrak{r}+1)}}{100}$, then $\parallel w\parallel =\frac{1}{100}$, $m=0.5$ $k=\frac{1}{2}$, $k^{*}=\frac{1}{2}$, $k_3=\frac{1}{4}$, $G=\frac{1}{2},$ and $\theta =\frac{1}{2}$. By using the above relations, we get $\mathrm{\Lambda }\simeq 1.9468834$. Hence, the above data satisfy the condition of Theorem 2:

$\begin{array}{ccc}\hfill \mathrm{\Lambda }[\parallel \omega \parallel (m+k_1(\parallel \delta \parallel +k_{2}r\left)\right)& +& (\parallel \omega \parallel r+G)k_1{k}_2]+\frac{1}{{\lambda }_2}[k_3\left[\parallel \omega \parallel r+G\right]\hfill \\ & +& \parallel \omega \parallel \left[k_{3}r+\mathrm{\Theta }\right]]\simeq 0.5826955222<1.\hfill \end{array}$

Using Theorem 5, then the problem (18) and (19) with multi-valued feedback (20) has a unique solution.

4. Conclusions

Many works in the literature and monographs have treated and developed mathematical models that appear in various real-world applications, for example, thermostats or heat controllers. One approach is to develop very complex iterations of popular models from real-world issues which can be described using inclusions or fractional differential equations [2,7,12,16].

In this work, we provide a comprehensive investigation of a class of hybrid fractional models of thermostats via nonlocal multi-valued boundary conditions (3) and (4) which satisfy multi-valued feedback control. The main tool of our study is applying Dhage’s hybrid fixed point theorem [19]. The use of various approaches for certain differential and integral problems, including constraints or control variables, has recently been developed by several scholars, for example, in refs. [27,28,29,30,31,32,33]. This feedback control may be in an implicit form as in [27,28,29,30], multi-valued feedback control as in [32], or fractal feedback control [33].

We have established the continuous dependence of the unique solution of our problem on the control variable and on the set $S_{\mathrm{\Phi }}.$ In this study, we have investigated some qualitative properties of the solution of this problem, which encourages us to investigate and discuss additional singular dynamical systems that appear in a variety of natural and engineering phenomena.

Author Contributions

Methodology, S.M.A.-I., A.M.A.E.-S. and H.H.G.H.; Formal analysis, S.M.A.-I., A.M.A.E.-S. and H.H.G.H.; Investigation, S.M.A.-I., A.M.A.E.-S. and H.H.G.H.; Writing—review & editing, S.M.A.-I., A.M.A.E.-S. and H.H.G.H. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the referees for their remarks and comments that helped to improve our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

References

1. Al-Issa, S.M.; El-Sayed, A.M.A.; Hashem, H.H.G. On Chandrasekhar hybrid Caputo fractional modeling for thermostat via hybrid boundary value conditions in Banach algebra. Proceedings of the 2023 International Conference on Fractional Differentiation and Its Applications (ICFDA); Roma, Italy, 11–15 September 2023.

2. Cahlon, B.; Schmidt, D.; Shillor, M.; Zou, X. Analysis of thermostat models. Eur. J. Appl. Math.; 1997; 8, pp. 437-457. [DOI: https://dx.doi.org/10.1017/S0956792597003240]

3. Zou, X.; Jordan, J.A.; Shillor, M. A dynamic model for a thermostat. J. Eng. Math.; 1999; 36, pp. 291-310. [DOI: https://dx.doi.org/10.1023/A:1004587425961]

4. Webb, J.R.L. Multiple positive solutions of some nonlinear heat flow problems. Conf. Publ.; 2005; 2005, pp. 895-903.

5. Webb, J.R.L. Existence of positive solutions for a thermostat model. Nonlinear Anal. Real World Appl.; 2012; 13, pp. 923-938. [DOI: https://dx.doi.org/10.1016/j.nonrwa.2011.08.027]

6. Shen, C.; Zhou, H.; Yang, L. Existence and nonexistence of positive solutions of a fractional thermostat model with a parameter. Math. Methods Appl. Sci.; 2016; 39, pp. 4504-4511. [DOI: https://dx.doi.org/10.1002/mma.3878]

7. Liang, Y.; Levine, D.I.; Shen, Z.J. Thermostats for the smart grid: Models, benchmarks and insights. Energy J.; 2012; 33, pp. 61-95. [DOI: https://dx.doi.org/10.5547/01956574.33.4.4]

8. Aydi, H.; Jleli, M.; Samet, B. On Positive Solutions for a Fractional Thermostat Model with a Convex–Concave Source Term via ψ-Caputo Fractional Derivative. Mediterr. J. Math.; 2020; 17, 16. [DOI: https://dx.doi.org/10.1007/s00009-019-1450-7]

9. El-Sayed, A.M.A.; Hashem, H.H.G.; Al-Issa, S.M. Analysis of a hybrid integro-differential inclusion. Bound. Value Probl.; 2022; 2022, 68. [DOI: https://dx.doi.org/10.1186/s13661-022-01650-w]

10. Rezapour, S.; Etemad, S.; Agarwal, R.P.; Nonlaopon, K. On a Lyapunov-Type Inequality for Control of a y-Model Thermostat and the Existence of Its Solutions. Mathematics; 2022; 10, 4023. [DOI: https://dx.doi.org/10.3390/math10214023]

11. Baleanu, D.; Etemad, S.; Rezapour, S. A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions. Bound. Value Probl.; 2020; 2020, 64. [DOI: https://dx.doi.org/10.1186/s13661-020-01361-0]

12. Mohammadi, H.; Kumar, S.; Rezapour, S.; Etemad, S. A theoretical study of the Caputo–Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control. Chaos Solitons Fractals; 2021; 144, 110668. [DOI: https://dx.doi.org/10.1016/j.chaos.2021.110668]

13. Kontes, G.D.; Giannakis, G.I.; Horn, P.; Steiger, S.; Rovas, D.V. Using thermostats for indoor climate control in office buildings: The effect on thermal comfort. Energies; 2017; 10, 1368. [DOI: https://dx.doi.org/10.3390/en10091368]

14. Dhage, B.C.; Lakshmikantham, V. Basic results on hybrid differential equation. Nonlinear Anal. Hybrid Syst.; 2010; 4, pp. 414-424. [DOI: https://dx.doi.org/10.1016/j.nahs.2009.10.005]

15. Baleanu, D.; Etemad, S.; Rezapour, S. On a fractional hybrid integro-differential equation with mixed hybrid integral boundary value conditions by using three operators. Alex. Eng. J.; 2020; 59, pp. 3019-3027.

16. Caballero, J.; Mingarelli, A.B.; Sadarangani, K. Existence of solutions of an integral equation of Chandrasekhar type in the theory of radiative transfer. Electron. J. Differ. Equ.; 2006; 2006, pp. 1-11.

17. Busbridge, W. On Solutions of Chandrasekhar’s Integral Equation. J. Am. Math. Soc.; 1962; 105, pp. 112-117.

18. Dhage, B.C. A fixed point theorem in Banach algebras involving three operators with applications. Kyungpook Math. J.; 2004; 44, pp. 145-155.

19. Dhage, B.C. Nonlinear Functional Boundary Value Problems in Banach Algebras Involving Carathéodories. Kyungpook Math. J.; 2006; 46, pp. 427-441.

20. Baleanu, D.; Etemad, S.; Pourrazi, S.; Rezapour, S. On the new fractional hybrid boundary value problems with three-point integral hybrid conditions. Adv. Differ. Equ.; 2019; 2019, 473. [DOI: https://dx.doi.org/10.1186/s13662-019-2407-7]

21. Matar, M.M.; Abbas, M.I.; Alzabut, J.; Kaabar, M.K.A.; Etemad, S.; Rezapour, S. Investigation of the p-Laplacian nonperiodic nonlinear boundary value problem via generalized Caputo fractional derivatives. Adv. Differ. Equ.; 2021; 2021, 68. [DOI: https://dx.doi.org/10.1186/s13662-021-03228-9]

22. Piazza, L.D.; Marraffa, V.; Satco, B. Approximating the solutions of differential inclusions driven by measures. Ann. Mat. Pur Appl. 1923; 2019; 198, pp. 2123-2140. [DOI: https://dx.doi.org/10.1007/s10231-019-00857-6]

23. Satco, B. Continuous dependence results for set-valued measure differential problems. Electron. Qual. Theory Differ. Equ.; 2015; 79, pp. 1-15. [DOI: https://dx.doi.org/10.14232/ejqtde.2015.1.79]

24. Satco, B. Existence results for Urysohn integral inclusions involving the Henstock integral. J. Math. Anal. Appl.; 2007; 336, pp. 44-53. [DOI: https://dx.doi.org/10.1016/j.jmaa.2007.02.050]

25. Baleanu, D.; Hedayati, V.; Rezapour, S.; Al Qurash, M. On two fractional differential inclusions. SpringerPlus; 2016; 5, 882. [DOI: https://dx.doi.org/10.1186/s40064-016-2564-z]

26. Aubin, J.P.; Cellina, A. Differential Inclusions: Set-Valued Maps and Viability Theory; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 264.

27. Chen, F. The permanence and global attractivity of Lotka–Volterra competition system with feedback controls. Nonlinear Anal. Real World Appl.; 2006; 7, pp. 133-143. [DOI: https://dx.doi.org/10.1016/j.nonrwa.2005.01.006]

28. Nasertayoob, P. Solvability and asymptotic stability of a class of nonlinear functional-integral equation with feedback control. Commun. Nonlinear Anal.; 2018; 5, pp. 19-27.

29. Nasertayoob, P.; Vaezpour, S.M. Positive periodic solution for a nonlinear neutral delay population equation with feedback control. J. Nonlinear Sci. Appl.; 2013; 6, pp. 152-161. [DOI: https://dx.doi.org/10.22436/jnsa.007.03.08]

30. El-Sayed, A.M.A.; Hashem, H.H.G.; Al-Issa, S.M. Analytical Contribution to a Cubic Functional Integral Equation with Feedback Control on the Real Half Axis. Mathematics; 2023; 11, 1133. [DOI: https://dx.doi.org/10.3390/math11051133]

31. El-Sayed, A.M.A.; Hamdallah, E.A.; Ahmed, R.G. On a nonlinear constrained problem of a nonlinear functional integral equation. Appl. Anal. Optim.; 2022; 6, pp. 95-107.

32. El-Sayed, A.M.A.; Hashem, H.H.G.; Al-Issa, S.M. New Aspects on the Solvability of a Multidimensional Functional Integral Equation with Multivalued Feedback Control. Axioms; 2023; 12, 653. [DOI: https://dx.doi.org/10.3390/axioms12070653]

33. Hashem, H.H.G.; El-Sayed, A.M.A.; Al-Issa, S.M. Investigating Asymptotic Stability for Hybrid Cubic Integral Inclusion with Fractal Feedback Control on the Real Half-Axis. Fractal Fract.; 2023; 7, 449. [DOI: https://dx.doi.org/10.3390/fractalfract7060449]