Skip to main content

Theory and Modern Applications

On nonlinear pantograph fractional differential equations with Atangana–Baleanu–Caputo derivative

Abstract

In this paper, we obtain sufficient conditions for the existence and uniqueness results of the pantograph fractional differential equations (FDEs) with nonlocal conditions involving Atangana–Baleanu–Caputo (ABC) derivative operator with fractional orders. Our approach is based on the reduction of FDEs to fractional integral equations and on some fixed point theorems such as Banach’s contraction principle and the fixed point theorem of Krasnoselskii. Further, Gronwall’s inequality in the frame of the Atangana–Baleanu fractional integral operator is applied to develop adequate results for different kinds of Ulam–Hyers stabilities. Lastly, the paper includes an example to substantiate the validity of the results.

1 Introduction

Fractional calculus (FC) has been growing quicker during the most recent few years, and numerous phenomena having the power-law impact have been described precisely with fractional models [19]. Numerous outstanding results of the fractional models have been acquired in different fields of science and engineering. One of the specificities of the FC is that we have numerous fractional derivatives (FDs) that offer the authors the chance to pick the specific FD which coincides better with a given real-world problem. The description of phenomena with memory effect is as yet a major test for the specialists. Along these lines, new tools and methods ought to be made to have the option to show a better improvement description of real-world phenomena and the existing models. In this regard, it appears that there is a need for new FDs with the nonsingular kernel. For the nonlocal FDs with the nonsingular exponential kernel, we allude to [10, 11], and for other local approaches of the FDs, we allude to the recent works [12, 13]. Probably the best competitor among the current kernels is the one dependent on Mittag-Leffler functions (MLF) [14]. In view of this, very lately a novel FD [14] (ABC fractional operators) was structured and applied to sundry real-world problems [15, 16]. Then, in [17, 18], the authors deliberated the discrete versions of those new operators. For the modeling and important applications in the frame of ABC fractional operator, see [1926]. Recent investigations of the existence and uniqueness of solutions for fractional differential equations (FDEs) of the impulsive, evolution, and functional problems with initial or boundary conditions can be found within the following research series [2730] and the references therein. Recent contributions on FDEs involving ABC-FDs can be found in the articles [13, 3138].

On the other hand, the pantograph is an apparatus employed in electric trains to collect electric currents from the overload lines. This type of equation was designed by Ockendon and Tayler [39]. Pantograph equations play a pivotal role in pure and applied mathematics and physics. Motivated by their significance, a ton of scientists generalized these equations into different types and presented the solvability aspect of such problems both numerically and theoretically; for additional subtleties, see [4046] and the references therein. Besides, some authors applied various kinds of fractional derivatives and studied the existence and stability of Ulam–Hyers, which can be found in [4751]. However, not many works have been proposed for pantograph FDEs, especially those involving ABC fractional operator and nonlocal conditions.

Motivated by the above argumentations, the intent of this work is to investigate the ABC- type pantograph FDEs with nonlocal conditions described by

$$ \textstyle\begin{cases} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\varsigma (\mathfrak{r})=f( \mathfrak{r},\varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r})),& \mathfrak{r}\in [ a,T], 0< \vartheta \leq 1, \\ \varsigma (a)=\sum_{k=1}^{m}c_{k}\varsigma (\mathfrak{r}_{k}),& \mathfrak{r}_{k}\in (a,T), \end{cases} $$
(1.1)

where \(0<\gamma <1\), \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\) is the AB-Caputo FD of order ϑ, \(f:[a,T]\times \mathbb{R} \times \mathbb{R} \rightarrow \mathbb{R} \) is a continuous function with \(f(a,\varsigma (a),\varsigma (\gamma a))=0\), and the constant \(c_{k}\) satisfies the condition \(\sum_{k=1}^{m}c_{k}\neq 1\). Note that, if \(\gamma =1\), then our problem reduces to

$$ \textstyle\begin{cases} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\varsigma (\mathfrak{r})=f( \mathfrak{r},\varsigma (\mathfrak{r})),& \mathfrak{r}\in [ a,T],0< \vartheta \leq 1, \\ \varsigma (a)=\sum_{k=1}^{m}c_{k}\varsigma (\mathfrak{r}_{k}),& \mathfrak{r}_{k}\in (a,T). \end{cases} $$
(1.2)

Therefore, if \(\gamma =1\), the results acquired in the present paper are also true for ABC-type pantograph FDEs (1.2).

Some fixed point theorems are applied to establish the existence and uniqueness of solution. The Ulam–Hyers stabilities are proved via Gronwall’s inequality in the frame of AB fractional integral operator. The proposed problems are more generalized, also the obtained results are recent studies and an extension of the development of FDEs involving this new operator. Moreover, the analysis of the results was limited to the minimum assumptions.

This paper is formatted as follows. Section 2 provides the background materials and preliminaries required for our analysis. Section 3 is devoted to obtaining a formula of solution to the ABC type pantograph FDEs (1.1). In Sect. 4, we prove the existence and uniqueness of solution to problems at hand by means of some techniques of FPTs. In Sect. 5, the Ulam–Hyers and generalized Ulam–Hyers stability of the pantograph ABC-FDEs (1.1) is discussed via Gronwall’s inequality in the frame of AB fractional integral operator. Finally, an illustrative example is offered in Sect. 6.

2 Background materials and preliminaries

Here, we recollect some requisite definitions and preliminary concepts related to our work.

Let \(\mathfrak{Z}=[a,T]\), \(\mathfrak{Z}^{\prime }=(a,T)\subset \mathbb{R} \), and \(\mathfrak{D}=C(\mathfrak{Z},\mathbb{R} )\) be the space of continuous functions \(\varsigma :\mathfrak{Z}\rightarrow \mathbb{R} \) with the norm

$$ \Vert \varsigma \Vert =\max \bigl\{ \bigl\vert \varsigma ( \mathfrak{r}) \bigr\vert :\mathfrak{r}\in \mathfrak{Z}\bigr\} . $$

Clearly, \(\mathfrak{D}\) is a Banach space with the norm \(\Vert \varsigma \Vert \).

Definition 2.1

([14])

Let \(\vartheta \in (0,1]\) and \(\mathfrak{p}\in H^{1}(\mathfrak{Z})\). Then the AB-Caputo and AB-Riemann–Liouville FDs of order ϑ for a function \(\mathfrak{p}\) are described by

$$ {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\mathfrak{p}( \mathfrak{r})= \frac{\mathfrak{N}(\vartheta )}{1-\vartheta } \int _{a}^{\mathfrak{r}} \mathbb{E}_{\vartheta } \biggl( \frac{-\vartheta }{\vartheta -1}(\mathfrak{r}- \mathfrak{s})^{\vartheta } \biggr) \mathfrak{p}^{\prime }(\mathfrak{s})\,d \mathfrak{s}, \quad \mathfrak{r}>a, $$

and

$$ {}^{ABR}\mathbb{D}_{a^{+}}^{\vartheta }\mathfrak{p}( \mathfrak{r})= \frac{\mathfrak{N}(\vartheta )}{1-\vartheta }\frac{d}{d\mathfrak{r}} \int _{a}^{\mathfrak{r}}\mathbb{E}_{\vartheta } \biggl( \frac{-\vartheta }{\vartheta -1}(\mathfrak{r}-\mathfrak{s})^{\vartheta } \biggr) \mathfrak{p}( \mathfrak{s})\,d\mathfrak{s}, \quad \mathfrak{r}>a, $$

respectively, where \(\mathfrak{N}(\vartheta )>0\) is a normalization function complying with \(\mathfrak{N}(0)=\mathfrak{N}(1)=1\), and \(\mathbb{E}_{\vartheta }\) is called the Mittag-Leffler function described by

$$ \mathbb{E}_{\vartheta } ( \mathfrak{p} ) =\sum_{k=0}^{ \infty } \frac{\mathfrak{p}^{k}}{\Gamma ( k \vartheta +1 ) }, \quad Re ( \vartheta ) >0, \mathfrak{p}\in \mathbb{C} . $$

The associated AB fractional integral is specified by

$$ {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }\mathfrak{p}( \mathfrak{r})= \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\mathfrak{p}(\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \mathfrak{p}(\mathfrak{r}), $$

where

$$ I_{a^{+}}^{\vartheta }\mathfrak{p}(\mathfrak{r})= \frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}}(\mathfrak{r}-\mathfrak{s})^{\vartheta -1} \mathfrak{p}(\mathfrak{s})\,d\mathfrak{s}. $$

Lemma 2.1

([17])

Let \(\vartheta \in (0,1]\) and \(\mathfrak{p}\in H^{1}(\mathfrak{Z})\), if AB-Caputo FD exists, then we have

$$ {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta } {}^{ABC} \mathbb{D}_{a^{+}}^{\vartheta }\mathfrak{p}(\mathfrak{r})= \mathfrak{p}( \mathfrak{r})-\mathfrak{p}(a). $$

Definition 2.2

([31])

The relation between the AB-Riemann–Liouville and AB-Caputo FDs is given by

$$ {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\mathfrak{p}( \mathfrak{r})= {}^{ABR}\mathbb{D}_{a^{+}}^{\vartheta } \mathfrak{p}(\mathfrak{r})- \frac{\mathfrak{N}(\vartheta )}{1-\vartheta }\mathfrak{p}(a)\mathbb{E}_{ \vartheta } \biggl( \frac{-\vartheta }{\vartheta -1}(\mathfrak{r}-a)^{ \vartheta } \biggr) . $$
(2.1)

Remark 2.1

Replacing \(\mathfrak{p}(\mathfrak{r})\) with \({}^{AB}\mathbb{I}_{a^{+}}^{\vartheta } {\mathfrak{p}(\mathfrak{r})}\) in Definition 2.2 and using Lemma 2.1, it can be shown that

$$ {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } {}^{AB} \mathbb{I}_{a^{+}}^{ \vartheta }\mathfrak{p}(\mathfrak{r})= \mathfrak{p}(\mathfrak{r})-\mathfrak{p}(a) \mathbb{E}_{\vartheta } \biggl( \frac{-\vartheta }{\vartheta -1}( \mathfrak{r}-a)^{\vartheta } \biggr). $$

Hence, under the condition that \(\mathfrak{p}(a)=0\), we get the identity

$$ {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } {}^{AB} \mathbb{I}_{a^{+}}^{ \vartheta }\mathfrak{p}(\mathfrak{r})= \mathfrak{p}(\mathfrak{r}). $$

Lemma 2.2

([14])

Let \(\vartheta >0\). Then \({}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }\) is bounded from \(\mathfrak{D}\) into \(\mathfrak{D}\).

Lemma 2.3

([13, 17])

Let \(\vartheta \in (0,1]\) and \(\varpi \in \mathfrak{D}\) with \(\varpi (a)=0\). Then the solution of the following problem

$$\begin{aligned}& {}^{ABC}D_{a^{+}}^{\vartheta }\mathfrak{p}(\mathfrak{r}) = \varpi ( \mathfrak{r}), \quad \mathfrak{r}\in \mathfrak{Z}, \\& \mathfrak{p}(a) = \mathfrak{p}_{a}, \end{aligned}$$

is given by

$$ \mathfrak{p}(\mathfrak{r})=\mathfrak{p}_{a}+ \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}}(\mathfrak{r}- \mathfrak{s})^{\vartheta -1} \varpi (\mathfrak{s})\,d\mathfrak{s}. $$

Theorem 2.1

([52], Banach’s contraction principle)

Let \(\mathcal{J}\) be a Banach space, and \(\mathfrak{K}\) be a nonempty closed subset of \(\mathcal{J}\). If \(\mathfrak{B}:\mathfrak{K}\longrightarrow \mathfrak{K}\) is a contraction, then there exists a unique fixed point of \(\mathfrak{B}\).

Theorem 2.2

([52], Krasnoselskii’s fixed point theorem)

Let \(\mathfrak{K}\) be a nonempty, closed, convex subset of a Banach space \(\mathcal{J}\). Let \(\mathfrak{B}_{1}\), \(\mathfrak{B}_{2}\) be two operators such that (i) \(\mathfrak{B}_{1}u+\mathfrak{B}_{2}v\in \mathfrak{K}\), \(\forall u,v\in \mathfrak{K}\); (ii) \(\mathfrak{B}_{1}\) is compact and continuous; (iii) \(\mathfrak{B}_{2}\) is a contraction mapping. Then there exists \(w\in \mathfrak{K}\) such that \(\mathfrak{B}_{1}w+\mathfrak{B}_{2}w=w\).

Theorem 2.3

([13], Generalized Gronwall’s inequality)

Suppose that \(0<\vartheta \leq 1\), \(a(\mathfrak{r}) ( 1- \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}b(\mathfrak{r}) ) ^{-1}\) is a nonnegative, nondecreasing, and locally integrable function on \([c,d)\), \(\frac{\vartheta b(t)}{\mathfrak{N}(\vartheta )} ( 1- \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}b(\mathfrak{r}) ) ^{-1}\) is nonnegative and bounded on \([c,d)\), and \(\sigma (\mathfrak{r})\) is nonnegative and locally integrable on \([c,d)\) with

$$ \sigma (\mathfrak{r})\leq a(\mathfrak{r})+b(\mathfrak{r}) \bigl( {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }\sigma \bigr) (\mathfrak{r}). $$

Then

$$ \sigma (\mathfrak{r})\leq \frac{a(\mathfrak{r})\mathfrak{N}(\vartheta )}{\mathfrak{N}(\vartheta )- ( 1-\vartheta ) b(\mathfrak{r})} \mathbb{E}_{\vartheta } \biggl( \frac{\vartheta b(\mathfrak{r})(\mathfrak{r}-a)^{\vartheta }}{\mathfrak{N}(\vartheta )- ( 1-\vartheta ) b(\mathfrak{r})} \biggr) . $$

3 Formulas of solution

This section is devoted to obtaining formulas of solution to linear problems corresponding to (1.1).

Theorem 3.1

Let \(0<\vartheta \leq 1\), \(\sum_{k=1}^{m}c_{k}\neq 1\), and let \(\varpi \in \mathfrak{D}\) with \(\varpi (a)=0\). A function \(\varsigma \in \mathfrak{D}\) is a solution of the fractional integral equation (FIE)

$$\begin{aligned} \varsigma (\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \varpi ( \mathfrak{r}) \end{aligned}$$
(3.1)

if and only if ς is a solution of the ABC-problem

$$ \begin{gathered} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } \varsigma (\mathfrak{r})=\varpi (\mathfrak{r}),\quad \mathfrak{r}\in \mathfrak{Z}, \\ \varsigma (a)=\sum_{k=1}^{m}c_{k} \varsigma (\mathfrak{r}_{k}),\quad \mathfrak{r}_{k} \in \mathfrak{Z}^{\prime }.\end{gathered} $$
(3.2)

Proof

Assume that ς satisfies the first equation of (3.2). From Lemma 2.3, we have

$$ \varsigma (\mathfrak{r})=\varsigma (a)+ \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}}(\mathfrak{r}- \mathfrak{s})^{\vartheta -1} \varpi (\mathfrak{s})\,d\mathfrak{s}. $$
(3.3)

Now, if we replace \(\mathfrak{r}=\mathfrak{r}_{k}\) and multiply both sides by \(c_{k}\) in (3.3), we get

$$ c_{k}\varsigma (\mathfrak{r}_{k})=c_{k} \varsigma (a)+ \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}c_{k}\varpi (\mathfrak{r}_{k})+ \frac{\vartheta c_{k}}{\mathfrak{N}(\vartheta )}\frac{1}{\Gamma (\vartheta )} \int _{a}^{ \mathfrak{r}_{k}}(\mathfrak{r}_{k}- \mathfrak{s})^{\vartheta -1}\varpi ( \mathfrak{s})d\mathfrak{s}. $$
(3.4)

From the nonlocal condition, we get

$$\begin{aligned} \varsigma (a) =&\sum_{k=1}^{m}c_{k} \varsigma (\mathfrak{r}_{k}) \\ =&\sum_{k=1}^{m}c_{k}\varsigma (a)+ \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k} \varpi (\mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}), \end{aligned}$$

which implies

$$ \varsigma (a)=\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{ \vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) . $$
(3.5)

By matching the two equations (3.3) and (3.5), we get

$$\begin{aligned} \varsigma (\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \varpi ( \mathfrak{r}). \end{aligned}$$

Thus (3.1) is satisfied.

Conversely, suppose that ς satisfies equation (3.1). Applying \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\) on both sides of (3.1), then using Remark 2.1, and from the fact \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }(k)=0\), for k= constant, we find that

$$\begin{aligned} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\varsigma (\mathfrak{r}) =&{}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } \Biggl[ \frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \Biggr] \\ &{}+{}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } \biggl[ \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (\mathfrak{r})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}) \biggr] \\ =&{}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }\varpi (\mathfrak{r}) \\ =&\varpi (\mathfrak{r}). \end{aligned}$$

On the other hand, by taking \(\mathfrak{r}\rightarrow a\) on both sides of (3.1), then using the fact that \(\varpi (a)=0\) and \(\lim_{\mathfrak{r}\rightarrow a}I_{a^{+}}^{\vartheta }\varpi (\mathfrak{r})=I_{a^{+}}^{\vartheta }\varpi (a)=0\), we get

$$\begin{aligned} \varsigma (a) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{ \vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\varpi (a)+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \varpi (a) \\ =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) . \end{aligned}$$
(3.6)

Substitute \(\mathfrak{r}=\mathfrak{r}_{k}\) and multiply by \(c_{k}\) in (3.1). Then we derive

$$\begin{aligned} \sum_{k=1}^{m}c_{k}\varsigma ( \mathfrak{r}_{k}) =& \Biggl( \frac{\sum_{k=1}^{m}c_{k}}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k} \varpi (\mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{ \vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k} \varpi (\mathfrak{r}_{k})+\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta }c_{k} \varpi (\mathfrak{r}_{k}) \\ =& \biggl( \frac{\sum_{k=1}^{m}c_{k}}{1-\sum_{k=1}^{m}c_{k}}+1 \biggr) \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k} \varpi (\mathfrak{r}_{k})+\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) \\ =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}\varpi ( \mathfrak{r}_{k})+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \varpi (\mathfrak{r}_{k}) \Biggr) . \end{aligned}$$
(3.7)

It follows from (3.6) and (3.7) that

$$ \varsigma (a)=\sum_{k=1}^{m}c_{k} \varsigma (\mathfrak{r}_{k}). $$

 □

As a result of Theorem 3.1, we have the subsequent theorem.

Theorem 3.2

Let \(0<\vartheta \leq 1\), \(\sum_{k=1}^{m}c_{k}\neq 1\), and \(f:\mathfrak{Z}\times \mathbb{R} \times \mathbb{R} \rightarrow \mathbb{R} \) be a continuous function along with \(f(a,\varsigma (a),\varsigma (\gamma a))=0\), \(\varsigma \in \mathfrak{D}\). Then the ABC- type pantograph FDEs (1.1) are equivalent to the following FIE:

$$\begin{aligned} \varsigma (\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}f\bigl( \mathfrak{r}_{k},\varsigma (\mathfrak{r}_{k}), \varsigma \bigl(\gamma ( \mathfrak{r}_{k} ) \bigr)\bigr) \\ &{} +\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta }f \bigl(\mathfrak{r}_{k}, \varsigma (\mathfrak{r}_{k}), \varsigma \bigl(\gamma ( \mathfrak{r}_{k} ) \bigr)\bigr) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr)+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta }f \bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr),\quad \mathfrak{r}\in \mathfrak{Z}. \end{aligned}$$

4 Existence and uniqueness theorems

This section is devoted to proving the existence and uniqueness theorems for the ABC-type pantograph FDEs (1.1). Before proceeding with the main findings, we are obligated to provide the following assumption:

\((A_{1})\):

The function \(f:\mathfrak{Z} \times \mathbb{R}\times \mathbb{R} \to \mathbb{R}\) is continuous, and there exists \(L_{f}>0\) such that

$$ \bigl\vert f(\mathfrak{r},u,\overline{u})-f(\mathfrak{r},v, \overline{v}) \bigr\vert \leq L_{f} \bigl( \vert u-v \vert + \vert \overline{u}-\overline{v} \vert \bigr) , \quad \mathfrak{r}\in \mathfrak{Z},\text{ and }u,v,\overline{u},\overline{v} \in \mathbb{R} . $$

Theorem 4.1

Suppose that hypothesis \((A_{1})\) holds. Then the ABC-type pantograph FDEs (1.1) have a unique solution, provided that

$$ \begin{aligned}[b] \mathfrak{P}_{1}&:= \Biggl( \frac{(1-\vartheta )}{1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl[ \frac{1}{ ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta }+(T-a)^{ \vartheta } \Biggr] \Biggr) \frac{2L_{f}}{\mathfrak{N}(\vartheta )}\\ &< 1. \end{aligned} $$
(4.1)

Proof

Define the operator \(\mathbb{T}:\mathfrak{D}\rightarrow \mathfrak{D}\) by \(\mathbb{T}\varsigma =\varsigma \), \(\varsigma \in \mathfrak{D}\), i.e.,

$$\begin{aligned} ( \mathbb{T}\varsigma ) (\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}f \bigl(\mathfrak{r}_{k},\varsigma (\mathfrak{r}_{k}), \varsigma \bigl(\gamma ( \mathfrak{r}_{k} ) \bigr)\bigr) \\ &{} +\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta }f \bigl(\mathfrak{r}_{k}, \varsigma (\mathfrak{r}_{k}), \varsigma \bigl(\gamma ( \mathfrak{r}_{k} ) \bigr)\bigr) \Biggr) \\ &{}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr)+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta }f \bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr). \end{aligned}$$

The operator \(\mathbb{T}\) is well defined, that is, \(\mathbb{T}(\mathfrak{D})\subseteq \mathfrak{D}\). Indeed, for any \(\varsigma \in \mathfrak{D}\), \(f(\mathfrak{\cdot },\varsigma (\mathfrak{\cdot }),\varsigma (\gamma ( \mathfrak{\cdot })))\) is continuous. Besides, by Lemma 2.2, \(\mathbb{T}\varsigma \in \mathfrak{D}\). Also, by Lemma 2.1 with Remark 2.1, we end up at

$$\begin{aligned} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } ( \mathbb{T}\varsigma ) (\mathfrak{r}) =& {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } \varsigma (a)+ {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl( \mathfrak{r},\varsigma (\mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr) \\ =&f\bigl(\mathfrak{r},\varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r}) \bigr). \end{aligned}$$

Since \(f(\mathfrak{r},\cdot ,\cdot )\) is continuous on \([a,T]\), then \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } ( \mathbb{T}\varsigma ) ( \mathfrak{r})\in \mathfrak{D}\).

Now, we need to prove that \(\mathbb{T}\) is a condensing map. Let \(\varsigma ,\overline{\varsigma }\in \mathfrak{D}\) and \(\mathfrak{r}\in \mathfrak{Z}\). Then

$$\begin{aligned}& \bigl\vert (\mathbb{T}\varsigma ) (\mathfrak{r})-(\mathbb{T} \overline{\varsigma }) (\mathfrak{r}) \bigr\vert \\& \quad \leq \frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k} \bigl\vert f\bigl( \mathfrak{r}_{k}, \varsigma (\mathfrak{r}_{k}), \varsigma (\gamma \mathfrak{r}_{k})\bigr)-f\bigl(\mathfrak{r}_{k},\overline{\varsigma }(\mathfrak{r}_{k}),\overline{\varsigma }(\gamma \mathfrak{r}_{k})\bigr) \bigr\vert \\& \qquad {}+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl( \mathfrak{r}_{k},\varsigma ( \mathfrak{r}_{k}),\varsigma (\gamma \mathfrak{r}_{k}) \bigr)-f\bigl(\mathfrak{r}_{k},\overline{\varsigma }( \mathfrak{r}_{k}), \overline{\varsigma }(\gamma \mathfrak{r}_{k}) \bigr) \bigr\vert \Biggr) \\& \qquad {}+\frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\bigl\vert f\bigl( \mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr)-f( \mathfrak{r},\overline{\varsigma }(\mathfrak{r}),\overline{\varsigma }(\gamma \mathfrak{r}) \bigr\vert \\& \qquad {}+\frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl( \mathfrak{r},\varsigma (\mathfrak{r}),\varsigma ( \gamma \mathfrak{r})\bigr)-f \bigl(\mathfrak{r},\overline{\varsigma }( \mathfrak{r}),\overline{\varsigma }(\gamma \mathfrak{r})\bigr) \bigr\vert . \end{aligned}$$

By assumption (\({A}_{1}\)), we obtain

$$\begin{aligned}& I_{a^{+}}^{\vartheta } \bigl\vert f\bigl(\mathfrak{r}_{k}, \varsigma ( \mathfrak{r}_{k}),\varsigma \bigl(\gamma ( \mathfrak{r}_{k} ) \bigr)\bigr)-f\bigl( \mathfrak{r}_{k},\overline{\varsigma }(\mathfrak{r}_{k}),\overline{\varsigma }\bigl( \gamma ( \mathfrak{r}_{k} ) \bigr)\bigr) \bigr\vert \\& \quad \leq \frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}_{k}}( \mathfrak{r}_{k}- \mathfrak{s})^{\vartheta -1}\bigl\vert f\bigl(\mathfrak{s}, \varsigma (\mathfrak{s}),\varsigma \bigl(\gamma (\mathfrak{s)}\bigr)\bigr)-f\bigl( \mathfrak{s}, \overline{\varsigma }(\mathfrak{s}),\overline{\varsigma } \bigl(\gamma ( \mathfrak{s)\bigr)}\bigr)\bigr\vert \,d\mathfrak{s} \\& \quad = \frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}_{k}}( \mathfrak{r}_{k}-\mathfrak{s})^{\vartheta -1}L_{f} \bigl( \bigl\vert \varsigma ( \mathfrak{s})-\overline{\varsigma }(\mathfrak{s}) \bigr\vert + \bigl\vert \varsigma \bigl(\gamma (\mathfrak{s)}\bigr)-\overline{\varsigma }\bigl( \gamma ( \mathfrak{s)}\bigr) \bigr\vert \bigr) \,d\mathfrak{s} \\& \leq \frac{2L_{f}(\mathfrak{r}_{k}-a)^{\vartheta }}{\Gamma (\vartheta +1)} \Vert \varsigma -\overline{\varsigma } \Vert . \end{aligned}$$

Similarly,

$$\begin{aligned}& I_{a^{+}}^{\vartheta } \bigl\vert f\bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr)-f\bigl( \mathfrak{r}, \overline{\varsigma }(\mathfrak{r}),\overline{\varsigma }(\gamma \mathfrak{r})\bigr) \bigr\vert \\& \quad \leq \frac{2L_{f}(T-a)^{\vartheta }}{\Gamma (\vartheta +1)} \Vert \varsigma -\overline{\varsigma } \Vert . \end{aligned}$$

Therefore,

$$\begin{aligned}& \bigl\Vert (\mathbb{T}\varsigma )-(\mathbb{T} \overline{\varsigma }) \bigr\Vert \\& \quad = \max_{\mathfrak{r\in Z}} \bigl\vert (\mathbb{T}\varsigma ) ( \mathfrak{r})-(\mathbb{T}\overline{\varsigma }) (\mathfrak{r}) \bigr\vert \\& \quad \leq \frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{2L_{f}(1-\vartheta )}{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k}+ \frac{2L_{f}}{\mathfrak{N}(\vartheta )\Gamma (\vartheta )}\sum _{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta } \Biggr) \Vert \varsigma - \overline{\varsigma } \Vert \\& \qquad {}+ \biggl( \frac{2L_{f}(1-\vartheta )}{\mathfrak{N}(\vartheta )}+ \frac{2L_{f}(T-a)^{\vartheta }}{\mathfrak{N}(\vartheta )\Gamma (\vartheta )} \biggr) \Vert \varsigma -\overline{\varsigma } \Vert \\& \quad = \Biggl( \frac{4(1-\vartheta )\sum_{k=1}^{m}c_{k}}{1-\sum_{k=1}^{m}c_{k}}+\frac{2}{\Gamma (\vartheta )} \Biggl[ \frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta }+(T-a)^{ \vartheta } \Biggr] \Biggr)\\& \qquad {}\times \frac{L_{f}}{\mathfrak{N}(\vartheta )} \Vert \varsigma -\overline{ \varsigma } \Vert \\& \quad = \mathfrak{P}_{1} \Vert \varsigma -\overline{\varsigma } \Vert . \end{aligned}$$

Condition (4.1) shows that \(\mathbb{T}\) is a condensing operator. Hence, by Theorem 2.1, \(\mathbb{T}\) has a unique fixed point. □

Theorem 4.2

Suppose that hypothesis \((A_{1})\) holds. Then there exists at least one solution of the ABC-type pantograph FDEs (1.1), provided that condition (4.1) is satisfied.

Proof

Consider the operator \(\mathbb{T}:\mathfrak{D}\rightarrow \mathfrak{D}\) defined by

$$ (\mathbb{T}\varsigma ) (\mathfrak{r})=(\mathbb{T}_{1}\varsigma ) ( \mathfrak{r})+(\mathbb{T}_{2}\varsigma ) (\mathfrak{r}),\quad \varsigma \in \mathfrak{D}, \mathfrak{r}\in \mathfrak{Z}, $$

where

$$\begin{aligned} (\mathbb{T}_{1}\varsigma ) (\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}f \bigl(\mathfrak{r}_{k},\varsigma (\mathfrak{r}_{k}), \varsigma (\gamma \mathfrak{r}_{k})\bigr) \\ &{} +\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta }f \bigl(\mathfrak{r}_{k}, \varsigma (\mathfrak{r}_{k}), \varsigma (\gamma \mathfrak{r}_{k})\bigr) \Biggr) \end{aligned}$$
(4.2)

and

$$ (\mathbb{T}_{2}\varsigma ) (\mathfrak{r})= \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}f\bigl( \mathfrak{r},\varsigma (\mathfrak{r}),\varsigma ( \gamma \mathfrak{r})\bigr)+ \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta }f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr). $$
(4.3)

Since \(f:\mathfrak{Z}\times \mathbb{R}\times \mathbb{R}\rightarrow \mathbb{R} \) is continuous, \(\mu _{f}:=\max \{ \vert f(\mathfrak{r},0,0) \vert :\mathfrak{r}\in \mathfrak{Z} \} \) exists. Let

$$ B_{\xi }=\bigl\{ \varsigma \in \mathfrak{D}: \Vert \varsigma \Vert \leq \xi \bigr\} $$
(4.4)

with the radius

$$ \xi \geq \frac{\mathfrak{P}_{2}}{1-\mathfrak{P}_{1}}, $$
(4.5)

where

$$ \mathfrak{P}_{2}:= \Biggl( \frac{1-\vartheta }{1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl[ \frac{1}{ ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta }+(T-a)^{ \vartheta } \Biggr] \Biggr) \frac{\mu _{f}}{\mathfrak{N}(\vartheta )}. $$
(4.6)

We will complete the proof in the following several steps.

Step1: We show that \(\mathbb{T}_{1}\varsigma +\mathbb{T}_{2}\upsilon \in B_{\xi }\) for all \(\varsigma ,\upsilon \in B_{\xi }\).

By (4.2), we have

$$\begin{aligned} \bigl\vert (\mathbb{T}_{1}\varsigma ) (\mathfrak{r}) \bigr\vert \leq &\frac{1}{1-\sum_{k=1}^{m}c_{k}} \Biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )}\sum _{k=1}^{m}c_{k} \bigl\vert f\bigl( \mathfrak{r}_{k},\varsigma ( \mathfrak{r}_{k}), \varsigma (\gamma \mathfrak{r}_{k})\bigr) \bigr\vert \\ &{} +\frac{\vartheta }{\mathfrak{N}(\vartheta )}\sum_{k=1}^{m}c_{k}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl( \mathfrak{r}_{k},\varsigma ( \mathfrak{r}_{k}),\varsigma (\gamma \mathfrak{r}_{k}) \bigr) \bigr\vert \Biggr) . \end{aligned}$$
(4.7)

Using hypothesis \((A_{1})\), for \(\varsigma \in B_{\xi }\) and for any \(\mathfrak{r}\in \mathfrak{Z}\), we have

$$\begin{aligned} \bigl\vert f\bigl(\mathfrak{r},\varsigma (\mathfrak{r}),\varsigma ( \gamma \mathfrak{r})\bigr) \bigr\vert \leq & \bigl\vert f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr)-f( \mathfrak{r},0,0) \bigr\vert + \bigl\vert f(\mathfrak{r},0,0) \bigr\vert \\ \leq &L_{f} \bigl\vert \varsigma (\mathfrak{r}) \bigr\vert +L_{f} \bigl\vert \varsigma (\gamma \mathfrak{r}) \bigr\vert +\mu _{f} \\ \leq &2L_{f}\xi +\mu _{f}. \end{aligned}$$
(4.8)

Further, by using (4.8), for any \(\mathfrak{r}\in \mathfrak{Z}\), we have

$$\begin{aligned} I_{a^{+}}^{\vartheta } \bigl\vert f\bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert \leq &I_{a^{+}}^{ \vartheta } ( 2L_{f}\xi +\mu _{f} ) \\ =& ( 2L_{f}\xi +\mu _{f} ) \frac{(\mathfrak{r}-a)^{\vartheta }}{\Gamma (\vartheta +1)}. \end{aligned}$$
(4.9)

Using inequalities (4.8) and (4.9) into inequality (4.7), we obtain

$$\begin{aligned} \bigl\Vert (\mathbb{T}_{1}\varsigma ) \bigr\Vert =&\max _{ \mathfrak{r}\in \mathfrak{Z}} \bigl\vert (\mathbb{T}_{1}\varsigma ) ( \mathfrak{r}) \bigr\vert \\ \leq &\frac{ ( 2L_{f}\xi +\mu _{f} ) }{\mathfrak{N}(\vartheta ) ( 1-\sum_{k=1}^{m}c_{k} ) } \Biggl( (1-\vartheta )\sum _{k=1}^{m}c_{k}+ \frac{1}{\Gamma (\vartheta )}\sum _{k=1}^{m}c_{k}(\mathfrak{r}_{k}-a)^{\vartheta } \Biggr) \\ =& \Biggl( (1-\vartheta ) \frac{\sum_{k=1}^{m}c_{k}}{1-\sum_{k=1}^{m}c_{k}}+\frac{1}{\Gamma (\vartheta ) ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta } \Biggr) \frac{2L_{f}}{\mathfrak{N}(\vartheta )}\xi \\ &{}+ \Biggl( (1-\vartheta ) \frac{\sum_{k=1}^{m}c_{k}}{1-\sum_{k=1}^{m}c_{k}}+\frac{1}{\Gamma (\vartheta ) ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta } \Biggr) \frac{\mu _{f}}{\mathfrak{N}(\vartheta )}. \end{aligned}$$
(4.10)

Also, for \(\upsilon \in B_{\xi }\),

$$\begin{aligned} \bigl\Vert (\mathbb{T}_{2}\upsilon ) \bigr\Vert =&\max _{ \mathfrak{r}\in \mathfrak{Z}} \bigl\vert (\mathbb{T}_{2}\upsilon ) ( \mathfrak{r}) \bigr\vert \\ \leq &\max_{\mathfrak{r}\in \mathfrak{Z}} \biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f \bigl(\mathfrak{r},\upsilon ( \mathfrak{r}),\upsilon (\gamma \mathfrak{r})\bigr) \bigr\vert + \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl(\mathfrak{r}, \upsilon (\mathfrak{r}),\upsilon ( \gamma \mathfrak{r})\bigr) \bigr\vert \biggr) \\ \leq &\frac{ ( 1-\vartheta ) }{\mathfrak{N}(\vartheta )} ( 2L_{f}\xi +\mu _{f} ) + \frac{ ( 2L_{f}\xi +\mu _{f} ) }{\mathfrak{N}(\vartheta )} \frac{(T-a)^{\vartheta }}{\Gamma (\vartheta )} \\ =& \biggl( ( 1-\vartheta ) + \frac{(T-a)^{\vartheta }}{\Gamma (\vartheta )} \biggr) \frac{2L_{f}}{\mathfrak{N}(\vartheta )}\xi \\ &{}+ \biggl( ( 1-\vartheta ) + \frac{(T-a)^{\vartheta }}{\Gamma (\vartheta )} \biggr) \frac{\mu _{f}}{\mathfrak{N}(\vartheta )}. \end{aligned}$$
(4.11)

Inequalities (4.10) and (4.11) give

$$\begin{aligned} & \bigl\Vert (\mathbb{T}_{1}\varsigma )+(\mathbb{T}_{2} \upsilon ) \bigr\Vert \\ &\quad \leq \bigl\Vert (\mathbb{T}_{1}\varsigma ) \bigr\Vert + \bigl\Vert (\mathbb{T}_{2}\upsilon ) \bigr\Vert \\ &\quad \leq \Biggl( \frac{1-\vartheta }{1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl[ \frac{1}{ ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta }+(T-a)^{ \vartheta } \Biggr] \Biggr) \frac{2L_{f}}{\mathfrak{N}(\vartheta )}\xi \\ &\qquad {} + \Biggl( \frac{1-\vartheta }{1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl[ \frac{1}{ ( 1-\sum_{k=1}^{m}c_{k} ) }\sum_{k=1}^{m}c_{k}( \mathfrak{r}_{k}-a)^{\vartheta }+(T-a)^{ \vartheta } \Biggr] \Biggr) \frac{\mu _{f}}{\mathfrak{N}(\vartheta )} \\ &\quad =\mathfrak{P}_{1}\xi +\mathfrak{P}_{2}. \end{aligned}$$

Using (4.1) and (4.5), we get

$$ \Vert \mathbb{T}_{1}\varsigma +\mathbb{T}_{2}\upsilon \Vert \leq \xi . $$

Thus, \(\mathbb{T}_{1}\varsigma +\mathbb{T}_{2}\upsilon \in B_{\xi }\) for all \(\varsigma ,\upsilon \in B_{\xi }\).

Step 2. \(\mathbb{T}_{1}\) is a condensing map. This is evident due to \(\mathbb{T}\) is a contraction map.

Step 3: \(\mathbb{T}_{2}\) is continuous and compact.

\(\mathbb{T}_{2}:B_{\xi }\rightarrow B_{\xi }\) is continuous due to f is continuous. Indeed, let \(\varsigma _{n}\) be a sequence such that \(\varsigma _{n}\rightarrow \varsigma \) in \(\mathfrak{D}\). Then, for all \(\mathfrak{r}\in \mathfrak{Z}\), one has

$$\begin{aligned} \bigl\vert \bigl(\mathbb{T}_{2}\varsigma _{n}( \mathfrak{r})\bigr)-\bigl( \mathbb{T}_{2}\varsigma ( \mathfrak{r})\bigr) \bigr\vert \leq &\frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f\bigl( \mathfrak{r},\varsigma _{n}( \mathfrak{r}),\varsigma _{n}(\gamma \mathfrak{r})\bigr)-f\bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert \\ &{}+\frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl( \mathfrak{r},\varsigma _{n}(\mathfrak{r}),\varsigma _{n}( \gamma \mathfrak{r})\bigr)-f\bigl(\mathfrak{r},\varsigma (\mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert \\ \leq &\frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\Vert f\bigl( \mathfrak{\cdot }, \varsigma _{n}(\mathfrak{\cdot }),\varsigma _{n}\bigl(\gamma ( \mathfrak{\cdot })\bigr)\bigr)-f\bigl(\mathfrak{\cdot },\varsigma ( \mathfrak{\cdot }), \varsigma \bigl(\gamma (\mathfrak{\cdot })\bigr)\bigr) \bigr\Vert \\ &{}+ \frac{(T-a)^{\vartheta }}{\mathfrak{N}(\vartheta )\Gamma (\vartheta )} \bigl\Vert f\bigl(\mathfrak{\cdot }, \varsigma _{n}(\mathfrak{\cdot }), \varsigma _{n}\bigl( \gamma (\mathfrak{\cdot })\bigr)\bigr)-f\bigl(\mathfrak{\cdot }, \varsigma ( \mathfrak{\cdot }),\varsigma \bigl(\gamma (\mathfrak{\cdot })\bigr) \bigr) \bigr\Vert . \end{aligned}$$

Since f is continuous, the operator \(\mathbb{T}_{2}\) is also continuous. Thus, we have

$$ \bigl\Vert (\mathbb{T}_{2}\varsigma _{n})-( \mathbb{T}_{2}\varsigma ) \bigr\Vert \rightarrow 0\quad \text{as }n \rightarrow \infty . $$

Next, \(\mathbb{T}_{2}\) is uniformly bounded on \(B_{\xi }\). For any \(\varsigma \in B_{\xi }\) and \(\mathfrak{r}\in \mathfrak{Z}\), we have

$$\begin{aligned} \bigl\vert (\mathbb{T}_{2}\varsigma ) \bigr\vert =&\max _{ \mathfrak{r}\in \mathfrak{Z}} \bigl\vert (\mathbb{T}_{2}\varsigma ) ( \mathfrak{r}) \bigr\vert \\ \leq &\max_{\mathfrak{r}\in \mathfrak{Z}} \biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f \bigl(\mathfrak{r},\varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert + \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}),\varsigma ( \gamma \mathfrak{r})\bigr) \bigr\vert \biggr) \\ \leq &\max_{\mathfrak{r}\in \mathfrak{Z}} \biggl( \frac{ ( 1-\vartheta ) }{\mathfrak{N}(\vartheta )} ( 2L_{f} \xi +\mu _{f} ) + \frac{ ( 2L_{f}\xi +\mu _{f} ) }{\mathfrak{N}(\vartheta )} \frac{(\mathfrak{r}-a)^{\vartheta }}{\Gamma (\vartheta )} \biggr) \\ \leq & \frac{2L_{f}\xi +\mu _{f} }{\mathfrak{N}(\vartheta )} \biggl( (1- \vartheta )+\frac{(T-a)^{\vartheta }}{\Gamma (\vartheta )} \biggr). \end{aligned}$$

This leads to \(\mathbb{T}_{2}\) is uniformly bounded on \(B_{\xi }\).

Now, we show that \(\mathbb{T}_{2}(B_{\xi })\) is equicontinuous. For that, let \(\varsigma \in B_{\xi }\) and \(a\leq \mathfrak{r}_{1}<\mathfrak{r}_{2}\leq T\). Then, by using (4.8), we have

$$\begin{aligned}& \bigl\vert (\mathbb{T}_{2}\varsigma ) (\mathfrak{r}_{2})-( \mathbb{T}_{2}\varsigma ) (\mathfrak{r}_{1}) \bigr\vert \\& \quad = \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f\bigl( \mathfrak{r}_{2}, \varsigma (\mathfrak{r}_{2}),\varsigma \bigl(\gamma ( \mathfrak{r}_{2})\bigr)\bigr)-f\bigl(\mathfrak{r}_{1}, \varsigma (\mathfrak{r}_{1}),\varsigma \bigl(\gamma ( \mathfrak{r}_{1})\bigr)\bigr) \bigr\vert \\& \qquad {}+ \biggl\vert \frac{\vartheta }{\mathfrak{N}(\vartheta )} \frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}_{2}}( \mathfrak{r}_{2}- \mathfrak{s})^{\vartheta -1}f \bigl( \mathfrak{s},\varsigma ( \mathfrak{s}), \varsigma \bigl(\gamma (\mathfrak{s})\bigr) \bigr) \,d\mathfrak{s} \\& \qquad {} -\frac{\vartheta }{\mathfrak{N}(\vartheta )} \frac{1}{\Gamma (\vartheta )} \int _{a}^{\mathfrak{r}_{1}}( \mathfrak{r}_{1}- \mathfrak{s})^{\vartheta -1}f \bigl( \mathfrak{s},\varsigma ( \mathfrak{s}), \varsigma \bigl(\gamma (\mathfrak{s})\bigr) \bigr) \,d\mathfrak{s} \biggr\vert \\& \quad \leq \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f\bigl( \mathfrak{r}_{2}, \varsigma (\mathfrak{r}_{2}),\varsigma \bigl(\gamma ( \mathfrak{r}_{2})\bigr)\bigr)-f\bigl(\mathfrak{r}_{1}, \varsigma (\mathfrak{r}_{1}),\varsigma \bigl(\gamma ( \mathfrak{r}_{1})\bigr)\bigr) \bigr\vert \\& \qquad {}+\frac{\vartheta }{\mathfrak{N}(\vartheta )} \frac{1}{\Gamma (\vartheta )}\int _{a}^{\mathfrak{r}_{1}} \bigl\vert (\mathfrak{r}_{1}- \mathfrak{s})^{\vartheta -1}-(\mathfrak{r}_{2}- \mathfrak{s})^{\vartheta -1} \bigr\vert \bigl\vert f \bigl( \mathfrak{s}, \varsigma ( \mathfrak{s}),\varsigma \bigl(\gamma (\mathfrak{s})\bigr) \bigr) \bigr\vert \,d\mathfrak{s} \\& \qquad {}+\frac{\vartheta }{\mathfrak{N}(\vartheta )} \frac{1}{\Gamma (\vartheta )}\int _{\mathfrak{r}_{1}}^{\mathfrak{r}_{2}}(\mathfrak{r}_{2}- \mathfrak{s})^{\vartheta -1} \bigl\vert f \bigl( \mathfrak{s}, \varsigma ( \mathfrak{s}),\varsigma \bigl(\gamma (\mathfrak{s})\bigr) \bigr) \bigr\vert \,d \mathfrak{s} \\& \quad \leq \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert f\bigl( \mathfrak{r}_{2}, \varsigma (\mathfrak{r}_{2}),\varsigma \bigl(\gamma ( \mathfrak{r}_{2})\bigr)\bigr)-f\bigl(\mathfrak{r}_{1}, \varsigma (\mathfrak{r}_{1}),\varsigma \bigl(\gamma ( \mathfrak{r}_{1})\bigr)\bigr) \bigr\vert \\& \qquad {}+ \frac{2 (2L_{f}\xi +\mu _{f}) }{\mathfrak{N}(\vartheta )\Gamma (\vartheta )}( \mathfrak{r}_{2}-\mathfrak{r}_{1})^{\vartheta }. \end{aligned}$$

Since \(f(\mathfrak{\cdot },\varsigma (\mathfrak{\cdot }),\varsigma (\gamma (\mathfrak{\cdot })))\) is continuous, \(\vert (\mathbb{T}_{2}\varsigma )(\mathfrak{r}_{2})-(\mathbb{T}_{2}\varsigma )(\mathfrak{r}_{1}) \vert \rightarrow 0\) as \(\mathfrak{r}_{2}\rightarrow \mathfrak{r}_{1}\). In light of the former steps with Arzela–Ascoli theorem, we derive that \((\mathbb{T}_{2}B_{\xi })\) is relatively compact, and hence \(\mathbb{T}_{2}\) is completely continuous. So, Theorem 2.2 shows that (1.1) has at least one solution. □

5 Ulam–Hyers stability

In this section, we discuss two types of stability for (1.1), namely Ulam–Hyers and generalized Ulam–Hyers stabilities. For \(\varepsilon >0\), we consider the following inequations:

$$ \bigl\vert {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta } \widetilde{ \varsigma }(\mathfrak{r})-f\bigl(\mathfrak{r},\widetilde{\varsigma }( \mathfrak{r}), \widetilde{\varsigma }(\gamma \mathfrak{r})\bigr) \bigr\vert \leq \varepsilon , \quad \mathfrak{r}\in \mathfrak{Z}. $$
(5.1)

Definition 5.1

The pantograph ABC-FDEs (1.1) are Ulam–Hyers stable if there exists a real number \(C_{f}>0\) such that, for each \(\varepsilon >0\) and for each solution \(\widetilde{\varsigma }\in \mathfrak{D}\) of inequality (5.1), there exists a unique solution \(\varsigma \in \mathfrak{D}\) of (1.1) with

$$ \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\varsigma ( \mathfrak{r}) \bigr\vert \leq C_{f}\varepsilon , \quad \mathfrak{r}\in \mathfrak{Z}. $$

And the pantograph ABC-FDEs (1.1) are generalized Ulam–Hyers stable if we can find \(\phi _{f}:\mathbb{R} ^{+}\rightarrow \mathbb{R} ^{+}\) with \(\phi _{f}(0)=0\) such that

$$ \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\varsigma ( \mathfrak{r}) \bigr\vert \leq \phi _{f}(\varepsilon ), \quad \mathfrak{r}\in \mathfrak{Z}. $$

Remark 5.1

Let \(\widetilde{\varsigma }\in \mathfrak{D}\) be the solution of inequality (5.1) if and only if we have a function \(h\in \mathfrak{D}\) which depends on ς̃ such that

  1. i)

    \(\vert h(\mathfrak{r}) \vert \leq \varepsilon \) for all \(\mathfrak{r}\in \mathfrak{Z}\),

  2. ii)

    \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\widetilde{\varsigma }(\mathfrak{r})=f(\mathfrak{r},\widetilde{\varsigma }(\mathfrak{r}), \widetilde{\varsigma }(\gamma \mathfrak{r}))+h(\mathfrak{r})\), \(\mathfrak{r}\in \mathfrak{Z}\).

