/ / Article

Research Article | Open Access

Volume 2020 |Article ID 1010382 | https://doi.org/10.1155/2020/1010382

Ayed Al e’damat, "Approximate Analytical Solution for a Coupled System of Fractional Nonlinear Integrodifferential Equations by the RPS Method", Complexity, vol. 2020, Article ID 1010382, 8 pages, 2020. https://doi.org/10.1155/2020/1010382

# Approximate Analytical Solution for a Coupled System of Fractional Nonlinear Integrodifferential Equations by the RPS Method

Accepted18 Jul 2020
Published13 Aug 2020

#### Abstract

In this work, a modified residual power series method is implemented for providing efficient analytical and approximate solutions for a class of coupled system of nonlinear fractional integrodifferential equations. The proposed algorithm is based on the concept of residual error functions and generalized power series formula. The fractional derivative is described under the Caputo concept. To illustrate the potential, accuracy, and efficiency of the proposed method, two numerical applications of the coupled system of nonlinear fractional integrodifferential equations are tested. The numerical results confirm the theoretical predictions and depict that the suggested scheme is highly convenient, is quite effective, and practically simplifies computational time. Consequently, the proposed method is simple, accurate, and convenient in handling different types of fractional models arising in the engineering and physical systems.

#### 1. Introduction

The theory of fractional calculus is indeed a generalization of the standard calculus that deals with differentiation and integration of a noninteger order, which is utilized to describe various real-world phenomena arising in natural sciences, applied mathematics, and engineering fields with great applications for these tools, including nonlinear oscillation of earthquakes, fractional fluid-dynamic traffic, economics, solid mechanics, viscoelasticity, and control theory . The major cause behind this is that modeling of a specific phenomenon does not depend only at the time instant but also the historical state, so the differential and integral operators for integer and fractional cases are found to be a superb tool in describing the hereditary and memory properties for different engineering and physical phenomena . However, several mathematical forms of the abovementioned issues contain nonlinear fractional integrodifferential equations (FIDEs), and other nonlinear models can be found in . Since most fractional differential and integrodifferential equations cannot be solved analytically, it is necessary to find an accurate numerical and analytical method to deal with the complexity of fractional operators involving such equations .

This paper aims to introduce a recent analytical as well as numerical method based on the use of fractional residual power series (RPS) technique for obtaining the approximate solution for a class of coupled system of fractional integrodifferential equations in the following form:

This is subject to the following initial conditions:where .

The residual power series method (RPSM) has a wide range of applications, especially in simulating nonlinear issues in a fractional meaning, which has been developed and modified over recent years as a powerful mathematical treatment indispensable in dealing with the emerging realistic system in physics, engineering, and natural sciences . More specifically, it is a modern analytical and approximation technique that relies on the expansion of the fractional power series and residual error functions, which was first proposed in 2013 to provide analytical series solutions to fuzzy differential equations of the first and second orders and minimize the residual errors. This method has many advantages and properties as follows: it is an accurate alternative instrument, it requires less effort to achieve results, it provides a rapid convergence rate to the exact solution, it deals directly with different types of nonlinear terms and complex functions, and it has the ability to choose any point in the integration domain, making the approximate solution applicable. It has also excellent estimating characteristics that reflect high reliability. Furthermore, it is a systematic tool to find sequential solutions for several types of nonlinear differential equations and integrodifferential equations of fractional order without having to discretize, linearize, or even expose to perturbation. Therefore, it attracted the attention of many researchers.

Freihet et al.  have used the fractional power series for solving the fractional stiff system and introduced some basic theorems related to RPS generalization in the sense of Caputo fractional derivative. The (2 + 1)-dimensional time-fractional Burgers–Kadomtsev–Petviashvili equation has been solved by the RPS method . In , analytic-approximate solutions for nonlinear coupled fractional Jaulent–Miodek equations with energy-dependent Schrödinger potential have been obtained using the RPS and q-homotopy analysis methods. Moreover, this method was successfully applied for solving both linear and nonlinear ordinary, partial, and fuzzy differential equations . Therefore, such adaptives can be used as an alternative technique in solving several nonlinear problems arising in engineering and sciences.

The outline of this paper is organized as follows: In the next section, we review some basic definitions and theories related to fractional differentiation and fractional power series representations. In Section 3, the solution by the RPS technique is provided. In Section 4, numerical application is performed to show accuracy and efficiency of the RPS method. Finally, we give concluding remarks in Section 5.

#### 2. Basic Mathematical Concepts

In this section, basic definitions and results related to fractional calculus are given and fractional power series concept is also represented.

Definition 1. The Riemann–Liouville fractional integral operator of order , over the interval for a function , is defined byFor and , the operator has the following basic properties:(1)(2)(3)

Definition 2. For . The Caputo fractional derivative of order is given byFor In case and then .
The following are some interesting properties of the operator :(1)For any constant , then .(2)(3)(4)

Definition 3. A fractional power series (FPS) representation at has the following form:where and and is the coefficient of the series.

Theorem 1. Suppose that has the following FPS representation at :where , and for . Then, the coefficients will be in the form such that (-times).

#### 3. Fractional RPS Method for the Coupled System of IDEs

The purpose of this section is to construct FPS solution for the coupled system of nonlinear fractional integrodifferential equations (1) and (2) by substituting the FPS expansion among the truncated residual functions.

The RPS algorithm proposing the solution of equations (1) and (2) about gives the following FPS expansion:

The truncated series form of equation (7) can be given by the following th-FPS approximate solution:

Clearly, if , then expansion (8) can be written as

Define the residual function for equations (1) and (2) as follows:

Also, define the th-residual function as

According to the RPS algorithm , we have the following relations:(i), for each (ii), for each

For obtaining the coefficients , , one can solve the solution of the following relation:

Algorithm 1. To find the coefficients in equation (9), perform the following steps:(i)Step 1: substitute expansion (9) function , into -th residual function (11) such that (ii)Step 2: find the relation of fractional formula of at (iii)Step 3: for obtain the relation through the fact . For obtain the relation through the fact . For obtain the relation through the fact , obtain the relation through the fact (iv)Step 4: solve the obtained algebraic fractional system (v)Step 5: substitute the values of back into equation (8) and then stop.

#### 4. Numerical Applications and Simulation

This section aims to test two applications of the system of nonlinear fractional IDEs to show the efficiency, accuracy, and applicability of the proposed method. In this section, all calculations are preformed using Wolfram-Mathematica 10.

Example 1. Consider the following nonlinear fractional integrodifferential equation:This is subject to the following initial conditions:The exact solution of this coupled system is and .
Using the RPS algorithm, the -th residual functions are given bywhere has the following form:Consequently,Numerical simulation specializes in advanced numeric or approximate methods for finding digital or approximate solutions and estimating errors of these approximations. For such purpose, a few decimal numbers are often recorded to calculate the absolute error, which produces some vague estimates related to significance and units, as well as does not provide clear and explicit evidence regarding the subject matter of the study. Therefore, the comparison of absolute error with the exact value leads to the determination of the relative error as a ratio between the value of absolute error and the exact, which gives some importance and reduces ambiguity for a deeper understanding of the behavior of the approximate solutions. By using the RPS method of Example 1, the numerical results of and are shown in Tables 1 and 2 at and The results obtained in Tables 1 and 2 show that the error estimate using the proposed method is very small and that the solutions correspond well to each other. In general, it should be noted that increasing the number of iterations will lead to an improvement in numerical solutions and approaching the exact value.
Tables 3 and 4 show the sixth approximate solutions of Example 1 at different values of such that with step size and . From these tables, it can be concluded that the RPS algorithm and the approximate solutions are consistent with each other and with the exact solutions for all values of in . Here, it is worth noting that the closer the value of the fractional derivative approaching the integer case , the closer the approximate solution is to the exact solution.

 0 1.0 1.0 0.0 0.1 1.1051709 1.1051708 0.2 1.2214028 1.2214000 0.3 1.3498588 1.3498375 0.4 1.4918247 1.4917333 0.5 1.6487213 1.6484375 0.6 1.8221188 1.8214000 0.7 2.0137527 2.0121708 0.8 2.2255409 2.2224000 0.9 2.4596031 2.4538375 1.0 2.7182818 2.7083333
 0 1.0 1.0 0.0 0.1 0.9048374 0.9048375 0.2 0.8187308 0.8187333 0.3 0.7408182 0.7408375 0.4 0.67033201 0.6704000 0.5 0.6065307 0.6067708 0.6 0.5488116 0.5494000 0.7 0.4965853 0.4978375 0.8 0.4493290 0.4517333 0.9 0.4065697 0.4108375 1.0 0.3678794 0.3750000
 0.0 1.0 1.0 1.0 1.0 0.1 1.0510833 1.0733961 1.1073884 1.1601510 0.2 1.1246667 1.1716823 1.2392129 1.3365098 0.3 1.2212500 1.2956381 1.3974229 1.5359285 0.4 1.3413333 1.4446154 1.5801578 1.7560644 0.5 1.4854167 1.6181748 1.7860776 1.9951767 0.6 1.6540000 1.8160345 2.0142021 2.2519736 0.7 1.8475833 2.0380217 2.2637874 2.5254605 0.8 2.0666667 2.2840401 2.5342526 2.8148458 0.9 2.3117500 2.5540495 2.8251333 3.1194830 1.0 2.5833333 2.8480518 3.1360514 3.4388329
 0.0 1.0 1.0 1.0 1.0 0.1 0.9048333 0.8780779 0.8460694 0.8088467 0.2 0.8186666 0.7855478 0.7499395 0.7123611 0.3 0.7404999 0.7071736 0.6734520 0.6391004 0.4 0.6693333 0.6386210 0.6084174 0.5772768 0.5 0.6041666 0.5772055 0.5505285 0.5214466 0.6 0.5439999 0.5209163 0.4969643 0.4684158 0.7 0.4878333 0.4681027 0.4456563 0.4160595 0.8 0.4346666 0.4173375 0.3949800 0.3628482 0.9 0.3835000 0.3673449 0.3435991 0.3076204 1.0 0.3333333 0.3169592 0.2903784 0.2494581

Example 2. Consider the following fractional integrodifferential equation:This is subject to the following initial conditions:The exact solution of this coupled system is and . Using the RPS algorithm, the -th residual functions and are given bywhere and have the following form:The absolute errors are listed in Tables 5 and 6. The results obtained by the RPS method show that the exact solutions are in good agreement with approximate solutions at , , and step size . Tables 7 and 8 show approximate solutions at different values of such that and with step size . From these tables, one can find that the RPS method provides us with an accurate approximate solution, which is in good agreement with the exact solutions for all values of in . Also, it is worth noting that the closer the value of the fractional derivative approaching the integer case , the closer the approximate solution is to the exact solution.

 0.1 0.0998334 0.0998257 0.2 0.1986693 0.1986624 0.3 0.2955202 0.2955216 0.4 0.3894183 0.3894152 0.5 0.4794255 0.4794232 0.6 0.5646425 0.5646492 0.7 0.6442177 0.6442199 0.8 0.7173561 0.7173533 0.9 0.7833269 0.7833286 1.0 0.8414710 0.8414475
 0.1 0.9950042 0.9950103 0.2 0.9800666 0.9800597 0.3 0.9553365 0.9553270 0.4 0.9210610 0.9210410 0.5 0.8775826 0.8775771 0.6 0.8253356 0.8253269 0.7 0.7648422 0.7645622 0.8 0.6967067 0.6966984 0.9 0.6216100 0.6215986 1.0 0.5403023 0.5402669
 0.0 0.0 0.0 0.0 0.1 0.0749501 0.1027865 0.1422048 0.2 0.1553618 0.2021055 0.2644892 0.3 0.2451462 0.3086629 0.3897369 0.4 0.3445577 0.4231060 0.5195960 0.5 0.4535802 0.5454302 0.6544017 0.6 0.5721434 0.6754941 0.7941689 0.7 0.7001695 0.8131291 0.9388129 0.8 0.8375854 0.9581693 1.0882176 0.9 0.9843269 1.1104609 1.2422602 1.0 1.1403382 1.2698635 1.4008207
 0.0 1.0 1.0 1.0 0.1 1.0036839 0.9995798 0.9902709 0.2 0.9918177 0.9770693 0.9526430 0.3 0.9678131 0.9409453 0.9018037 0.4 0.9329689 0.8938412 0.8413882 0.5 0.8881139 0.8372777 0.7733064 0.6 0.8338592 0.7722922 0.6987749 0.7 0.7706901 0.6996578 0.6186510 0.8 0.6990108 0.6199832 0.5335787 0.9 0.6191685 0.5337664 0.4440635 1.0 0.5314684 0.4414262 0.3505153

