Fractional physical models based on falling body problem

Bahar Acay; Ramazan Ozarslan; Erdal Bas; Bahar Acay; Ramazan Ozarslan; Erdal Bas

doi:10.3934/math.2020170

AIMS Mathematics

2020, Volume 5, Issue 3: 2608-2628. doi: 10.3934/math.2020170

Previous Article Next Article

Research article

Fractional physical models based on falling body problem

Department of Mathematics, Science Faculty, Firat University, 23119 Elazig, Turkey

Received: 19 December 2019 Accepted: 26 February 2020 Published: 12 March 2020
MSC : 26A33, 97M50, 70E17

This article is devoted to investigate the fractional falling body problem relied on Newton's second law. We analyze this physical model by means of Atangana-Baleanu fractional derivative in the sense of Caputo (ABC), generalized fractional derivative introduced by Katugampola and generalized ABC containing the Mittag-Leffler function with three parameters

$\mathbb{E}_{\alpha, \mu}^{\gamma}(.)$ . For that purpose, the Laplace transform (LT) is utilized to obtain fractional solutions. In order to maintain the dimensionality of the physical parameter in the model, we employ an auxiliary parameter

$\sigma$ having a relation with the order of fractional operator. Moreover, simulation analysis is carried out by comparing the underlying fractional derivatives with traditional one to grasp the virtue of the results.

Keywords:

Citation: Bahar Acay, Ramazan Ozarslan, Erdal Bas. Fractional physical models based on falling body problem[J]. AIMS Mathematics, 2020, 5(3): 2608-2628. doi: 10.3934/math.2020170

Related Papers:

[1]	Sitalakshmi Venkatraman, Ramanathan Venkatraman . Big data security challenges and strategies. AIMS Mathematics, 2019, 4(3): 860-879. doi: 10.3934/math.2019.3.860
[2]	Atiqe Ur Rahman, Muhammad Saeed, Mazin Abed Mohammed, Alaa S Al-Waisy, Seifedine Kadry, Jungeun Kim . An innovative fuzzy parameterized MADM approach to site selection for dam construction based on sv-complex neutrosophic hypersoft set. AIMS Mathematics, 2023, 8(2): 4907-4929. doi: 10.3934/math.2023245
[3]	Pengyu Liu, Jie Jian . Effects of network topology and trait distribution on collective decision making. AIMS Mathematics, 2023, 8(5): 12287-12320. doi: 10.3934/math.2023619
[4]	Shahryar Sorooshian, Seyedh Mahboobeh Jamali, Nader Ale Ebrahim . Performance of the decision-making trial and evaluation laboratory. AIMS Mathematics, 2023, 8(3): 7490-7514. doi: 10.3934/math.2023376
[5]	Erdal Bayram, Gülşah Çelik, Mustafa Gezek . An advanced encryption system based on soft sets. AIMS Mathematics, 2024, 9(11): 32232-32256. doi: 10.3934/math.20241547
[6]	Esmail Hassan Abdullatif Al-Sabri, Muhammad Rahim, Fazli Amin, Rashad Ismail, Salma Khan, Agaeb Mahal Alanzi, Hamiden Abd El-Wahed Khalifa . Multi-criteria decision-making based on Pythagorean cubic fuzzy Einstein aggregation operators for investment management. AIMS Mathematics, 2023, 8(7): 16961-16988. doi: 10.3934/math.2023866
[7]	Muhammad Arshad, Muhammad Saeed, Khuram Ali Khan, Nehad Ali Shah, Wajaree Weera, Jae Dong Chung . A robust MADM-approach to recruitment-based pattern recognition by using similarity measures of interval-valued fuzzy hypersoft set. AIMS Mathematics, 2023, 8(5): 12321-12341. doi: 10.3934/math.2023620
[8]	Muhammad Ihsan, Muhammad Saeed, Atiqe Ur Rahman, Hüseyin Kamacı, Nehad Ali Shah, Wajaree Weera . An MADM-based fuzzy parameterized framework for solar panels evaluation in a fuzzy hypersoft expert set environment. AIMS Mathematics, 2023, 8(2): 3403-3427. doi: 10.3934/math.2023175
[9]	Esra Korkmaz . An innovative algorithm based on weighted fuzzy soft multisets and its application in selecting optimal construction materials. AIMS Mathematics, 2024, 9(10): 27512-27534. doi: 10.3934/math.20241336
[10]	Jia-Bao Liu, Rashad Ismail, Muhammad Kamran, Esmail Hassan Abdullatif Al-Sabri, Shahzaib Ashraf, Ismail Naci Cangul . An optimization strategy with SV-neutrosophic quaternion information and probabilistic hesitant fuzzy rough Einstein aggregation operator. AIMS Mathematics, 2023, 8(9): 20612-20653. doi: 10.3934/math.20231051

Abstract

1. Introduction

Fractional differential equations with various types of fractional derivatives arise in modeling some dynamical processes (see, for example, ^[15] for the globally projective synchronization of Caputo fractional-order quaternion-valued neural networks with discrete and distributed delays, ^[18] for the quasi-uniform synchronization issue for fractional-order neural networks with leakage and discrete delays and ^[11] for Mittag-Leffler stability and adaptive impulsive synchronization of fractional order neural networks in quaternion field). In contrast to the classical derivative the fractional derivative is nonlocal and it depends significantly on its lower limit. As it is mentioned in ^[13], this leads to some obstacles for studying impulsive fractional differential equations.

Since many phenomena are characterized by abrupt changes at certain moments it is important to consider differential equations with impulses. In the literature there are two main approaches used to introduce impulses to fractional equations:

(i) With a fixed lower limit of the fractional derivative at the initial time- the fractional derivative of the unknown function has a lower limit equal to the initial time point over the whole interval of study;

(ii) With a changeable lower limit of the fractional derivative at each time of impulse- the fractional derivative on each interval between two consecutive impulses is changed because the lower limit of the fractional derivative is equal to the time of impulse.

Both interpretations of impulses are based on corresponding interpretations of impulses in ordinary differential equation, which coincide in the case of integer derivatives. However this is not the case for fractional derivatives. In the literature both types of interpretations are discussed and studied for Caputo fractional differential equations of order ${\alpha}\in(0, 1)$ . We refer the reader to the papers ^{[6,7,12,13,16]} as well as the monograph ^[3].

We note in the case of the Caputo fractional derivative there is a similarity between both the initial conditions and the impulsive condition between fractional equations and ordinary equations (see, for example, ^[10] concerning the impulsive control law for the Caputo delay fractional-order neural network model). However for Riemann-Liouville fractional differential equations both the initial condition and impulsive conditions have to be appropriately given (which is different in the case of ordinary derivatives as well as the case of Caputo fractional derivatives). Riemann-Liouville fractional differential equations are considered, for example, in ^[1,2] for integral presentation of solutions in the case of the fractional order ${\alpha}\in(0, 1)$ , ^[5,8,17] for the case of the fractional order ${\alpha}\in(1, 2)$ .

In ^[14] the authors studied the following coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives of the form

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} \mathcal{D}^{{\alpha}}u(t)-\phi_1(t, \mathcal{I}^{\alpha} u(t), \mathcal{I}^{\beta} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_i, \ \ i = 1, 2, \dots, p, \\ \Delta u(t_j)-\mathcal{E}_j(u(t_j)) = 0, \ \ \ \ \ \ \Delta u'(t_j)-\mathcal{E}_j^*(u(t_j)) = 0, \ \ j = 1, 2, \dots, p, \\ \nu_1\mathcal{D}^{{\alpha}-2}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_2 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} \mathcal{D}^{{\beta}}v(t)-\phi_2(t, \mathcal{I}^{\alpha} u(t), \mathcal{I}^{\beta} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_k, \ \ k = 1, 2, \dots, q, \\ \Delta v(t_k)-\mathcal{E}_k(v(t_k)) = 0, \ \ \ \ \ \ \Delta v'(t_k)-\mathcal{E}_j^*(v(t_k)) = 0, \ \ k = 1, 2, \dots, q, \\ \nu_3\mathcal{D}^{{\beta}-2}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_4 \mathcal{I}^{{\beta}-1} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(1.1)

where ${\alpha}, {\beta}\in (1, 2]$ , $I = [0, T]$ , $\phi_1, \phi_2:I\times {\mathbb{R}} \times {\mathbb{R}}\to {\mathbb{R}}$ are continuous functions, $\Delta u(t_j) = u(t_j^+)-u(t_j^-), \ \ \Delta u'(t_j) = u'(t_j^+)-u'(t_j^-)$ , $\Delta v(t_k) = v(t_k^+)-v(t_k^-), \ \ \Delta v'(t_j) = v'(t_k^+)-v'(t_k^-)$ , where $u(t_j^+), v(t_k^+)$ and $v(t_j^-), v(t_k^-)$ are the right limits and left limits respectively, $\mathcal{E}_j, \mathcal{E}^*_j, \mathcal{E}_k, \mathcal{E}^*_k:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, and $\mathcal{D}^{{\alpha}}, \mathcal{I}^{{\alpha}}$ are the ${\alpha}$ -order Riemann-Liouville fractional derivative and integral operators respectively and $\mathcal{D}^{{\beta}-2} = \mathcal{I}^{2-{\beta}}$ .

Since fractional integrals and derivatives have memories, and their lower limits are very important we will use the notations $^{RL}D^{{\alpha}}_{a, t}$ and $I^{{\beta}}_{a, t}$ , respectively, instead of $\mathcal{D}^{{\alpha}}$ and $\mathcal{I}^{{\beta}}$ , i.e. the Riemann-Liouville fractional derivative is defined by (see, for example ^[9])

$\begin{equation} ^{RL}D^{{\alpha}}_{a, t}u(t) = \frac{1}{\Gamma(2-{\alpha})}\Big(\frac{d}{dt}\Big)^2\int_a^t\frac{u(s)}{(t-s)^{{\alpha}-1}}ds, \ \ \ t > a, \ \ {\alpha}\in(1, 2), \end{equation}$

(1.2)

and the Riemann-Liouville fractional integral $I^{{\beta}}_{a, t}$ of order ${\alpha} > 0$ is defined by (see, for example, ^[9])

$\begin{equation} I^{{\beta}}_{a, t}u(t) = \frac{1}{\Gamma({\beta})}\int _a^t(t-s)^{{\beta}-1}u(s)ds, t > a, \end{equation}$

(1.3)

where $a\geq 0$ , ${\beta} > 0$ are given numbers.

Note there are some unclear parts in the statement of coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives (1.1), such as:

a). The presence of two different integers $p$ and $q$ in (1.1) leads to different domains of both the unknown functions $u$ and $v$ . For example, in Corollary 1 ^[14] the solutions $u(t)$ and $v(t)$ are defined on $[0, t_{p+1}]$ and $[0, t_{q+1}]$ respectively, which causes some problems in the definitions of formulas (3.7) or (3.8) (^[14]);

b). The impulsive functions $\mathcal{E}_j, \mathcal{E}_k, \ j = 1, 2, \dots, p, \ k = 1, 2, \dots, q$ are assumed different but they are not (it is clear for example, for $j = k = 1$ ). The same is about the functions $\mathcal{E}^*_j, \mathcal{E}^*_k, \ j = 1, 2, \dots, p, \ k = 1, 2, \dots, q$ .

In this paper we sort out the above mentioned points by setting up the cleared statement of the boundary value problem with the Riemann-Liouville (RL) fractional integral for the impulsive Riemann-Liouville fractional differential equation studied in ^[14], and we prove a new the integral presentations of the solutions. To be more general, we study two different interpretations for the presence of impulses in fractional differential equations. The first one is the case of the fixed lower limit of the RL fractional derivative at the initial time 0 and the second one is the case of the changed lower limit of the fractional RL derivative at any point of impulse. In both cases the integral presentation of the solution is provided.

2. RL-fractional derivatives with fixed lower limit

2.1. Statement of the problem

Define two different sequences of points of impulses

$0 = t_0 < t_1 < t_2 < \dots < t_p < t_{p+1} = T\ \ \mbox{and}\ \ \ 0 = {\tau}_0 < {\tau}_1 < {\tau}_2 < \dots < {\tau}_q < {\tau}_{q+1} = T,$

where $p, q$ are given natural numbers.

We will consider the following nonlinear boundary value problem for the coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives with a lower limit at $0$

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} ^{RL}D^{{\alpha}}_{0, t}u(t)-\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = t_i, \ \ i = 1, 2, \dots, p, \\ \Delta u(t_j)-\mathcal{E}_j(u(t_j)) = 0, \ \ \ \ \ \ \Delta u'(t_j)-\mathcal{E}_j^*(u(t_j)) = 0, \ \ j = 1, 2, \dots, p, \\ - \mathcal{I}^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1}_{0, t} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} ^{RL}D^{{\beta}}_{0, t}v(t)-\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t}v(t)) = 0\ \ \ \mbox{for}\ t\in I, \ t\not = {\tau}_k, \ \ k = 1, 2, \dots, q, \\ \Delta v({\tau}_k)-\mathcal{S}_k(v({\tau}_k)) = 0, \ \ \ \ \ \ \Delta v'({\tau}_k)-\mathcal{S}_k^*(v({\tau}_k)) = 0, \ \ k = 1, 2, \dots, q, \\ \mathcal{I}^{2-{\beta}}_{0, t}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_2 \mathcal{I}^{{\beta}-1}_{0, t} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(2.1)

where ${\alpha}, {\beta}\in (1, 2]$ , $I = [0, T]$ , $\phi_1, \phi_2:I\times{\mathbb{R}}\times{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $\Delta u(t_j) = u(t_j^+)-u(t_j^-), \ \ \Delta u'(t_j) = u'(t_j^+)-u'(t_j^-)$ , $\Delta v({\tau}_k) = v({\tau}_k^+)-v({\tau}_k^-), \ \ \Delta v'({\tau}_k) = v'({\tau}_k^+)-v'({\tau}_k^-)$ , where $u(t_j^+), v({\tau}_k^+), u'(t_j^+), v'({\tau}_k^+)$ and $u(t_j^-)$ , $v(t_k^-)$ , $u'(t_j^-)$ , $v'(t_k^-)$ are the right limits and left limits respectively, $\mathcal{E}_j, \mathcal{E}^*_j, \mathcal{S}_k, \mathcal{S}^*_k:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, and $^{RL}D^{{\beta}}_{0, t}, \mathcal{I}^{{\alpha}}_{0, t}$ are the ${\alpha}$ -order Riemann-Liouville fractional derivative and integral operators, respectively, $\mu_i, \nu_i, u_k, v_k, \ i = 1, 2,$ are given constants.

In the statement of the problem (2.1) some parts of (1.1) are cleared: there are two different points of impulses; the lower limits of the fractional integrals and fractional integrals are written; different functions at different points of impulses are used.

In our proofs we will use the following well known properties for fractional integrals (see, for example ^[9]).

$\begin{equation} \label{{7}}\begin{split} &I^{{\alpha}}_{a, t}I^{{\beta}}_{a, t}u(t) = I^{{\alpha}+{\beta}}_{a, t}u(t), \ \ \ {\alpha}, {\beta} > 0, \\ &I^{{\alpha}}_{a, t}(t-a)^q = \frac{\Gamma(q+1)}{\Gamma(q+{\alpha}+1)}(t-a)^{q+{\alpha}}, \ \ {\alpha} > 0, \ q > -1.\ \ t > a.\end{split} \end{equation}$

(2.2)

We will apply the following auxiliary result which is a generalization of the result in ^[4] for an arbitrary lower limit of the fractional derivative:

Lemma 1. (^[4]). The general solution of the Riemann-Liouville fractional differential equation

$\begin{equation} \label{{9}}^{RL}D^{{\alpha}}_{a, t}w(t) = g(t), \ \ \ t\in (a, T], \ \ \ {\alpha}\in(1, 2) \end{equation}$

(2.3)

is given by

$\begin{equation} \label{{400}}\begin{split} w(t)& = c_1(t-a)^{{\alpha}-1}+c_0(t-a)^{{\alpha}-2}+I^{{\alpha}}_{a, t}\ g(t)\\ & = c_1(t-a)^{{\alpha}-1}+c_0(t-a)^{{\alpha}-2}+ \frac{1}{\Gamma({\alpha})}\int_a^t(t-s)^{{\alpha}-1}g(s)ds, \ t\in(a, T], \end{split} \end{equation}$

(2.4)

where $c_0, c_1, a\geq 0$ are arbitrary real constants.

2.2. Integral presentation of the linear problem

We will consider an appropriate boundary value problem for a scalar impulsive linear equation, we will prove a formula for its solution and later we will apply it to obatin the main result.

Consider the following boundary value problem for the linear impulsive fractional differential equation with Riemann-Liouville derivatives of the form

$\begin{equation} \label{{999}}\begin{split}&^{RL}D^{{\alpha}}_{0, t}u(t) = f(t), \ \ \ t\in (0, T], \ \ t\not = t_j, \ j = 1, 2, \dots p, \ \ {\alpha}\in(1, 2), \\ &\Delta u(t_j) = \mathcal{E}_j(u(t_j)), \ \ \ \ \ \ \Delta u'(t_j) = \mathcal{E}_j^*(u(t_j))\ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(2.5)

where $f:\ [0, T]\to {\mathbb{R}}$ is a continuous function, $\mathcal{E}_j, \mathcal{E}^*_j:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $u_1, u_2\in{\mathbb{R}}$ .

Lemma 2. The solution of (2.5) satisfies the integral equation

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds, & t\in (0, t_1], \\ c_0t^{{\alpha}-1}+\Big(\sum _{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum _{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum _{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum _{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, & t\in (t_j, t_{j+1}], \\ \end{array} \right.$

where

$\begin{equation} \label{{302}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}f(s)ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}.\end{split} \end{equation}$

(2.6)

Proof. We will use induction.

For $t\in (0, t_1]$ we apply Lemma 1 with $a = 0$ and we get

$\begin{equation} \label{{200}} u(t) = c_0t^{{\alpha}-1}+c_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(2.7)

and

$\begin{equation} \label{{209}} u'(t) = c_0({\alpha}-1)t^{{\alpha}-2}+c_1({\alpha}-2)t^{{\alpha}-3}+\frac{{\alpha}-1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.8)

From the initial condition $I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1$ and equalities (2.2), (2.7) we get

$I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = c_0\frac{\Gamma({\alpha})}{\Gamma(1)}t|_{t = 0}+c_1\frac{\Gamma({\alpha}-1)}{\Gamma(0)}|_{t = 0}+I^{2}_{0, t}f(t)|_{t = 0} = c_1\Gamma({\alpha}-1),$

i.e. $c_1 = \frac{u_1}{\Gamma({\alpha}-1)}.$

For $t\in (t_1, t_2]$ by Lemma 1 with $a = 0$ we get

$\begin{equation} \label{{210}} u(t) = b_0t^{{\alpha}-1}+b_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(2.9)

and

$\begin{equation} \label{{219}} u'(t) = b_0({\alpha}-1)t^{{\alpha}-2}+b_1({\alpha}-2)t^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.10)

From the impulsive condition $u(t_1+0)-u(t_1-0) = \mathcal{E}_1(u(t_1))$ we obtain

$\begin{equation} \label{{300}} \begin{split} u(t_1^+) = &b_0t_1^{{\alpha}-1}+b_1t_1^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^{t_1}(t_1-s)^{{\alpha}-1}f(s)ds\\ &-c_0t_1^{{\alpha}-1}-\frac{u_1}{\Gamma({\alpha}-1)}t_1^{{\alpha}-2}-\frac{1}{\Gamma({\alpha})}\int_0^{t_1}(t_1-s)^{{\alpha}-1}f(s)ds\\ = &(b_0-c_0)t_1^{{\alpha}-1}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})t_1^{{\alpha}-2} = \mathcal{E}_1(u(t_1)). \end{split} \end{equation}$

(2.11)

and from the impulsive condition $u'(t_1+0)- u'(t_1) = \mathcal{E}_1^*(u(t_1))$ we get

$\begin{equation} \label{{310}} \begin{split} u'(t_1^+) = &b_0({\alpha}-1)t_1^{{\alpha}-2}+b_1({\alpha}-2)t_1^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^{t_1}(t_1-s)^{{\alpha}-2}f(s)ds\\ &-c_0({\alpha}-1)t_1^{{\alpha}-2}-\frac{u_1}{\Gamma({\alpha}-1)}({\alpha}-2)t_1^{{\alpha}-3}-\frac{1}{\Gamma({\alpha}-1)}\int_0^{t_1}(t_1-s)^{{\alpha}-2}f(s)ds\\ = &(b_0-c_0)({\alpha}-1)t_1^{{\alpha}-2}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})({\alpha}-2)t_1^{{\alpha}-3} = \mathcal{E}_1^*(u(t_1)).\end{split} \end{equation}$

(2.12)

Thus we get the linear system w.r.t. $b_0$ and $b_1$

$\begin{equation} \label{{8}}\begin{split} &(b_0-c_0)({\alpha}-1)t_1^{{\alpha}-2}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})({\alpha}-2)t_1^{{\alpha}-3} = \mathcal{E}_1^*(u(t_1))\\ &(b_0-c_0)t_1^{{\alpha}-1}+(b_1-\frac{u_1}{\Gamma({\alpha}-1)})t_1^{{\alpha}-2} = \mathcal{E}_1(u(t_1))\nonumber\end{split} \end{equation}$

$b_0 = c_0+t_1^{2-{\alpha}}\mathcal{E}_1^*(u(t_1))-({\alpha}-2)t_1^{1-{\alpha}}\mathcal{E}_1(u(t_1))$

$b_1 = ({\alpha}-1)t_1^{2-{\alpha}}\mathcal{E}_1(u(t_1))-t_1^{3-{\alpha}}\mathcal{E}_1^*(u(t_1))+\frac{u_1}{\Gamma({\alpha}-1)}.$

Therefore,

$\begin{equation} \label{{220}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(t_1^{2-{\alpha}}\mathcal{E}_1^*(u(t_1))-({\alpha}-2)t_1^{1-{\alpha}}\mathcal{E}_1(u(t_1))\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(({\alpha}-1)t_1^{2-{\alpha}}\mathcal{E}_1(u(t_1))-t_1^{3-{\alpha}}\mathcal{E}_1^*(u(t_1))+\frac{u_1}{\Gamma({\alpha}-1)}\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_1, t_2].\end{split} \end{equation}$

(2.13)

Assume the integral presentation of $u(t)$ is correct on $(t_{j-1}, t_j]$ , i.e

$\begin{equation} \label{{2209}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))-({\alpha}-2)\sum\limits_{k = 1}^{j-1}t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(({\alpha}-1)\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{j-1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))+\frac{u_1}{\Gamma({\alpha}-1)}\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_{j-1}, t_{j}].\end{split} \end{equation}$

(2.14)

Denote

$m_0 = c_0+\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))-({\alpha}-2)\sum\limits_{k = 1}^{j-1}t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))$

and

$m_1 = ({\alpha}-1)\sum\limits_{k = 1}^{j-1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{j-1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))+\frac{u_1}{\Gamma({\alpha}-1)}.$

Let $t\in (t_j, t_{j+1}], \ j = 2, \dots, p,$ . By Lemma 1 with $a = 0$ we get

$\begin{equation} \label{{1210}} u(t) = k_0t^{{\alpha}-1}+k_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \end{equation}$

(2.15)

and

$\begin{equation} \label{{1219}} u'(t) = k_0({\alpha}-1)t^{{\alpha}-2}+k_1({\alpha}-2)t^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{0}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(2.16)

From the impulsive conditions and the equality (2.14) we obtain the linear system w.r.t. $k_0$ and $k_1$

$\begin{equation} \label{{18}}\begin{split} &(k_0-m_0)({\alpha}-1)t_j^{{\alpha}-2}+(k_1-m_1)({\alpha}-2)t_j^{{\alpha}-3} = \mathcal{E}_2^*(u(t_2))\\ &(k_0-m_0)t_j^{{\alpha}-1}+(k_1-m_1)t_j^{{\alpha}-2} = \mathcal{E}_j(u(t_j))\nonumber\end{split} \end{equation}$

$k_0 = c_0+\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]$

$k_1 = \frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k)).$

Therefore,

$\begin{equation} \label{{1220}} \begin{split}u(t)& = c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ &\ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ &\ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}f(s)ds, \ \ \ t\in(t_j, t_{j+1}], \ j = 1, 2, \dots, p.\end{split} \end{equation}$

(2.17)

From the boundary condition $\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2$ we get

$\mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = \frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds$

and

$\begin{equation} \label{{888}}\begin{split} &\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T}\\ & = \mu_1 c_0T^{{\alpha}-1}+\mu_1 \Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)T^{{\alpha}-1}\\ &\ \ \ \ \ +\mu_1 \Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{{\alpha}-2}\\ &\ \ \ \ \ +\mu_1 \frac{1}{\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ +\nu_1\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds\\ & = u_2. \end{split} \end{equation}$

(2.18)

From (2.18) we have (2.6).

We will give an example to illustrate the claim of Lemma 2.

Example 1. Consider the following boundary value problem for the scalar RL fractional differential equation with an impulse at $t = 1$

$\begin{equation} \label{{99987}}\begin{split}&^{RL}D^{1.5}_{0, t}u(t) = t, \ \ \ t\in (0, 1]\cup (1, 2], \\ &\Delta u(1) = 1, \ \ \ \ \ \ \Delta u'(1) = 0, \\ &I^{0.5}_{0, t}u(t)|_{t = 0} = 0, \ \ \ \ \ u(t)|_{t = 2}+ I^{0.5}_{0, t}u(t)|_{t = 2} = 1.\end{split} \end{equation}$

(2.19)

The solution of (2.19) satisfies the integral equation

$u(t) = \left\{ \begin{array}{ll} c_0t^{0.5}+\frac{t^{-0.5}}{\Gamma(0.5)}+\frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (0, 1] \\ c_0t^{0.5}+0.5t^{0.5} +\Big(\frac{1}{\Gamma(0.5)}+0.5\Big)t^{-0.5}+\frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (1, 2], \\ \end{array} \right.$

where

$\begin{equation} \label{{30298}} \begin{split} c_0 & = 0.25-\frac{1.50849}{2^{0.5}\Gamma(1.5)}-\frac{ 1}{2^{0.5}\Gamma(0.5)}\int _0^2(2-s)^{-0.5}u(s)ds.\end{split} \end{equation}$

(2.20)

Consider the boundary value problem for the nonlinear impulsive fractional integro-differential equations with Riemann-Liouville derivatives of the form

$\begin{equation} \label{{88}}\begin{split}&^{RL}D^{{\alpha}}_{0, t}u(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t)), \ \ \ t\in (0, T], \ \ t\not = t_j, \ j = 1, 2, \dots p, \ \ {\alpha}\in(1, 2)\\ &\Delta u(t_j) = \mathcal{E}_j(u(t_j)), \ \ \ \ \ \ \Delta u'(t_j) = \mathcal{E}_j^*(u(t_j))\ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(2.21)

where $\phi_1:\ [0, T]\times {\mathbb{R}}\to {\mathbb{R}}$ is a continuous function, $\mathcal{E}_j, \mathcal{E}^*_j:{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $u_1, u_2\in{\mathbb{R}}$ .

Corollary 1. The solution of (2.21) satisfies the integral equation

$\begin{equation} \label{{66}} \begin{split} u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), )ds , \ \ t\in (0, t_1] \\ c_0t^{{\alpha}-1}+\Big(\sum _{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum _{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum _{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum _{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s))ds, \ \ \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.22)

where

$\begin{equation} \label{{87}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}.\end{split} \end{equation}$

(2.23)

The proof of Corollary 1 is similar to the one of Lemma 2 with $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t))$ and we omit it.

Remark 1. Corollary 1 and the integral presentation (2.22) correct Theorem 3.1 and the formula (3.2) ^[14]. The main mistake in the proof of formula (3.2) ^[14] is the incorrect application of Lemma 1 with $a = 0$ on $(t_1, t_2]$ and taking the lower limit of the integral in (3.5) ^[14] incorrectly at $t_1 (\sigma_1)$ instead of 0. A similar comment applies to all the other intervals $(t_j, t_{j+1}]$ .

2.3. Integral presentation of the problem (2.1)

Following the proof of Lemma 2 and the integral presentation (2.22) of problem (2.21), we have the following result:

Theorem 1. The solution of (2.1) satisfies the integral equations

$\begin{equation} \label{{44}}\begin{split} u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ t\in (0, t_1], \\ c_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.24)

and