Lemma 5.1

If \(\widetilde{\varsigma }\in \mathfrak{D}\) is a solution of inequality (5.1), then ς̃ is a solution of the following inequality:

$$ \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\mathcal{R}_{ \widetilde{\varsigma }}-{}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl( \mathfrak{r},\widetilde{ \varsigma }(\mathfrak{r}),\widetilde{\varsigma }(\gamma \mathfrak{r}) \bigr) \bigr\vert \leq \Lambda \varepsilon , $$
(5.2)

where

$$ \mathcal{R}_{\widetilde{\varsigma }}:= \frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum _{k=1}^{m}c_{k} {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl( \mathfrak{r}_{k}, \widetilde{\varsigma }(\mathfrak{r}_{k}),\widetilde{\varsigma }( \gamma \mathfrak{r}_{k})\bigr) $$

and

$$ \Lambda := \Biggl( \frac{ 1-\vartheta }{1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl( \frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum_{k=1}^{m}c_{k} (\mathfrak{r}_{k}-a)\Biggr)^{\vartheta }+(T-a)\Biggr)^{ \vartheta } \frac{1}{\mathfrak{N}(\vartheta )}. $$

Proof

In view of Remark 5.1, we have

$$\begin{aligned}& {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\widetilde{\varsigma }( \mathfrak{r}) = f\bigl(\mathfrak{r},\widetilde{\varsigma }(\mathfrak{r}), \widetilde{\varsigma }(\gamma \mathfrak{r})\bigr)+h(\mathfrak{r}), \quad \mathfrak{r}\in \mathfrak{Z}, \\& \widetilde{\varsigma }(a) = \sum_{k=1}^{m}c_{k} \widetilde{\varsigma }(\mathfrak{r}_{k}), \quad \mathfrak{r}_{k}\in \mathfrak{Z}^{ \prime }. \end{aligned}$$

Then, by Theorem 3.1, we get

$$\begin{aligned} \widetilde{\varsigma }(\mathfrak{r}) =&\frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum _{k=1}^{m}c_{k} {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta } \bigl( f\bigl(\mathfrak{r}_{k},\widetilde{\varsigma }(\mathfrak{r}_{k}), \widetilde{\varsigma }(\gamma \mathfrak{r}_{k})\bigr)+h( \mathfrak{r}_{k}) \bigr) \\ &{}+{}^{AB}\mathbb{I}_{a^{+}}^{\vartheta } \bigl( f\bigl( \mathfrak{r}, \widetilde{\varsigma }(\mathfrak{r}),\widetilde{\varsigma }(\gamma \mathfrak{r})\bigr)+h(\mathfrak{r}) \bigr) \\ =&\mathcal{R}_{\widetilde{\varsigma }}+ \frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum _{k=1}^{m}c_{k} {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }h( \mathfrak{r}_{k})+{}^{AB} \mathbb{I}_{a^{+}}^{\vartheta } \bigl( f\bigl(\mathfrak{r}, \widetilde{\varsigma }(\mathfrak{r}),\widetilde{\varsigma }(\gamma \mathfrak{r})\bigr)+h(\mathfrak{r}) \bigr) . \end{aligned}$$

From this it follows that

$$\begin{aligned} & \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\mathcal{R}_{ \widetilde{\varsigma }} - {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }f\bigl( \mathfrak{r},\widetilde{\varsigma }(\mathfrak{r}),\widetilde{\varsigma }(\gamma \mathfrak{r})\bigr) \bigr\vert \\ &\quad \leq \frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum_{k=1}^{m}c_{k} {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta } \bigl\vert h(\mathfrak{r}_{k}) \bigr\vert + {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta } \bigl\vert h( \mathfrak{r}) \bigr\vert \\ &\quad =\frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum_{k=1}^{m}c_{k} \biggl( \frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert h(\mathfrak{r}_{k}) \bigr\vert +\frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert h(\mathfrak{r}_{k}) \bigr\vert \biggr) \\ &\qquad {} +\frac{1-\vartheta }{\mathfrak{N}(\vartheta )} \bigl\vert h( \mathfrak{r}) \bigr\vert + \frac{\vartheta }{\mathfrak{N}(\vartheta )}I_{a^{+}}^{\vartheta } \bigl\vert h(\mathfrak{r}) \bigr\vert \\ &\quad \leq \Biggl( \frac{ 1-\vartheta }{ 1-\sum_{k=1}^{m}c_{k}}+ \frac{1}{\Gamma (\vartheta )} \Biggl( \frac{1}{1-\sum_{k=1}^{m}c_{k}} \sum_{k=1}^{m}c_{k} (\mathfrak{r}_{k}-a)\Biggr)^{\vartheta }+(T-a)\Biggr)^{\vartheta } \frac{\varepsilon }{\mathfrak{N}(\vartheta )} \\ &\quad =\Lambda \varepsilon . \end{aligned}$$

 □

Theorem 5.1

Suppose that hypothesis \((A_{1})\) holds. If \(\mathfrak{N}(\vartheta )- ( 1-\vartheta ) 2L_{f}<1\), then the pantograph ABC-FDEs (1.1) are Ulam–Hyers stable.

Proof

Let \(\varepsilon >0\) and \(\widetilde{\varsigma }\in \mathfrak{D}\) be a function which satisfies inequality (5.1), and let \(\varsigma \in \mathfrak{D}\) be the unique solution of the following problem:

$$ \textstyle\begin{cases} {}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\varsigma (\mathfrak{r})=f( \mathfrak{r},\varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r})),& \mathfrak{r}\in \mathfrak{Z}, \\ \varsigma (a)=\sum_{k=1}^{m}c_{k}\varsigma (\mathfrak{r}_{k})=\sum_{k=1}^{m}c_{k}\widetilde{\varsigma }(\mathfrak{r}_{k})= \widetilde{\varsigma }(a),&\mathfrak{r}_{k}\in \mathfrak{Z}^{\prime }.\end{cases} $$
(5.3)

Using Theorem 3.1, we obtain

$$ \varsigma (\mathfrak{r})=\mathcal{R}_{\varsigma }+{}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl(\mathfrak{r},\varsigma ( \mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr), $$

where

$$ \mathcal{R}_{\varsigma }:=\frac{1}{1-\sum_{k=1}^{m}c_{k}}\sum_{k=1}^{m}c_{k} {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }f\bigl( \mathfrak{r}_{k}, \varsigma (\mathfrak{r}_{k}), \varsigma (\gamma \mathfrak{r}_{k})\bigr). $$

Since \(\varsigma (a)=\widetilde{\varsigma }(a)\) and \(\sum_{k=1}^{m}c_{k}\varsigma (\mathfrak{r}_{k})=\sum_{k=1}^{m}c_{k} \widetilde{\varsigma }(\mathfrak{r}_{k})\), then \(\mathcal{R}_{\varsigma }=\mathcal{R}_{\widetilde{\varsigma }}\). Hence

$$ \varsigma (\mathfrak{r})=\mathcal{R}_{\widetilde{\varsigma }}+{}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl(\mathfrak{r},\varsigma ( \mathfrak{r}), \varsigma (\gamma \mathfrak{r})\bigr). $$

It follows from Lemma 5.1 and (\({A}_{1}\)) that

$$\begin{aligned}& \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\varsigma ( \mathfrak{r}) \bigr\vert \\& \quad \leq \bigl\vert \widetilde{\varsigma }(\mathfrak{r})- \mathcal{R}_{\widetilde{\varsigma }}-{}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }f \bigl( \mathfrak{r},\widetilde{\varsigma }(\mathfrak{r}),\widetilde{ \varsigma }(\gamma \mathfrak{r})\bigr) \bigr\vert \\& \qquad {}+ \bigl\vert {}^{AB}\mathbb{I}_{a^{+}}^{\vartheta }f \bigl( \mathfrak{r},\widetilde{\varsigma }(\mathfrak{r}),\widetilde{ \varsigma }(\gamma \mathfrak{r})\bigr)- {}^{AB} \mathbb{I}_{a^{+}}^{\vartheta }f\bigl(\mathfrak{r}, \varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert \\& \quad \leq \Lambda \varepsilon + {}^{AB}\mathbb{I}_{a^{+}}^{ \vartheta } \bigl\vert f\bigl(\mathfrak{r},\widetilde{\varsigma }( \mathfrak{r}), \widetilde{\varsigma }(\gamma \mathfrak{r})\bigr)-f\bigl(\mathfrak{r}, \varsigma ( \mathfrak{r}),\varsigma (\gamma \mathfrak{r})\bigr) \bigr\vert \\& \quad \leq \Lambda \varepsilon +2L_{f} {}^{AB} \mathbb{I}_{a^{+}}^{ \vartheta } \bigl\vert \widetilde{\varsigma }( \mathfrak{r})- \varsigma (\mathfrak{r}) \bigr\vert . \end{aligned}$$

Using Theorem 2.3 with \(\sigma (\mathfrak{r})= \vert \widetilde{\varsigma }(t)-\varsigma (t) \vert \), \(a(\mathfrak{r})=\Lambda \varepsilon \), and \(b(\mathfrak{r})=2L_{f}\), we get

$$\begin{aligned} \bigl\vert \widetilde{\varsigma }(\mathfrak{r})-\varsigma ( \mathfrak{r}) \bigr\vert \leq &\frac{\Lambda \varepsilon \mathfrak{N}(\vartheta )}{\mathfrak{N}(\vartheta )- ( 1-\vartheta ) 2L_{f}} \mathbb{E}_{\vartheta } \biggl( \frac{2\vartheta L_{f}(\mathfrak{r}-a)^{\vartheta }}{\mathfrak{N}(\vartheta )-2 ( 1-\vartheta ) L_{f}} \biggr) \\ \leq &C_{f}\varepsilon , \end{aligned}$$

where

$$ C_{f}:= \frac{\Lambda \mathfrak{N}(\vartheta )}{\mathfrak{N}(\vartheta )- ( 1-\vartheta ) 2L_{f}} \mathbb{E}_{\vartheta } \biggl( \frac{2\vartheta L_{f}(T-a)^{\vartheta }}{\mathfrak{N}(\vartheta )-2 ( 1-\vartheta ) L_{f}} \biggr) . $$

 □

Corollary 5.1

Under the hypotheses of Theorem 5.1, if there exists \(\phi _{f}:\mathbb{R} ^{+}\rightarrow \mathbb{R} ^{+}\) with \(\phi _{f}(0)=0\), then the pantograph ABC problem (1.1) has generalized Ulam–Hyers stability.

Proof

Choosing \(\phi _{f}(\varepsilon )=C_{f}\varepsilon \) and \(\phi _{f}(0)=0\), from Theorem 5.1 we obtain \(\vert \widetilde{\varsigma }(\mathfrak{r})-\varsigma (\mathfrak{r}) \vert \leq \phi _{f}( \varepsilon )\). □

6 Examples

In this section, we justify the validity of Theorems 4.1, 4.2, and 5.1 through an example.

Example 6.1

For \(\vartheta \in (0,1]\), we consider the following ABC fractional problem:

$$ \textstyle\begin{cases} {}^{ABC}\mathbb{D}_{0^{+}}^{\frac{1}{2}}\varsigma (\mathfrak{r})= \frac{\mathfrak{r}^{2}}{10} ( e^{-\mathfrak{r}}+ \frac{ \vert \varsigma (\mathfrak{r}) \vert }{10 ( 1+ \vert \varsigma (\mathfrak{r}) \vert ) }+ \frac{ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert }{10 ( 1+ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert ) } ) ,& \mathfrak{r}\in [ 0,1], \\ \varsigma (0)=\frac{1}{4}\varsigma (\frac{1}{3}), \end{cases} $$
(6.1)

where \(\vartheta =\frac{1}{2}\), \(c_{1}=\frac{1}{4}\), \(\mathfrak{r}_{1}=\frac{1}{3}\) \((m=1)\), \(f(\mathfrak{r},\varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r}))=\frac{\mathfrak{r}^{2}}{10} ( e^{-\mathfrak{r}}+ \frac{ \vert \varsigma (\mathfrak{r}) \vert }{10 ( 1+ \vert \varsigma (\mathfrak{r}) \vert ) }+ \frac{ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert }{10 ( 1+ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert ) } ) \), and \(\gamma =\frac{1}{2}\).

Clearly, \(f(0,\varsigma (0),\overline{\varsigma }(0))=0\). Moreover,

$$ \mu _{f}=\max_{\mathfrak{r}\in [ 0,1]} \bigl\vert f( \mathfrak{r} ,0,0) \bigr\vert =\max_{\mathfrak{r}\in [ 0,1]} \frac{\mathfrak{r}^{2}}{ 10}e^{-\mathfrak{r}}= \frac{1}{10e}. $$

Let \(\mathfrak{r}\in [ 0,1]\) and \(\varsigma ,\upsilon \in \mathbb{R} \). Then

$$\begin{aligned}& \biggl\vert f(\mathfrak{r},\varsigma (\mathfrak{r}),\varsigma \biggl( \frac{1}{2}\mathfrak{r}\biggr)-f(\mathfrak{r},\upsilon ( \mathfrak{r}),\upsilon \biggl( \frac{1}{2}\mathfrak{r}\biggr) \biggr\vert \\& \quad \leq \biggl\vert \frac{\mathfrak{r}^{2}}{10} \biggl( e^{- \mathfrak{r}}+ \frac{ \vert \varsigma (\mathfrak{r}) \vert }{10 ( 1+ \vert \varsigma (\mathfrak{r}) \vert ) }+ \frac{ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert }{10 ( 1+ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert ) } \biggr) \\& \qquad {}- \frac{\mathfrak{r}^{2}}{10} \biggl( e^{-\mathfrak{r}}+ \frac{ \vert \upsilon (\mathfrak{r}) \vert }{10 ( 1+ \vert \upsilon (\mathfrak{r}) \vert ) }+ \frac{ \vert \upsilon (\frac{\mathfrak{r}}{2}) \vert }{10 ( 1+ \vert \upsilon (\frac{\mathfrak{r}}{2}) \vert ) } \biggr) \biggr\vert \\& \quad \leq \frac{1}{10} \biggl( \frac{10 \vert \varsigma (\mathfrak{r})-\upsilon (\mathfrak{r}) \vert }{100 ( 1+ \vert \varsigma (\mathfrak{r}) \vert ) ( 1+ \vert \upsilon (\mathfrak{r}) \vert ) }+ \frac{10 \vert \varsigma (\frac{\mathfrak{r}}{2})-\upsilon (\frac{\mathfrak{r}}{2}) \vert }{100 ( 1+ \vert \varsigma (\frac{\mathfrak{r}}{2}) \vert ) ( 1+ \vert \upsilon (\frac{\mathfrak{r}}{2}) \vert ) } \biggr) \\& \quad \leq \frac{1}{10} \biggl( \bigl\vert \varsigma (\mathfrak{r})- \upsilon (\mathfrak{r}) \bigr\vert + \biggl\vert \varsigma \biggl( \frac{\mathfrak{r}}{2}\biggr)-\upsilon \biggl(\frac{\mathfrak{r}}{2}\biggr) \biggr\vert \biggr) . \end{aligned}$$

Therefore, hypothesis \((A_{1})\) holds with \(L_{f}=\frac{1}{10}\). We shall examine that condition (4.1) is satisfied with \(\mathfrak{N}(\vartheta )=1\) and \(m=1\). Hence, by some simple calculations, we find that

$$ \mathfrak{P}_{1}= \biggl( \frac{1}{15}+ \frac{9+\sqrt{3}}{45\sqrt{\pi }} \biggr) \approx 0.2< 1. $$

Thus Theorems 4.1 and 4.2 guarantee the existence and uniqueness of solution on \([0,1]\) problem (6.1).

Finally, since \(\mathfrak{N}(\vartheta )- ( 1-\vartheta ) 2L_{f}=\frac{9}{10}<1\), problem (6.1) is Ulam–Hyers and generalized Ulam–Hyers stable with

$$ C_{f}=\frac{10}{9}\Lambda \mathbb{E}_{\frac{1}{2}} \biggl( \frac{1}{\mathfrak{9}} \biggr) \quad \text{and}\quad \Lambda = \biggl( \frac{1}{3}+ \frac{9+\sqrt{3}}{9\sqrt{\pi }} \biggr) . $$

Remark 6.1

  1. (1)

    If \(\gamma =1\), then our problem (1.1) reduces to problem (1.2). Therefore, all the results mentioned in this work are also valid for problem (1.2).

  2. (2)

    If we replace the nonlocal condition \(\varsigma (a)=\sum_{k=1}^{m}c_{k}\varsigma (\mathfrak{r}_{k})\) with the initial condition \(\varsigma (0)=\varsigma _{0}\) \((a=0)\) and use \({}^{CF}\mathbb{D}_{a^{+}}^{\vartheta }\) (Caputo–Fabrizio FD) instead of \({}^{ABC}\mathbb{D}_{a^{+}}^{\vartheta }\) (ABC derivative), then our problem (1.1) reduces to the following problem:

    $$ \textstyle\begin{cases} {}^{CF}\mathbb{D}_{0^{+}}^{\vartheta }\varsigma (\mathfrak{r})=f( \mathfrak{r},\varsigma (\mathfrak{r}),\varsigma (\gamma \mathfrak{r})),& \mathfrak{r}\in [ a,T], 0< \vartheta \leq 1, \\ \varsigma (0)=\varsigma _{0}. \end{cases} $$

7 Conclusions

The theory of fractional operators including nonsingular kernels is novel and of significant recent interest, thus there is a need to study the qualitative properties of FDEs involving such operators. In this paper, we have obtained the existence and uniqueness of solutions for the pantograph FDEs with nonlocal conditions involving ABC fractional derivative. Our approach is based on the reduction of ABC-type pantograph FDEs into FIE and some fixed point theorems of Banach and Krasnoselskii. Further, we have applied Gronwall’s inequality in the frame of AB fractional integral operator to develop adequate results for various types of Ulam–Hyers stability. Pertinent examples are provided to justify the results. The problems scrutinized are also valid for some special cases, in other words, they are reduced to the corresponding problems that contain Caputo–Fabrizio fractional derivative operator. Besides, the analysis of the obtained results was restricted to a minimum of assumptions.

Availability of data and materials

Not applicable.

References

  1. Podlubny, I.: Fractional Differential Equations. Academic Press, San Diego (1999)

    MATH  Google Scholar 

  2. Samko, S.G., Kilbas, A.A., Marichev, O.I., et al.: Fractional Integrals and Derivatives. Gordon & Breach, Yverdon (1993)

    MATH  Google Scholar 

  3. Diethelm, K.: The analysis of fractional differential equations. 2004 of lecture notes in mathematics, (2010)

  4. Magin, R.L.: Fractional Calculus in Bioengineering. 2. Begell House, Redding (2006)

    Google Scholar 

  5. Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies. Elsevier, Amsterdam (2006)

    MATH  Google Scholar 

  6. Machado, J.T., Kiryakova, V., Mainardi, F.: Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 16(3), 1140–1153 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  7. Hilfer, R., et al.: Applications of Fractional Calculus in Physics, vol. 35. World Scientific, Singapore (2000)

    Book  MATH  Google Scholar 

  8. Lorenzo, C.F., Hartley, T.T.: Variable order and distributed order fractional operators. Nonlinear Dyn. 29(1–4), 57–98 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  9. Agrawal, O.P.: Formulation of Euler-Lagrange equations for fractional variational problems. J. Math. Anal. Appl. 272(1), 368–379 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  10. Losada, J., Nieto, J.J.: Properties of a new fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 1(2), 87–92 (2015)

    Google Scholar 

  11. Caputo, M., Fabrizio, M.: A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 1(2), 1–13 (2015)

    Google Scholar 

  12. Khalil, R., Al Horani, M., Yousef, A., Sababheh, M.: A new definition of fractional derivative. J. Comput. Appl. Math. 264, 65–70 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  13. Jarad, F., Abdeljawad, T., Hammouch, Z.: On a class of ordinary differential equations in the frame of Atangana-Baleanu fractional derivative. Chaos Solitons Fractals 117, 16–20 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  14. Atangana, A., Baleanu, D.: New fractional derivatives with nonlocal and nonsingular kernel: theory and application to heat transfer model (2016). arXiv:1602.03408. arXiv preprint

  15. Atangana, A., Koca, I.: Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order. Chaos Solitons Fractals 89, 447–454 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  16. Alkahtani, B.S.T.: Chua’s circuit model with Atangana-Baleanu derivative with fractional order. Chaos Solitons Fractals 89, 547–551 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  17. Abdeljawad, T., Baleanu, D.: Discrete fractional differences with nonsingular discrete Mittag-Leffler kernels. Adv. Differ. Equ. 2016(1), 232 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  18. Abdeljawad, T., Baleanu, D.: On fractional derivatives with exponential kernel and their discrete versions. Rep. Math. Phys. 80(1), 11–27 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  19. Singh, J., Kumar, D., Hammouch, Z., Atangana, A.: A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl. Math. Comput. 316, 504–515 (2018)

    MathSciNet  MATH  Google Scholar 

  20. Uïcar, S., Uïcar, E., Ozdemir, N., Hammouch, Z.: Mathematical analysis and numerical simulation for a smoking model with Atangana-Baleanu derivative. Chaos Solitons Fractals 118, 300–306 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  21. Atangana, A.: Modelling the spread of Covid-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination. Chaos Solitons Fractals 136, 109860 (2020)

    Article  MathSciNet  Google Scholar 

  22. Zhang, Z.: A novel COVID-19 mathematical model with fractional derivatives: singular and nonsingular kernels. Chaos Solitons Fractals 139, 110060 (2020)

    Article  MathSciNet  Google Scholar 

  23. Behzad, G., Atangana, A.: A new application of fractional Atangana-Baleanu derivatives: designing ABC-fractional masks in image processing. Phys. A, Stat. Mech. Appl. 542, 123516 (2020)

    Article  MathSciNet  Google Scholar 

  24. Atangana, A., Gómez-Aguilar, J.F.: Fractional derivatives with no-index law property: application to chaos and statistics. Chaos Solitons Fractals 114, 516–535 (2018)

    Article  MathSciNet  MATH  Google Scholar 

  25. Atangana, A.: Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system. Chaos Solitons Fractals 102, 396–406 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  26. Toufik, M., Atangana, A.: New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models. Eur. Phys. J. Plus 132, 444 (2017). https://doi.org/10.1140/epjp/i2017-11717-0

    Article  Google Scholar 

  27. Abdo, M.S., Panchal, S.K.: Fractional integro-differential equations involving ψ-Hilfer fractional derivative. Adv. Appl. Math. Mech. 11(1), 1–22 (2019)

    MathSciNet  Google Scholar 

  28. Abdo, M.S., Panchal, S.K., Saeed, A.M.: Fractional boundary value problem with ψ-Caputo fractional derivative. Proc. Math. Sci. 129(5), 65 (2019)

    Article  MathSciNet  MATH  Google Scholar 

  29. Ahmad, B., Ntouyas, S.K., Alsaedi, A.: Sequential fractional differential equations and inclusions with semi-periodic and nonlocal integro-multipoint boundary conditions. J. King Saud Univ., Sci. 31(2), 184–193 (2019)

    Article  Google Scholar 

  30. Gu, H., Trujillo, J.J.: Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput. 257, 344–354 (2015)

    MathSciNet  MATH  Google Scholar 

  31. Abdeljawad, T.: A Lyapunov type inequality for fractional operators with nonsingular Mittag-Leffler kernel. J. Inequal. Appl. 2017(1), 130 (2017)

    Article  MathSciNet  MATH  Google Scholar 

  32. Abdo, M.S., Shah, K., Wahash, H.A., Panchal, S.K.: On a comprehensive model of the novel coronavirus (Covid-19) under Mittag-Leffler derivative. Chaos Solitons Fractals 2020, 109867 (2020). https://doi.org/10.1016/j.chaos.2020.109867

    Article  MathSciNet  Google Scholar 

  33. Abdo, M.S., Abdeljawad, T., Shah, K., Jarad, F.: Study of impulsive problems under Mittag-Leffler power law. Heliyon 6(10), e05109 (2020)

    Article  Google Scholar 

  34. Redhwan, S.S., Abdo, M.S., Shah, K., Abdeljawad, T., Dawood, S., Abdo, H.A., Shaikhh, S.L.: Mathematical modeling for the outbreak of the coronavirus (COVID-19) under fractional nonlocal operator. Results Phys. 19, 103610 (2020). https://doi.org/10.1016/j.rinp.2020.103610

    Article  Google Scholar 

  35. Thabet, S.T.M., Abdo, M.S., Shah, K., Abdeljawad, T.: Study of transmission dynamics of COVID-19 mathematical model under ABC fractional order derivative. Results Phys. 19, 103507 (2020). https://doi.org/10.1016/j.rinp.2020.103507

    Article  Google Scholar 

  36. Yavuz, M., Ozdemir, N., Baskonus, H.M.: Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel. Eur. Phys. J. Plus 133(6), 215 (2018)

    Article  Google Scholar 

  37. Kucche, K.D., Sutar, S.T.: Analysis of nonlinear fractional differential equations involving Atangana-Baleanu-Caputo derivative. Chaos Solitons Fractals 143, 110556 (2020)

    Article  MathSciNet  Google Scholar 

  38. Sutar, S.T., Kucche, K.D.: On nonlinear hybrid fractional differential equations with Atangana-Baleanu-Caputo derivative (2020). arXiv:2007.11034. arXiv preprint

  39. Ockendon, J.R., Tayler, A.B.: The dynamics of a current collection system for an electric locomotive. Proc. R. Soc. Lond. Ser. A, Math. Phys. Sci. 322(1551), 447–468 (1971)

    Google Scholar 

  40. Balachandran, K., Kiruthika, S., Trujillo, J.J.: Existence of solutions of nonlinear fractional pantograph equations. Acta Math. Sci. 33(3), 712–720 (2013)

    Article  MathSciNet  MATH  Google Scholar 

  41. Vivek, D., Kanagarajan, K., Sivasundaram, S.: Theory and analysis of nonlinear neutral pantograph equations via Hilfer fractional derivative. Nonlinear Stud. 24(3), 699–712 (2017)

    MathSciNet  MATH  Google Scholar 

  42. Anguraj, A., Vinodkumar, A., Malar, K.: Existence and stability results for random impulsive fractional pantograph equations. Filomat 30(14), 3839–3854 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  43. Vivek, D., Kanagarajan, K., Harikrishnan, S.: Dynamics and stability of Hilfer-Hadamard type fractional pantograph equations with boundary conditions. J. Nonlinear Anal. Appl. 2018(1), 1–13 (2018)

    MATH  Google Scholar 

  44. Elsayed, E., Harikrishnan, S., Kanagarajan, K.: Analysis of nonlinear neutral pantograph differential equations with Hilfer fractional derivative. MathLAB 1, 231–240 (2018)

    MATH  Google Scholar 

  45. Harikrishnan, S., Ibrahim, R., Kanagarajan, K.: Establishing the existence of Hilfer fractional pantograph equations with impulses. Fundam. J. Math. Appl. 1(1), 36–42 (2018)

    Google Scholar 

  46. Ahmed, I., Kumam, P., Shah, K., et al.: Stability results for implicit fractional pantograph differential equations via Hilfer fractional derivative with a nonlocal Riemann-Liouville fractional integral condition. Mathematics 8(1), 94 (2020)

    Article  Google Scholar 

  47. Borisut, P., Kumam, P., Ahmed, I., Sitthithakerngkiet, K.: Nonlinear Caputo fractional derivative with nonlocal Riemann-Liouville fractional integral condition via fixed point theorems. Symmetry 11(6), 829 (2019)

    Article  Google Scholar 

  48. Ahmed, I., Kumam, P., Shah, K., Borisut, P., Sitthithakerngkiet, K., Ahmed, D.: Stability results for implicit fractional pantograph differential equations via ϕ-Hilfer fractional derivative with a nonlocal Riemann-Liouville fractional integral condition. Mathematics 8(1), 94 (2020)

    Article  Google Scholar 

  49. Borisut, P., Kumam, P., Ahmed, I., Jirakitpuwapat, W.: Existence and uniqueness for ψ-Hilfer fractional differential equation with nonlocal multi-point condition. Math. Methods Appl. Sci. (2020). https://doi.org/10.1002/mma.6092

    Article  Google Scholar 

  50. Ahmed, I., Kumam, P., Jarad, F., Borisut, P., Sitthithakerngkiet, K., Ibrahim, A.: Stability analysis for boundary value problems with generalized nonlocal condition via Hilfer-Katugampola fractional derivative. Adv. Differ. Equ. 2020(1), 225 (2020)

    Article  MathSciNet  Google Scholar 

  51. Borisut, P., Kumam, P., Ahmed, I., Sitthithakerngkiet, K.: Positive solution for nonlinear fractional differential equation with nonlocal multi-point condition. Fixed Point Theory 21(2), 427–440 (2020)

    Article  MATH  Google Scholar 

  52. Zhou, Y., Wang, J., Zhang, L.: Basic Theory of Fractional Differential Equations. World Scientific, Singapore (2016)

    Book  Google Scholar 

Download references

Acknowledgements

All authors have read and approved the last draft of the paper.

Funding

This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed equally to this work. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Thabet Abdeljawad or Manar A. Alqudah.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Abdo, M.S., Abdeljawad, T., Kucche, K.D. et al. On nonlinear pantograph fractional differential equations with Atangana–Baleanu–Caputo derivative. Adv Differ Equ 2021, 65 (2021). https://doi.org/10.1186/s13662-021-03229-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13662-021-03229-8

MSC

Keywords