#### 5. Concluding Remarks

In this work, a class of a coupled system of nonlinear fractional integrodifferential equations of fractional order has been discussed by using the RPS method under the Caputo fractional derivative. The RPS algorithm has been given to optimize the approximate solution by minimizing a residual error function with the help of generalized Taylor formula. To demonstrate the consistency with the theoretical framework, two illustrative examples have been provided. The obtained results indicated that the approximate solutions are coinciding with the exact solution and with each other for different values of the fractional order over the selected nods and parameters. From our results, we can conclude that the RPS method is a systematic and suitable scheme to address many fractional initial value problems with great potential in scientific applications. The calculations have been performed by using Wolfram-Mathematica 10.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The author declares that there are no conflicts of interest.

1. Y. A. Rossikhin and M. V. Shitikova, “Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids,” Applied Mechanics Reviews, vol. 50, no. 1, pp. 15–67, 1997. View at: Publisher Site | Google Scholar
2. R. T. Baillie, “Long memory processes and fractional integration in econometrics,” Journal of Econometrics, vol. 73, no. 1, pp. 5–59, 1996. View at: Publisher Site | Google Scholar
3. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley and Sons, Hoboken, NY, USA, 1993.
4. X.-J. Yang, General Fractional Derivatives: Theory, Methods and Applications, CRC Press, New York, NY, USA, 2019.
5. X.-J. Yang, F. Gao, and Y. Ju, General Fractional Derivatives with Applications in Viscoelasticity, Academic Press, Cambridge, MA, USA, 2020.
6. X. J. J. A. T. Yang and J. A. Tenreiro Machado, “A new fractal nonlinear Burgers’ equation arising in the acoustic signals propagation,” Mathematical Methods in the Applied Sciences, vol. 42, no. 18, pp. 7539–7544, 2019. View at: Publisher Site | Google Scholar
7. L. Akinyemi and O. S. Iyiola, “A reliable technique to study nonlinear time-fractional coupled Korteweg-de Vries equations,” Advances in Difference Equations, vol. 2020, no. 1, Article ID 169, 27 pages, 2020. View at: Publisher Site | Google Scholar
8. M. Al-Smadi, “Simplified iterative reproducing kernel method for handling time-fractional BVPs with error estimation,” Ain Shams Engineering Journal, vol. 9, no. 4, pp. 2517–2525, 2018. View at: Publisher Site | Google Scholar
9. S. Momani, A. Freihat, and M. Al-Smadi, “Analytical study of fractional-order multiple chaotic Fitzhugh-Nagumo neurons model using multistep generalized differential transform method,” Abstract and Applied Analysis, vol. 2014, Article ID 276279, 10 pages, 2014. View at: Publisher Site | Google Scholar
10. Z. Altawallbeh, M. Al-Smadi, I. Komashynska, and A. Ateiwi, “Numerical solutions of fractional systems of two-point BVPs by using the iterative reproducing kernel algorithm,” Ukrainian Mathematical Journal, vol. 70, no. 5, pp. 687–701, 2018. View at: Publisher Site | Google Scholar
11. O. Abu Arqub and M. Al-Smadi, “An adaptive numerical approach for the solutions of fractional advection-diffusion and dispersion equations in singular case under Riesz’s derivative operator,” Physica A: Statistical Mechanics and Its Applications, vol. 540, Article ID 123257, 13 pages, 2020. View at: Publisher Site | Google Scholar
12. M. Al-Smadi and O. A. Arqub, “Computational algorithm for solving fredholm time-fractional partial integrodifferential equations of dirichlet functions type with error estimates,” Applied Mathematics and Computation, vol. 342, pp. 280–294, 2019. View at: Publisher Site | Google Scholar
13. M. Al-Smadi, “Solving fractional system of partial differential equations with parameters derivative by combining the GDTM and RDTM,” Nonlinear Studies, vol. 26, no. 3, pp. 587–601, 2019. View at: Google Scholar
14. M. Al-Smadi, O. Abu Arqub, and S. Momani, “A computational method for two-point boundary value problems of fourth-order mixed integrodifferential equations,” Mathematical Problems in Engineering, vol. 2013, Article ID 832074, 10 pages, 2013. View at: Publisher Site | Google Scholar
15. D. Zeidan, A. Slaouti, E. Romenski, and E. F. Toro, “Numerical solution for hyperbolic conservative two-phase flow equations,” International Journal of Computational Methods, vol. 4, no. 2, pp. 299–333, 2017. View at: Publisher Site | Google Scholar
16. M. Al-Smadi, “Reliable Numerical algorithm for handling fuzzy integral equations of second kind in Hilbert spaces,” Filomat, vol. 33, no. 2, pp. 583–597, 2019. View at: Publisher Site | Google Scholar
17. G. N. Gumah, M. F. M. Naser, M. Al-Smadi, and S. K. Al-Omari, “Application of reproducing kernel Hilbert space method for solving second-order fuzzy Volterra integro-differential equations,” Advances in Difference Equations, vol. 2018, no. 1, Article ID 457, 15 pages, 2018. View at: Publisher Site | Google Scholar
18. M. Al-Smadi, A. Freihat, O. Abu Arqub, and N. Shawagfeh, “A novel multistep generalized differential transform method for solving fractional-order Lü chaotic and hyperchaotic systems,” Journal of Computational Analysis and Applications, vol. 19, no. 4, pp. 713–724, 2015. View at: Google Scholar
19. M. Al-Smadi, O. A. Arqub, N. Shawagfeh, and S. Momani, “Numerical investigations for systems of second-order periodic boundary value problems using reproducing kernel method,” Applied Mathematics and Computation, vol. 291, pp. 137–148, 2016. View at: Publisher Site | Google Scholar
20. E. Goncalvès and D. Zeidan, “Numerical simulation of unsteady cavitation in liquid hydrogen flows,” International Journal of Engineering Systems Modelling and Simulation, vol. 9, no. 1, pp. 41–51, 2017. View at: Publisher Site | Google Scholar
21. G. Gumah, M. F. M. Naser, M. Al-Smadi, S. K. Q. Al-Omari, and D. Baleanu, “Numerical solutions of hybrid fuzzy differential equations in a hilbert space,” Applied Numerical Mathematics, vol. 151, pp. 402–412, 2020. View at: Publisher Site | Google Scholar
22. M. Al-Smadi, A. Freihat, H. Khalil, S. Momani, and R. A. Khan, “Numerical multistep approach for solving fractional partial differential equations,” International Journal of Computational Methods, vol. 14, no. 2, Article ID 1750029, 15 pages, 2017. View at: Publisher Site | Google Scholar
23. M. H. AL-Smadi and G. N. Gumah, “On the homotopy analysis method for fractional SEIR epidemic model,” Research Journal of Applied Sciences, Engineering and Technology, vol. 7, no. 18, pp. 3809–3820, 2014. View at: Publisher Site | Google Scholar
24. I. Komashynska and M. Al-Smadi, “Iterative reproducing kernel method for solving second-order integrodifferential equations of fredholm type,” Journal of Applied Mathematics, vol. 2014, Article ID 459509, 11 pages, 2014. View at: Publisher Site | Google Scholar
25. I. Komashynska, M. Al-Smadi, O. Abu Arqub, and S. Momani, “An efficient analytical method for solving singular initial value problems of nonlinear systems,” Applied Mathematics & Information Sciences, vol. 10, no. 2, pp. 647–656, 2016. View at: Publisher Site | Google Scholar
26. A. El-Ajou, M. Al-Smadi, M. a. N. Oqielat, S. Momani, and S. Hadid, “Smooth expansion to solve high-order linear conformable fractional PDEs via residual power series method: applications to physical and engineering equations,” Ain Shams Engineering Journal, 2020, In Press. View at: Publisher Site | Google Scholar
27. I. Komashynska, M. Al-Smadi, A. Ateiwi, and S. Al-Obaidy, “Approximate analytical solution by residual power series method for system of Fredholm integral equations,” Applied Mathematics & Information Sciences, vol. 10, no. 3, pp. 975–985, 2016. View at: Publisher Site | Google Scholar
28. M. Al-Smadi, O. Abu Arqub, and S. Momani, “Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense,” Physica Scripta, vol. 95, no. 7, Article ID 075218, 2020. View at: Publisher Site | Google Scholar
29. A. Freihet, S. Hasan, M. Al-Smadi, M. Gaith, and S. Momani, “Construction of fractional power series solutions to fractional stiff system using residual functions algorithm,” Advances in Difference Equations, vol. 95, 2019. View at: Publisher Site | Google Scholar
30. M. Senol, “Analytical and approximate solutions of (2 + 1)-dimensional time-fractional Burgers-Kadomtsev-Petviashvili equation,” Communications in Theoretical Physics, vol. 72, no. 5, Article ID 055003, 2020. View at: Publisher Site | Google Scholar
31. M. Şenol, O. S. Iyiola, H. D. Kasmaei, and L. Akinyemi, “Efficient analytical techniques for solving time-fractional nonlinear coupled Jaulent–Miodek system with energy-dependent Schrödinger potential,” Advances in Difference Equations, vol. 2019, Article ID 459509, 11 pages, 2019. View at: Publisher Site | Google Scholar
32. M. Al-Smadi, O. Abu Arqub, and S. Hadid, “An attractive analytical technique for coupled system of fractional partial differential equations in shallow water waves with conformable derivative,” Communications in Theoretical Physics, vol. 72, no. 8, Article ID 085001, 2020. View at: Publisher Site | Google Scholar
33. S. Hasan, M. Al-Smadi, A. Freihet, and S. Momani, “Two computational approaches for solving a fractional obstacle system in Hilbert space,” Advances in Difference Equations, vol. 2019, no. 1, Article ID 55, 17 pages, 2019. View at: Publisher Site | Google Scholar
34. R. Saadeh, M. Alaroud, M. Al-Smadi, R. Ahmad, and U. Salma Din, “Application of fractional residual power series algorithm to solve newell-whitehead-segel equation of fractional order,” Symmetry, vol. 11, no. 12, p. 1431, 2019. View at: Publisher Site | Google Scholar
35. A. Freihet, H. Shatha, M. Alaroud, M. Al-Smadi, R. R. Ahmad, and U. K. Salma Din, “Toward computational algorithm for time-fractional Fokker–Planck models,” Advances in Mechanical Engineering, vol. 11, no. 10, pp. 1–11, 2019. View at: Publisher Site | Google Scholar
36. S. Alshammari, M. Al-Smadi, M. Al Shammari, I. Hashim, and M. A. Alias, “Advanced analytical treatment of fractional logistic equations based on residual error functions,” International Journal of Differential Equations, vol. 2019, Article ID 7609879, pp. 1–11, 2019. View at: Publisher Site | Google Scholar
37. S. Alshammari, M. Al-Smadi, I. Hashim, and M. A. Alias, “Residual power series technique for simulating fractional Bagley-Torvik problems emerging in applied physics,” Applied Sciences, vol. 9, no. 23, p. 5029, 2019. View at: Publisher Site | Google Scholar