Research article

On Cerone's and Bellman's generalization of Steffensen's integral inequality via conformable sense

  • Received: 12 September 2022 Revised: 30 September 2022 Accepted: 10 October 2022 Published: 27 October 2022
  • MSC : 26D15, 26A33

  • By making use of the conformable integrals, we establish some new results on Cerone's and Bellman's generalization of Steffensen's integral inequality. In fact, we provide a variety of generalizations of Steffensen's integral inequality by using conformable calculus.

    Citation: Mohammed S. El-Khatib, Atta A. K. Abu Hany, Mohammed M. Matar, Manar A. Alqudah, Thabet Abdeljawad. On Cerone's and Bellman's generalization of Steffensen's integral inequality via conformable sense[J]. AIMS Mathematics, 2023, 8(1): 2062-2082. doi: 10.3934/math.2023106

    Related Papers:

    [1] Lv Zhang, Qingbiao Wu . Modified Newton-EHS method for solving nonlinear problems with complex symmetric Jacobian matrices. AIMS Mathematics, 2023, 8(10): 24233-24253. doi: 10.3934/math.20231236
    [2] Lu-Chuan Ceng, Shih-Hsin Chen, Yeong-Cheng Liou, Tzu-Chien Yin . Modified inertial subgradient extragradient algorithms for generalized equilibria systems with constraints of variational inequalities and fixed points. AIMS Mathematics, 2024, 9(6): 13819-13842. doi: 10.3934/math.2024672
    [3] Xiao Guo, Chuanpei Xu, Zhibin Zhu, Benxin Zhang . Nonmonotone variable metric Barzilai-Borwein method for composite minimization problem. AIMS Mathematics, 2024, 9(6): 16335-16353. doi: 10.3934/math.2024791
    [4] Fouad Mohammad Salama, Nur Nadiah Abd Hamid, Umair Ali, Norhashidah Hj. Mohd Ali . Fast hybrid explicit group methods for solving 2D fractional advection-diffusion equation. AIMS Mathematics, 2022, 7(9): 15854-15880. doi: 10.3934/math.2022868
    [5] Mashail M. Al Sobhi . The modified Kies-Fréchet distribution: Properties, inference and application. AIMS Mathematics, 2021, 6(5): 4691-4714. doi: 10.3934/math.2021276
    [6] Zui-Cha Deng, Fan-Li Liu, Liu Yang . Numerical simulations for initial value inversion problem in a two-dimensional degenerate parabolic equation. AIMS Mathematics, 2021, 6(4): 3080-3104. doi: 10.3934/math.2021187
    [7] Xin-Hui Shao, Wan-Chen Zhao . Relaxed modified Newton-based iteration method for generalized absolute value equations. AIMS Mathematics, 2023, 8(2): 4714-4725. doi: 10.3934/math.2023233
    [8] Krunal B. Kachhia, Jyotindra C. Prajapati . Generalized iterative method for the solution of linear and nonlinear fractional differential equations with composite fractional derivative operator. AIMS Mathematics, 2020, 5(4): 2888-2898. doi: 10.3934/math.2020186
    [9] Huiqin Zhang, Yanping Chen . Two-grid finite element method with an H2N2 interpolation for two-dimensional nonlinear fractional multi-term mixed sub-diffusion and diffusion wave equation. AIMS Mathematics, 2024, 9(1): 160-177. doi: 10.3934/math.2024010
    [10] Cui-Xia Li, Long-Quan Yong . Modified BAS iteration method for absolute value equation. AIMS Mathematics, 2022, 7(1): 606-616. doi: 10.3934/math.2022038
  • By making use of the conformable integrals, we establish some new results on Cerone's and Bellman's generalization of Steffensen's integral inequality. In fact, we provide a variety of generalizations of Steffensen's integral inequality by using conformable calculus.



    Fractional-order calculus has received considerable attention in the engineering and physical sciences over the last few decades to model a number of diverse phenomena in robotic-technology, bio-engineering, control theory, viscoelasticity diffusion model, relaxation processes and signal processing [1,2]. The order of derivatives, as well as integrals in the fractional-order calculus, is arbitrary. Therefore, fractional-order NPDEs have developed a fundamental interest in generalising integer-order NPDEs to model complex systems in thermodynamics, engineering, fluid dynamics and optical physics [3].

    The enormous advantage of using fractional differential equations (FDEs) in modeling real-world problems is their global behavior together with preserving memory [4] which is not present in integer-order differential equations. It has also been noted that FDEs fastly converge to ordinary differential equations (ODEs) in a case when fractional-order is equal to one. Moreover, fractional calculus can clarify the basic features of various models and processes them more precisely than integer-order [5]. Several techniques have been applied to study analytical as well as numerical solutions of FNPDEs such as, Variational Iteration Method (VIM) [6], Laplace transforms Method [7], the double Laplace transform [8], invariant subspace method [9], Integral transform [10], Sumudu transform Method (STM) [11], natural transform [12] and Adomian decomposition method (ADM) [13].

    The Klein-Gordon equation (KGE) considered herein is a basic non-linear evolution equation that arises in relativistic quantum Mechanics. It was formulated by Erwin Schrödinger for the non-relativistic wave equation in quantum physics, while precisely studied by the famous physicists O. Klein and W. Gordon (as it is named after their work) in 1926 [14,15]. The KGE has an extensive variety of applications in classical field theory [16] as well as in quantum field theory [17]. It has also been extensively used in numerous areas of physical phenomena such as in solid-state physics, dispersive wave-phenomena, nonlinear optics, elementary particles behavior, dislocations propagation in crystals, and different class of soliton solutions [18]. Here, we investigate equation of the form [19]

    ωψtω2ψx2+pψ+qg(ψ)=r(x,t),  1<ω2, (1.1)

    together with

    ψ(x,0)=F(x),ψt(x,0)=G(x),

    where ψ=ψ(x,t), g(ψ) and r(x,t) represent nonlinear term and external function respectively.

    The nonlinear differential equations involve numerous fractional differential operators, such as, Caputo, Hilfer, Riemann-Liouville (R-L), Atangana-Baleanu in Caputo's sense, and Caputo-Fabrizio, [20]. The above fractional operators are very useful in FC due to the complexities of fractional-PDEs/ODEs because standard operators cannot handle some equations to obtain explicit solutions. The Caputo fractional derivative is the basic idea of fractional derivatives. All the fractional derivatives will reduce in Caputo or Riemann-Liouville fractional derivatives after some parametric replacement. One can assume that the fractional derivative could provide a power-law of the local behavior of non-differentiable functions. The Caputo fractional derivative was introduced by Michele Caputo in 1967 [21] to study initial/boundary value problems in many areas of real-world phenomena. The Caputo's derivative has many advantages as it is the most important tool for dealing with integer order models in a fractional sense with suitable subsidiary conditions [23]. Most of the problems have been handled precisely using Caputo operator [24].

    The integer order KGE has broadly studied by using a variety of methods [25]. Time-fractional Klein Gordon equations with Caputo's fractional operator have also been extensively studied using a variety of numerical and analytical techniques [26]. Here, we apply double Laplace transform with decomposition method to study the general solutions of the governing model with power law. Some particular examples are also studied numerically with some physical analysis. For preliminaries and some basic definitions of Caputo's derivative, see [31] and the reference therein.

    Definition 1. Let us suppose ψ(x,t) lies in xtplane, the double Laplace transform (DLT) of ψ(x,t) is defined by [32]

    LxLt[ψ(x,t)]=0erx0estψdxdt,

    where, r,s(C).

    Definition 2. Application of DLT on fractional-order operator in Caputo's sense gives

    LxLt{CDωxψ(x,t)}=rω¯ψ(r,s)n1k=0rω1kLt{kψ(0,t)xk},

    and

    LxLt{CDβtψ(x,t)}=sβ¯ψ(p,s)m1k=0sβ1kLx{kψ(x,0)tk},

    where, m=[β]+1 and n=[ω]+1. Hence, we infer that

    LxLtψ(x)v(t)=¯ψ(p)¯v(s)=Lxψ(x)Ltv(t).

    The inverse DLT L1xL1t{¯ψ}=ψ, is represented by

    L1xL1t{¯ψ(x,t)}=12πic+iciestd+idiepx¯ψ(p,s)dpds,

    where Re(p)c and Re(s)d, and c,dR to be chosen appropriately.

    It is often more challenging to find the closed form of series solution to a nonlinear FDEs due to their complexity. Therefore question arises about the existence of the solution to such FDEs. For this, we utilize the applications of fixed point theory to study whether the solution of our considered system exists. So far, in the literature there exists no such theory for the existence of our considered system. We use here for the first time the βl-Geraghty type contraction to show that there exists a solution to the considered model. So we progress as follows

    CaDωtψψx+pψ+qg(ψ)=r,1<ω2, (2.1)

    with

    ψ(x,0)=F(x),ψt(x,0)=G(x). (2.2)

    The above equation can also be expressed in the form

    CaDωtψ=H(x,t,ψ),1<ω2, (2.3)

    where

    H(x,t,ψ)=ψxpψqg(ψ)+r. (2.4)

    For the existence of the above model, we use the following notions.

    Let Ω be the family of continuous and increasing functions defined as l:[0,)[0,) satisfying

    l(qx)ql(x)qx,q>1,

    and the elements of Θ are non-decreasing functions, such that

    ε:[0,)[0,1ρ21),whereρ11.

    Definition 3. Suppose that (M,d) be a complete b-metric space: Let T:MM also consider that F:M×M[0,) with F(m,n)l(ρ31d(Tm,Tn))ε(ld(m,n))l(d(m,n)), for m,nM,whereρ11,εΘandlΩ. Then T is known a generalized Fl-Gergaghty type contraction mapping.

    Definition 4. Consider T:MM, where M is non-empty and F:M×M[0,), where β(m,n)1β(Mm,Mn)m,nM, then T is called β-admissible mapping.

    First we show that there exists the fixed point for the considered model Eq (2.3), for this we apply the following theorem.

    Theorem 1. [33] Let T:MM be a generalized Fl-Gergaghty type contraction such that

    (1) T is β-admissible.

    (2) There exists υ0Mwithβ(υ0,Tυ0)1.

    (3) {υn}M,limnυn=υ, where υM and β(υn,υn+1)1β(υn,υ)1,

    then a fixed point for T. Let M=C(π,R) and d:M×M[0,), where π=[0,1]×[0,1] given by

    d(u,v)=||(uv)2||=supm[0,M]t[0,T](uv)2,

    thus (M,d) be a complete b-metric space. The following theorem shows the existence of solution of the considered model Eq (2.3).

    Theorem 2. Suppose that J:R2R such that

    (1) |H(x,t,ψ(x,t))H(x,t,ϕ(x,t))|α+133ε(l|uv|2)l(|uv|2), for x[0,X],t[0,T], and u,vM with J(u,v)0.

    (2) there exists u1M with J(u1,Tu1)0, where T:MM is defined by Tuj=u0u0+u1ut+IωtH(x,t,ψ(x,t)).

    (3) for u,vM,J(u,v)0J(Tun,Tv).

    (4) {un}M, unu where uM and mathcalJ(un,un+1)mathcalJ(un,u)0, for nN.Then there exists a solution of the model Eq (2.3).

    Proof. Applying the fractional integral to Eq (2.3), we obtain

    ψ(x,t)=C0ψ(x,0)+C1ψt(x,0)+IωtH(x,t,ψ(x,t))=Tψ(x,t).

    Here we prove that T has a fixed point using the above technique, thus

    |Tψ(x,t)Tϕ(x,t)|2=|IωtH(x,t,ψ(x,t))IωtH(x,t,ϕ(x,t))|2Iωt{|H(x,t,ψ(x,t))H(x,t,ϕ(x,t))|}2={1ω|H(x,t,ψ(x,t))H(x,t,ϕ(x,t))|}2{ωω33ωTωt0(ts)ω1ε(l(|uv|2)l(|uv|2))|}2{ωTωt0(ts)ω1ε(l(supx[0,X]t[0,T]|uv|2)l(supx[0,X]t[0,T]|uv|2))|}2133ε(l(d(u,v))l(d(u,v))).

    Hence for u,vC(π with J(u,v)0 we have 27||(TuTv)2||ε(l(d(u,v))l(d(u,v))). Now F:C([0,X]×[0,T],R)×C([0,X]×[0,T],R)[0,) by

    β(u,v)={1ifJ(u,v)0,0else,

    and

    β(u,v)l(27d(Tu,Tv))27d(Tu,Tv)ε(ld(u,v)ld(u,v)).

    Thus, T is an βl-contraction. Now to show that T is β-admissible, we have from condition (iii)

    β(u,v)1J(u,v)0J(Tu,Tv)0β(Tu,Tv)1.

    For u,vC(π,R), from condition (ii) we have uC(π,R). Such that β(u0,Tu0)1. Similarly from (iv) and Theorem 1, there exists uC(π,R), such that u=Tu. Therefore we proved that the model Eq (2.3) has a solution.

    Here, we study the above technique, which is a composition of DLT with the decomposition method. This method can be applied to find the general series solutions for various PDEs/ODEs. This is an efficient technique to study the analytical solutions of several nonlinear systems [34]. Let us consider the general non-linear system

    Lψ+Rψ+Nψ=r(x,t). (3.1)

    Here, L and R is linear and nonlinear operators, r(x,t) is some particular external function and N is nonlinearity in the system. The convergence analysis of the considered technique can be seen in [35].

    General solution of proposed model in Caputo's sense

    Using the technique defined above and expressing Eq (1.1) in the form

    CDωtψ(x,t)2ψx2+pψ+qg(ψ)=r(x,t),1<ω2, (3.2)

    with

    ψ(x,0)=F(x),ψt(x,t)=G(x). (3.3)

    Applying DLT to above equation, we obtain

    LxLt{CDωtψ}LxLt{2x2ψ}+pLxLt{ψ}+qLxLt{gψ}=LxLt{r(x,t)}. (3.4)

    Applying DLT on fractional order, gives

    LxLt{ψ}=1sLx{ψ(x,0)}+1s2Lx{ψt(x,0)}+1sωLxLt{2x2ψ}+p1sωLxLt{ψ}+q1sωLxLt{g(ψ)}+LxLt{r(x,t)}. (3.5)

    Similarly, applying Laplace transform on Eq (3.3), gives

    Lx{ψ(x,0)}=¯F(p),Lx{ψt(x,0)}=t¯G(p). (3.6)

    Now consider

    ψ=n=0ψn, (3.7)

    where the non-linear term can be degraded as

    g(ψ)=i=0An, (3.8)

    where An, is given by [36]

    An=1n!dndλn[nk=0λkg(ψk)]λ=0. (3.9)

    Finally, applying inverse DLT to Eq (3.2), using Eq (3.6) and Eq (3.9), gives

    ψ0=L1xL1t[1s¯F(p,0)]+tL1xL1t[1s2¯G(p,0)]=ψ(x,0),ψ1=L1xL1t[1sωLxLt{ψ0xx}]pL1xL1t[1sωLxLt{ψ0}]qL1xL1t[1sωLxLt{A0}]+[1sωLxLt{r(x,t)}],ψ2=L1xL1t[1sωLxLt{ψ1xx}]pL1xL1t[1sωLxLt{ψ1}]qL1xL1t[1sωLxLt{A1}],ψ3=L1xL1t[1sωLxLt{ψ2xx}]pL1xL1t[1sωLxLt{ψ2}]qL1xL1t[1sωLxLt{A2}].

    In a similar manner, other terms can be computed. Final result can be obtained as

    ψ(x,t)=n=0ψn(x,t). (3.10)

    which is the general solution of Eq (3.2) in series form by using the proposed method as discussed above.

    Here, we present numerical examples on the TFKG equation in Caputo's sense given as Eq (3.2) and discuss the behaviour of each example. We apply the aforesaid technique discussed in Section 3, to obtain the approximate solution of the problems.

    Example 1. Consider the nonlinear TFKG equation

    Dωtψ2ψx2+34ψ32ψ3=0,1<ω2,p=34,q=32, (4.1)
    g(ψ)=ψ3,r(x,t)=0, (4.2)

    with

    ψ(x,0)=sech(x),ψt(x,0)=12sech(x)tanh(x). (4.3)

    For α=2, the exact solution of Eq (4.1) can be obtained in the form [19]

    ψ(x,t)=sech(x+t2). (4.4)

    Consider TFKG Eq (4.1) in Caputo's sense

    CDωtψ2ψx2+34ψ32ψ3=0,1<ω2. (4.5)

    Applying MDLDM scheme discussed in Section 3, we obtain

    ψ0=sech(x)+t2sech(x)tanh(x),ψ1=tωΓ(ω+1)[1432sech2(x)]sech(x)+tω+1Γ(ω+2)[11874sech2(x)]sech(x)tanh(x)2!tω+2Γ(ω+3)[98sech3(x)tanh2(x)]+3!tω+3Γ(ω+4)[316sech3(x)tanh3(x)].

    In a similar manner, other terms can be computed. Final result can be obtained as

    ψ(x,t)=n=0ψn(x,t). (4.6)

    Discussion

    The error analysis between series solution Eq (4.6) and the exact solution Eq (4.4) is shown in Table 1, while the surface behaviour is shown in Figure 1 reveals that Eq (4.1) is depends on time (t). It should be noted that when the time (t) is small enough, there is less extent of error exist amongst the approximate and exact solutions obtained by the MDLDM method. Figure 2 [left panel] shows the absolute of wave solution Eq (4.6) with deviations in (α) with t=0.65 in comparison with exact solution Eq (4.4). Notice that numerical result, Eq (4.6) exactly matches to the exact solution Eq (4.4). This shows that the governing equation admits a soliton solution. Figure 2 [right panel] represents Eq (4.6) reveals that the amplitude of the solitary potentials blow-up as t increases. The 3D profiles for Eq (4.6) is shown in Figure 3 versus x for t=0.65. It reveals the progression of localized mode in the governing system. The solution obtained as Eq (4.6) versus x for t=1(dashedline),0.8(solidcurve),0.6(dottedcurve) Figure 4, when ω=2and1.7 are also depicted. Clearly, one can see the wave amplitude enhancement with variations in t that concludes that coefficient (ω) considerably increases the wave amplitudes.

    Table 1.  Comparison between the exact solution (4.4) with approximate solution obtained in the form (4.6).
    (x, t) Exact ψ Exactψ (x, t) Exact ψ Exactψ
    (-6, 0.6) 0.0067 0.0067 2.4083 ×105 (-4, 0.6) 0.0494 0.0492 1.7609×104
    (-2, 0.6) 0.3536 0.3529 6.2194×104 (0, 0.6) 0.9566 0.9550 1.6000×103
    (2, 0.6) 0.1985 0.1992 6.7883×104 (4, 0.6) 0.0271 0.0273 1.5198×104
    (6, 0.6) 0.0037 0.0037 2.0728×105 (-6, 0.2) 0.0055 0.0055 8.4718 ×107
    (-4, 0.2) 0.0405 0.0405 6.2013×106 (-2, 0.2) 0.2926 0.2926 2.4219×105
    (0, 0.2) 0.9950 0.9950 2.0749×105 (2, 0.2) 0.2413 0.2413 2.4874×105
    (4, 0.2) 0.0331 0.0331 5.9043×106 (6, 0.2) 0.0045 0.0045 8.0587×107
    (-6, 0.05) 0.0051 0.0051 1.2989×108 (-4, 0.05) 0.0375 0.0375 9.5114×108
    (-2, 0.05) 0.2723 0.2723 3.8310×107 (0, 0.05) 0.9997 0.9997 8.1360×108
    (2, 0.05) 0.2595 0.2595 3.8564×107 (4, 0.05) 0.0357 0.0357 9.3954×108
    (6, 0.05) 0.0048 0.0048 1.2828×108

     | Show Table
    DownLoad: CSV
    Figure 1.  The surface plot of the error analysis given in Table 1.
    Figure 2.  The left plot portrays comparison between Eq (4.4) and Eq (4.6) for various values of ω, while the right panel portrays solution profiles of ψ(x,t) vs t for various values of ω.
    Figure 3.  The surface plot of Eq (4.6) for the parameters used in the Figure 2's left panel.
    Figure 4.  The solution profiles of Eq (4.6) for desperate values of ω with desperate values of t.

    Example 2. Consider the time-fractional nonlinear KGE in the form

    Dωtϕ2ϕx2+qϕ3=0,1<ω2,p=0q=1, (4.7)
    y(ϕ)=ϕ3,r(x,t)=0, (4.8)

    with

    ϕ(x,0)=Rtan(λx),ϕt(x,0)=Rηλsec2(λx), (4.9)

    where

    R=ρκandλ=ρ2(σ+η2).

    The parameters ρ, κ, σ, and ηR. It should be noted that, for α=2, an exact solution of Eq (4.7) can be obtained in the form [37]

    ϕ(x,t)=Rtan[λ(x+ηt)]. (4.10)

    Writing Eq (4.7) in Caputo's sense gives

    CDωtϕ2ϕx2+ϕ3=0,1<ω2. (4.11)

    The series solution of Eq (4.11) with conditions (4.9) gives

    ϕ0=Rtan(λx)+tRηλsec2(λx),ϕ1=(tωΓ(ω+1))[2Rλ2sec2(λx)tan(λx)R2tan3(λx)]+Rηλ(tω+1Γ(ω+2))[4λ2sec2(λx)tan2(λx)+2λ2sec4(λx)3R2tan2(λx)sec2(λx)]3R3λ2η2(2!tω+2Γ(ω+3))[sec4(λx)tan(λx)]R3λ3η3(3!tω+3Γ(ω+4))[sech6(λx)].

    In a similar manner, other terms can be computed. The final result can be written in the form

    ϕ=n=0ϕn. (4.12)

    Discussion

    The parameters as κ=1,σ=8.5,η=0.05, and ρ=1 are considered for numerical illustration. The error amongst the approximate and exact solutions of Eq (4.7) is shown in Table 2 and its corresponding surface plot is presented in Figure 5. The numerical solution, Eq (4.12) and exact solution Eq (4.10) is depicted in Figure 6 [left panel], for t=7 and for different values of time-fractional coefficient (ω). It is noted that TFKG Eq (4.12) may admits the excitation in the system. This amount enrichment in ω overturned the wave amplitude as it interrupt the dispersion/nonlinearity effects. To see the effect of a temporal variable (t) on the wave solution, Eq (4.12) is displayed in Figure 6 [right panel] which shows that ϕ(x,t) increases with time. Further, the 3D profiles for Eq (4.12) is shown versus x with t=7 in Figure 7 for (ω=2), which represents the physical behaviour of Eq (4.12). It shows the advancement of localized mode excitations in the governing equation. The solution of Eq (4.12) versus x with t=7 (dashed curve), 6 (solid green curve), 5 (dotted curve), in Figure 8, with ω=2and1.5 respectively is depicted. Clearly, the wave amplitude increases with deviations in t. It infers that the fractional order (ω) significantly increases the wave amplitudes.

    Table 2.  Comparison between approximate solution obtained in the form (4.12) with exact solution (4.7).
    (x, t) Exact ϕ Exactϕ (x, t) Exact ϕ Exactϕ
    (-10, 0.02) 0.9845 0.9845 4×107 (-8, 0.02) 0.9596 0.9596 4×106
    (-6, 0.02) 0.8968 0.8965 3×105 (-4, 0.02) 0.7488 0.7486 2×105
    (-2, 0.02) 0.4504 0.4503 1×105 (0, 0.02) 0 0 0
    (2, 0.02) 0.4504 0.4503 1×105 (4, 0.02) 0.7488 0.7486 2×105
    (6, 0.02) 0.8968 0.8965 3×105 (8, 0.02) 0.9596 0.9596 4×106
    (10, 0.02) 0.9845 0.9845 4×107 (-10, 0.05) 0.9845 0.9845 1.8751×105
    (-8, 0.05) 0.9595 0.9596 4.8543×105 (-6, 0.05) 0.8966 0.8968 1.2196×104
    (-4, 0.05) 0.7486 0.7488 2.8545 ×104 (-2, 0.05) 0.4499 0.4499 5.8734×104
    (0, 0.05) 0.0012 0.0012 1.1300×106 (2, 0.05) 0.4508 0.4505 3.7952×104
    (4, 0.05) 0.7491 0.7488 2.4730×104 (6, 0.05) 0.8969 0.8968 1.1554×104
    (8, 0.05) 0.9596 0.9596 4.7525×105 (10, 0.05) 0.9845 0.9845 1.8587×105

     | Show Table
    DownLoad: CSV
    Figure 5.  The surface plot for Table 2.
    Figure 6.  Comparison between Eq (4.10) and Eq (4.12) for different values of ω [left panel]. The solution profiles of ψ(x,t) against time (t) interval with various values of ω [right panel].
    Figure 7.  The surface plot for the parameters used for the left panel of Figure 6.
    Figure 8.  The solution profiles of Eq (4.12) for different ω with different values of time(t).

    We have studied the TFKG equation using double Laplace transforms with the decomposition method. The general solution of proposed system is obtained as a class of general series solution. It is relevant to note that following only two iterations, fairly precise results are obtained that converges to the exact solution of the governing equation. The proposed method offers perfect numerical results without any alteration and complicated numerical methods for the governing equation in the fractional case. The numerical results obtained for particular examples are compared with the exact solutions at the classical order. The result profiles with physical interpretations for different fraction orders were revealed explicitly.

    Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R8). Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

    It is declared that all the authors have no conflict of interest regarding this manuscript.



    [1] A. A. K. Abu Hany, M. S. El-Khatib, K. A. K. Shourrab, On some fractional Hardy-Hilbert's integral inequalities, Gen. Lett. Math., 6 (2019), 35–44. https://doi.org/10.31559/glm2019.6.1.5 doi: 10.31559/glm2019.6.1.5
    [2] A. A. K. Abu Hany, M. Al Agha, Some extensions on Cerone's generalizations of Steffensen's inequality, Gen. Lett. Math., 3 (2017), 112–120. https://doi.org/10.31559/glm2016.3.2.4 doi: 10.31559/glm2016.3.2.4
    [3] A. A. K. Abu Hany, M. S. El-Khatib, K. A. K. Shourrab, Some new Hilbert type fractional inequalities, J. Nat. Stud., 27 (2019), 44–52.
    [4] D. S. Mitrinvic, The Steffensen inequality, Publ. Elektroteh. Fak. Univ. Beogradu Ser. Mat. Fiz., 1969,247–273. Available from: https://www.jstor.org/stable/43667360.
    [5] J. F. Steffensen, On certain inequalities between mean values, and their application to actuarial problems, Scand. Actuar. J., 1918 (1918), 82–97. https://doi.org/10.1080/03461238.1918.10405302 doi: 10.1080/03461238.1918.10405302
    [6] P. Cerone, On some generalizations of Steffensen's inequality and related results, J. Inequal. Pure Appl. Math., 2 (2001), 28.
    [7] U. N. Katugampola, A new fractional derivative with classical properties, arXiv, 2014. https://doi.org/10.48550/arXiv.1410.6535
    [8] D. R. Anderson, Taylor's formula and integral inequalities for conformable fractional derivatives, Contrib. Math. Eng., 2016, 25–43. https://doi.org/10.1007/978-3-319-31317-7-2 doi: 10.1007/978-3-319-31317-7-2
    [9] M. Z. Sarikaya, H. Yaldiz, H. Budak, Steffensen's integral inequality for conformable fractional integrals, Int. J. Anal. Appl., 15 (2016), 23–30.
    [10] R. Khalil, M. Al Horani, A. Yousef, M. Sababheh, A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65–70. https://doi.org/10.1016/j.cam.2014.01.002 doi: 10.1016/j.cam.2014.01.002
    [11] D. S. Mitrinović, J. E. Pečarić, A. M. Fink, Classical and new inequalities in analysis, Kluwer Dordrecht, 1993. https://doi.org/10.1007/978-94-017-1043-5
    [12] T. Abdeljawad, On conformable fractional calculus, J. Comput. Appl. Math., 279 (2015), 57–66. https://doi.org/10.1016/j.cam.2014.10.016 doi: 10.1016/j.cam.2014.10.016
    [13] K. S. Miller, B. Ross, An introdsction to fractional calculus and fractional differential equations, New York Wiley, 1993.
    [14] I. Podlubny, Fractional differential equations, Academic Press, 1998.
    [15] I. E. Abo Amra, M. M. Matar, Coupled system of fractional differential equations with impulsive and nonlocal coupled boundary conditions, Bol. Soc. Matematica Mex., 26 (2020), 477–497. https://doi.org/10.1007/s40590-019-00254-2 doi: 10.1007/s40590-019-00254-2
    [16] N. Tabouche, A. Berhail, M. M. Matar, J. Alzabut, A. G. M Selvam, D. Vignesh, Existence and stability analysis of solution for Mathieu fractional differential equations with applications on some physical phenomena, Iran. J. Sci. Technol. Trans. A Sci., 45 (2021), 973–982. https://doi.org/10.1007/s40995-021-01076-6 doi: 10.1007/s40995-021-01076-6
    [17] I. Suwan, M. S. Abdo, T. Abdeljawad, M. M. Matar, A. Boutiara, M. Almalahi, Existence theorems for ψ-fractional hybrid systems with periodic boundary conditions, AIMS Math., 7 (2021), 171–186. https://doi.org/10.3934/math.2022010 doi: 10.3934/math.2022010
    [18] S. Rezapour, S. T. M. Thabet, M. M. Matar, J. Alzabut, S. Etemad, Some existence and stability criteria to a generalized FBVP having fractional Composite-Laplacian operator, J. Funct. Spaces, 2021 (2021), 9554076. https://doi.org/10.1155/2021/9554076 doi: 10.1155/2021/9554076
    [19] M. M. Matar, M. Abu Jarad, M. Ahmad, A. Zada, S. Etemad, S. Rezapour, On the existence and stability of two positive solutions of a hybrid differential system of arbitrary fractional order via Avery–Anderson–Henderson criterion on cones, Adv. Differ. Equ., 2021 (2021), 423. https://doi.org/10.1186/s13662-021-03576-6 doi: 10.1186/s13662-021-03576-6
    [20] M. E. Samei, M. M. Matar, S. Etemad, S. Rezapour, On the generalized fractional snap boundary problems via G-Caputo operators: Existence and stability analysis, Adv. Differ. Equ., 2021 (2021), 498. https://doi.org/10.1186/s13662-021-03654-9 doi: 10.1186/s13662-021-03654-9
    [21] M. M. Matar, O. M. Al-Salmy, Existence and uniqueness of solution for conformable sequential differential equations, J. Nat. Sci., 19 (2017), 41–56.
    [22] S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional integrals and derivatives: Theory and applications, Gordon and Breach Science Publishers, 1993.
  • This article has been cited by:

    1. Yuan Xu, Siu-Long Lei, Hai-Wei Sun, An efficient red–black skewed extrapolation cascadic multigrid method for two-dimensional Poisson equation, 2024, 43, 2238-3603, 10.1007/s40314-023-02514-4
  • Reader Comments
  • © 2023 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1483) PDF downloads(64) Cited by(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog