Investigation on integro-differential equations with fractional boundary conditions by Atangana-Baleanu-Caputo derivative

Abstract

We establish, the existence and uniqueness of solutions to a class of Atangana-Baleanu (AB) derivative-based nonlinear fractional integro-differential equations with fractional boundary conditions by using special type of operators over general Banach and Hilbert spaces with bounded approximation numbers. The Leray-Schauder alternative theorem guarantees the existence solution and the Banach contraction principle is used to derive uniqueness solutions. Furthermore, we present an implicit numerical scheme based on the trapezoidal method for obtaining the numerical approximation to the solution. To illustrate our analytical and numerical findings, an example is provided and concluded in the final section.

1 Introduction

Fractional differential equations (FDEs) seemed as an excellent mathematical tool for modeling of many physical phenomena appearing in various branches, such as viscoelasticity, self-similar protein dynamics, continuum and statistical mechanics, dynamics of particles etc [1, 2]. A number of non-integer derivatives and their associated integrals have been developed and a systematic classification of fractional integrals with certain generalized Leibniz (type) rules was presented in [3]. Consequently, fractional calculus allows us to choose what kind of derivative to be used to solve a model, resulting in more accurate solutions [4–7].

Recently many authors have studied mathematical models involving AB fractional derivative. AB derivative in Caputo sense has also received considerable attention from researchers in the theoretical area of nonlinear FDEs [8, 9]. For instance, [10–12] addressed the existence of solutions to nonlinear Cauchy problems containing AB derivative.

Researchers in the theoretical domain of nonlinear fractional differential equations have shown significant interest in the Caputo sense AB derivative as well [13–16]. In this situation, a fractional derivative with a nonsingular kernel that relied on the M-L function was introduced by Atangana and Baleanu(A-B). This new dimension demonstrates the close connection between fractional calculus and the M-L function. It helps us to more effectively handle computing needs, nonlocal dynamics, and distinguish between various traits. [17–19].

In addition, the Gronwall inequality was established in the context of AB integral in order to analyze the Ulam-Hyers stability of the solution [20, 21]. The existence and qualitative behaviors of solutions of fractional integro-differential equations(FIDE) were investigated and which involved an extension model of [15, 22–24]. In [26], the authors proposed a new explicit numerical scheme based on two-step Lagrange polynomial for solving nonlinear FDEs of AB derivative and also presented the error analysis.

[25–28] has defined higher-order AB fractional operators and established some useful relations among the operators. Some recent progress on the existence of solutions to various higher-order nonlinear FDEs based on AB derivative with classical boundary conditions can be found in a series of papers [29–37].

Now, we consider the nonlinear FIDE with fractional boundary conditions of the following from:(1.1) (1.2)Where and denote the AB fractional derivatives (FD) in Caputo sense; is any real-valued absolutely continuous function and w₀, w₁ are real numbers. Here, are given function which fulfills certain conditions to be defined later on.

For our convenient, we consider and .

Then (1.1) becomes,(1.3)

In this work, we derive the necessary preliminaries in Section 2. In Section 3, some basic properties of AB operators, along with an auxiliary lemma for the linear version of (1.1), are derived. In Section 4, we study the existence of solution to (1.1). Section 5 contains an implicit numerical scheme for solving the proposed nonlinear FDE (1.1). In Section 6, illustrative examples for our obtained results are provided with their numerical simulations.

2 Preliminaries

Here, we recall basic concepts and useful lemmas related to AB fractional derivatives and AB fractional integral. Let C[a, b], a < b, be the Banach space of all real-valued continuous functions on [a, b] with the norm ||y||_∞ = sup_t∈[a,b]{|y(t)|}. Let AC[a, b] be the space of all real value absolutely continuous functions on [ab] and

Definition 2.1 [19, 26] The classical & two-parametric Mittag-Leffler functions are,respectively, where E_p,1(z) = E_p(z).

Definition 2.2 [2, 26, 38] The RL fractional integral of order q > 0 is, where w ∈ L¹(a, b).

Remark 2.3 [26] For 0 < q < 1, maps AC[a, b] into AC[a, b] and maps C[a, b] into C[a, b].

Definition 2.4 [8, 10, 18, 26, 32] Let 0 < q < 1. The AB FD in Caputo sense is, where w ∈ AC[a, b].

The AB FD in RL sense is, where w ∈ L¹(a, b).

The AB fractional integral is, where w ∈ L¹(a, b). Now, we consider the normalization function to be real-valued strictly positive such that .

Lemma 2.5 [8, 26] Let 0 < q < 1. Then (2.1)where w ∈ AC[a, b].

The higher-order AB fractional derivatives and integral defined in [3, 26, 29–37] based on Definition 2.4 are as follow:

Definition 2.6 [3, 26] Let m < q < m + 1(m ∈ N₀) and p = q − m and N₀ is denoting 0.

The AB FD in caputo sense is,where w ∈ AC^m+1[a, b].

The AB FD in RL sense is, where w ∈ L¹(a, b).

The AB fractional integral is,where w ∈ L¹(a, b).

Lemma 2.7 [3, 26] Let m < q < m + 1(m ∈ N₀). Then (2.2)where w ∈ AC^m+1[a, b].

3 Equivalent linear system solutions

This section is devoted to obtaining several important properties and results that will be helpful in our forthcoming discussion. We first define the Laplace transform of higher-order AB operators as follows:

Proposition 3.1. [26] Let . Then, in particular a = 0, we getprovided the Laplace transform of w exists.

Lemma 3.2 [26] Let m < q < m + 1(m ∈ N₀). Then, in particular a = 0 yields the following result,where w ∈ AC^m+1[a, b].

Proposition 3.3 [26] Let m < q < m + 1(m ∈ N₀). Then,where w ∈ AC[a, b].

Proof: From [26] and Definition 2.4, by setting p = q − m and using AB fractional derivative in Caputo sense and AB fractional integral operators, we get

Hence the final conclusion follows from Eq 2.1 and [[3], Proposition 3.1].

Lemma 3.4 [26] Let 0 < α ≤ 1 and λ > 0. Then, for all x > 0, the function xE_α,2(−λx^α) is strictly positive and monotonically increasing with respect to ‘x’.

Proof. From [2, 26], we have

In particular, by setting β = 1, the above equation reduces to

From[[39] in Lemma 2.2], xE_α,2(−λx^α), x > 0, is strictly positive and monotonically increasing with respect to ‘x’.

Remark 3.5 [26] Let, , 1 < q < 2, from Lemma 3.4, for all x > 0,

Lemma 3.6 Let with and 1 < q < 2. Then the solution to the linear system(3.1)is(3.2)where,

Proof. Use of Atangana-Baleanu integral operator on both sides of the equation in (3.1) with the help of Eq 2.2, results in(3.3)where . The condition w(0) = w₀ yields c₁ = w₀. Moreover, by Definition 2.6, we have

In view of Proposition 3.3, we obtain

By using boundary condition at t = b given in (3.1), it follows that

Substitution of the above values of c₁ and c₂ into (3.3) leads to the desired solution given in (3.2). This Lemma derive the integral equation for the problem (1.1).

Remark 3.7 [26] From 3.2, one can obtain

Therefore, our assumption together with Remark 2.3, that y ∈ AC²(J). By Definition 2.6, we can define the AB fractional derivative of the function y in Caputo sense on J. Moreover, it follows from Proposition 3.3 that w satisfies the linear problem (3.1).

For the elementary results of the Eq (1.3), we need some following hypotheses:

1. (H1) , for all (t, w) ∈ Ω and .

2. (H2) for a.e. and .

3. (H3) Let w ∈ C[0, T] & function is continuous function & a positive constant ζ₁, ζ₂ & ζ s.t., w₁, w₂, w₃, φ₁, φ₂, φ₃ in w. .

4 Existence result

We first convert the FIDE (1.1) to an integral equation and then consider a fixed-points problem for an integral operator. On the basis of Lemma 3.6, and solution operator as follows:(4.1)

We now only need to seek the fixed-points of in , because the fixed-points of are solutions of (1.1). For brevity, we use the notations

By considering the growth condition on nonlinearity, we now establish the result based on Leray-Schauder alternative theorem [40].

Theorem 4.1 Let and |w| ≤ λ for fixed λ > 0}. Suppose that is absolutely continuous and there exist φ_i ∈ C(J, [0, ∞)), i = 1, 2 and if , then (1.1) has a solution in C(T), provided

Proof. The operator is completely continuous.

Suppose is an arbitrary bounded set. Then, by using (H1), there is a number RL > 0, such that for

Now, , we obtain is bounded.

Next, we show that is equicontinuous. For all and 0 ≤ τ₁ < τ₂ ≤ b, we have

Since the function τ^q, 1 < q < 2, is uniformly continuous on is a family of equicontinuous functions. Therefore, following the Arzela-Ascoli theorem, is relatively compact set in C(T). Hence, we conclude that is completely continuous.

Now, we prove and λ ∈ (0, 1)} is a bounded set.

Let , then for t ∈ T. By (H1), for t ∈ T, one has

Hence, for all , we have is bounded. Therefore, the Leray-Schauder alternative theorem ensures us that has at least one fixed-point in C(T), as required. Hence the theorem follows. By considering the generalized Lipschitz condition on nonlinearity.

5 Uniqueness result

Theorem 5.1 Let and |w| ≤ λ for fixed λ > 0}. Suppose that is absolutely continuous and there exist ϕ ∈ L¹(J, [0, ∞)). If , then (1.1) has a unique solution in C(T), provided(5.1)

Proof. The operator is defined in (4.1). Let us first denote M = and fix a real number ρ, 0 < ρ ≤ λ, s.t.,

We now show that , where . By (H2), for all , one has

Therefore, for all , we deduce that maps .

Next, we prove that is a operator in C(T). Let . By (H2), for each t ∈ T, one has

Hence, for all we deduce that

The assumption (5.1) ensures that is contractive. Hence the Banach fixed point theorem allows us to conclude that has a unique fixed-point in C(T), as required. This fixed-point is the aspired solution to (1.1).

6 Numerical result

This section is devoted to describe a numerical scheme based on the classical and fractional trapezoidal methods, which are used to interpolate the vector field by the first-degree polynomials [26, 41].

The fractional trapezoidal method is a numerical technique used for approximating the definite integral of a function. It is a modification of the traditional trapezoidal rule, which divides the area under a curve into trapezoids and sums their areas to estimate the integral [42–46]. It is particularly useful when dealing with functions that may have singularities or other complexities and also can provide more accurate results compared to the traditional trapezoidal rule, especially when dealing with functions that have rapidly changing behavior or singularities. It’s a simple numerical integration method that can be computationally efficient for many problems but may require more intervals to achieve high accuracy [47–49].

Now, we can directly discretize the integral equation of (1.3) to derive the numerical scheme.

Let h be the step size, and choose n + 1 equi-spaced grid points on [0, b] such that 0 = t₀ < t₁ < t₂ < … < t_n = b with some positive integers n and .

We first rewrite the integral equation of (1.3) at any point t = t_k in the following piecewise way,

In order to formulate the numerical scheme for (1.3), the nonlinear function is approximated, on each subinterval [t_τ, t_τ + 1], by the following piecewise linear function,

After integrating, one obtains the following unknown approximations(6.1)

Due to the nonlinear nature of the function and the presence of all unknown quantities w(t_k) on the RHS of each equation of (6.1), we cannot solve (6.1) for w(t_k) directly, in general.

Then, the Newton-Raphson method is a numerical technique for finding approximate solutions to equations of the form f(x) = 0. It is a powerful iterative method commonly used for root finding and optimization problems.

The method starts with an initial guess and then refines that guess in each iteration until it converges to a root of the equation. It’s important to note that the Newton-Raphson method may not always converge to a solution, and it can even diverge if certain conditions are not met.

Therefore, it’s essential to choose a suitable initial guess and be aware of the characteristics of the function being analyzed. So Newton-Raphson method is used to determine each w(t_k);k = 0, 1, 2, …, n from the above implicit numerical scheme (6.1). This implicit method is seldomly encountered in the numerical simulation for solving fractional differential equations. However, this method is useful in dealing with the problem (1.3).

7 Example

Our theoretical results are illustrated here with the example.

Let the normalization function as We also present the trajectory of the numerical solution to the nonlinear problem (1.1) by using the proposed numerical scheme (6.1).(7.1)

Now, Consider the linear problem(7.2)where λ > 0 is a real number. Its exact solution is,

Now we consider t = 1, α = 0.5, Δ = 2, comparison of the exact and numerical solutions to Eq (6.1) for q = 1.5, λ = 1.5, b = 2, w₀ = 1, w₁ = 17 and h = 0.05.

Hence, t = 1; w(t) = 12.4692.

The numerical solution for t = 1 has been derived above. The derivation of numerical solutions w(t) for t ranging from 0 to 2 is similar and the values are presented in Table 1 below:

[Figure omitted. See PDF.]

The obtained results are demonstrated as a graph in Fig 1. The numerical solution shows agreement with the solution in the entire interval.

[Figure omitted. See PDF.]

8 Conclusion

In this paper, the essential outcomes of the existence and uniqueness solutions are gained by using the Leray-Schauder alternative theorem and Banach fixed point theorem respectively. We also discussed a numerical scheme for (1.1), which helps to study the numerical solution when is nonlinear. Our theoretical and numerical findings were verified with examples. We have presented some important properties of higher-order AB operators. However, the stability analysis of the solution may be considered as an extension work for the problem (1.1).

Acknowledgments

The authors extend their appreciation to the Deputy ship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IF2/PSAU/2022/01/22564).

Citation: Harisa SA, Faried N, Vijayaraj V, Ravichandran C, Morsy A (2024) Investigation on integro-differential equations with fractional boundary conditions by Atangana-Baleanu-Caputo derivative. PLoS ONE 19(5): e0301338. https://doi.org/10.1371/journal.pone.0301338

About the Authors:

Samy A. Harisa

Roles: Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

¶☯ All these authors are contributed equally to this work.

Affiliations: Department of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Aldawaser, Saudi Arabia, Department of Mathematics, Faculty of Science, Ain Shams University, Cairo, Egypt

Nashat Faried

Roles: Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

¶☯ All these authors are contributed equally to this work.

Affiliation: Department of Mathematics, Faculty of Science, Ain Shams University, Cairo, Egypt

V. Vijayaraj

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

¶☯ All these authors are contributed equally to this work.

Affiliation: Department of Mathematics, Kongunadu Arts and Science College, Coimbatore, Tamil Nadu, India

ORICD: https://orcid.org/0000-0002-9549-962X

C. Ravichandran

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

¶☯ All these authors are contributed equally to this work.

Affiliation: Department of Mathematics, Kongunadu Arts and Science College, Coimbatore, Tamil Nadu, India

ORICD: https://orcid.org/0000-0003-0214-1280

Ahmed Morsy

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

E-mail: [email protected]

¶☯ All these authors are contributed equally to this work.

Affiliation: Department of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Aldawaser, Saudi Arabia

ORICD: https://orcid.org/0000-0003-0375-6295

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Miller KS, Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley and Sons, New York. 1993.

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

3. Abdeljawad T. A Lyapunov type inequality for fractional operators with nonsingular Mittag-Leffler kernel. Journal of Inequalities and Applications. 2017;130: 1–11. pmid:28680233

4. Kamalesh K, Pramod CP, Pratibhamoy D, Higinio R. A graded mesh refinement approach for boundary layer originated singularly perturbed time-delayed parabolic convection diffusion problems. Mathematical Methods in the Applied Sciences. 2021;44(21): 12332–12350.

5. Kamal S, Arfan M, Ullah A, Qasem AM, Khursheed JA. T. Abdeljawad, Computational study on the dynamics of fractional order differential equations with applications. Chaos, Solitons & Fractals. 2022;157: 111955.

6. Pallavi B, Anoop K, Abdeljawad T, Aziz K. Study of Hilfer fractional evolution equations by the properties of controllability and stability. Alexandria Engineering Journal. 2021; 60(4): 3741–3749.

7. Pallavi B, Anoop K, Abdeljawad T, Khan Zareen A., Aziz K. Existence and approximate controllability of Hilfer fractional evolution equations with almost sectorial operators. Advances in Difference Equations. 2020; (2020): 615.

8. Atangana A, Baleanu D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Thermal Science. 2016;20(2): 763–769.

9. Al Nuwairan M, Gamal Ibrahim A. Nonlocal impulsive differential equations and inclusions involving Atangana-Baleanu fractional derivative in infinite dimensional spaces. AIMS Mathematics. 2023;8(5): 11752–11780. https://www.aimspress.com/article/doi/10.3934/math.2023595

10. Fernandez A, Joel ER, Suragan D. Linear differential equations with variable coefficients and Mittag-Leffler kernels. Alexandria Engineering Journal. 2021;61(2).

11. Hussain Alharthi N, Atangana A, Alkahtani Badr S. Study of a cauchy problem of fractional order derivative with variable order fractal dimension. Results in Physics. 2023;49: 106524.

12. Sombir D, Anoop K, Aziz K, Abdeljawad T, Alqudah Manar A. Existence Results of Langevin Equations with Caputo-Hadamard Fractional Operator. Journal of Mathematics. 2023;2023: 12.

13. Atangana A, Igret Araz S. Step forward on nonlinear differential equations with the Atangana-Baleanu derivative: Inequalities, existence, uniqueness and method. Chaos, Solitons & Fractals. 2023;173: 113700.

14. Atangana A, Seda IA. Analysis of a new partial integro-differential equation with mixed fractional operators. Chaos, Solitons & Fractals. 2019;127: 257–271.

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

16. Karthikeyan K, Trujillo JJ. Existence and uniqueness results for fractional integro-differential equations with boundary value conditions. Communications in Nonlinear Science and Numerical Simulation. 2012;17(11): 4037–4043.

17. Atangana A, Mishra J. Analysis of nonlinear ordinary differential equations with the generalized Mittag-Leffler kernel. Mathematical Biosciences & Engineering. 2023;20(11): 19763–19780. https://www.aimspress.com/article/doi/10.3934/mbe.2023875 pmid:38052623

18. Baleanu D, Fernandez A. On some new properties of fractional derivatives with Mittag-Leffler kernel. Communications in Nonlinear Science and Numerical Simulation. 2018;59: 444–462.

19. Gorenflo R, Kilbas AA, Mainardi F, Rogosin SV. Mittag-Leffler Functions, Related Topics and Applications. Vol. 2. Berlin: Springer. 2014.

20. Alb Lupas A, Catas A. Applications of the Atangana-Baleanu Fractional Integral Operator. Symmetry. 2022; 14(3): 630.

21. Srivastava HM, Nain AK, Vats RK, Pratibhamoy D. A theoretical study of the fractional-order p-Laplacian nonlinear Hadamard type turbulent flow models having the Ulam-Hyers stability. Revista de la Real Academia de Ciencias Exactas, Fisicas y Naturales. Serie A. Matematicas. Springer, 2023; 117:160. https://link.springer.com/article/10.1007/s13398-023-01488-6

22. Agarwal RP, Zhou Y, He Y. Existence of fractional neutral functional differential equations. Computer and Mathematics with Application. 2010;59(3): 1095–1100.

23. Kirti K, Anoop K, Aziz K, Abdeljawad T. Existence of solutions by fixed point theorem of general delay fractional differential equation with p-Laplacian operator. AIMS Mathematics. 2023; 8(5): 10160–10176. https://www.aimspress.com/article/doi/10.3934/math.2023514

24. Ravichandran C, Logeswari K, Jarad F. New results on existence in the framework of Atangana-Baleanu derivative for fractional integro-differential equations. Chaos, Solitons & Fractals. 2019;125: 194–200.

25. Abdeljawad T, Sabri TMT, Kedim I, Iadh AM, Aziz K. A higher-order extension of Atangana-Baleanu fractional operators with respect to another function and a Gronwall-type inequality. Boundary Value Problems. 2023;2023: 49.

26. Kumar SK, Sukavanam N, Pan S. Existence and uniqueness of solutions to fractional differential equations with fractional boundary conditions. Alexandria Engineering Journal. 2023;72(1): 147–155. https://doi.org/10.1016/j.aej.2023.03.076

27. Pratibhamoy D, Subrata R, Higinio R. Homotopy perturbation method for solving Caputo-type fractional-order Volterra-Fredholm integro-differential equations. Computational and Mathematical Methods. 2019;1(5): e1047. https://doi.org/10.1002/cmm4.1047

28. Sudarshan S, Jugal M, Pratibhamoy D, Debajyoti C. Higher order approximations for fractional order integro-parabolic partial differential equations on an adaptive mesh with error analysis. Computers & Mathematics with Applications. 2023;150: 87–101.

29. Pratibhamoy D, Subrata R, Higinio R. A Perturbation Based Approach for Solving Fractional Order Volterra-Fredholm Integro Differential Equations and its Convergence Analysis. International Journal of Computer Mathematics. 2019;97(10): 1–18. https://doi.org/10.1080/00207160.2019.1673892

30. Pratibhamoy D, Subrata R. Theoretical prospects of fractional order weakly singular Volterra Integro differential equations and their approximations with convergence analysis. Mathematical Methods in the Applied Sciences. 2021;44: 9419–9440.

31. Pratibhamoy D, Subrata R, Higinio R. On the approximate solutions of a class of fractional order nonlinear Volterra integro-differential initial value problems and boundary value problems of first kind and their convergence analysis. Journal of Computational and Applied Mathematics. 2022;404: 113116.

32. Pratibhamoy D, Subrata R, Jesus VA. Higher order accurate approximations on equidistributed meshes for boundary layer originated mixed type reaction diffusion systems with multiple scale nature. Applied Numerical Mathematics. 2020;148: 79–97.

33. Pratibhamoy D. A higher order difference method for singularly perturbed parabolic partial differential equations. Journal of Difference Equations and Applications. 2018;24(3): 452–477.

34. Pratibhamoy D. Comparison of a priori and a posteriori meshes for singularly perturbed nonlinear parameterized problems. Journal of Computational and Applied Mathematics. 2015;290: 16–25. https://doi.org/10.1016/j.cam.2015.04.034

35. Sumit S, Pratibhamoy D, Sunil K. Computatioal cost reduction for coupled system of multiple scale reaction diffusion problems with mixed type boundary conditions having boundary layers. Revista de la Real Academia de Ciencias Exactas, Físicas y Naturales. Serie A. Matematicas. 2023;117: article number 66. https://doi.org/s13398-023-01397-8

36. Ram S, Vembu S, Pratibhamoy D. A higher order hybrid-numerical approximation for a class of singularly perturbed two-dimensional convection-diffusion elliptic problem with non-smooth convection and source terms. Computers & Mathematics with Applications. 2023;142: 9–30.

37. Renu C, Satpal S, Pratibhamoy D, Devendra K. A higher order stable numerical approximation for time-fractional non-linear Kuramoto-Sivashinsky equation based on quintic B-spline. Mathematical Methods in the Applied Sciences. 2024;1–23. https://doi.org/10.1002/mma.9778

38. Kilbas AA, Srivastava HM, Trujillo JJ. Theory and Applications of Fractional Differential Equations, North-Holland Mathematics Studies. Vol. 204. Elsevier Science B. V., Amsterdam, 2006.

39. Wei Z, Li Q, Che J. Initial value problems for fractional differential equations involving Riemann-Liouville sequential fractional derivative. Journal of Mathematical Analysis and Applications. 2010;367(1): 260–272.

40. Granas A, Dugundji J. Fixed Point Theory. Vol. 14. New York: Springer. 2003.

41. Rohul A, Shah K, Asif M, Khan I, Ullah F. An efficient algorithm for numerical solution of fractional integro-differential equations via Haar wavelet. Journal of Computational and Applied Mathematics. 2021;381: 113028.

42. Deepti S, Jugal M, Pratibhamoy D, Jesus VA. A moving mesh refinement based optimal accurate uniformly convergent computational method for a parabolic system of boundary layer originated reaction-diffusion problems with arbitrary small diffusion terms. Journal of Computational and Applied Mathematics. 2022;404: 113167.

43. Pratibhamoy D. An a posteriori based convergence analysis for a nonlinear singularly perturbed system of delay differential equation equations on an adaptive mesh. Numerical Algorithms. 2019;81(1): 465–487. https://doi.org/10.1007/s11075-018-0557-4

44. Pratibhamoy D, Srinivasan N. Numerical solution of a system of singularly perturbed convection-diffusion boundary-value problems using mesh equidistribution technique. The Australian Journal of Mathematical Analysis and Applications. 2013;10(1): 1–17.

45. Seda IA. Numerical analysis of a new volterra integro-differential equation involving fractal-fractional operators. Chaos, Solitons & Fractals. 2020;130: 109396.

46. Sumit S, Pratibhamoy D, Sunil K, Parameter uniform higher order numerical treatment for singularly perturbed Robin type parabolic reaction diffusion multiple scale problems with large delay in time. Applied Numerical Mathematics. 2024;196: 1–21.

47. Li C, Zeng F. Numerical Methods for Fractional Calculus. Vol. 24. CRC Press. 2015.

48. Garrappa R. Numerical solution of fractional differential equations: A survey and a software tutorial. Mathematics. 2018;6(2): 16.

49. Diethelm K. Fundamental approaches for the numerical handling of fractional operators and time-fractional differential equations. In: Karniadakis, G.E. (ed.) Handbook of Fractional Calculus with Applications: Numerical Methods, vol. 3. De Gruyter, Berlin. 2019. https://doi.org/10.1515/9783110571684-001