Existence of mild solutions for fractional non-instantaneous impulsive integro-differential equations with nonlocal conditions

Arshi Meraj, Dwijendra N. Pandey

Arab Journal of Mathematical Sciences

ISSN: 1319-5166

Open Access. Article publication date: 27 November 2018

Issue publication date: 31 August 2020

648

Abstract

This paper is concerned with the existence of mild solutions for a class of fractional semilinear integro-differential equations having non-instantaneous impulses. The result is obtained by using noncompact semigroup theory and fixed point theorem. The obtained result is illustrated by an example at the end.

Keywords

Citation

Meraj, A. and Pandey, D.N. (2020), "Existence of mild solutions for fractional non-instantaneous impulsive integro-differential equations with nonlocal conditions", Arab Journal of Mathematical Sciences, Vol. 26 No. 1/2, pp. 3-13. https://doi.org/10.1016/j.ajmsc.2018.11.002

Publisher

:

Emerald Publishing Limited

Copyright © 2019, Arshi Meraj and Dwijendra N. Pandey

License

Published in the Arab Journal of Mathematical Sciences. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) license. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this license may be seen at http://creativecommons.org/licences/by/4.0/legalcode


1. Introduction

The objective of this paper is to study the existence of mild solutions to the following abstract integro-differential equations of fractional order with non-instantaneous impulses and nonlocal conditions in a Banach space X:

msupDqu(t)+Au(t)=f(t,u(t),0tK(t,s)u(s)ds),tk=0m(sk,tk+1],
u(t)=γk(t,u(t)),tk=1m(tk,sk],
(1.1)u(0)+g(u)=u0,
where msupDq is the Caputo fractional derivative of order q(0<q<1),A:D(A)XX is closed linear operator, A is the infinitesimal generator of an equicontinuous and uniformly bounded C0 semigroup T(t)(t0) on X,J=[0,a],a>0 is a constant, 0<t1<t2<<tm<tm+1:=a,s0:=0 and sk(tk,tk+1) for each k=1,2,,m,f:J×X×XX,g:PC(J,X)X are given functions satisfying certain assumptions, γk:(tk,sk]×XX are non-instantaneous impulsive functions for all k=1,2,,m and KC(D,+) where D:={(t,s):0s<ta} and u0X.

In the past decades, many researchers paid attention to study the differential equations with instantaneous impulses, which have been used to describe abrupt changes such as shocks, harvesting and natural disasters. Particularly, the theory of instantaneous impulsive equations have wide applications in control, mechanics, electrical engineering, biological and medical fields. For more details on the differential equations with instantaneous impulses one may see [,,,,].

It seems that models with instantaneous impulses could not explain the certain dynamics of evolution process in pharmacotherapy. For example, one considers the hemodynamic equilibrium of a person, the introduction of the drugs in bloodstream and the consequent absorption for the body are gradual and continuous process. Hernández and O’Regan [] and Pierri et al. [], initially studied Cauchy problems for first order evolution equations with non-instantaneous impulses. The recent results for evolution equations with non-instantaneous impulses can be found in [,,,] and the references therein.

The nonlocal problem was motivated by physical problems. Indeed it is demonstrated that the nonlocal problems have better effects in applications than the classical Cauchy problems. For example it is used to represent mathematical models for evolution of various phenomena such as nonlocal neutral networks, nonlocal pharmacokinetics, nonlocal pollution and nonlocal combustion (see []). The existence results to evolution equations with nonlocal conditions in Banach space were first studied by Byszewski []. Deng [] used the nonlocal condition to describe the diffusion phenomenon of a small amount of gas in a transparent tube.

To the best of our knowledge, there is no work yet reported on fractional non-instantaneous impulsive integro-differential equations with nonlocal conditions when the corresponding semigroup T(t)(t0) is noncompact. Therefore inspired by the previous works, we will study the existence of PC-mild solutions for under the assumption that the corresponding C0 semigroup is noncompact, by using the properties of Kuratowski measure of noncompactness, and ρ-set contraction mapping fixed point theorem (see ). We conclude this section by summarizing the contents of this paper. In the next section, we will introduce some basic definitions, notations and preliminary lemmas. In Section , we will prove existence of mild solutions for the problem also we will give an example to illustrate the feasibility of our abstract result.

2. Preliminaries

Let X be a Banach space with norm ·, we use θ to denote the zero function in PC(J,X) and J=[0,a] for any constant a>0. Let C(J,X) be a Banach space of all continuous functions from J into X endowed with supremum norm uC=suptJu(t). Consider the space PC(J,X)={u:JX:uis continuous at ttk,u(tk)=u(tk)and u(tk+)exists for all k=1,2,.m}, which is a Banach space endowed with supremum norm uPC=suptJu(t). For each finite constant r>0, let Ωr={uPC(J,X):u(t)r,tJ}. Let Lp(J,X)(1p<) be the Banach space of all X-valued Bochner integrable functions defined on J with norm uLp(J,X)=(0au(t)pdt)1p. Denote Gu(t):=0tK(t,s)u(s)ds, and let G=suptJ0tK(t,s)ds<. Let M=suptJT(t)L(X), where L(X) stands for the Banach space of all linear and bounded operators on X, note that M1. A C0-semigroup T(t)(t0) is called equicontinuous if the operator T(t) is continuous by the operator norm for every t>0.

Lemma 2.1

([]). If h satisfies a uniform Hölder continuity with exponent β(0,1], then the unique solution of the following linear Cauchy problem:

(2.1)msupDqu(t)+Au(t)=h(t),tJ,
u(0)=x0X,
is given by
(2.2)u(t)=U(t)x0+0t(t-s)q-1V(ts)h(s)ds,
where
(2.3)U(t)=0ζq(θ)T(tqθ)dθ,V(t)=q0θζq(θ)T(tqθ)dθ,
(2.4)ζq(θ)=1qθ11qρq(θ1q),ρq(θ)=1πn=0(1)n1θqn1Γ(nq+1)n!sin(nπq),θ(0,),
ζq(θ) is a probability density function defined on (0,).
Remark 2.2.

ζq(θ)0,θ(0,),0ζq(θ)dθ=1,0θζq(θ)dθ=1Γ(1+q).

Lemma 2.3

([]). The operators U(t)(t0) and V(t)(t0) have the following properties:

  1. For any fixed t0, U(t) and V(t) are strongly continuous.

  2. For any fixed t0, U(t) and V(t) are linear bounded operators, moreover for any uX,

    U(t)uMu,V(t)uMΓ(q)u.

  3. If T(t)(t0) is an equicontinuous semigroup, then U(t) and V(t) are continuous for t>0 by the operator norm, which means that for 0<t<ta, we have

    U(t)U(t)0andV(t)V(t)0 as tt.

Definition 2.4

([]). A function uPC(J,X) is said to be a mild solution of the problem if u(0)=u0g(u), u(t)=γk(t,u(t)) for all tk=1m(tk,sk], and

u(t)={U(t)(u0-g(u)) +0t(t - s)q-1V(t-s) f(s, u(s), Gu(s))ds,t(0, t1];U(t-sk)γk(sk, u(sk))+skt(t - s)q-1V(t-s) f(s, u(s), Gu(s))ds,t(sk, tk+1], k = 1, 2, . . . ,m.

Now, we recall some properties of measure of noncompactness which are useful to prove our main result. For the details about measure of noncompactness, one may see [,]. Let α(·) denotes the Kuratowski measure of noncompactness of the bounded set.

Lemma 2.5

([]). Let X be a Banach space, and UC(J,X), U(t)={u(t):uU}(tJ). If U is bounded and equicontinuous in C(J,X), then α(U(t)) is continuous on J, and α(U)=maxtJα(U(t)).

Lemma 2.6

([]). If X be a Banach space and D={un}n=1PC(J,X) be a bounded and countable set, then α(D(t)) is Lebesgue integrable on J, and

α({0tun(s)ds}n=1)20tα({un(s)}n=1)ds.
Lemma 2.7

([]). Let X be a Banach space and U is bounded subset of X, then there exists a countable set D={un}n=1U such that α(U)2α(D).

Lemma 2.8

([]). Let X and E be Banach spaces and Q:D(Q)EX is Lipschitz continuous with constant L, then α(Q(V))Lα(V) for any bounded subset VD(Q).

Definition 2.9

([]). Let X be a Banach space, and S be a nonempty subset of X. A continuous map Q:SX is called ρ-set contractive if there exists a constant ρ[0,1) such that for every bounded set ΩS,

α(Q(Ω))ρα(Ω).
Lemma 2.10

([]). Let X be a Banach space, ΩX be a closed bounded and convex subset, and the operator Q:ΩΩ is ρ-set contractive, then Q has at least one fixed point in Ω.

3. Main result and example

In this section, we will discuss the existence of mild solutions for the system , then we will present an example to illustrate our proved result. Let us introduce the required assumptions which are needed to prove our main result:

(H1)

For each tJ, the function f(t,,):X×XX is continuous and for all (x,y)X×X, the function f(,x,y):JX is Lebesgue measurable.

(H2)

There exist a continuous nondecreasing function ψ:[0,)(0,), a constant q1(0,q), and a function φL1q1(J,+) such that

f(t,x,y)φ(t)ψ(x),x,yX;tJ.
(H3)

g:PC(J,X)X is continuous and there exists a constant α>0 such that

g(x)g(y)αxy,x,yPC(J,X).
(H4)

γk:[tk,sk]×XX are continuous and there exist constants Kγk>0,k=1,2,,m such that

γk(t,x)γk(t,y)Kγkxy,x,yX;t[tk,sk].
(H5)

There exist positive constants Lk and Nk,k=0,1,2,,m such that for any countable sets D1,D2X,

α(f(t,D1,D2))Lkα(D1)+Nkα(D2),t(sk,tk+1],k=0,1,2,,m.

Let us denote:

K=maxk=1,2,,mKγk,K=max{K,α},
(3.1)L=maxk=0,1,2,,m(Lk+NkG)(tk+1sk)q.
Theorem 3.1.

Assume that the semigroup T(t)(t0) generated by A is equicontinuous, the functions g(θ) and γk(·,θ) are bounded for k=1,2,,m, and the assumptions (H1) (H5) are satisfied, then the system has at least one PC- mild solution provided that

(3.2)max{Λ1,Λ2}<1,
where Λ1=M(α+K) and Λ2=M(K+4LΓ(q+1)).
Proof.

Define the operator F:PC(J,X)PC(J,X) as

(3.3)(Fu)(t)=(F1u)(t)+(F2u)(t),
where
(3.4)(F1u)(t)={U(t)(u0g(u)),t[0,t1];γk(t,u(t)),t(tk,sk],k=1,2,,m,U(tsk)γk(sk,u(sk)),t(sk,tk+1],k=1,2,,m.
(3.5)(F2u)(t)={sktt(ts)q1V(ts)f(s,u(s),Gu(s))ds,t(sk,tk+1],k=0,1,2,,m,0,otherwise

It is easy to see that F is well defined. From , one can easily see that the PC-mild solution of the system is equivalent to a fixed point of the operator F defined by . Now, we will prove that the operator F has a fixed point.

Let uΩR for some R>0, q2=q11q1(1,0) and M1=ψ(R)φL1q1(J,R+), by using Hölder inequality and (H2), we obtain

(3.6)0t(ts)q1f(s,u(s),Gu(s))ds(0t(ts)q2ds)1q1ψ(R)φL1q1(J,R+)M1(1+q2)1q1a(1+q2)(1q1).

Now, we divide the proof into the following steps:

Step I:

We prove that there exists a constant R>0 such that F(ΩR)ΩR.

If this is not true, then for each r>0, there will exist urΩr and trJ such that (Fur(tr))>r. If tr[0,t1], then by , , and (H3) we have

(3.7)(Fur)(tr)M(u0+αurθ+g(θ)+MM1Γ(q)(1+q2)1q1a(1+q2)(1q1)M(αr+u0+g(θ)+MM1Γ(q)(1+q2)1q1a(1+q2)(1q1).

If tr(tk,sk],k=1,2,,m, then by and (H4), we obtain

(3.8)(Fur)(tr)=(γk(tr,ur(tr)))Kγkur(tr)+γk(tr,θ)Kγkr+β,
where β=maxk=1,2,,m{suptJγk(t,θ)}. If tr(sk,tk+1],k=1,2,,m, then by , , and (H4) we have
(3.9)(Fur)(tr)M(Kγkr+β)+Msktr(trs)q1f(s,ur(s),Gur(s))dsM(Kγkr+β)+MM1Γ(q)(1+q2)1q1a(1+q2)(1q1)

Combining with the fact r<(Fur)(tr), we obtain

(3.10)r<(Fur)(tr)M(αr+u0+g(θ)+M(Kr+β)+MM1Γ(q)(1+q2)1q1a(1+q2)(1q1).

Dividing both sides of by r and taking limit as r, we have

(3.11)1M(α+K),
which contradicts .
Step II:

We prove that the operator F1:ΩRΩR is Lipschitz continuous.

For t[0,t1] and u,vΩR, using and (H3) we have

(3.12)(F1u)(t)(F1v)(t)Mg(u)g(v)Mαuv.

For t(tk,sk],k=1,2,,m and u,vΩR, by and the assumption (H4), we obtain

(3.13)(F1u)(t)(F1v)(t)Kγku(t)v(t)MKuv.

For t(sk,tk+1],k=1,2,,m and u,vΩR, using (H4), we have

(3.14)(F1u)(t)(F1v)(t)Mγk(sk,u(sk))γk(sk,v(sk))MKuv.

From , we obtain

(3.15)F1uF1vMKuv,
where K:=m{K,α}.
Step III:

In this step, we prove that F2 is continuous on ΩR.

Let {un} be a sequence in ΩR such that limnun=u in ΩR. By the continuity of nonlinear term f with respect to second and third variables, for each sJ, we have

(3.16)limnf(s,un(s),Gun(s))=f(s,u(s),Gu(s)).

So, we can conclude that

(3.17)suptJf(s,un(s),Gun(s))f(s,u(s),Gu(s))0asn.

For s[sk,t] and t(sk,tk+1],k=0,1,2,,m, by and , we obtain

(3.18)(F2un)(t)(F2u)(t)MΓ(q)skt(ts)q1f(s,un(s),Gun(s))f(s,u(s),Gu(s))dsMaqΓ(q+1)suptJf(s,un(s),Gun(s))f(s,u(s),Gu(s))
(3.19)0asn.
Hence,
(3.20)F2unF2u0asn,
which means that F2 is continuous on ΩR.
Step IV:

Now, we show F2:ΩRΩR is equicontinuous.

For any uΩR and skt<ttk+1 for k=0,1,2,,m, we have

(F2u)(t)(F2u)(t)=skt(ts)q1V(ts)f(s,u(s),Gu(s))dsskt(ts)q1V(ts)f(s,u(s),Gu(s))dstt(ts)q1V(ts)f(s,u(s),Gu(s))ds+skt[(ts)q1(ts)q1]V(ts)f(s,u(s),Gu(s))ds+skt(ts)q1[V(ts)V(ts)]f(s,u(s),Gu(s))ds=I1+I2+I3,
where,
I1=tt(ts)q1V(ts)f(s,u(s),Gu(s))ds,
I2=skt[(ts)q1(ts)q1]V(ts)f(s,u(s),Gu(s))ds,
I3=skt(ts)q1[V(ts)V(ts)]f(s,u(s),Gu(s))ds.

Now, we only need to check that I1,I2 and I3 tend to 0 independently of uΩR when tt. By , we have

I1M1MΓ(q)(1+q2)1q1(tt)(1+q2)(1q1)0astt.

For I2, by (H2), , Hölder inequality, and [], we get that

I2MΓ(q)(skt[(ts)q1(ts)q1]11q1ds)1q1ψ(R)φL1q1(J,)M1MΓ(q)(skt[(ts)q2(ts)q2]ds)1q1M1MΓ(q)(1+q2)1q1[(t)1+q2(t)1+q2+(tt)1+q2]1q1M1MΓ(q)(1+q2)1q1(tt)(1+q2)(1q1)0astt.

For t=sk, it is easy to see that I3=0. For t>sk and ϵ>0 small enough, by (H2), , and the equicontinuity of T(t) , we estimate

I3sktϵ(ts)q1[V(ts)V(ts)]f(s,u(s),Gu(s))ds+tϵt(ts)q1[V(ts)V(ts)]f(s,u(s),Gu(s))dssktϵ(ts)q1f(s,u(s),Gu(s))dssups[sk,tϵ]V(ts)V(ts)+2MΓ(q)tϵt(ts)q1f(s,u(s),Gu(s))dsM1(1+q2)1q1((t)1+q2ϵ1+q2)1q1sups[sk,tϵ]V(ts)V(ts)+2M1MΓ(q)(1+q2)1q1ϵ(1+q2)(1q1)0astt.

As a result, (F2u)(t)(F2u)(t)0 independently of uΩR as tt, which means that F2:ΩRΩR is equicontinuous.

Step V:

We show that F:ΩRΩR is a ρ-set contractive map.

For any bounded set DΩR, by , we know that there exists a countable set D0={un}D such that

(3.21)α(F2(D))2α(F2(D0)).

Since F2(D0)F2(ΩR) is bounded and equicontinuous, by , we get

(3.22)α(F2(D0))=maxt[sk,tk+1],k=0,1,2,,mα(F2(D0)(t)).

For every t[sk,tk+1],k=0,1,2,,m, by , the assumption (H5) and , we have

(3.23)α(F2(D0)(t))=α({skt(ts)q1V(ts)f(s,un(s),Gun(s))ds})2MΓ(q)skt(ts)q1α({f(s,un(s),Gun(s))})ds2MΓ(q)skt(ts)q1[Lkα(D0(s))+Nkα(GD0(s))]ds.

Meanwhile, we have

(3.24)α(GD0(s))α(GD0)Gα(D0)Gα(D0)Gα(D).

Therefore,

(3.25)α(F2(D0)(t))2MΓ(q+1)(Lk+NkG)(tk+1sk)qα(D)2MLΓ(q+1)α(D).

From and , we obtain

(3.26)α(F2(D))4MLΓ(q+1)α(D).

From and , we know that for any bounded set DΩR,

(3.27)α(F1(D))MKα(D).

Therefore, by and , we obtain

(3.28)α(F(D))α(F1(D))+α(F2(D))M(K+4LΓ(q+1))α(D)=Λ2α(D).

Now combining with and , we get that F:ΩRΩR is a ρ-set-contractive map with ρ=Λ2. Hence implies that F has at least one fixed point uΩR, which is a PC-mild solution of .

Next, we present an example to illustrate our main result.

Example.

Consider the following fractional partial differential system with non-instantaneous impulses and nonlocal conditions:

(3.29){msupD12u(t,x)+2x2u(t,x)=125et1+etu(t,x)+0t150esu(s,x)ds,x(0,1),t(0,13](23,1];u(t,0)=u(t,1)=0,t[0,1];u(t,x)=e(t13)4|u(t,x)|1+|u(t,x)|,x(0,1),t(13,23];u(0,x)+i=1213iu(1i,x)=u0(x),x[0,1].

Let X=L2[0,1] and Au=u with D(A)={uX:u,uare absolutely continuous anduX,u(0)=u(1)=0}. It is well known by [], that A generates an equicontinuous C0-semigroup T(t)(t0) on X, and T(t)1, for any t0. Let a=t2=1,t0=s0=0,t1=13,s1=23. By putting

u(t)=u(t,),
f(t,u(t),Gu(t))=125et1+etu(t,)+0t150esu(s,)ds,
Gu(t)=0t150esu(s,)ds,
γ1(t,u(t))=e(t13)4|u(t,)|1+|u(t,)|,
g(u)=i=1213iu(1i,),
the parabolic partial differential equation can be rewritten into the abstract form of for m=1. It is easy to verify that the assumptions (H1)–(H5) and condition hold with
q=12,M=1,φ(t)=125et1+et+150,ψ;(r)=r,
α*=49,K=Kγ1=14,L=0.02,Λ1=0.69<1,Λ2=0.53<1.

Therefore, is applicable, so the system has at least one PC-mild solution.

References

[1]L. Bai, J.J. Nieto, Variational approach to differential equations with not instantaneous impulses, Appl. Math. Lett. 73 (2017) 4448.

[2]D.D. Bainov, V. Lakshmikantham, P.S. Simeonov, Theory of Impulsive Differential Equations, in: Series in Modern Applied Mathematics, World Scientific, Singapore, 1989.

[3]J. Banas, K. Goebel, Measure of Noncompactness in Banach Space, Marcal Dekker Inc., New York, 1980.

[4]M. Benchohra, J. Henderson, S. Ntouyas, Impulsive Differential Equations and Inclusions, in: Contemp. Math. Appl., Hindawi Publ. Corp., New York, 2006.

[5]D. Bothe, Multivalued perturbations of m-accretive differential inclusions, Israel J. Math. 108 (1998) 109138.

[6]L. Byszewski, Theorem about the existence and uniqueness of solutions of a semilinear evolution nonlocal Cauchy problem, J. Math. Appl. Anal. 162 (1991) 494505.

[7]P. Chen, Y. Li, Mixed monotone iterative technique for a class of semilinear impulsive evolution equations in Banach spaces, Nonlinear Anal. 74 (2011) 35783588.

[8]P. Chen, X. Zhang, Y. Li, Existence of mild solutions to partial differential equations with non-instantaneous impulses, Electron. J. Differential Equations 241 (2016) 111.

[9]K. Deng, Exponential decay of solutions of semilinear parabolic equations with nonlocal initial conditions, J. Math. Anal. Appl. 179 (1993) 630637.

[10]M.M. El-Borai, Some probability densities and fundamental solutions of fractional evolution equations, Chaos Solitons Fractals 14 (2002) 433440.

[11]H. Heinz, On the behaviour of measures of noncompactness with respect to differentiation and integration of vector valued functions, Nonlinear Anal. 7 (1983) 13511371.

[12]E. Hernández, D. O’Regan, On a new class of abstract impulsive differential equations, Proc. Amer. Math. Soc. 141 (2013) 16411649.

[13]P. Kumar, D.N. Pandey, D. Bahuguna, On a new class of abstract impulsive functional differential equations of fractional order, J. Nonlinear Sci. Appl. 7 (2014) 102114.

[14]J. Liang, J.H. Liu, T.J. Xiao, Nonlocal impulsive problems for integrodifferential equations, Math. Comput. Modelling 49 (2009) 789804.

[15]S. Liang, R. Mei, Existence of mild solutions for fractional impulsive neutral evolution equations with nonlocal conditions, Adv. Differential Equations (2014) http://dx.doi.org/10.1186/1687-1847-2014-101.

[16]M. McKibben, Discovering Evolution Equations with Applications, Chapman and Hall/CRC, Boca Raton, 2011.

[17]A. Pazy, Semigroup of Linear Operators and Applications to Partial Differential Equations, in: Applied Mathematical Sciences, Springer-Verlag, Berlin, 1983.

[18]M. Pierri, D. O’Regan, V. Rolnik, Existence of solutions for semi-linear abstract differential equations with non instantaneous impulses, Appl. Math. Comput. 219 (2013) 67436749.

[19]J. Wang, X. Li, Periodic BVP for integer/fractional order nonlinear differential equations with non-instantaneous impulses, J. Appl. Math. Comput. 46 (2014) 321334.

[20]D. Yang, J. Wang, Integral boundary value problems for nonlinear non-instataneous impulsive differential equations, J. Appl. Math. Comput. (2016) http://dx.doi.org/10.1007/s12190-016-1025-8.

[21]X. Yu, J. Wang, Periodic boundary value problems for nonlinear impulsive evolution equations on Banach spaces, Commun. Nonlinear Sci. Numer. Simul. 22 (2015) 980989.

[22]X. Zhang, P. Chen, Fractional evolution equation nonlocal problems with noncompact semigroups, Opuscula Math. 36 (2016) 123137.

Acknowledgements

The authors would like to express thanks to the editor and referees for their careful reading of the manuscript and valuable comments. The work of first author is supported by the “Ministry of Human Resource Development, India under Grant Number: MHR-01-23-200-428”. The publisher wishes to inform readers that the article “Existence of mild solutions for fractional non-instantaneous impulsive integro-differential equations with nonlocal conditions” was originally published by the previous publisher of the Arab Journal of Mathematical Sciences and the pagination of this article has been subsequently changed. There has been no change to the content of the article. This change was necessary for the journal to transition from the previous publisher to the new one. The publisher sincerely apologises for any inconvenience caused. To access and cite this article, please use Meraj, A., Pandey, D.N. (2018), “Existence of mild solutions for fractional non-instantaneous impulsive integro-differential equations with nonlocal conditions”, Arab Journal of Mathematical Sciences, Vol. 26 No. 1/2, pp. 3-13. The original publication date for this paper was 27/11/2018.

Corresponding author

Arshi Meraj can be contacted at: arshimeraj@gmail.com

Related articles