$\begin{equation} \label{{45}}\begin{split} v(t) = \left\{ \begin{array}{ll} b_0t^{{\alpha}-1}+\frac{v_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ t\in (0, t_1] \\ b_0t^{{\alpha}-1}+\Big(\sum\limits_{k = 1}^j[t_k^{2-{\alpha}}\mathcal{S}_k^*(v(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^j[t_k^{1-{\alpha}}\mathcal{S}_k(v(t_k))]\Big)t^{{\alpha}-1}\\ \ \ \ \ \ +\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^jt_k^{2-{\alpha}}\mathcal{S}_k(v(t_k))-\sum\limits_{k = 1}^jt_k^{3-{\alpha}}\mathcal{S}_k^*(v(t_k))\Big)t^{{\alpha}-2}\\ \ \ \ \ \ +\frac{1}{\Gamma({\alpha})}\int_{0}^t(t-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds, \ \ \ \ t\in (t_j, t_{j+1}], \\ \end{array} \right. \end{split} \end{equation}$

(2.25)

where

$\begin{equation} \label{{5302}} \begin{split} c_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{E}_k^*(u(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{E}_k(u(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{u_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{E}_k(u(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{E}_k^*(u(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds+\frac{u_2}{\mu_1}\\ b_0 & = -\Big(\sum\limits_{k = 1}^{p+1}[t_k^{2-{\alpha}}\mathcal{S}_k^*(v(t_k))]-({\alpha}-2)\sum\limits_{k = 1}^{p+1}[t_k^{1-{\alpha}}\mathcal{S}_k(v(t_k))]\Big)\\ &\ \ \ \ \ -\Big(\frac{v_1}{\Gamma({\alpha}-1)}+({\alpha}-1)\sum\limits_{k = 1}^{p+1}t_k^{2-{\alpha}}\mathcal{S}_k(v(t_k))-\sum\limits_{k = 1}^{p+1}t_k^{3-{\alpha}}\mathcal{S}_k^*(v(t_k))\Big)T^{-1}\\ &\ \ \ \ \ -\frac{1}{T^{{\alpha}-1}\Gamma({\alpha})}\int_{0}^T(T-s)^{{\alpha}-1}\phi_2(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds\\ &\ \ \ \ \ \ -\frac{\nu_2}{\mu_2T^{{\alpha}-1}\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}v(s)ds+\frac{v_2}{\mu_2}.\end{split} \end{equation}$

(2.26)

The proof is similar to the one of Lemma 2 applied twice to each of the both components $u$ and $v$ of the coupled system (2.1) for impulsive points $t_i, \ i = 1, 2, \dots, p$ and ${\tau}_i, \ i = 1, 2, \dots, p$ and the functions $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ and $f(t) = \phi_2(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ respectively.

Remark 2. Note the integral presentation (2.24), (2.25) of the solutions of the coupled system is the correction of Corollary 1 and integral presentation (3.7), (3.8) in ^[14].

3. RL fractional derivatives with changeable lower limits

3.1. Statement of the problem with changeable lower limit of the RL fractional derivatives

Consider the following nonlinear boundary value problem for the coupled impulsive fractional integro-differential equations with Riemann-Liouville derivatives with lower limits at impulsive points $t_i, \ i = 0, 1, 2, \dots, p-1$ and ${\tau}_k, \ k = 0, 1, 2, \dots, q-1$ , respectively,

$\begin{align} \left\{ \begin{array}{ll} \left\{ \begin{array}{ll} ^{RL}D^{{\alpha}}_{t_i, t}u(t)-\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)) = 0\ \ \ \mbox{for}\ t\in (t_i, t_{i+1}]\ \ i = 0, 1, 2, \dots, p, \\ I^{2-{\alpha}}_{t_j, t}u'(t)|_{t = t_j} = \mathcal{P}_ju(t_j)+\mathcal{Q}_j, \ \ j = 1, 2, \dots, p, \\ I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t} u(t)|_{t = T} = u_2, \end{array} \right.\\ \left\{ \begin{array}{ll} ^{RL}D^{{\beta}}_{{\tau}_k, t}v(t)-\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t}v(t)) = 0\ \ \ \mbox{for}\ t\in ({\tau}_k, {\tau}_{k+1}]\ \ k = 0, 1, 2, \dots, q, \\ I^{2-{\alpha}}_{{\tau}_k, t}u'(t)|_{t = {\tau}_k} = \mathcal{P}_k^*v({\tau}_k)+\mathcal{Q}_jk^*, \ \ k = 1, 2, \dots, q, \\ I^{2-{\beta}}_{0, t}v(t)|_{t = 0} = v_1, \ \ \ \ \ \mu_2 v(t)|_{t = T}+\nu_2 I^{{\beta}-1}_{0, t} v(t)|_{t = T} = v_2, \end{array} \right.\end{array} \right. \end{align}$

(3.1)

where ${\alpha}, {\beta}\in (1, 2]$ , $\phi_1, \phi_2:I\times{\mathbb{R}}\times{\mathbb{R}}\to{\mathbb{R}}$ are continuous functions, $\mathcal{P}_j, \mathcal{Q}_j, \ j = 1, 2, \dots, p,$ and $\mathcal{P}_j^*, \mathcal{Q}_j^*, \ j = 1, 2, \dots, q,$ are real numbers, $^{RL}D^{{\alpha}}_{t_i, t}$ and $^{RL}D^{{\beta}}_{{\tau}_k, t}$ are the ${\alpha}$ -order Riemann-Liouville fractional derivatives with lower limits at $t_i$ and ${\tau}_k$ , respectively, $\mu_k, \nu_k, u_k, v_k, \ i = 1, 2,$ are given constants.

Remark 3. Note problem (3.1) differs from problem (2.1):

$-$ The lower limits of the RL fractional derivatives $^{RL}D^{{\alpha}}_{t_j, t}$ and $^{RL}D^{{\beta}}_{{\tau}_k, t}$ in (3.1) are changed at any time of impulse $t_j$ and ${\tau}_k$ , respectively.

$-$ The impulsive conditions are changed in (3.1). This is because the values of the unknown functions after the impulse, $u(t_j+0)$ and $v({\tau}_k+0)$ , respectively, are considered as initial values at that point. But the RL fractional derivative has a singularity at its lower limit. It requires the chang of the impulsive conditions for the unknown functions.

3.2. Integral presentation of the linear problem

Consider the following boundary value problem for the scalar linear impulsive fractional equation with Riemann-Liouville derivatives of the form

$\begin{equation} \label{{1999}}\begin{split}&^{RL}D^{{\alpha}}_{t_i, t}u(t) = f(t), \ \ \ \mbox{for}\ t\in (t_i, t_{i+1}], \ \ i = 0, 1, 2, \dots, p, \ \ {\alpha}\in(1, 2)\\ &I^{2-{\alpha}}_{t_j, t}u'(t)|_{t = t_j} = \mathcal{P}_ju(t_j)+\mathcal{Q}_j, \ \ j = 1, 2, \dots, p, \\ &I^{2-{\alpha}}_{0, t}u(t)|_{t = 0} = u_1, \ \ \ \ \ \mu_1 u(t)|_{t = T}+\nu_1 I^{{\alpha}-1}_{0, t}u(t)|_{t = T} = u_2, \end{split} \end{equation}$

(3.2)

where the function $f:\ [0, T]\to {\mathbb{R}}$ is a continuous function, $\mathcal{P}_j, \mathcal{Q}_j, \ j = 1, 2, \dots, p$ are real numbers, $u_1, u_2\in {\mathbb{R}}$ .

Now we will provide an integral presentation of the solution of (3.2).

Lemma 3. The solution of (3.2) satisfies the integral equation

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}f(t), \ \ \ t\in (0, t_1] \\ \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(t-t_m)^{{\alpha}-1}\prod\limits_{k = 1}^m\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{t_m, t}f(t)+(t-t_m)^{{\alpha}-1}\sum\limits_{k = 1}^m\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in (t_m, t_{m+1}], m = 1, 2, \dots, p, \end{array} \right.$

where

$\begin{equation} \label{{1302}}\begin{split} c_0& = \frac{u_2}{\mu_1}(T-t_p)^{1-{\alpha}}\mathcal{M} -\frac{u_1}{t_1\Gamma({\alpha}-1)} -(T-t_p)^{1-{\alpha}}\mathcal{M} I^{{\alpha}}_{t_p, t}f(t)|_{t = T}\\ &-\mathcal{M}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1} (T-t_p)^{1-{\alpha}}\mathcal{M}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{M}& = ({\alpha}-1)^p\Gamma^p({\alpha}-1)\prod\limits_{k = 1}^p\mathcal{P}_k(t_k-t_{k-1})^{1-{\alpha}}. \end{split} \end{equation}$

(3.3)

Proof. We will use induction.

For $t\in (0, t_1]$ similar to Lemma 2 we get

$\begin{equation} \label{{1200}} u(t) = c_0t^{{\alpha}-1}+c_1t^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_0^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(3.4)

where $c_1 = \frac{u_1}{\Gamma({\alpha}-1)}.$

For $t\in (t_1, t_2]$ by Lemma 1 with $a = t_1$ we get

$\begin{equation} \label{{12105}} u(t) = b_0(t-t_1)^{{\alpha}-1}+b_1(t-t_1)^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{t_1}^t(t-s)^{{\alpha}-1}f(s)ds \end{equation}$

(3.5)

and

$\begin{equation} \label{{12195}} u'(t) = b_0({\alpha}-1)(t-t_1)^{{\alpha}-2}+b_1({\alpha}-2)(t-t_1)^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{t_1}^t(t-s)^{{\alpha}-2}f(s)ds. \end{equation}$

(3.6)

From the impulsive condition $I^{2-{\alpha}}_{t_1, t}u'(t)|_{t = t_1} = \mathcal{P}_1u(t_1)+\mathcal{Q}_1$ , equalities (2.2) and $I^{{\alpha}}_{t_1, t}f(t))|_{t = t_1} = 0$ we obtain

$\begin{equation} \label{{1309}} \begin{split}& I^{2-{\alpha}}_{t_1, t}u'(t) = b_0({\alpha}-1)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-2}+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3} + I^{2-{\alpha}}_{t_1, t} I^{{\alpha}-1}_{t_1, t}f(t)\\ & = b_0({\alpha}-1)\Gamma({\alpha} -1)+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3} + I^{1}_{t_1, t}f(t)\end{split} \end{equation}$

(3.7)

and

$\begin{equation} \label{{6309}} \begin{split}& I^{2-{\alpha}}_{t_1, t}u'(t)|_{t = t_1} = b_0({\alpha}-1)\Gamma({\alpha} -1)+b_1({\alpha}-2)I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3}|_{t = t_1} = \mathcal{P}_1u(t_1)+\mathcal{Q}_1 < \infty.\end{split} \end{equation}$

(3.8)

Since the integral $I^{2-{\alpha}}_{t_1, t}(t-t_1)^{{\alpha}-3}$ does not converge, it follows that $b_1 = 0$ and

$b_0 = \frac{\mathcal{P}_1}{({\alpha}-1)\Gamma({\alpha}-1)}(c_0(t_1-t_0)^{{\alpha}-1}+\frac{u_1}{t_1\Gamma({\alpha}-1)}(t_1-t_0)^{{\alpha}-1}+I^{{\alpha}}_{t_0, t}f(t)|_{t = t_1})+\frac{\mathcal{Q}_1}{({\alpha}-1)\Gamma({\alpha}-1)}.$

Thus,

$\begin{equation} \label{{12196}} \begin{split}u(t)& = \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)\frac{\mathcal{P}_1(t_1-t_0)^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}(t-t_1)^{{\alpha}-1}\\ &\ \ \ +\frac{\mathcal{P}_1}{({\alpha}-1)\Gamma({\alpha}-1)}I^{{\alpha}}_{t_0, t}f(t)|_{t = t_1}(t-t_1)^{{\alpha}-1}\\ &\ \ \ +\frac{\mathcal{Q}_1}{({\alpha}-1)\Gamma({\alpha}-1)}(t-t_1)^{{\alpha}-1} +I^{{\alpha}}_{t_1, t}f(t).\end{split} \end{equation}$

(3.9)

Similarly, for $t\in (t_j, t_{j+1}], \ j = 1, 2, \dots, p,$ by Lemma 1 with $a = t_j$ we get

$\begin{equation} \label{{12106}} u(t) = k_0(t-t_j)^{{\alpha}-1}+k_1(t-t_j)^{{\alpha}-2}+\frac{1}{\Gamma({\alpha})}\int_{t_j}^t(t-s)^{{\alpha}-1}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds \end{equation}$

(3.10)

and

$\begin{equation} \label{{12198}} u'(t) = k_0({\alpha}-1)(t-t_j)^{{\alpha}-2}+k_1({\alpha}-2)(t-t_j)^{{\alpha}-3}+\frac{1}{\Gamma({\alpha}-1)}\int_{t_j}^t(t-s)^{{\alpha}-2}\phi_1(s, I^{{\alpha}}_{0, s}u(s), I^{{\beta}}_{0, s} v(s))ds. \end{equation}$

(3.11)

From the impulsive conditions we obtain $k_1 = 0$ and

$k_0 = \frac{\mathcal{P}_j}{({\alpha}-1)\Gamma({\alpha}-1)}u(t_j)+\frac{\mathcal{Q}_j}{({\alpha}-1)\Gamma({\alpha}-1)}.$

From the boundary condition $\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = u_2$ we get

$\mathcal{I}^{{\alpha}-1} u(t)|_{t = T} = \frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds$

and

$\begin{equation} \label{{6888}}\begin{split} &\mu_1 u(t)|_{t = T}+\nu_1 \mathcal{I}^{{\alpha}-1} u(t)|_{t = T}\\ & = \mu_1 \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(T-t_p)^{{\alpha}-1}\prod\limits_{k = 1}^p\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ +\mu_1 I^{{\alpha}}_{t_p, t}f(t)|_{t = T}+\mu_1(T-t_p)^{{\alpha}-1}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}f(t)|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ +\nu_1\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds\\ & = u_2 \end{split} \end{equation}$

(3.12)

and we obtain (3.3).

Example 2. Consider the equation

$\begin{equation} \label{{99988}}\begin{split}&^{RL}D^{1.5}_{0, t}u(t) = t, \ \ \ t\in (0, 1], \ \ \ \ \ \ \ ^{RL}D^{1.5}_{1, t}u(t) = t, \ t\in(1, 2], \\ &I^{0.5}_{1, t}u'(t)|_{t = 1} = 1, \\ &I^{0.5}_{0, t}u(t)|_{t = 0} = 0, \ \ \ \ \ u(t)|_{t = 2}+ I^{0.5}_{0, t}u(t)|_{t = 2} = 1.\end{split} \end{equation}$

(3.13)

Note the Eq (3.13) is similar to (2.19) but the lower limit of the fractional derivative is changed at the point of the impulse. The solution of (3.13) satisfies the integral equation

$u(t) = \left\{ \begin{array}{ll} \frac{0.266667 t^{2.5}}{\Gamma(1.5)}, & t\in (0, 1], \\ \frac{1}{\Gamma(1.5)}(-0.4 (t-1)^{2.5} + \frac{2 (t-1)^{1.5} t}{3})+\frac{(t-1)^{0.5}}{(0.5)\Gamma(0.5)}, & t\in (1, 2]. \end{array} \right.$

It is clear the change of the lower limits of the fractional derivatives has a huge influence on the solution of the equation.

3.3. Integral presentation of the problem (3.1)

Based on the integral presentation of the linear problem (3.2) and Lemma 3, we obtain the following result:

Theorem 2. The solution of (3.1) satisfies the integral equations

$u(t) = \left\{ \begin{array}{ll} c_0t^{{\alpha}-1}+\frac{u_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)), \ \ \ t\in (0, t_1] \\ \Big(c_0+\frac{u_1}{t_1\Gamma({\alpha}-1)}\Big)(t-t_m)^{{\alpha}-1}\prod _{k = 1}^m\frac{\mathcal{P}_k(t_k-t_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{t_m, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))\\ \ \ \ \ +(t-t_m)^{{\alpha}-1}\sum _{k = 1}^m\frac{\mathcal{P}_{k}\ I^{{\alpha}}_{t_{k-1}, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod _{j = k+1}^{m}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in (t_m, t_{m+1}], m = 1, 2, \dots, p, \end{array} \right.$

and

$v(t) = \left\{ \begin{array}{ll} b_0t^{{\alpha}-1}+\frac{v_1}{\Gamma({\alpha}-1)}t^{{\alpha}-2}+I^{{\alpha}}_{0, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t)), \ \ \ t\in (0, {\tau}_1] \\ \Big(b_0+\frac{v_1}{{\tau}_1\Gamma({\alpha}-1)}\Big)(t-{\tau}_m)^{{\alpha}-1}\prod _{k = 1}^m\frac{\mathcal{P}_k^*({\tau}_k-{\tau}_{k-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ \ \ \ \ +I^{{\alpha}}_{{\tau}_m, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))\\ \ \ \ \ +(t-{\tau}_m)^{{\alpha}-1}\sum _{k = 1}^m\frac{\mathcal{P}_{k}^*\ I^{{\alpha}}_{{\tau}_{k-1}, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = {\tau}_{k}}+\mathcal{Q}_{k}^*}{({\alpha}-1)\Gamma({\alpha}-1)}\prod _{j = k+1}^{m}\frac{\mathcal{P}_j^*({\tau}_j-{\tau}_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}, \\ \ \ \ \ \ \mathit{\mbox{for}}\ \ \ \ t\in ({\tau}_m, {\tau}_{m+1}], m = 1, 2, \dots, q, \end{array} \right.$

where

$\begin{equation} \label{{53029}} \begin{split} c_0& = \frac{u_2}{\mu_1}(T-t_p)^{1-{\alpha}}\mathcal{M} -\frac{u_1}{t_1\Gamma({\alpha}-1)} \\ &\ \ \ -(T-t_p)^{1-{\alpha}}\mathcal{M} I^{{\alpha}}_{t_q, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = T}\\ &\ \ \ -\mathcal{M}\sum\limits_{k = 1}^p\frac{\mathcal{P}_{k}I^{{\alpha}}_{t_{k-1}, t}\phi_1(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = t_{k}}+\mathcal{Q}_{k}}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{p}\frac{\mathcal{P}_j(t_j-t_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_1}{\mu_1} (T-t_p)^{1-{\alpha}}\mathcal{M}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{M}& = ({\alpha}-1)^p\Gamma^p({\alpha}-1)\prod\limits_{k = 1}^p\mathcal{P}_k(t_k-t_{k-1})^{1-{\alpha}}, \nonumber\end{split} \end{equation}$

$\begin{equation} \label{{53022}} \begin{split}b_0& = \frac{v_2}{\mu_2}(T-{\tau}_q)^{1-{\alpha}}\mathcal{C} -\frac{v_1}{{\tau}_1\Gamma({\alpha}-1)} \\ &\ \ \ -(T-{\tau}_q)^{1-{\alpha}}\mathcal{C} \ I^{{\alpha}}_{{\tau}_q, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))_{t = T}\\ &\ \ \ -\mathcal{C}\sum\limits_{k = 1}^q\frac{\mathcal{P}_{k}^*\ I^{{\alpha}}_{{\tau}_{k-1}, t}\phi_2(t, I^{{\alpha}}_{0, t} u(t), I^{{\beta}}_{0, t} v(t))|_{t = {\tau}_{k}}+\mathcal{Q}_{k}^*}{({\alpha}-1)\Gamma({\alpha}-1)}\prod\limits_{j = k+1}^{q}\frac{\mathcal{P}_j^*({\tau}_j-{\tau}_{j-1})^{{\alpha}-1}}{({\alpha}-1)\Gamma({\alpha}-1)}\\ &\ \ \ \ \ \ -\frac{\nu_2}{\mu_2} (T-t_q)^{1-{\alpha}}\mathcal{C}\frac{1}{\Gamma({\alpha}-1)}\int _0^T(T-s)^{{\alpha}-2}u(s)ds, \\ \mathcal{C}& = ({\alpha}-1)^q\Gamma^q({\alpha}-1)\prod\limits_{k = 1}^q\mathcal{P}_k^*({\tau}_k-{\tau}_{k-1})^{1-{\alpha}}.\nonumber\end{split} \end{equation}$

The proof is similar to the one of Lemma 3 applied twice to each of the both components $u$ and $v$ of the coupled system (3.1) for impulsive points $t_i, \ i = 1, 2, \dots, p$ and ${\tau}_i, \ i = 1, 2, \dots, p$ and the functions $f(t) = \phi_1(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ and $f(t) = \phi_2(t, I^{{\alpha}}_{0, t}u(t), I^{{\beta}}_{0, t} v(t))$ respectively.

4. Conclusions

In this paper we set up and study a scalar nonlinear integro-differential equation with Riemann-Liouville fractional derivative and impulses. We consider a boundary value problem for the studied equation with Riemann-Liouville fractional derivative of order in $(1, 2)$ . Note for Riemann-Liouville fractional differential equations both the initial condition and impulsive conditions have to be appropriately given (which is different in the case of ordinary derivatives as well as the case of Caputo fractional derivatives). We consider both interpretations in the literature on the presence of impulses in fractional differential equations: With fixed lower limit of the fractional derivative at the initial time point and with lower limits changeable at each impulsive time point. In both cases we set up in an appropriate way impulsive conditions which are dependent on the Riemann-Liouville fractional derivative. We obtain integral presentations of the solutions in both cases. These presentations could be successfully used for furure studies of existence, stability and other qualitative properties of the solutions of the integro-differential equations with Riemann-Liouville fractional derivative of order in $(1, 2)$ and impulses.

Acknowledgments

S. H. is partially supported by the Bulgarian National Science Fund under Project KP-06-N32/7.

Conflict of interest

The authors declare that they have no competing interests.

References

[1]	A. A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier Science Limited, 204 (2006).
[2]	F. Jarad, T. Abdeljawad, Generalized fractional derivatives and Laplace transform, Discrete Contin. Dyn. Syst.-S, 13 (2020), 709-722. doi: 10.3934/dcdss.2020039
[3]	F. Jarad, T. Abdeljawad, A modified Laplace transform for certain generalized fractional operators, Results Nonlinear Anal., 1 (2018), 88-98.
[4]	E. Bonyah, A. Atangana, M. Chand, Analysis of 3D IS-LM macroeconomic system model within the scope of fractional calculus, Chaos, Solitons Fractals: X, 2 (2019), 100007.
[5]	M. Yavuz, N. Ozdemir, European vanilla option pricing model of fractional order without singular kernel, Fractal Fractional, 2 (2018), 1-11.
[6]	A. Atangana, S. I. Araz, Fractional stochastic modelling illustration with modified Chua attractor, Eur. Phys. J. Plus, 134 (2019), 1-23. doi: 10.1140/epjp/i2019-12286-x
[7]	S. Qureshi, A. Yusuf, Modeling chickenpox disease with fractional derivatives: From caputo to atangana-baleanu Chaos, Solitons Fractals, 122 (2019), 111-118.
[8]	S. Qureshi, A. Yusuf, Mathematical modeling for the impacts of deforestation on wildlife species using Caputo differential operator, Chaos, Solitons Fractals, 126 (2019), 32-40. doi: 10.1016/j.chaos.2019.05.037
[9]	M. Yavuz, E. Bonyah, New approaches to the fractional dynamics of schistosomiasis disease model, Phys. A: Stat. Mech. Appl., 525 (2019), 373-393. doi: 10.1016/j.physa.2019.03.069
[10]	F. A. Rihan, Q. M. Al-Mdallal, H. J. AlSakaji, et al. A fractional-order epidemic model with time-delay and nonlinear incidence rate, Solitons Fractals, 126 (2019), 97-105. doi: 10.1016/j.chaos.2019.05.039
[11]	E. Bas, R. Ozarslan, Real world applications of fractional models by Atangana-Baleanu fractional derivative, Chaos, Solitons Fractals, 116 (2018), 121-125. doi: 10.1016/j.chaos.2018.09.019
[12]	E. Bas, B. Acay, R. Ozarslan, Fractional models with singular and non-singular kernels for energy efficient buildings, Chaos: Interdiscip. J. Nonlinear Sci., 29 (2019), 023110.
[13]	A. Atangana, E. Bonyah, Fractional stochastic modeling: New approach to capture more heterogeneity, Chaos: Interdiscip. J. Nonlinear Sci., 29 (2019), 013118.
[14]	T. Abdeljawad, Fractional operators with boundary points dependent kernels and integration by parts, Discrete Contin. Dyn. Syst.-S, 13 (2019), 1098-1107.
[15]	Q. M. Al-Mdallal, On fractional-Legendre spectral Galerkin method for fractional Sturm-Liouville problems, Chaos, Solitons Fractals, 116 (2018), 261-267. doi: 10.1016/j.chaos.2018.09.032
[16]	B. Acay, E. Bas, T. Abdeljawad, Non-local fractional calculus from different viewpoint generated by truncated M-derivative, J. Comput. Appl. Math., 366 (2019), 112410.
[17]	E. Bas, B. Acay, R. Ozarslan, The price adjustment equation with different types of conformable derivatives in market equilibrium, AIMS Math., 4 (2019), 805-820. doi: 10.3934/math.2019.3.805
[18]	N. Sene, Stability analysis of the generalized fractional differential equations with and without exogenous inputs, J. Nonlinear Sci. Appl., 12 (2019), 562-572. doi: 10.22436/jnsa.012.09.01
[19]	M. Al-Refai, M. A. Hajji, Analysis of a fractional eigenvalue problem involving Atangana-Baleanu fractional derivative: A maximum principle and applications, Chaos: Interdiscip. J. Nonlinear Sci., 29 (2019), 013135.
[20]	E. Bas, R. Yilmazer, E. Panakhov, Fractional Solutions of Bessel Equation with N-Method, Sci. World J., 2013 (2013), 1-9.
[21]	M. G. Sakar, O. Saldır, Improving variational iteration method with auxiliary parameter for nonlinear time-fractional partial differential equations, J. Optim. Theory Appl., 174 (2017), 530-549. doi: 10.1007/s10957-017-1127-y
[22]	M. G. Sakar, Numerical solution of neutral functional-differential equations with proportional delays, Int. J. Optim. Control: Theory Appl. (IJOCTA), 7 (2017), 186-194. doi: 10.11121/ijocta.01.2017.00360
[23]	E. Bas, R. Ozarslan, D. Baleanu, et al. Comparative simulations for solutions of fractional Sturm-Liouville problems with non-singular operators, Adv. Differ. Equations, 2018 (2018), 1-19. doi: 10.1186/s13662-017-1452-3
[24]	S. A. A. Shah, M. A. Khan, M. Farooq, et al., A fractional order model for Hepatitis B virus with treatment via Atangana-Baleanu derivative, Phys. A: Stat. Mech. Appl., 538 (2020), 122636.
[25]	E. O. Alzahrani, M. A. Khan, Comparison of numerical techniques for the solution of a fractional epidemic model, Eur. Phys. J. Plus, 135 (2020), 1-28. doi: 10.1140/epjp/s13360-019-00059-2
[26]	M. A. Khan, O. Kolebaje, A. Yildirim, et al. Fractional investigations of zoonotic visceral leishmaniasis disease with singular and non-singular kernel, Eur. Phys. J. Plus, 134 (2019), 1-29. doi: 10.1140/epjp/i2019-12286-x
[27]	M. A. Khan, A. Khan, A. Elsonbaty, et al. Modeling and simulation results of a fractional dengue model, Eur. Phy. J. Plus, 134 (2019), 1-15. doi: 10.1140/epjp/i2019-12286-x
[28]	S. Ullah, M. A. Khan, M. Farooq, et al. A fractional model for the dynamics of tuberculosis (TB) using Atangana-Baleanu derivative, Discrete Contin. Dyn. Syst.-S, 13 (2019), 937-956.
[29]	M. A. Khan, S. Ullah, M. Farhan, The dynamics of Zika virus with Caputo fractional derivative, AIMS Math., 4 (2019), 134-146. doi: 10.3934/Math.2019.1.134
[30]	J. R. Garcia, M. G. Calderon, J. M. Ortiz, et al. Motion of a particle in a resisting medium using fractional calculus approach, Proc. Romanian Acad. A, 14 (2013), 42-47.
[31]	S. Salahshour, A. Ahmadian, F. Ismail, et al. A fractional derivative with non-singular kernel for interval-valued functions under uncertainty, Optik, 130 (2017), 273-286. doi: 10.1016/j.ijleo.2016.10.044
[32]	J. Losada, J. J. Nieto, Properties of a new fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 87-92.
[33]	T. Abdeljawad, Fractional operators with generalized Mittag-Leffler kernels and their iterated differintegrals, Chaos: Interdiscip. J. Nonlinear Sci., 29 (2019), 023102.
[34]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1 (2015), 1-13.
[35]	A. Atangana, D. Baleanu, New fractional derivatives with non-local and nonsingular kernel: Theory and application to heat transfer model, Therm. Sci., 20 (2016), 763-769. doi: 10.2298/TSCI160111018A
[36]	U. N. Katugampola, New approach to a generalized fractional integral, Appl. Math. Comput., 218 (2011), 860-865.
[37]	F. Jarad, T. Abdeljawad, D. Baleanu, On the generalized fractional derivatives and their Caputo modification, J. Nonlinear Sci. Appl., 10 (2017), 2607-2619. doi: 10.22436/jnsa.010.05.27
[38]	S. Liang, R. Wu, L. Chen, Laplace transform of fractional order differential equations. Electron. J. Differ. Equations, 139 (2015), 1-15.
[39]	A. Atangana, D. Baleanu, Caputo-Fabrizio derivative applied to groundwater flow within confined aquifer, J. Eng. Mech., 143 (2017), 1-5.
[40]	T. Abdeljawad, Fractional difference operators with discrete generalized Mittag-Leffler kernels, Chaos, Solitons Fractals, 126 (2019), 315-324. doi: 10.1016/j.chaos.2019.06.012

This article has been cited by:

1.	Vitor Miguel Ribeiro, Pioneering paradigms: unraveling niche opportunities in green finance through bibliometric analysis of nation brands and brand culture, 2024, 6, 2643-1092, 287, 10.3934/GF.2024012
2.	Manuel Mojica, Pedro R. Palos-Sanchez, Edgar Cabanas, Is there innovation management of emotions or just the commodification of happiness? A sentiment analysis of happiness apps, 2024, 1460-1060, 10.1108/EJIM-11-2023-0963
3.	Pedro R. Palos-Sanchez, Raúl J. Chang-Tam, José A. Folgado-Fernández, The Role of Neobanks and FinTech in Sustainable Finance and Technology. The Customer/User Perspective for Entrepreneurs, 2025, 27730328, 100109, 10.1016/j.stae.2025.100109

Reader Comments

Your name:*

Email:*
© 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)