Next Article in Journal
Editorial Conclusion for the Special Issue “Advances in Computational Mechanics for Symmetrical Engineering Systems”
Next Article in Special Issue
On the Solution of Fractional Biswas–Milovic Model via Analytical Method
Previous Article in Journal
Solitons Solution of Riemann Wave Equation via Modified Exp Function Method
Previous Article in Special Issue
Fractional Series Solution Construction for Nonlinear Fractional Reaction-Diffusion Brusselator Model Utilizing Laplace Residual Power Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulations of the Fractional Systems of Volterra Integral Equations within the Chebyshev Pseudo-Spectral Method

1
Department of Mathematics and Computer Science, Faculty of Science and Technology, Rajamangala University of Technology Thanyaburi (RMUTT), Pathumthani 12110, Thailand
2
Department of Mathematics, Deanship of Applied Sciences, Umm Al-Qura University, Makkah 517, Saudi Arabia
3
Department of Mechanical Engineering, Sejong University, Seoul 05006, Republic of Korea
4
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work and are co-first authors.
Symmetry 2022, 14(12), 2575; https://doi.org/10.3390/sym14122575
Submission received: 17 November 2022 / Revised: 28 November 2022 / Accepted: 2 December 2022 / Published: 6 December 2022
(This article belongs to the Special Issue Functional Analysis, Fractional Operators and Symmetry/Asymmetry)

Abstract

:
In this article, we find the solutions to fractional Volterra-type integral equation nonlinear systems through a Chebyshev pseudo-spectral method (CPM). The fractional derivative is described in the Caputo manner. The suggested method’s accuracy and reliability are confirmed by the results. The proposed method is implemented for solving various nonlinear systems; the results we obtained were compared with the exact solution and other method solutions. The graphical representation and tables show that our method’s error quickly converges as compared to other methods. By comparing the proposed method’s solution with the actual solution and other methods, we can confirm that CPM is more accurate and closer to the exact solution. We display the pointwise solution in the tables, which verifies the proposed method’s accuracy at each point and aids in a better comprehension of the suggested approach. Moreover, the results of using the suggested method at different fractional orders are examined, showing that when a value moves from a fractional order to an integer order, the result is closer to the precise solution. Furthermore, the proposed technique for handling fractional-order linear and non-linear physical problems in science and engineering is straightforward to implement.

1. Introduction

Fractional-order calculus has been around for as long as integer-order calculus. It can be seen in a letter by Leibniz to L’Hopital on 30 September 1695. To this day, the question about D n ψ D ψ n , the Leibniz notation of the nth derivative of the linear function ξ ( ψ ) = ψ , appears in a letter from L’Hopital. L’Hopital strangely inquired, “What will be the result if n = 1 2 ? ”, to which Leibniz responded, “An seeming paradox, one day for which a useful result will be drawn”. The integer-order calculus is generalized to the fractional order calculus. Fractional calculus applications include anomalous transports in disordered systems [1], the time-fractional Belousov–Zhabotinsky reaction [2], dielectric relaxation phenomena in polymeric materials [3], long-time memory in a financial time series [4], and transport dynamics in a system governed by anomalous diffusion [5].
Many somatic problems in real life may be demonstrated using mathematical formulations, which convert physical occurrences into sophisticated mathematical formulae. Differential equations are used to simulate various physical phenomena, such as population growth or decay models [6,7,8]. However, some physical phenomena cannot be adequately represented using integer-order differential equations. As a result, the researchers created a new branch of mathematics known as fractional differential equations (FDEs). In comparison to integer-order differential equations, FDEs are employed to accurately simulate a variety of physical phenomena. In recent years, FDEs have gained prominence in the modeling of real-world physical problems, such as colored noise [9], economics [10], earthquake oscillation [11], and bioengineering [12]. The other applications are control theory [13], rheology [14], signal processing [15], damping method [16], polymers [17], and so on [18,19,20,21,22]. Mathematicians are interested in the solutions to these FDEs in order to obtain the objective solution of mathematical models. Analytical and numerical solutions are the most common forms of solutions. However, the analytical solutions we obtain are complex and final for all real-world problems. As a result, mathematicians solve these issues with approximate solutions [23,24]. Symmetry analysis is lovely to study when studying differential equations, more specifically when studying equations from the mathematics of finance. The secret to nature is symmetry, but most observations in the natural world lack it. A powerful technique for disguising symmetry is the occurrence of spontaneous symmetry-breaking. Finite and infinitesimal are two types of symmetries. Finite symmetries may be continuous or discrete. While parity and temporal reversal are discrete natural symmetries, space undergoes continuous modifications. Mathematicians have always been fascinated by patterns. Classifications of spatial and planar patterns made significant achievements in the eighteenth century.
Every branch of engineering and science uses fractional integral- and integro-differential equations. When a physical phenomenon is described using differential equations, the result is a differential equation, integro-differential equation, or integral equation. With these types of equations, some applications are glass-forming processes [25], nanohydrodynamics [26], drop-wise condensation [27], or wind ripples in the desert [28]. In most circumstances, there is no analytical solution to integral- or integro-differential equations. It is difficult to find, even if it exists in some circumstances. For approximating the solutions of integral- and integro-differential equations, various numerical approaches have been developed. To solve these challenges, a variety of approaches have been used, such as the nonstandard difference method (NDM) [29], Adams–Bashforth–Moulton method (ABMM) [30], homotopy analysis method (HAM) [31], generalized differential transform method (GDTM) [32], mixed interpolation collocation method (MICM) [33], collocation method (CM) [34], iterated Galerkin method (IGM) [35], spline collocation method (SCM) [36], Legendre wavelet method (LWM) [37], Newton polynomial (NP) [38], predictor–corrector (PC) [39], and Galerkin method (GM) [40].
To calculate the solutions of FIDEs in this study, we applied a new technique called the Chebyshev pseudo-spectral method (CPM). Our method reduces the proposed models to linear/non-linear systems of algebraic equations, considerably simplifying the problems; an appropriate method is then utilized to solve the resulting system. When compared to other methods, the suggested method has higher accuracy and better convergence. The numerical results show that the suggested technique is effective and reliable. In addition, CPM can be used in a variety of other physical models.

2. Definitions and Preliminaries Concept

This unit introduces the fundamental concepts related to fractional calculus.
Definition 1.
A real function ξ ( ψ ) , ψ > 0 , will be in the space C ν , ν R if a real number exists p > ν , with ξ ( ψ ) = ψ p ξ 1 ( ψ ) where ξ 1 ( ψ ) [ 0 , ) , and will be in the space C ν m if and only if ξ ( m ) C ν , m N .
Definition 2.
The fractional Caputo derivative of order α is stated as [41,42]
D α ξ ( ψ ) = 1 Γ ( n α ) 0 ψ ( ψ t ) n α 1 ξ ( n ) ( t ) d t ,
for n 1 < α n , n N , ψ > 0 , ξ C 1 m .
Definition 3.
Jin-Hunan’s fractional derivatives are referred to as [42]
D α ξ ( ψ ) D ψ α = Γ ( 1 + α ) lim Δ ψ = ψ 1 ψ 2 L ξ ( ψ 1 ) ξ ( ψ 2 ) ( ψ 1 ψ 2 ) α ,
where Δ ψ does not approach zero.
Definition 4.
Xiao-Jun defines fractional order derivatives as [42]
D ψ α ξ ( ψ 0 ) = ξ α ( ψ 0 ) = d α ξ ( ψ ) d ψ α ψ = ψ 0 = lim ψ ψ 0 Δ α ( ξ ( ψ ) ξ ( ψ 0 ) ) ( ψ ψ 0 ) α ,
where
Δ α ( ξ ( ψ ) ξ ( ψ 0 ) ) Γ ( 1 + α ) Δ ( ξ ( ψ ) ξ ( ψ 0 ) )
Definition 5.
The integral operator in the Riemann–Liouville sense is stated as [41,42]
I α ξ ( ψ ) = 1 Γ ( α ) 0 ψ ( ψ t ) α 1 ξ ( t ) d t ,
with the following properties
D α I α ξ ( ψ ) = ξ ( ψ ) ,
I α D α ξ ( ψ ) = ξ ( ψ ) k = 0 n 1 ξ ( k ) ( 0 + ) k ! ψ k , ψ 0 n 1 < α < n .

3. Chebyshev Pseudo-Spectral Method (CPM)

Chebyshev polynomials are defined as the interval [ 1 , 1 ] and are demonstrated using recurrence equations [43,44].
T n + 1 ( t ) = 2 u T n ( ψ ) T n 1 ( ψ ) , n = 1 , 2 , . . .
where
T 0 ( ψ ) = 1 , T 1 ( ψ ) = ψ .
In order to apply the Chebyshev polynomials in the interval [ 0 , 1 ] , Chebyshev’s shifted polynomials are described as T ^ n ( ψ ) , which explains (in the same way) the Chebyshev polynomials T n ( ψ ) by relation
T ^ n ( ψ ) = T n ( 2 ψ 1 ) .
The recurrence formula is as follows
T ^ n + 1 ( ψ ) = 2 ( 2 ψ 1 ) T n ^ ( ψ ) T ^ n 1 ( ψ ) , n = 1 , 2 , . . .
where
T ^ 0 ( ψ ) = 1 , T ^ 1 ( ψ ) = 2 ψ 1 .
In terms of Chebyshev’s shifted polynomials, a function ξ ( ψ ) L 2 [ 0 , 1 ] is described as
ξ ( ψ ) = n = 1 c n T ^ n ( ψ ) .
The first (m + 1) terms of Chebyshev’s shifted polynomials are considered as
ξ m ( ψ ) = n = 0 m c n T ^ n ( ψ ) ,
0 D x α n = 0 m c n T ^ n ( ψ ) + n = 0 m c n T ^ n ( ψ ) = g ( ψ , ξ ) .
We have the ability to find a system of equations as
0 D x α n = 0 m c n T ^ n ( ψ i ) + n = 0 m c n T ^ n ( ψ i ) = g ( ψ i , ξ ) .
Whereas
ψ i = i 0.5 2 k 1 M .
We used maple software to solve the resultant system, which provides a CPM solution.

4. Applications

Example 1.
Consider the following fractional integral equation nonlinear system
0 D ψ α ( ( 1 ψ 2 + t ) + ( ξ ( t ) + ζ 3 ( t ) ) d t = 1 12 ψ 6 2 15 ψ 5 + 1 4 ψ 4 + 1 3 ψ 3 , 0 D ψ α ( ( 5 + ψ t ) ( ξ 3 ( t ) ζ ( t ) ) d t = 1 56 ψ 8 + 5 7 ψ 7 1 6 ψ 3 5 2 ψ 2 , 0 < α 1
with the exact solution ξ ( ψ ) = ψ 2 , ζ ( ψ ) = ψ ,
In Table 1, we present the accurate and numerical results obtained by implementing the suggested approach and the results obtained by CWM while Table 2 shows the absolute error comparison of the suggested approach and results obtained through CWM at m = 6 and the radial basis function network (RBFN) at m = 8 . Figure 1 and Figure 2 show the behavior of the exact solution (E.S) and approximate solution (our technique) of this case when α = 1 , while Figure 3 and Figure 4 show the error comparison of CPM, CWM, and RBFN. Moreover, Figure 5 demonstrate the graphical behaviors of the solutions for various fractional orders, showing that as the value of α moves toward the integer-order from the fractional order, the solution converges to the exact.
Example 2.
Consider the following fractional integral equation nonlinear system
0 D ψ α ( ξ 2 ( t ) ζ 2 ( t ) ) d t = ξ ( ψ ) sin ( ψ ) + ψ , 0 D ψ α ( ξ ( t ) ζ ( t ) ) d t = ζ ( ψ ) + 1 2 sin 2 ψ + cos ( ψ ) , 0 < α 1
with the exact solution ξ ( ψ ) = sin ( ψ ) , ζ ( ψ ) = cos ( ψ ) ,
The accurate and numerical results obtained by implementing the suggested approach and results obtained by CWM for m = 6 are shown in Table 3. In Table 4, the errors implemented by the current approach are compared to those acquired by CWM. In Figure 6 and Figure 7, we compare the accurate results with our method’s results, showing that they are very near to each other. Moreover, Figure 8 and Figure 9 illustrate the CPM and CWM error comparisons, showing that CPM is in good agreement with the exact results, whereas Figure 10 and Figure 11 show graphical representations of the solutions for various fractional orders.
Example 3.
Consider the following fractional integral equation nonlinear system:
0 D ψ α ( ( 5 + ψ t ) ξ ( t ) + ( ψ 2 2 + t ) ζ ( t ) ϕ ( t ) ) = 1 48 ψ 6 + 19 270 ψ 5 + 19 72 ψ 4 + 7 6 ψ 3 + ψ 2 + 5 ψ , 0 D ψ α ( ( ψ 2 2 + t ) ξ ( t ) + ( 3 + ψ t ) ζ ( t ) + 1 4 ( ψ 2 t 2 ) ϕ ( t ) ) = 1 24 ψ 5 + 35 288 ψ 4 + 17 18 ψ 3 + 5 4 ψ 2 9 2 ψ , 0 D ψ α ( t ξ ( t ) ζ ( t ) ψ t ζ 2 ( t ) 5 ϕ ( t ) ) = 1 72 x 6 1 4 ψ 5 1 54 ψ 7 + 17 96 ψ 4 9 8 ψ 3 1 2 ψ 2 10 3 ψ , 0 < α 1
with the exact solution ξ ( ψ ) = 1 + 1 4 ψ 2 , ζ ( ψ ) = 3 2 + 1 2 ψ 2 ,   ϕ ( ψ ) = 2 3 + 1 2 ψ We use the technique described in Section 2 for α = 1 with m = 6 to solve this problem. Table 5 shows the exact and CPM estimated solutions. Table 6 compares the absolute error of our technique to that of CWM. The results of the proposed technique are in good agreement with the exact results, as shown in Figure 12 and Figure 13. The error comparison between CPM and CWM is provided in Figure 14, Figure 15 and Figure 16 to demonstrate the usefulness of CPM. In addition, in Figure 17, Figure 18 and Figure 19 we can see the estimated solutions for the various fractional orders.

5. Conclusions

In this article, we find the solutions to nonlinear fractional equation systems by implementing the Chebyshev pseudo-spectral method. This type of problem is reduced to the solution of a system of linear and nonlinear algebraic equations using the method proposed. The solutions obtained by utilizing the suggested approach are in good agreement with the actual results and are more accurate than those of other methods. Moreover, it can be confirmed from the figures and tables that our method’s error quickly converges compared to other techniques. Maple was used to perform the calculations in this article.

Author Contributions

Conceptualization, J.D.C.; Methodology, M.N.; Software, P.S.; Formal analysis, M.N. and N.A.S.; Investigation, J.D.C.; Data curation, N.A.S.; Writing—original draft, P.S. and R.S.; Writing—review & editing, J.D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work under Grant Code number: 22UQU4310396DSR37.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bisquert, J.; Garcia-Belmonte, G.; Fabregat-Santiago, F.; Compte, A. Anomalous transport effects in the impedance of porous film electrodes. Electrochem. Commun. 1999, 1, 429–435. [Google Scholar] [CrossRef]
  2. Alaoui, M.K.; Fayyaz, R.; Khan, A.; Shah, R.; Abdo, M.S. Analytical investigation of Noyes-Field model for time-fractional Belousov-Zhabotinsky reaction. Complexity 2021, 2021, 1–21. [Google Scholar] [CrossRef]
  3. Reyes-Melo, E.; Martinez-Vega, J.; Guerrero-Salazar, C.; Ortiz-Mendez, U. Application of fractional calculus to the modeling of dielectirc relaxtion phenomena in polymetric materials. J. Appl. Polym. Sci. 2005, 98, 923. [Google Scholar] [CrossRef]
  4. Picozzi, S.; West, B.J. Fractional Langevin model of memory in financial market. Phys. Rev. E 2002, 66, 046118. [Google Scholar] [CrossRef]
  5. Henry, B.I.; Wearne, S.L. Existance of Turing instabilities in a two-species fractional reaction-diffusion system. SIAM J. Appl. Math. 2002, 62, 870. [Google Scholar] [CrossRef] [Green Version]
  6. Xie, Z.; Feng, X.; Chen, X. Partial least trimmed squares regression. Chemom. Intell. Lab. Syst. 2022, 221, 104486. [Google Scholar] [CrossRef]
  7. Kovalnogov, V.N.; Fedorov, R.V.; Karpukhina, T.V.; Simos, T.E.; Tsitouras, C. Sixth Order Numerov-Type Methods with Coefficients Trained to Perform Best on Problems with Oscillating Solutions. Mathematics 2021, 9, 2756. [Google Scholar] [CrossRef]
  8. Shao, Z.; Zhai, Q.; Han, Z.; Guan, X. A linear AC unit commitment formulation: An application of data-driven linear power flow model. Int. J. Electr. Power Energy Syst. 2023, 145, 108673. [Google Scholar] [CrossRef]
  9. Mandelbrot, B. Some noises with i/f spectrum, a bridge between direct current and white noise. IEEE Trans. Inf. Theory 1967, 13, 289–298. [Google Scholar] [CrossRef]
  10. Baillie, R.T. Long memory processes and fractional integration in econometrics. J. Econom. 1996, 73, 5–59. [Google Scholar] [CrossRef]
  11. He, J.H. Nonlinear oscillation with fractional derivative and its applications. In Proceedings of the International Conference on Vibrating Engineering, Dalian, China, 6–9 August 1998; Volume 98, pp. 288–291. [Google Scholar]
  12. Magin, R.L. Fractional calculus in bioengineering, part 1. Crit. Rev. Biomed. Eng. 2004, 32, 1–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Vinagre, B.M.; Podlubny, I.; Hernandez, A.; Feliu, V. Some approximations of fractional order operators used in control theory and applications. Fract. Calc. Appl. Anal. 2000, 3, 231–248. [Google Scholar]
  14. Smit, W.; De Vries, H. Rheological models containing fractional derivatives. Rheol. Acta 1970, 9, 525–534. [Google Scholar] [CrossRef]
  15. Cuesta, E.; Fernandez-Manso, A.; Quintano, C. Generalized fractional integrals in advanced remote sensing. In Proceedings of the 2016 12th IEEE/ASME International Conference on Mechatronic and Embedded Systems and Applications (MESA), Auckland, New Zealand, 29–31 August 2016; pp. 1–4. [Google Scholar]
  16. Salo, M.; Tuusa, H. A vector controlled current-source pwm rectifier with a novel current damping method. IEEE Trans. Power Electron. 2000, 15, 464–470. [Google Scholar] [CrossRef]
  17. Metzler, R.; Schick, W.; Kilian, H.-G.; Nonnenmacher, T.F. Relaxation in filled polymers: A fractional calculus approach. J. Chem. Phys. 1995, 103, 7180–7186. [Google Scholar] [CrossRef]
  18. Areshi, M.; Khan, A.; Shah, R.; Nonlaopon, K. Analytical investigation of fractional-order Newell-Whitehead-Segel equations via a novel transform. AIMS Math. 2022, 7, 6936–6958. [Google Scholar] [CrossRef]
  19. Zidan, A.M.; Khan, A.; Shah, R.; Alaoui, M.K.; Weera, W. Evaluation of time-fractional Fishers equations with the help of analytical methods. AIMS Math. 2022, 7, 18746–18766. [Google Scholar] [CrossRef]
  20. Alyobi, S.; Shah, R.; Khan, A.; Shah, N.A.; Nonlaopon, K. Fractional Analysis of Nonlinear Boussinesq Equation under Atangana-Baleanu-Caputo Operator. Symmetry 2022, 14, 2417. [Google Scholar] [CrossRef]
  21. Sunthrayuth, P.; Alyousef, H.A.; El-Tantawy, S.A.; Khan, A.; Wyal, N. Solving Fractional-Order Diffusion Equations in a Plasma and Fluids via a Novel Transform. J. Funct. Spaces 2022, 2022, 1899130. [Google Scholar] [CrossRef]
  22. Shah, N.A.; El-Zahar, E.R.; Akgul, A.; Khan, A.; Kafle, J. Analysis of Fractional-Order Regularized Long-Wave Models via a Novel Transform. J. Funct. Spaces 2022, 2022, 2754507. [Google Scholar] [CrossRef]
  23. Liu, K.; Yang, Z.; Wei, W.; Gao, B.; Xin, D.; Sun, C.; Wu, G. Novel detection approach for thermal defects: Study on its feasibility and application to vehicle cables. High Volt. 2022, 1–10. [Google Scholar] [CrossRef]
  24. Liu, P.; Shi, J.; Wang, Z.-H. Pattern formation of the attraction-repulsion Keller-Segel system. Discret. Contin. Dyn. Syst.-B 2013, 18, 2597–2625. [Google Scholar] [CrossRef]
  25. Wang, H.; Fu, H.M.; Zhang, H.F.; Hu, Z.Q. A practical thermodynamic method to calculate the best glass-forming composition for bulk metallic glasses. Int. J. Nonlinear Sci. Numer. Simul. 2007, 8, 171–178. [Google Scholar] [CrossRef] [Green Version]
  26. Xu, L.; He, J.H.; Liu, Y. Electrospun nanoporous spheres with Chinese drug. Int. J. Nonlinear Sci. Numer. Simul. 2007, 8, 199–202. [Google Scholar] [CrossRef]
  27. Sun, F.Z.; Gao, M.; Lei, S.H.; Zhao, Y.; Wang, K.; Shi, Y.T.; Wang, N.H. The fractal dimension of the fractal model of dropwise condensation and its experimental study. Int. J. Nonlinear Sci. Numer. Simul. 2007, 8, 211–222. [Google Scholar] [CrossRef]
  28. Bo, T.L.; Xie, L.; Zheng, X.J. Numerical approach to wind ripple in desert. Int. J. Nonlinear Sci. Numer. Simul. 2007, 8, 223–228. [Google Scholar] [CrossRef]
  29. Diethelm, K.; Ford, N.J.; Freed, A.D. Detailed error analysis for a fractional adams method. Numer. Algorithms 2004, 36, 31–52. [Google Scholar] [CrossRef] [Green Version]
  30. Diethelm, K.; Ford, N.J. Analysis of fractional differential equations. J. Math. Anal. Appl. 2002, 265, 229–248. [Google Scholar] [CrossRef] [Green Version]
  31. Hashim, I.; Abdulaziz, O.; Momani, S. Homotopy analysis method for fractional ivps. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 674–684. [Google Scholar] [CrossRef]
  32. Odibat, Z.; Momani, S.; Erturk, V.S. Generalized differential transform method: Application to differential equations of fractional order. Appl. Math. Comput. 2008, 197, 467–477. [Google Scholar] [CrossRef]
  33. Brunner, H.; Makroglou, A.; Miller, R.K. Mixed interpolation collocation methods for first and second order volterra integro-differential equations with periodic solution. Appl. Numer. Math. 1997, 23, 381–402. [Google Scholar] [CrossRef]
  34. Rawashdeh, E.A. Numerical solution of fractional integro-differential equations by collocation method. Appl. Math. Comput. 2006, 176, 1–6. [Google Scholar] [CrossRef]
  35. Volk, W. The iterated galerkin method for linear integro-differential equations. J. Comput. Appl. Math. 1988, 21, 63–74. [Google Scholar] [CrossRef]
  36. Pedas, A.; Tamme, E. Spline collocation method for integro-differential equations with weakly singular kernels. J. Comput. Appl. Math. 2006, 197, 253–269. [Google Scholar] [CrossRef] [Green Version]
  37. Yousefi, S.; Razzaghi, M. Legendre wavelets method for the nonlinear volterra-fredholm integral equations. Math. Comput. Simul. 2005, 70, 1–8. [Google Scholar] [CrossRef]
  38. Atangana, A.; Araz, S.I. New numerical approximation for Chua attractor with fractional and fractal-fractional operators. Alex. Eng. J. 2020, 59, 3275–3296. [Google Scholar] [CrossRef]
  39. Mekkaoui, T.; Atangana, A.; Araz, S.I. Predictor-corrector for non-linear differential and integral equation with fractal-fractional operators. Eng. Comput. 2021, 37, 2359–2368. [Google Scholar] [CrossRef]
  40. Subasi, M.; Araz, S.I.; Gungor, H. On the numerical solution of two dimensional Schrodinger equation. Int. J. Math. Res. 2017, 6, 1–12. [Google Scholar] [CrossRef] [Green Version]
  41. Nonlaopon, K.; Alsharif, A.M.; Zidan, A.M.; Khan, A.; Hamed, Y.S.; Shah, R. Numerical investigation of fractional-order Swift-Hohenberg equations via a Novel transform. Symmetry 2021, 13, 1263. [Google Scholar] [CrossRef]
  42. He, J.H.; Li, Z.B.; Wang, Q.L. A new fractional derivative and its application to explanation of polar bear hairs. J. King Saud-Univ.-Sci. 2016, 28, 190–192. [Google Scholar] [CrossRef] [Green Version]
  43. Sunthrayuth, P.; Ullah, R.; Khan, A.; Shah, R.; Kafle, J.; Mahariq, I.; Jarad, F. Numerical analysis of the fractional-order nonlinear system of Volterra integro-differential equations. J. Funct. Spaces 2021, 2021, 1537958. [Google Scholar] [CrossRef]
  44. Sweilam, N.H.; Khader, M. A Chebyshev pseudo-spectral method for solving fractional-order integro-differential equations. ANZIAM J. 2010, 51, 464–475. [Google Scholar] [CrossRef]
Figure 1. Example 1 of the solution graph, CPM solution and exact solution for ξ ( ψ ) .
Figure 1. Example 1 of the solution graph, CPM solution and exact solution for ξ ( ψ ) .
Symmetry 14 02575 g001
Figure 2. Example 1 of the solution graph, CPM solution and exact solution for ζ ( ψ ) .
Figure 2. Example 1 of the solution graph, CPM solution and exact solution for ζ ( ψ ) .
Symmetry 14 02575 g002
Figure 3. Our method, CWM, and RBFN error graph for Example 1.
Figure 3. Our method, CWM, and RBFN error graph for Example 1.
Symmetry 14 02575 g003
Figure 4. Our method, CWM, and RBFN error graph for Example 1.
Figure 4. Our method, CWM, and RBFN error graph for Example 1.
Symmetry 14 02575 g004
Figure 5. At different fractional orders, the absolute error graph of Example 1.
Figure 5. At different fractional orders, the absolute error graph of Example 1.
Symmetry 14 02575 g005
Figure 6. Example 2: CPM solution and exact solution for ξ ( ψ ) .
Figure 6. Example 2: CPM solution and exact solution for ξ ( ψ ) .
Symmetry 14 02575 g006
Figure 7. Problem 2: CPM solution and exact solution for ζ ( ψ ) .
Figure 7. Problem 2: CPM solution and exact solution for ζ ( ψ ) .
Symmetry 14 02575 g007
Figure 8. Our method and the CWM error graph for Example 2.
Figure 8. Our method and the CWM error graph for Example 2.
Symmetry 14 02575 g008
Figure 9. Our method and the CWM error graph for Example 2.
Figure 9. Our method and the CWM error graph for Example 2.
Symmetry 14 02575 g009
Figure 10. At different fractional orders, the absolute error graph of Example 2.
Figure 10. At different fractional orders, the absolute error graph of Example 2.
Symmetry 14 02575 g010
Figure 11. At different fractional orders, the absolute error graph of Example 2.
Figure 11. At different fractional orders, the absolute error graph of Example 2.
Symmetry 14 02575 g011
Figure 12. Example 3: CPM solution and exact solution for ξ ( ψ ) .
Figure 12. Example 3: CPM solution and exact solution for ξ ( ψ ) .
Symmetry 14 02575 g012
Figure 13. Example 3: CPM solution and exact solution for ϕ ( ψ ) .
Figure 13. Example 3: CPM solution and exact solution for ϕ ( ψ ) .
Symmetry 14 02575 g013
Figure 14. Our method and the CWM error graph for Example 3.
Figure 14. Our method and the CWM error graph for Example 3.
Symmetry 14 02575 g014
Figure 15. Our method and the CWM error graph for Example 3.
Figure 15. Our method and the CWM error graph for Example 3.
Symmetry 14 02575 g015
Figure 16. Our method and the CWM error graph for Example 3.
Figure 16. Our method and the CWM error graph for Example 3.
Symmetry 14 02575 g016
Figure 17. At different fractional orders, the absolute error graph of Example 3.
Figure 17. At different fractional orders, the absolute error graph of Example 3.
Symmetry 14 02575 g017
Figure 18. At different fractional orders, the absolute error graph of problem 3.
Figure 18. At different fractional orders, the absolute error graph of problem 3.
Symmetry 14 02575 g018
Figure 19. At different fractional orders, the absolute error graph of Example 3.
Figure 19. At different fractional orders, the absolute error graph of Example 3.
Symmetry 14 02575 g019
Table 1. Example 1: exact vs. CPM, CWM solution at m = 6.
Table 1. Example 1: exact vs. CPM, CWM solution at m = 6.
ψ E.S ξ ( ψ ) E.S ζ ( ψ ) CPM ξ ( ψ ) CPM ζ ( ψ ) CWM ξ ( ψ ) CWM ζ ( ψ )
00.00000000.00000000.00000000000.0000000000−0.0000260700.000222296
0.20.04000000.20000000.03999999990.20000000000.0399870350.199964603
0.40.16000000.40000000.15999999990.40000000000.1599826150.400041558
0.60.36000000.60000000.35999999990.60000000000.3600478970.599977336
0.80.64000000.80000000.64000000000.79999999990.6399664610.799997734
1.01.00000001.00000001.00000000000.99999999991.0001858560.999922349
Table 2. CPM versus the CWM error comparison at m = 6 of problem 1.
Table 2. CPM versus the CWM error comparison at m = 6 of problem 1.
ψ Error ( ξ CPM ) Error ( ζ CPM ) Error ( ξ CWM ) Error ( ζ CWM ) Error ( ξ RBFN ) Error ( ζ RBFN )
00.00000 × 10 + 00 0.00000 × 10 + 00 0.000026070.00022227.54 × 10 05 2.09 × 10 05
0.21.60000 × 10 15 1.00000 × 10 15 0.000012960.00003532.28 × 10 05 5.70 × 10 06
0.43.00000 × 10 15 1.00000 × 10 15 0.000017380.00004151.98 × 10 06 4.79 × 10 06
0.63.00000 × 10 15 0.00000 × 10 + 00 0.000047890.00002261.16 × 10 05 6.33 × 10 06
0.84.00000 × 10 15 2.00000 × 10 15 0.000033530.00000223.80 × 10 05 1.29 × 10 05
1.03.00000 × 10 14 1.30000 × 10 14 0.000185850.00007765.82 × 10 05 3.89 × 10 05
Table 3. Problem 2: exact vs. CPM; CWM solution at m = 6.
Table 3. Problem 2: exact vs. CPM; CWM solution at m = 6.
ψ E.S ξ ( ψ ) E.S ζ ( ψ ) CPM ξ ( ψ ) CPM ζ ( ψ ) CWM ξ ( ψ ) CWM ζ ( ψ )
00.000000000001.00000000000.000000000000.99999999999−0.0000014031.000000319
0.20.198669330790.98006657780.198669330720.980066577860.19866963930.9800669460
0.40.389418342300.92106099400.389418342020.921060994060.38941657870.9210604974
0.60.564642473390.82533561490.564642473400.825335614820.56464339640.8253349025
0.80.717356090890.69670670930.717356090790.696706709280.71735576410.6967067839
1.00.841470984800.54030230580.841470980760.540302307570.84147198420.5403020838
Table 4. CPM versus the CWM error comparison at m = 6 of problem 2.
Table 4. CPM versus the CWM error comparison at m = 6 of problem 2.
s Error ( ξ CPM ) Error ( ζ CPM ) Error ( ξ CWM ) Error ( ζ CWM )
00.0000000000 × 10 + 00 0.0000000000 × 10 + 00 0.00000140370.000000319
0.26.7443073256 × 10 11 1.9129717529 × 10 11 0.00000030850.0000003682
0.42.8634095988 × 10 10 6.3350567575 × 10 11 0.00000017630.0000004966
0.61.4540273872 × 10 11 8.3812136477 × 10 11 0.0000009230.0000007124
0.81.0840388114 × 10 10 6.5764647225 × 10 11 0.00000032680.0000000746
1.04.0385798791 × 10 9 1.7087357128 × 10 9 0.00000099940.0000002221
Table 5. Problem 3: exact vs. CPM solution at m = 6.
Table 5. Problem 3: exact vs. CPM solution at m = 6.
ψ E.S ξ ( ψ ) E.S ζ ( ψ ) E.S ϕ ( ψ ) CPM ξ ( ψ ) CPM ζ ( ψ ) CPM ϕ ( ψ )
01.0000000000001.500000000.66666666661.0000000001.500000000.6666666666
0.11.0025000000001.503333330.71666666661.0025000001.503333330.7166666666
0.21.0100000000001.513333330.76666666661.0099999991.513333330.7666666666
0.31.0225000000001.530000000.81666666661.0225000001.530000000.8166666666
0.41.0400000000001.553333330.86666666661.0400000001.553333330.8666666666
0.51.0625000000001.583333330.91666666661.0624999991.583333330.9166666666
0.61.0900000000001.620000000.96666666661.0900000001.620000000.9666666666
0.71.1225000000001.663333331.01666666661.1225000001.663333331.0166666666
0.81.1600000000001.713333330.06666666661.1600000001.713333331.0666666666
0.91.2025000000001.770000001.11666666661.2024999991.770000001.1166666666
1.01.2500000000001.833333331.16666666661.2499999991.833333331.1666666666
Table 6. CPM versus the CWM error comparison at m = 6 of problem 3.
Table 6. CPM versus the CWM error comparison at m = 6 of problem 3.
ψ ( ξ CPM ) ( ζ CPM ) ( ϕ CPM ) ( ξ CWM ) ( ζ CWM ) ( ϕ CWM )
00.00000 × 10 + 00 1.00000 × 10 14 3.100000 × 10 15 0.0000011580.0000006420.0000097282
0.21.00000 × 10 14 0.00000 × 10 + 00 4.100000 × 10 15 0.0000009590.0000004600.0000008219
0.41.00000 × 10 14 1.00000 × 10 14 2.100000 × 10 15 0.0000000240.0000000790.0000030059
0.60.00000 × 10 + 00 2.00000 × 10 14 3.100000 × 10 15 0.0000023010.0000010720.0000071362
0.84.00000 × 10 14 0.00000 × 10 + 00 2.100000 × 10 14 0.0000044580.0000024750.0000081324
1.06.40000 × 10 13 5.20000 × 10 13 8.100000 × 10 14 0.0000229450.000013260.0000439496
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sunthrayuth, P.; Naeem, M.; Shah, N.A.; Shah, R.; Chung, J.D. Numerical Simulations of the Fractional Systems of Volterra Integral Equations within the Chebyshev Pseudo-Spectral Method. Symmetry 2022, 14, 2575. https://doi.org/10.3390/sym14122575

AMA Style

Sunthrayuth P, Naeem M, Shah NA, Shah R, Chung JD. Numerical Simulations of the Fractional Systems of Volterra Integral Equations within the Chebyshev Pseudo-Spectral Method. Symmetry. 2022; 14(12):2575. https://doi.org/10.3390/sym14122575

Chicago/Turabian Style

Sunthrayuth, Pongsakorn, Muhammad Naeem, Nehad Ali Shah, Rasool Shah, and Jae Dong Chung. 2022. "Numerical Simulations of the Fractional Systems of Volterra Integral Equations within the Chebyshev Pseudo-Spectral Method" Symmetry 14, no. 12: 2575. https://doi.org/10.3390/sym14122575

APA Style

Sunthrayuth, P., Naeem, M., Shah, N. A., Shah, R., & Chung, J. D. (2022). Numerical Simulations of the Fractional Systems of Volterra Integral Equations within the Chebyshev Pseudo-Spectral Method. Symmetry, 14(12), 2575. https://doi.org/10.3390/sym14122575

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop