Loading [MathJax]/jax/output/SVG/jax.js
Research article

Numerical method for solving the subdiffusion differential equation with nonlocal boundary conditions

  • Received: 09 September 2024 Revised: 23 December 2024 Accepted: 26 December 2024 Published: 31 December 2024
  • MSC : 35K20, 35R11, 65M06

  • This work was devoted to the construction of a numerical algorithm for solving the initial boundary value problem for the subdiffusion equation with nonlocal boundary conditions. For the case of not strongly regular boundary conditions, the well-known methods cannot be used. We applied an algorithm that consists of reducing the nonlocal problem to a sequential solution of two subproblems with local boundary conditions. The solution to the original problem was summed up from the solutions of the subproblems. To solve the subproblems, we constructed implicit difference schemes on the basis of the L1 formula for approximating the Caputo fractional derivative and central difference for approximating the space derivatives. Stability and convergence of the schemes were established. The Thomas algorithm was used to solve systems of linear algebraic equations. Numerical experiments were conducted to study the constructed algorithm. In terms of accuracy and stability, the algorithm performs well. The results of experiments confirmed that the convergence order of the method coincides with the theoretical one, O(τ2α+h2).

    Citation: Murat A. Sultanov, Vladimir E. Misilov, Makhmud A. Sadybekov. Numerical method for solving the subdiffusion differential equation with nonlocal boundary conditions[J]. AIMS Mathematics, 2024, 9(12): 36385-36404. doi: 10.3934/math.20241726

    Related Papers:

    [1] Muhammad Nazam, Aftab Hussain, Asim Asiri . On a common fixed point theorem in vector-valued b-metric spaces: Its consequences and application. AIMS Mathematics, 2023, 8(11): 26021-26044. doi: 10.3934/math.20231326
    [2] Shaoyuan Xu, Yan Han, Suzana Aleksić, Stojan Radenović . Fixed point results for nonlinear contractions of Perov type in abstract metric spaces with applications. AIMS Mathematics, 2022, 7(8): 14895-14921. doi: 10.3934/math.2022817
    [3] Mustafa Mudhesh, Hasanen A. Hammad, Eskandar Ameer, Muhammad Arshad, Fahd Jarad . Novel results on fixed-point methodologies for hybrid contraction mappings in Mb-metric spaces with an application. AIMS Mathematics, 2023, 8(1): 1530-1549. doi: 10.3934/math.2023077
    [4] Amer Hassan Albargi, Jamshaid Ahmad . Fixed point results of fuzzy mappings with applications. AIMS Mathematics, 2023, 8(5): 11572-11588. doi: 10.3934/math.2023586
    [5] Aftab Hussain . Fractional convex type contraction with solution of fractional differential equation. AIMS Mathematics, 2020, 5(5): 5364-5380. doi: 10.3934/math.2020344
    [6] Gunaseelan Mani, Arul Joseph Gnanaprakasam, Choonkil Park, Sungsik Yun . Orthogonal F-contractions on O-complete b-metric space. AIMS Mathematics, 2021, 6(8): 8315-8330. doi: 10.3934/math.2021481
    [7] Moirangthem Pradeep Singh, Yumnam Rohen, Khairul Habib Alam, Junaid Ahmad, Walid Emam . On fixed point and an application of C-algebra valued (α,β)-Bianchini-Grandolfi gauge contractions. AIMS Mathematics, 2024, 9(6): 15172-15189. doi: 10.3934/math.2024736
    [8] Qing Yang, Chuanzhi Bai . Fixed point theorem for orthogonal contraction of Hardy-Rogers-type mapping on O-complete metric spaces. AIMS Mathematics, 2020, 5(6): 5734-5742. doi: 10.3934/math.2020368
    [9] Mustafa Mudhesh, Aftab Hussain, Muhammad Arshad, Hamed AL-Sulami, Amjad Ali . New techniques on fixed point theorems for symmetric contraction mappings with its application. AIMS Mathematics, 2023, 8(4): 9118-9145. doi: 10.3934/math.2023457
    [10] Fatima M. Azmi . New fixed point results in double controlled metric type spaces with applications. AIMS Mathematics, 2023, 8(1): 1592-1609. doi: 10.3934/math.2023080
  • This work was devoted to the construction of a numerical algorithm for solving the initial boundary value problem for the subdiffusion equation with nonlocal boundary conditions. For the case of not strongly regular boundary conditions, the well-known methods cannot be used. We applied an algorithm that consists of reducing the nonlocal problem to a sequential solution of two subproblems with local boundary conditions. The solution to the original problem was summed up from the solutions of the subproblems. To solve the subproblems, we constructed implicit difference schemes on the basis of the L1 formula for approximating the Caputo fractional derivative and central difference for approximating the space derivatives. Stability and convergence of the schemes were established. The Thomas algorithm was used to solve systems of linear algebraic equations. Numerical experiments were conducted to study the constructed algorithm. In terms of accuracy and stability, the algorithm performs well. The results of experiments confirmed that the convergence order of the method coincides with the theoretical one, O(τ2α+h2).



    The Banach contraction principle [1] is considered a landmark fixed point theorem in the metric fixed point theory. The main component of the metric fixed-point theory is the exploration of new contraction principles to give fresh and helpful fixed-point theorems. Rakotch [2] was the first mathematician who involved a function instead of a Lipschitz constant in the Banach contraction principle. Then, Boyd and Wong [3] generalized the Rakotch contraction principle. In an effort to give a new contraction principle, Kannan [4] introduced a contraction principle that characterizes the metric completeness, and gave a new direction in metric fixed point theory that led many mathematicians to introduce various contraction principles. Among the classical contraction principles, the most famous are: Meir and Keeler contraction principle [5], Chatterjea contraction principle [6], Riech contraction principle [7], Hardy and Rogers contraction principle [8], Círíc contraction principle [9] and Caristi contraction principle [10].

    Recently, Wardowski [11] generalized Banach contraction principle by using an auxiliary nonlinear function F:(0,)(,) that satisfied three conditions. In the literature, this new contraction principle is known as F-contraction principle. This idea proved another milestone in metric fixed point theory. The F-contraction principle has been revisited and generalized in many abstract spaces (see [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] and references therein).

    On the other hand, metric generalization also has a significant impact on metric fixed point theory. Many mathematicians have contributed in this direction, producing many generalizations of a metric space (see [28]). One of these generalizations was done by Perov [29]. By extending the co-domain of the metric function from R to Rn, Perov [29] gave a vector version of the metric and hence produced another generalization of the Banach contraction principle. By following the Perov, Boriceanu [30] introduced a vector-valued b-metric and hence extended Perov [29] fixed point result to vector-valued b-metric spaces.

    In 2020, Altun et al. [31] presented a vector version of F-contraction principle and obtained an extension of Wardowski [11] fixed point theorem in the vector-valued metric spaces as follow:

    Theorem 1.1. [31] Every self-mapping T on a complete vector-valued metric space (X,d) that satisfies the following inequality:

    d(T(q),T(h))0IF(d(T(q),T(h)))F(d(q,h))q,hX,

    admits a unique fixed point, provided F satisfies (AF1)(AF3) and I=(τi)mi=10.

    Where the operator F:PmRm satisfies the following conditions:

    (AF1) Q,WPm with QW, we have F(Q)F(W);

    (AF2) {vn:nN}Pm, we have

    limnv(i)n=0 if and only if limnu(i)n=, for each i;

    (AF3) κ(0,1) satisfying limvi0+(vi)κui=0.

    In this paper, we investigate the possible conditions on F and T for which the mapping T admits a unique fixed point in a vector-valued b-metric space and in this way, we generalize the Theorem 1.1. Moreover, we introduce the (β,F)-contraction principle and call it Perov's type (β,F)-contraction principle. We obtain fixed point theorems on Perov's type (β,F)-contraction in a vector-valued b-metric space. Illustrative example and an application of the obtained fixed point theorem are given.

    In this section we present a summary of prerequisites and notations to be considered in the sequel. Let Rm={v=(xi)mi=1=(x1,x2,,xm)|ixiR} represents all matrices of order m×1 (will be called vectors), then (Rm,,) is a linear space with respect and defined by

    vw=(xi+yi)mi=1 for all v=(xi)mi=1 and w=(yi)mi=1Rm.
    kv=(kxi)mi=1 for all v=(xi)mi=1Rm and kR.

    Note that + and represent usual addition and multiplication of scalars. By using above operations, we can define difference of vectors as: vw=v(1)w. Define the relations and on Rm by

    vwxiyi and vwxi<yi;i. (2.1)

    The relation defines a partial-order on Rm. Let Pm denotes the set of positive definite vectors, that is, if v=(xi)mi=10 (zero vector of order m×1) and v=(xi)mi=1Rm then v=(xi)mi=1Pm. Also let Rm0={v=(xi)mi=1=(x1,x2,,xm)|ixi[0,)}. The two vectors are considered equal if their corresponding coordinates are equal.

    Definition 2.1. [32] 1) Let V=[vij] be an m×m complex matrix having eigenvalues λi, 1in. Then, the spectral radius ρ(V) of matrix V is defined by ρ(V)=max1im|λi|.

    2) The matrix V converges to zero, if the sequence {Vn;nN} converges to zero matrix O.

    Theorem 2.2. [32] Let V be any complex matrix of order m×m, then V is convergent if and only if ρ(V)<1.

    Perov [29] applied Theorem 2.2 to obtain the following result in the vector-valued metric spaces. It states that:

    Theorem 2.3. [29] Every self-mapping J defined on a complete vector-valued metric space (X,d) satisfying the following inequality:

    d(J(g),J(h))Ad(g,h)(g,h)X×X,

    admits a unique fixed point provided ρ(A)<1; A is a positive square matrix of order m.

    By a vector-valued metric, we mean a mapping d:X×XRm obeying all the axioms of the metric. The object d(x,y) is an m-tuple. Let

    v=(vi)mi=1=(v1,v2,v3,,vm)Pm,
    vn=(v(i)n)mi=1=(v(1)n,v(2)n,v(3)n,,v(m)n)Pm,
    F(v)=(ui)mi=1=(u1,u2,u3,,um)Rm and
    F(vn)=(u(i)n)mi=1=(u(1)n,u(2)n,u(3)n,,u(m)n)Rm,
    (v)m1=(v,v,v,,v)Rm.

    We organize this paper as follows:

    Section 2 contains definition and related properties of the vector-valued b-metric space. Section 3 consists of necessary lemmas and Perov's type (β,F)-contraction principle, related fixed point theorem and an example explaining hypothesis of obtained result. Section 4 consists of an application of the main theorem.

    In light of the definitions of b-metric and vector-valued metric given by Czerwik [33] and Perov [31] respectively, we proceed with the following definition.

    Definition 3.1. (vector-valued b-metric)[30] Let G be a non-empty set. The operator A:G×GRm0 satisfying the axioms (A1)(A3) given below is known as a vector-valued b-metric. For all q,t,gG, we have

    (A1) q=t if and only if A(q,t)=0.

    (A2) A(q,t)=A(t,q).

    (A3) A(q,g)s[A(q,t)A(t,g)];s1.

    The pair (G,A,s) represents a vector-valued b-metric-space.

    For s=1, every vector-valued b-metric space is a vector-valued metric-space, but this is not true when s>1. Thus, it can be remarked that every vector-valued metric-space is a vector-valued b-metric-space but not conversely.

    Example 3.2. Let G=R and the operator A:G×GRm0 is defined by

    A(l,q)=(|H|2,|H|3,,|H|m+1)l,qG,

    where H=|lq|. Then (G,A,s=2m) is a vector-valued b-metric space. Note that it is not vector-valued-metric space.

    Example 3.3. Let G be a non-empty set. Let di:G×G[0,) be a b-metric for each i with respective constant si1 (1im) for each positive integer i. The mapping A:G×GRm0 defined by

    A(q,t)=(d1(q,t),d2(q,t),,dm(q,t)) for all q,tG

    defines a vector-valued b-metric on G.

    The axioms (A1) and (A2) hold trivially. For (A3), for all l,q,tG, consider

    A(l,t)=(d1(l,t),d2(l,t),,dm(l,t))(s1(d1(l,q)+d1(q,t)),s2(d2(l,q)+d2(q,t)),,sm(dm(l,q)+dm(q,t)))s(d1(l,q)+d1(q,t),d2(l,q)+d2(q,t),,dm(l,q)+dm(q,t))=s((d1(l,q),d2(l,q),,dm(l,q))(d1(q,t),d2(q,t),,dm(q,t)))=s(A(l,q)A(q,t));s=max{s1,s2,,sm}.

    In general, the vector-valued b-metric is discontinuous, so, it requires an auxiliary convergence result to establish fixed point theorems in the vector-valued b-metric spaces. For this purpose we give the following lemma (Lemma 3.4).

    Lemma 3.4. Let (G,A,s) be a vector-valued b-metric space. If l,gG and {ln}nN is such that limnln=l, then

    1sA(l,g)limninfA(ln,g)limnsupA(ln,g)sA(l,g).

    Proof. By (A3), we have

    1sA(l,g)A(l,ln)A(ln,g)1sA(l,g)A(ln,l)A(ln,g).

    Taking liminf, we have

    1sA(l,g)limninfA(ln,g). (3.1)

    Again by (A3), we get

    A(ln,g)s(A(ln,l)A(l,g)).

    This implies

    limnsupA(ln,g)sA(l,g). (3.2)

    But also, we know that

    limninfA(ln,g)limnsupA(ln,g). (3.3)

    Combining (3.1)–(3.3), we get required result.

    Apart from Lemma 3.4, to fulfill the objective of this paper, the following compatibility condition is required:

    (AF4): for every positive term sequence vn=(x(i)n)mi=1, I=(τi)mi=10 satisfying

    IF(svn)F(vn1) implies IF(snvn)F(sn1vn1).

    Our findings rely mostly on the class of vector-valued nonlinear functions satisfying (AF1), (AF3) and (AF4) denoted by Πbs.

    Remark 3.5. The collection of vector-valued nonlinear functions Πbs is non-empty.

    Let F:PmRm be defined by F((xi)mi=1)=(ln(xi+1))mi=1 for all vPm, then (AF1) and (AF3) are obvious.

    We establish (AF4):

    Let IF(svn)F(vn1), then for m-tuple I=(ln(sn1),ln(sn1),,ln(sn1))=(ln(sn1))m1, we have

    (ln(sn1))m1F(s(x(i)n)mi=1)F((x(i)n1))mi=1(ln(sn1))m1(ln(sx(i)n+1))mi=1(ln(x(i)n1+1))mi=1(ln(snx(i)n+sn1))mi=1(ln(x(i)n1+1))mi=1ln(snx(i)n+sn1)ln(x(i)n1+1) for each isnx(i)nx(i)n1+1sn1 for each i.

    Now consider

    IF(snvn)=(ln(sn1))m1(ln(snx(i)n+1))mi=1(ln(sn1))m1(ln(x(i)n1+1sn1+1))mi=1=(ln(sn1x(i)n1s2n2+2sn1))mi=1=(ln(sn1x(i)n1+sn1(2sn1)))mi=1(ln(sn1x(i)n1+1))mi=1=F((sn1x(i)n1)mi=1)=F(sn1vn1).

    Hence, FΠbs.

    Example 3.6. Let F:PmRm be defined by

    (a) F(v)=(ln(xi))mi=1;

    (b) F(v)=(xi+ln(xi))mi=1;

    (c) F(v)=ln(x2i+xi)mi=1;

    (d) F(v)=(1xi)mi=1;

    (e) F(v)=(xai)mi=1;a>0;

    (f) F(v)=(ln(xi+1))mi=1.

    Among these definitions, definition (a)–(d) satisfy (AF1)(AF3) and definition (e), (f) belong to the family Πbs.

    Define F:R2+R2 by F((g1,g2))=(gt1,ln(g2+1));t>0 then FΠbs.

    The following lemma explains the reasons to omit axiom (AF2).

    Lemma 3.7. Let F satisfies (AF1) and {vn}nNPm is a decreasing sequence satisfying limnu(i)n=, then limnv(i)n=0 for each i{1,2,,m}.

    Proof. We note that for each i{1,2,,m}, {v(i)n}nN is bounded below and decreasing sequence of real numbers, so it is convergent. Let limnv(i)n=ζ0 for each i. Suppose on contrary ζ>0. Since, v(i)nζ for each i, therefore, u(i)nF(ζ). Thus, F(ζ)limnu(i)n=, a contradiction. Hence, limnv(i)n=0.

    Recently, Altun et al. [31] obtained an existence theorem involving vector-valued nonlinear function and explained it through various nontrivial examples. We will introduce and investigate the notion of (β,F)-contractions where the function F is taken from Πbs and β is defined below.

    Definition 4.1. Let there exists FΠbs and I0, the mapping T:(G,A,s)(G,A,s) is said to be a (s,F)-contraction, if it satisfies the following inequality:

    A(T(l),T(q))0IF(sA(T(l),T(q))F(A(l,q)), for alll,qG. (4.1)

    Remark 4.2. Note that for s=1, Definition 4.1 is identical to Perov's type F-contraction introduced by Altun et al. [31]. Thus, class of (s,F)-contractions (defined in Definition 4.1) is more wider as compared to that of Perov's type F-contraction introduced by Altun et al. [31]. Now we explain inequality (4.1) by the following example (Example 4.3).

    Example 4.3. Let G={ln=2n2n|nN}. Define A:G×GPm by A(l,q)=(|lq|2)m1, then (G,A,s=2) is a vector-valued b-metric space. Define the mapping ϕ:GG by

    ϕ(l)={2n12(n1)ifl=ln;l0ifl=l0.

    Take (1)m1=I0 and define F:PmRm by F((gi)mi=1)=(gi)mi=1. Then for every l,qG such that ϕ(l)ϕ(q), we have

    F(2A(ϕ(l),ϕ(q)))F(A(l,q))I.

    Indeed for l=ln+k and q=ln, consider

    2A(ϕ(ln+k),ϕ(ln))A(ln+k,ln)=((2n+k2(n+k1)2n2(n1))2(2n+k2(n+k)2n2(n))2)m1=(2n(12k2)(2k2(2n+2k1)(2n1)))m1(1)m1=(1)m1.

    Also we see that FΠbs. Indeed, for F((gi)mi=1)=(gi)mi=1, axioms (AF1) and (AF3) hold. For axiom (AF4), we proceed as follow: let IF(s(g(i)n)mi=1)F((g(i)n1))mi=1 that is 1+g(i)ng(i)n1 for each i{1,2,,m}. Now consider

    IF(n(g(i)n)mi=1)=In(g(i)n)mi=1=In1(g(i)n)mi=1In1(g(i)n11)mi=1=1+n1g(i)n1n1=1n1+n1g(i)n1 for each i(n1g(i)n1)mi=1=F(n1(g(i)n1)mi=1).

    This shows that for I=(1)m1, ϕ is an F-contraction.

    Remark 4.4 We observe that the function αs (defined in [34]) is superficial because we can always have a function β:G×G[0,) defined by β(l,q)=αs(l,q)s2 with following properties:

    (1) (ϕ is β-admissible)

    β(l,q)1impliesβ(ϕ(l),ϕ(q))1 for all l,qG,

    (2) the αs-completeness implies β-completeness and vice versa.

    Definition 4.5. Let G be a non-empty set and β:G×G[0,). The function ϕ:GG is said to be β-admissible if

    β(l,q)1 implies β(ϕ(l),ϕ(l))1 for all l,qG and 

    triangular β-admissible if in addition β follows:

    β(l,j)1,β(j,q)1, imply β(l,q)1.

    Definition 4.6. Let (G,A,s) be a vector-valued b-metric space and let β:G×G[0,). Let lG and sequence {ln}G. A mapping q:GG is a β-continuous at l=l0, if whenever,

    limnA(ln,l)=0 and β(ln,ln+1)1 we have limnA(q(ln),q(l))=0.

    Example 4.7. Let G=[0,) and define A:G×GRm0 by

    A(l,q)=(|H|2,|H|3,,|H|m+1)l,qG,

    where H=|lq|, and let q:GG be defined by

    q(l)={sin(πl)ifl[0,1];cos(πl)+2ifl(1,), β(l,q)={l+q+1ifl,q[0,1];0otherwise.

    Obviously, q is not continuous at l0=1, however, q is a β-continuous mapping at this point. Indeed, the assumption limnA(ln,l0)=0 leads us to choose ln=11n[0,1] and β(ln,ln+1)1 directs to choose [0,1] as domain of mapping q. Thus,

    limn|q(ln)q(l)|i=limn(sin(π(11n)))i=0 for each i;2im+1.

    Hence limnA(q(ln),q(l))=0.

    Definition 4.8. If an arbitrary Cauchy sequence {ln}G satisfying β(ln,ln+1)1 converges in G, the space (G,A,s) is called β-complete.

    Remark 4.9. Every complete vector-valued b-metric space is a β-complete vector-valued b-metric space but not conversely.

    Look at the following example.

    Example 4.10. Let G=(0,) and define the vector-valued b-metric A:G×GRm0 by

    A(l,q)=(|H|2,|H|3,,|H|m+1) for all l,qG,

    where H=|lq|. Define β:G×G[0,) by

    β(l,q)={l2+q2if l,q[2,5];0if not in[2,5].

    We observe that the space (G,A,s) is not a complete but it satisfies β-completeness criteria. Indeed, if {ln} is a Cauchy sequence in G such that β(ln,ln+1)1, for all nN, then ln[2,5]. Since [2,5] is a closed subset of R, so, there exists l[2,5] such that A(ln,l)0 as n. Hence (G,A,2m) is a β-complete vector-valued b-metric space.

    Definition 4.11. If an arbitrary sequence {ln}G satisfies the condition:

    β(ln,ln+1)1 and A(ln,l)0β(ln,l)1,

    nN. Then the space (G,A,s) is known as β -regular space.

    Let G=[2,5] and define vector-valued b-metric as in Example 3.2 and β as in Example 4.10. Let ln=2+3n be nth term of a sequence in G. Then (G,A,s) is a β-regular space.

    Definition 4.12. Let there exists FΠbs and I0 (G,A,s), the mapping T:GG is said to be a (β,F)-contraction, if it satisfies the following inequality:

    A(T(l),T(q))0l,qG,β(l,q)1 imply 
    IF(sβ(l,q)A(T(l),T(q))F(A(l,q)). (4.2)

    Remark 4.13. Every F-contraction is (β,F)-contraction but not conversely. For β(l,q)=1, we have (s,F)-contraction

    Suzuki [35] established the following lemma.

    Lemma 4.14. [35] If there is a number C>0 such that the sequence {xn}(G,d) satisfies the inequality:

    d(xn,xn+1)Cnν for every ν>1+log2s.

    Then {xn} is a Cauchy sequence.

    Now we have another lemma that extends the Lemma 4.14.

    Lemma 4.15. If there is a number C>0 such that the sequence {xn}(G,A,s) satisfies the inequality:

    A(xn,xn+1)(Cnν)m1foreveryν>1+log2sandforeverypositiveintegern.

    Then {xn} is a Cauchy-sequence.

    Proof. Let G be any non-empty set and s=max{si:1im}. Let di:G×G[0,) be a b-metric for every i{1,2,3,,m}and si1. Define the vector-valued b-metric A by

    A(q,t)=(di(q,t))mi=1 for all q,tG.

    Let {xn} be a sequence in G and assume that

    A(xn,xn+1)(Cnν)m1 for every ν>1+log2s and for every nN.

    Then, by definition of partial order defined by (2.1), we have for each i

    di(xn,xn+1)Cnν for every ν>1+log2si and for every nN.

    Since Lemma 4.14 does not depend on a particular b-metric, therefore, Lemma 4.14 can be applied for each di (1im). Thus, {xn} is a Cauchy sequence with respect to every di (1im). Thus,

    di(xn,xm)0 as n,m for each i.

    This leads us to write that

    A(xn,xm)(0,0,,0)=O as n,m.

    Hence {xn} is a Cauchy-sequence in (G,A,s).

    Now we have an analogue of Lemma 4.15 subject to F contraction.

    Lemma 4.16. Let {Dn} be a sequence in Pm where Dn:=(j(i)n)mi=1. Assume that there exist a mapping F:PmRm, I=(τi)mi=1O and k(0,):=1/1+log2s satisfying (AF3) and the following:

    nIF(snDn)F(D0). (4.3)

    Then Dn(Cn1k)m1.

    Proof. The inequality (4.3) implies limnF(snDn)=()m1 and by Lemma 3.7, we get limnsnDn=0. By (AF3),

    limn(snj(i)n)kϑ(i)n=0 for eachi;F(snDn):=(ϑ(i)n)mi=1Rm.

    By (4.3), we also have the following information for each i.

    (snj(i)n)kϑ(i)n(snj(i)n)kϑ(i)0(snj(i)n)knτi0. (4.4)

    As n in (4.4), we have

    limnn(snj(i)n)k=0 for eachi.

    Equivalently there exists a positive integer N1 such that n(snj(i)n)k1 for nN1. It then follows for each i that

    snj(i)n1n1kj(i)n1snn1k1sn1k.

    This implies j(i)nCn1k for nN1 and for each i, where C=s1. Hence Dn(Cn1k)m1.

    Now we give main theorem of this paper.

    Theorem 4.17. Let there exists l0G such that β(l0,f(l0))1, then, every β-admissible (β,F)-contraction f:GG defined on the β-complete vector-valued b-metric space (G,A,s) admits a fixed point provided it is β-continuous on (G,A,s) or (G,A,s) is a β-regular space with additional assumption that F is continuous.

    Proof. (a) Let l0G be as assumed and construct a Picard iterative-sequence {ln} of points in G such that, l1=f(l0), l2=f(r1) and generally ln=f(ln1). Given β(l0,f(l0))=β(l0,l1)1 and f is β-admissible, so, β(f(l0),f(l1))1 i.e β(l1,l2)1. This leads to a general formula β(ln),ln+1)1 for all non-negative integers n. If A(ln,ln+1)=0, then ln=ln+1=f(ln) as required. (Step 1) Let An=A(ln,ln+1)0 and since β(ln,ln+1)1 for all non-negative integers n, so, by contractive condition (4.1), we get

    F(sA(f(ln1),f(ln)))F(sβ(ln1,ln)A(f(ln1),f(ln)))F(An1)I.

    This implies

    IF(sAn)F(An1). (4.5)

    Due to (AF4), inequality (4.5) implies

    IF(snAn)F(sn1An1). (4.6)

    (Step 2) Let An1=A(ln1,ln)0 and since β(ln1,ln)1 for all positive integers n, so, by contractive condition (4.1), we get

    F(sA(f(ln2),f(ln1)))F(sβ(ln2,ln1)A(f(ln2),f(ln1)))F(An2)I.

    This implies

    IF(sAn1)F(An2). (4.7)

    The condition (AF4) in association with (4.7) implies

    IF(sn1An1)F(sn2An2).

    Thus, inequality (4.6) leads to have a new inequality:

    F(snAn)F(sn2An2)2I.

    Finally, Step n provides the following inequality:

    F(snAn)F(A0)nI, for all positive integers n. (4.8)

    By Lemma 4.16, {An}(Cn1k)m1. Since 1k>1+log2s, so, by Lemma 4.14, {ln} is a Cauchy sequence in β-complete vector-valued b-metric space (G,A,s). Thus, there exists (say) lG such that limnA(ln,l)=0. As, β(ln,ln+1)1 for all non-negative integers, so, by the β-continuity of f, we have

    limnA(f(ln),f(l))=limnA(ln+1,f(l))=0f.

    Since (G,A,s) is not continuous in general, by Lemma 3.4, we have

    1sA(l,f(l))limninfA(ln+1,f(l))=limnA(ln+1,f(l))=0f.

    This implies that A(l,f(l))=0f and by axiom (A1) of vector-valued b-metric, we get l=f(l). Hence, f admits a fixed point l.

    (b) Case 1. If there exists a subsequence {lni} of {ln} such that lni=f(l) for all positive integers i, then, l=limilni=limif(l)=f(l). As required.

    Case 2. Let there is no such subsequence of {ln} as in Case 1. We have proved in part (a) that limnA(ln,l)=0, equivalently, K0N such that for all nK0, A(f(ln),l)A(l,f(l)). Since (G,A,s) is β-regular, thus, β(ln,l)1. By contractive condition (4.1) and monotonicity of F, we have

    IF(sA(ln+1,f(l)))IF(sβ(ln,l)A(f(ln),f(l)))F(A(ln,l))F(A(l,f(l))). (4.9)

    We are looking for A(l,f(l))=0, aiming at contradiction, suppose on contrary that A(l,f(l))0. By Lemma 3.4, monotonicity and continuity of F, we have

    IF(A(l,f(l)))IF(slimninfA(ln+1,f(l))=IlimninfF(sA(ln+1,f(l))F(A(l,f(l))) by (4.9).

    This is a contradiction to I0. Thus, A(l,f(l))=0. Finally, by axiom A1, we obtain l=f(l).

    Remark 4.18. Additionally, if g is also a fixed point of f such that β(l,g)1, then f admits a unique fixed point.

    Now we explain the hypothesis of the above theorem with an example.

    Example 4.19. Let G=[0,) and define the mapping A:G×GRm0 by

    A(l,q)=(|H|2,|H|3,,|H|m+1) for all l,qG,

    where H=|lq|. Define β:G×G[0,) by

    β(l,q)={Kifflq(K[1,));0otherwise.

    so that (G,A,s) is a β-complete vector-valued b-metric space with s=2m. Define the mapping f:GG, for all lG by

    f(l)=ln(1+l6). Then,

    f is β-continuous self-mapping: Indeed, consider the sequence ln=Kn2 for all positive integers n. As, Kn2K(n+1)2, so, β(ln,ln+1)1 and limnA(ln,l)=0 implies (l2,l3,,lm+1)=0. This is true for l=0.

    Now limnA(f(ln),f(l))=limn((ln(1+K6n2))i)m+1i=2=0f.

    Thus, whenever β(ln,ln+1)1 and limnA(ln,l)=0, we have limnA(f(ln),f(l))=0f.

    f is β-admissible mapping: Indeed, let β(l,q)1, then, lq, thus, we have ln(1+l6)ln(1+q6) i.e β(f(l),f(q))1. Also let l0=1, f(l0)=ln(76). As, l0f(l0), so, β(l0,f(l0))1.

    Finally, for each l,qG with lq and choosing K such that >sK6i (2im+1), we have for Δf=|f(l)f(q)|

    (sβ(l,q)|Δf|i+1)m+1i=2=(sK|ln(1+l6)ln(1+q6)|i+1)m+1i=2(sK(l6q6)i+1)m+1i=2=(sK6i(lq)i+1)m+1i=2.
    Thus, sβ(l,q)|Δf|i+1sK6i(lq)i+1 for eachi.
    This implies, ln(6isK)+ln(sβ(l,q)|Δf|i+1)ln((lq)i+1) for each i.
    So that, (ln(6isK))m+1i=2(ln(sβ(l,q)|Δf|i+1))m+1i=2(ln(lq)i+1)m+1i=2.

    Let F:Rm+Rm be defined by F(v)=(ln(xi+1))mi=1, for all v=(xi)m+1i=2Rm+, then FΠbs (as shown in Section 2). Thus, for all l,qG such that A(f(l),f(q))0, I=(ln(6isK))m+1i=2 we obtain

    IF(sβ(l,q)A(f(l),g(q)))F(A(l,q)).

    Note that f has a unique fixed point l=0.

    Corollary 4.20. Let A:G×GR40 and (G,A,s) be a β-complete vector-valued b-metric space. If f:GG be a continuous self-mapping satisfying the inequality:

    s3A(f(l),f(q)))QA(l,q), (4.10)

    l,qG, τ>0 be such that eτ(0,11+log2s). Then f has a unique fixed point in G.

    Proof. Define F:R4+R4 by F(v)=(ln(vi))4i=1 and mapping β by β(l,q)=s2 l,qG in the proof of Theorem 4.17 and Q=(eτ0000eτ0000eτ0000eτ). Then the inequality (4.10) reduces to (4.1) and f has a unique fixed point in G.

    Note that for s=1, Corollary 4.20 represents the Perov's fixed point theorem [29].

    Perov [29], presented some applications of his fixed point theorem to Cauchy problems for the system of ordinary differential equations and respectively, to the boundary value problems. Here we will see that Theorem 4.17 is also applicable to the following delay integro-differential problem:

    l(t)=ttLf(h,l(h),l(h))dh. (5.1)

    The Eq (5.1) generalizes the following delay-integral-equation:

    l(t)=ttLf(h,l(h))dh. (5.2)

    The model for the spread of a few infectious diseases with a seasonally variable contact rate is represented by the Eq (5.2), where

    (a) l(t): The prevalence of infection at time t in the population.

    (b) 0<L: The amount of time a person can still spread disease.

    (c) l(t): The current rate of infectivity.

    (d) f(t,l(t),l(t)): The rate of newly acquired infections per unit of time.

    Now we look for existence and uniqueness of the positive, periodic solution to (5.1) by the application of Theorem 4.17.

    Let p>0 and fC(R×R+×R) satisfying

    f(t+p,l,q)=f(t,l,q)(t,l,q)R×R+×R.

    Let us define the functional spaces by

    F(p)={lC1(R):l(t+p)=l(t)tR}.F+(p)={lF(p):l(t)0tR}.

    Let V=F+(p)×F(p), define a metric A:V×VR2 by

    A((l1,q1),(l2,q2))=(l1l22,q1q22),

    where, l=max{|l(t)|:t[0,p],lF(p)}. The function β:V×V[1,) defined by β(l,q)=K2 for all l,qV. Then (V,A,s) is a β-complete vector-valued b-metric space. Now we develop the structure to apply Theorem 4.17. Let g(t)=l(t), then, we have

    g(t)=f(t,l(t),g(t))f(tL,l(tL),g(tL)).

    Thus, Eq (5.1) can be written as:

    {l(t)=ttLf(h,l(h),g(h))dhg(t)=f(t,l(t),g(t))f(tL,l(tL),g(tL)).

    Let Y:VC(R)×C(R) be a mapping defined by

    Y(l,q)=(Y1(l,q),Y2(l,q) for all (l,q)V,

    where Y1 and Y2 are defined by the following matrix equation:

    (Y1(l,q)(t)Y2(l,q)(t))=(ttLf(h,l(h),q(h))dhf(t,l(t),q(t))f(tL,l(tL),q(tL))).

    We need to assume the following conditions about function f:

    (I1) fC(R×R+×R) and there exists N1,N20 such that

    N1f(t,l,q)N2 for all (t,l,q)R×R+×R;

    (I2) fC(R×R+×R) follows the equation:

    f(t+p,l,q)=f(t,l,q);p>0;

    (I3) there exist ω,K>0 so that tR,l(t),u(t)R+ and q(t),ϵ(t)R

    |f(t,l(t),q(t))f(t,u(t),ϵ(t))|eωsK|l(t)u(t)|,

    (I4) there exists l0V such that β(l0,Y(l0))=K2.

    By the assumption of I1, we infer that

    N1LY1(l,q)(t)N2L for all tR and (l,q)V.

    Theorem 5.1. Let the function fC(R×R+×R) satisfies conditions (I1)(I4) and e2ωs2K2<14+L2, then, Eq (5.1) admits a solution in F+(p).

    Proof. We note that the conditions (1) and (2) of Theorem 4.17 can be verified by using (I4) and continuity of f respectively. In the following, we prove the contractive condition (4.2). By definition,

    Y1(l,q)(t+p)=t+pt+pLf(h,l(h),q(h))dh=ttLf(up,l(up),q(up))du=ttLf(up+p,l(up+p),q(up+p))du=ttLf(u,l(u),q(u))du=Y1(l,q)(t) for all tR,(l,q)F.

    This shows Y1(V)F+(p). Similarly, we have Y2(V)F+(p). Let (l1,q1),(l2,q2)V and consider

    |Y1(l1,q1)(t)Y1(l2,q2)(t)|2=|ttLf(h,l1(h),q1(h))dhttLf(h,l2(h),q2(h))dh|2(ttL|f(h,l1(h),q1(h))f(h,l2(h),q2(h))|dh)2(ttL(eωsK(|l1(h)l2(h)|))dh)2e2ωL2s2K2l1l22

    and

    |Y2(l1,q1)(t)Y2(l2,q2)(t)|2=|f(t,l1(t),q1(t))f(tL,l1(tL),q1(tL))f(t,l2(t),q2(t))+f(tL,l2(tL),q2(tL))|2(|f(t,l1(t),q1(t))f(t,l2(t),q2(t))|+|f(tL,l1(tL),q1(tL))f(tL,l2(tL),q2(tL))|)2(eωsK(|l1(t)l2(t)|)+eωsK(|l1(tL)l2(tL)|))24e2ωs2K2l1l22.

    Consequently, we obtain the following matrix inequality:

    (Y1(l1,q1)Y1(l2,q2)2Y2(l1,q1)Y2(l2,q2)2)(e2ωL2s2K2l1l224e2ωs2K2l1l22)=(e2ωL2s2K2004e2ωs2K2)(l1l22l1l22).

    Now define the mappings F:R2+R2 by F(l,q)=(ln(l),ln(q)) and I=(se2ωL2,s4e2ω), we obtain

    IF(sβ((l1,q1),(l2,q2))A(Y(l1,q1),Y(l2,q2)))F(A((l1,q1),(l2,q2))).

    Finally, keeping in mind the definition of mapping β and above inequality, we say that the mapping Y satisfies all the requirements of Theorem 4.17 and hence, admits a fixed point in its domain. Consequently, the Eq (5.1) has a positive, periodic solution.

    The findings and analyses discussed here could inspire further investigation into this topic by interested academics. A fundamental finding in vector-valued b-metric space, the main result (Theorem 4.17), concerns F-contraction in vector-valued b-metric spaces. In order to demonstrate the presence of solutions to various linear and nonlinear equations reflecting models of the associated real-world issues, the application approach is also discussed.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This project is funded by King Saud University, Riyadh, Saudi Arabia. Research Supporting Project number (RSP2023R167), King Saud University, Riyadh, Saudi Arabia.

    The authors declare that they have no competing interests.



    [1] H. Sun, Y. Zhang, D. Baleanu, W. Chen, Y. Chen, A new collection of real world applications of fractional calculus in science and engineering, Commun. Nonlinear Sci., 64 (2018), 213–231. https://doi.org/10.1016/j.cnsns.2018.04.019 doi: 10.1016/j.cnsns.2018.04.019
    [2] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier, 2006.
    [3] J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
    [4] R. Metzler, J. H. Jeon, A. G. Cherstvy, E. Barkai, Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking, Phys. Chem. Chem. Phys., 16 (2014), 24128–24164. https://doi.org/10.1039/C4CP03465A doi: 10.1039/C4CP03465A
    [5] M. J. Buckingham, On pore-fluid viscosity and the wave properties of saturated granular materials including marine sediments, J. Acoust. Soc. Am., 122 (2007), 1486–1501. https://doi.org/10.1121/1.2759167 doi: 10.1121/1.2759167
    [6] W. Chen, S. Hu, W. Cai, A causal fractional derivative model for acoustic wave propagation in lossy media, Arch. Appl. Mech., 86 (2016), 529–539. https://doi.org/10.1007/s00419-015-1043-2 doi: 10.1007/s00419-015-1043-2
    [7] M. Du, Z. Wang, H. Hu, Measuring memory with the order of fractional derivative, Sci. Rep., 3 (2013), 3431. https://doi.org/10.1038/srep03431 doi: 10.1038/srep03431
    [8] J. Zhang, Z. Wei, L. Xiao, Adaptive fractional-order multi-scale method for image denoising, J. Math. Imaging Vis., 43 (2012), 39–49. https://doi.org/10.1007/s10851-011-0285-z doi: 10.1007/s10851-011-0285-z
    [9] Q. Yang, D. Chen, T. Zhao, Y. Chen, Fractional calculus in image processing: a review, Fract. Calc. Appl. Anal., 19 (2016), 1222–1249. https://doi.org/10.1515/fca-2016-0063 doi: 10.1515/fca-2016-0063
    [10] R. Magin, Fractional calculus in bioengineering, part 1, Crit. Rev. Biomed. Eng., 32 (2004), 1–104. https://doi.org/10.1615/CritRevBiomedEng.v32.i1.10 doi: 10.1615/CritRevBiomedEng.v32.i1.10
    [11] C. Pinto, A. Carvalho, Fractional complex-order model for HIV infection with drug resistance during therapy, J. Vib. Control, 22 (2016), 2222–2239. https://doi.org/10.1177/1077546315574964 doi: 10.1177/1077546315574964
    [12] Y. Zhang, B. Baeumer, L. Chen, D. Reeves, H. Sun, A fully subordinated linear flow model for hillslope subsurface stormflow, Water Resour. Res., 53 (2017), 3491–3504. https://doi.org/10.1002/2016WR020192 doi: 10.1002/2016WR020192
    [13] V. E. Tarasov, V. V. Tarasova, Long and short memory in economics: fractional-order difference and differentiation, IRA-International Journal of Management and Social Sciences, 5 (2016), 327–334. https://dx.doi.org/10.21013/jmss.v5.n2.p10 doi: 10.21013/jmss.v5.n2.p10
    [14] A. V. Bitsadze, A. A. Samarskii, Some elementary generalizations of linear elliptic boundary value problems, Dokl. Akad. Nauk SSSR, 185 (1969), 739–740.
    [15] V. A. Il'in, Existence of a reduced system of eigen- and associated functions for a nonselfadjoint ordinary differential operator, Trudy Mat. Inst. Steklov., 142 (1976), 148–155.
    [16] N. I. Ionkin, The solution of a certain boundary value problem of the theory of heat conduction with a nonclassical boundary condition, Differ. Uravn., 13 (1977), 294–304.
    [17] N. I. Ionkin, E. I. Moiseev, A problem for a heat equation with two-point boundary conditions, Differ. Uravn., 15 (1979), 1284–1295.
    [18] V. P. Mikhailov, On Riesz bases in L2(0,1), Dokl. Akad. Nauk SSSR, 144 (1962), 981–984.
    [19] M. A. Naimark, Linear differential operators, London: Harrap, 1967.
    [20] N. Dunford, J. Schwartz, Linear operators, part III: spectral operators, New York: John Wiley & Sons Inc, 1957.
    [21] A. Mokin, On a family of initial-boundary value problems for the heat equation, Diff. Equat., 45 (2009), 126–141. https://doi.org/10.1134/S0012266109010133 doi: 10.1134/S0012266109010133
    [22] M. A. Sadybekov, Initial-boundary value problem for a heat equation with not strongly regular boundary conditions, In: Functional analysis in interdisciplinary applications, Cham: Springer, 2017,330–348. https://doi.org/10.1007/978-3-319-67053-9_32
    [23] A. S. Makin, On spectral decompositions corresponding to non-self-adjoint Sturm-Liouville operators, Dokl. Math., 73 (2006), 15–18. https://doi.org/10.1134/S1064562406010042 doi: 10.1134/S1064562406010042
    [24] A. V. Gulin, On the spectral stability in subspaces for difference schemes with nonlocal boundary conditions, Diff. Equat., 49 (2013), 815–823. https://doi.org/10.1134/S0012266113070045 doi: 10.1134/S0012266113070045
    [25] R. Ashurov, Y. Fayziev, On the nonlocal problems in time for time-fractional subdiffusion equations, Fractal Fract., 6 (2022), 41. https://doi.org/10.3390/fractalfract6010041 doi: 10.3390/fractalfract6010041
    [26] E. Ozbilge, F. Kanca, E. Özbilge, Inverse problem for a time fractional parabolic equation with nonlocal boundary conditions, Mathematics, 10 (2022), 1479. https://doi.org/10.3390/math10091479 doi: 10.3390/math10091479
    [27] M. A. Almalahi, M. S. Abdo, S. K. Panchal, Periodic boundary value problems for fractional implicit differential equations involving Hilfer fractional derivative, Probl. Anal. Issues Anal., 9 (2020), 16–44. https://doi.org/10.15393/j3.art.2020.7410 doi: 10.15393/j3.art.2020.7410
    [28] M. A. Almalahi, A. B. Ibrahim, A. Almutairi, O. Bazighifan, T. A. Aljaaidi, J. A. Awrejcewicz, A qualitative study on second-order nonlinear fractional differential evolution equations with generalized ABC operator, Symmetry, 14 (2022), 207. https://doi.org/10.3390/sym14020207 doi: 10.3390/sym14020207
    [29] A. A. Alikhanov, Stability and convergence of difference schemes approximating a two-parameter nonlocal boundary value problem for time-fractional diffusion equation, Comput. Math. Model., 26 (2015), 252–272. https://doi.org/10.1007/s10598-015-9271-4 doi: 10.1007/s10598-015-9271-4
    [30] N. B. Kerimov, M. I. Ismailov, Direct and inverse problems for the heat equation with a dynamic-type boundary condition, IMA J. Appl. Math., 80 (2015), 1519–1533. https://doi.org/10.1093/imamat/hxv005 doi: 10.1093/imamat/hxv005
    [31] N. H. Tuan, N. A. Triet, N. H. Luc, N. D. Phuong, On a time fractional diffusion with nonlocal in time conditions, Adv. Differ. Equ., 2021 (2021), 204. https://doi.org/10.1186/s13662-021-03365-1 doi: 10.1186/s13662-021-03365-1
    [32] Z. Sun, G. Gao, Fractional differential equations: finite difference methods, Berlin: De Gruyter, 2020. https://doi.org/10.1515/9783110616064
    [33] M. A. Sadybekov, I. N. Pankratova, Correct and stable algorithm for numerical solving nonlocal heat conduction problems with not strongly regular boundary conditions, Mathematics, 10 (2022), 3780. https://doi.org/10.3390/math10203780 doi: 10.3390/math10203780
    [34] M. A. Sultanov, E. N. Akimova, V. E. Misilov, Y. Nurlanuly, Parallel direct and iterative methods for solving the time-fractional diffusion equation on multicore processors, Mathematics, 10 (2022), 323. https://doi.org/10.3390/math10030323 doi: 10.3390/math10030323
    [35] E. N. Akimova, M. A. Sultanov, V. E. Misilov, Y. Nurlanuly, Parallel algorithm for solving the inverse two-dimensional fractional diffusion problem of identifying the source term, Fractal Fract., 7 (2023), 801. https://doi.org/10.3390/fractalfract7110801 doi: 10.3390/fractalfract7110801
    [36] A. A. Alikhanov, A time-fractional diffusion equation with generalized memory kernel in differential and difference settings with smooth solutions, Comput. Meth. Appl. Math., 17 (2017), 647–660. https://doi.org/10.1515/cmam-2017-0035 doi: 10.1515/cmam-2017-0035
    [37] M. A. Sultanov, M. I. Akylbaev, R. Ibragimov, Conditional stability of a solution of a difference scheme for an ill-posed Cauchy problem, Electron. J. Differ. Eq., 2018 (2018), 33.
    [38] K. Cao, D. Lesnic, M. I. Ismailov, Determination of the time-dependent thermal grooving coefficient, J. Appl. Math. Comput., 65 (2021), 199–221. https://doi.org/10.1007/s12190-020-01388-7 doi: 10.1007/s12190-020-01388-7
    [39] N. H. Luc, D. Baleanu, R. P. Agarwal, L. D. Long, Identifying the source function for time fractional diffusion with non-local in time conditions, Comp. Appl. Math., 40 (2021), 159. https://doi.org/10.1007/10.1007/s40314-021-01538-y doi: 10.1007/10.1007/s40314-021-01538-y
    [40] M. Slodička, Uniqueness for an inverse source problem of determining a space-dependent source in a non-autonomous time-fractional diffusion equation, Fract. Calc. Appl. Anal., 23 (2020), 1702–1711. https://doi.org/10.1515/fca-2020-0084 doi: 10.1515/fca-2020-0084
    [41] E. N. Akimova, V. E. Misilov, M. A. Sultanov, Regularized gradient algorithms for solving the nonlinear gravimetry problem for the multilayered medium, Math. Method. Appl. Sci., 45 (2022), 8760–8768. https://doi.org/10.1002/mma.7012 doi: 10.1002/mma.7012
    [42] Y. Zhang, A finite difference method for fractional partial differential equation, Appl. Math. Comput., 215 (2009), 524–529. https://doi.org/10.1016/j.amc.2009.05.018 doi: 10.1016/j.amc.2009.05.018
    [43] I. Podlubny, Fractional differential equations, London: Academic Press, 1999.
    [44] A. Samarskii, E. Nikolaev, Numerical methods for grid equations, volume I: direct methods, Basel: Birkhäuser, 1989. https://doi.org/10.1007/978-3-0348-9272-8
    [45] Y. Dimitrov, R. Miryanov, V. Todorov, Asymptotic expansions and approximations for the Caputo derivative, Comp. Appl. Math., 37 (2018), 5476–5499. https://doi.org/10.1007/s40314-018-0641-3 doi: 10.1007/s40314-018-0641-3
    [46] G. Gao, Z. Sun, H. Zhang, A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications, J. Comput. Phys., 259 (2014), 33–50. https://doi.org/10.1016/j.jcp.2013.11.017 doi: 10.1016/j.jcp.2013.11.017
    [47] A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, J. Comput. Phys., 280 (2015), 424–438. https://doi.org/10.1016/j.jcp.2014.09.031 doi: 10.1016/j.jcp.2014.09.031
  • This article has been cited by:

    1. Aziza Souad Nouar, Khadidja Nisse, Said Beloul, Solving a system of Caputo-Hadamard fractional differential equations via Perov’s fixed point theorem, 2025, 61, 1841-3307, 1, 10.2478/awutm-2025-0001
  • Reader Comments
  • © 2024 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(826) PDF downloads(58) Cited by(0)

Figures and Tables

Figures(2)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog