Solving nonlinear Volterra integro-differential equations of fractional order by using Euler wavelet method

Advances in Difference Equations, Jan 2017

In this paper, a wavelet numerical method for solving nonlinear Volterra integro-differential equations of fractional order is presented. The method is based upon Euler wavelet approximations. The Euler wavelet is first presented and an operational matrix of fractional-order integration is derived. By using the operational matrix, the nonlinear fractional integro-differential equations are reduced to a system of algebraic equations which is solved through known numerical algorithms. Also, various types of solutions, with smooth, non-smooth, and even singular behavior have been considered. Illustrative examples are included to demonstrate the validity and applicability of the technique.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://www.advancesindifferenceequations.com/content/pdf/s13662-017-1085-6.pdf

Solving nonlinear Volterra integro-differential equations of fractional order by using Euler wavelet method

Wang and Zhu Advances in Difference Equations Solving nonlinear Volterra integro-differential equations of fractional order by using Euler wavelet method Yanxin Wang 0 Li Zhu 0 1 0 School of Science, Ningbo University of Technology , Ningbo, 315211 , China 1 School of Applied Mathematics, Xiamen University of Technology , Xiamen, 361024 , China In this paper, a wavelet numerical method for solving nonlinear Volterra integro-differential equations of fractional order is presented. The method is based upon Euler wavelet approximations. The Euler wavelet is first presented and an operational matrix of fractional-order integration is derived. By using the operational matrix, the nonlinear fractional integro-differential equations are reduced to a system of algebraic equations which is solved through known numerical algorithms. Also, various types of solutions, with smooth, non-smooth, and even singular behavior have been considered. Illustrative examples are included to demonstrate the validity and applicability of the technique. Volterra integro-differential equations; Euler wavelet; operational matrix; Caputo derivative; numerical solution 1 Introduction The fractional calculus is a mathematical discipline that is  years old, and it has developed progressively up to now. The concept of differentiation to fractional order was defined in the th century by Riemann and Liouville. In various problems of physics, mechanics, and engineering, fractional differential equations and fractional integral equations have been proved to be a valuable tool in modeling many phenomena [, ]. However, most fractional-order equations do not have analytic solutions. Therefore, there has been significant interest in developing numerical schemes for the solutions of fractional-order differential equations. In the past  years, the theory and applications of the fractional-order partial differential equations (FPDEs) have become of increasing interest for the researchers to generalize the integer-order differential equations. Conventionally various technologies, e.g. modified homotopy analysis transform method (MHATM) [], modified homotopy analysis Laplace transform method [], homotopy analysis transform method (HATM) [, ], fractional homotopy analysis transform method (FHATM) [], local fractional variational iteration algorithms [] were used for the solutions of the FPDEs. Meanwhile, local fractional similarity solution for the diffusion equation was discussed in []. The inverse problems for the fractal steady heat transfer described by the local fractional Volterra integrodifferential equations were considered in []. © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Recently, many effective methods for obtaining approximations or numerical solutions of fractional-order integro-differential equations have been presented. These methods include the variational iteration method [–], the adomian decomposition method [], the fractional differential transform method [], the reproducing kernel method [], the collocation method [, ], and the wavelet method [–]. Wavelet theory is a relatively new and an emerging area in the field of applied science and engineering. Wavelets permit the accurate representation of a variety of functions and operators. Moreover, wavelets establish a connection with fast numerical algorithms []. So the wavelet method is a new numerical method for solving the fractional equations and it needs a small amount of calculation. However, the method will produce a singularity in the case of certain increased resolutions. Using wavelet numerical method has several advantages: (a) the main advantage is that after discretizing the coefficient matrix of the algebraic equation shows sparsity; (b) the wavelet method is computer oriented, thus solving a higher-order equation becomes a matter of dimension increasing; (c) the solution is a multi-resolution type; (d) the solution is convergent, even the size of the increment may be large []. Many researchers started using various wavelets for analyzing problems of high computational complexity. It is proved that wavelets are powerful tools to explore new directions in solving differential equations and integral equations. In this paper, the main purpose is to introduce the Euler wavelet operational matrix method to solve the nonlinear Volterra integro-differential equations of fractional order. The Euler wavelet is first presented and it is constructed by Euler polynomials. The method is based on reducing the equation to a system of algebraic equations by expanding the solution as Euler wavelet with unknown coefficients. The characteristic of the operational method is to transform the integro-differential equations into the algebraic one. It not only simplifies the problem but also speeds up the computation. It is worth noting that the Euler polynomials are not based on orthogonal functions, nevertheless, they possess the operational matrix of integration. Also the Euler wavelet is superior to the Legendre wavelet and the Chebyshev wavelet for approximating an arbitrary function, which can be verified by numerical examples. The structure of this paper is as follows: In Section , we recall some basic definitions and properties of the fractional calculus theory. In Section , the Euler wavelets are constructed and the operational matrix of the fractional integration is derived. In Section , we summarize the application of the Euler wavelet operational matrix method to the solution of the fractional integro-differential equations. Some numerical examples are provided to clarify the approach in Section . The conclusion is given in Section . 2 Fractional calculus There are various definitions of fractional integration and derivatives. The widely used definition of a fractional integration is the Riemann-Liouville definition and the definition of a fractional derivative is the Caputo definition. Definition  The Rieman-Liouville fractional integral operator Itα of order α is defined as [] (α) t(t – τ )α–f (τ ) dτ , α > , t > , f (t), α = . For the Riemann-Liouville fractional integral we have Definition  The Caputo definition of fractional differential operator is given by  t f (n)(τ ) (n – α)  (t – τ )α+–n dτ , n –  < α ≤ n, n ∈ N , For the Caputo derivative we have the following two basic properties: If we restrict the parameters a and b to discrete values as a = a–k, b = nba–k, a > , b > , we have the following family of discrete wavelets: where ψkn form a wavelet basis for L(R). In particular, when a =  and b =  then ψkn(t) form an orthonormal basis. The Euler wavelet ψnm(t) = ψ (k, n, m, t) involves four arguments, n = , . . . , k–, k is assumed to be any positive integer, m is the degree of the Euler polynomials, and t is the normalized time. They are defined on the interval [, ) as 3 Euler wavelet operational matrix of the fractional integration 3.1 Wavelets and Euler wavelet Wavelets constitute a family of functions constructed from dilation and translation of a single function ψ (x) called the mother wavelet. When the dilation parameter a and the translation parameter b vary continuously we have the following family of continuous wavelets [, ]: k ψkn(t) = |a|  ψ akt – nb , k, n ∈ Z, In particular, the rational numbers Em = mEm(/) are called the classical Euler numbers. Also, the Euler polynomials of the first kind for k = , . . . , m can be constructed from the following relation: f (t) = n= m∈Z n= m= and we can approximate the function f (t) by the truncated series km Ek(t) + Em(t) = tm, where (km) is a binomial coefficient. Explicitly, the first basic polynomials are expressed by E(t) = , E(t) = t – t, These polynomials satisfy the following formula: m, n ≥ , and the Euler polynomials form a complete basis over the interval [, ]. Furthermore, when t = , we have E() = , E() = –  , E() =  , E() = –  , 3.2 Function approximation A function f (t), square integrable in [, ], may be expressed in terms of the Euler wavelet as where the coefficient vector C and the Euler function vector (t) are given by n= m= AT = CTD, Taking the collocation points as follows: ti = ik–M , i = , , . . . , k–M, we define the Euler wavelet matrix mˆ × mˆ = where mˆ = k–M. Notation: from now we define mˆ = k–M. To evaluate C, we let aij = Using equation () we obtain aij = n= m= where dnijm =  ψnm(t)ψij(t)f (t)dt and i = , , . . . , k–, j = , , . . . , M – . Therefore, D = dnijm , D = where D is a matrix of order k–M × k–M and is given by The matrix D in equation () can be calculated by using equation () in each interval n = , , . . . , k–. For example, with k =  and M = , D the identity matrix, and for k =  and M =  we have CT = AT D–. Hence, CT in equation () is given by – √      – √  ⎥⎦  k(s, t) = (s)TK (t), where K is a k–M × k–M matrix given by [] K = D– k(s, t) (s) (t)dt D–. 3.3 Convergence of Euler wavelets basis We first state some basic results as regards Euler polynomials approximations. The important properties will enable us to establish the convergence theorem of the Euler wavelets basis. The Euler polynomials of degree m are defined by []. Now we defined (t) = [E(t), E(t), . . . , EN (t)]T , so a function f (t) ∈ L[, ] can be expressed in terms of the Euler polynomials basis (t). Hence, eiEi(t) = ET (t), where E = [e, e, . . . , eN ]T . Lemma  Suppose that the function f : [, ] → R is m+ times continuously differentiable, and f ∈ Cm+[, ], Y = span{E, E, . . . , EN } is vector space. If ET (t) is the best approximation of f out of Y, then the mean error bound is presented as follows: Proof Consider the Taylor polynomials fˆ(t) = f (t) + f (t)(t – t) + f (t) where we have Since ET (t) is the best approximation of f (t), we have M˜ Sm+ ≤ [(m + )!](m + ) . Theorem  Suppose that the function f : [, ] → R is m +  times continuously differentiable and f ∈ Cm+[, ]. Then f˜(t) = CT (t) approximates f (t) with mean error bounded as follows: where M˜ = maxt∈[,] |f (m+)(t)|. Proof We divide the interval [, ] into subintervals Ik,n = [ nk–– , kn– ], n = , . . . , k– with the restriction that f˜(t) is a polynomial of degree less than m +  that approximates f with minimum mean error. The approximation approaches the exact solution as k approaches ∞. We use Lemma , to obtain f (t) – f˜(t)  = M˜  ≤ (k–)(m+)[(m + )!](m + ) , where M˜n = maxt∈Ik,n |f (m+)(t)|. By taking the square roots we arrive at the upper bound. The error of the approximation f˜(t) of f (t) therefore decays like –(m+)(k–). bi(t) = bi(t)bj(t) = f (t) ≈ fibi(t) = FTBm(t), fi = m bi(t)f (t) dt. where F = [f, f, . . . , fm–]T, Bm(t) = [b(t), b(t), . . . , bm–(t)]T. By using the orthogonality of BPFs, for i = , , . . . , m – , the coefficients fi can be obtained: By using the disjointness of the BPFs and the representation of Bm(t), we have The block pulse operational matrix of the fractional integration Fα has been given in [], ξ ⎥⎥⎦  ξκ = (κ + )α+ – κα+ + (κ – )α+. Figure 1 1/2 order integration of t. There is a relation between the block pulse functions and Euler wavelets, 4 Method of numerical solution Consider the nonlinear fractional-order integro-differential equation k(x, t) y(t) p dt + g(x), subject to the initial conditions y(i)() = δi, i = , , . . . , r – , r ∈ N , where y(i)(x) stands for the ith-order derivative of y(x), Dtα (r –  < α ≤ r) denotes the Caputo fractional-order derivative of order α, g(x) ∈ L[, ], k ∈ L([, ]) are given func(t) = where matrix Pα is called the Euler wavelet operational matrix of fractional integration. Using equations () and (), we have B mˆ(t) = Itα Bmˆ (t) ≈ Combining equation () and equation (), we can get We select the function t to verify the correctness of fractional integration operational matrix Pα . The fractional integration of order α for the function f (t) = t is given by k(x, t) = (x)TK (t) g(x) ≈ GT (x), where K = [kij], i, j = , , . . . , mˆ, and G = [g, g, . . . , gmˆ ]T . Using equations () and (), we have Bmˆ (x) + Y˜ T B mˆ(x). y(x) ≈ EB mˆ(x). In the above summation, we substitute the supplementary conditions () and approximate it with the Euler wavelet, we can get where Y˜ is an mˆ -vector. According to equation (), the above equation can be written as tions, y(x) is the solution to be determined, λ is a real constant, and p ∈ N . The given functions g, k are assumed to be sufficiently smooth. Now we approximate Dtαy(x), k(x, t), and g(x) in terms of Euler wavelets as follows: By using the disjointness property of the BPFs, we have = eb(x) + eb(x) + · · · + emˆ –bmˆ –(x) = e, e, . . . , emˆ – Bmˆ (x) = EB mˆ(x), y(x) p ≈ ep, ep, . . . , epmˆ – B mˆ(x) = EpBmˆ (x), where p is any positive integer. Using equations (), (), and () we will have k(x, t) y(t) p dt =  x T(x)K (t)BTmˆ(t)EpT dt B mˆ(t)BTmˆ(t)EpT dt = BTmˆ (x) TK = Q˜ TBmˆ (x), Bmˆ (t)BTmˆ(t)EpTdt diag(Ep)B mˆ(t) dt diag(Ep) diag(Ep)FB mˆ(x) diag(Ep)FB mˆ(x) where Q˜ is an mˆ -vector with elements equal to the diagonal entries of the following matrix: Q = diag(Eq)F. Y T B mˆ(x) = λ Q˜TBmˆ (x) + GT B mˆ(x). Substituting the above equations into equation (), we have Using Bmˆ (x) to multiply two sides of equation () and integration in the interval [, ], according to orthogonality of the BPFs we can get which is a nonlinear system of algebraic equations. By solving this system we can obtain the approximation of equation (), and we solve the nonlinear system by using the Newton iterative method. 5 Numerical examples In this section, six examples are given to demonstrate the applicability and accuracy of our method. Examples - have smooth solutions, while Example  has a non-smooth and singular solution. In all examples the package of Matlab . has been used to solve the test problems considered in this paper. Using equation () the absolute error function is defined as R mˆ(x) = Y T B mˆ(x) – λQ˜ TB mˆ(x) – GT B mˆ(x) , where mˆ = k–M; M is the degree of the Euler polynomials and usually takes small values in a computation. Since the truncated Euler wavelet series is an approximate solution of equation (), we must have Rmˆ (x) ≈ . In the following examples, we can find that when M is fixed, the larger the value of k, the more accurate the approximation solution of equation. So the optimum value of k is determined by the prescribed accuracy. To demonstrate the effectiveness of this method, we will adopt the same error definition as []. The approximate norm- of the absolute error is given by emˆ (x)  = y(x) – ymˆ (x)  ≈ N i= where y(x) is the exact solution and ymˆ (x) is the approximation solution obtained. k = 4 5.3744e-003 5.9284e-004 1.6221e-003 1.6909e-003 4.6755e-004 2.0212e-003 4.8121e-004 1.7295e-003 1.7247e-003 4.6700e-004 k = 4 6.0169e-003 1.2504e-003 5.1499e-005 1.1123e-004 1.0434e-003 3.4768e-003 9.2652e-004 3.6739e-004 4.0781e-004 8.0375e-004 k = 5 2.1912e-003 3.6354e-004 1.0167e-004 1.1436e-004 4.3188e-004 5.0356e-004 4.3422e-004 1.2103e-004 1.1954e-004 4.2982e-004 Table 2 Approximate norm-2 of absolute errors for some k of the Euler and SCW 8.3942e-007 9.4203e-007 7.2293e-008 5.9209e-008 7.1159e-009 3.7129e-009 k = 5 2.5058e-003 4.3350e-004 6.1143e-004 5.5204e-004 2.0180e-004 1.1110e-003 1.5066e-004 4.4277e-004 4.2369e-004 9.2562e-005 7.1538e-007 1.6350e-005 k = 6 9.7161e-004 1.6976e-005 1.0476e-004 1.0746e-004 3.0374e-005 1.2551e-004 3.0730e-005 1.0861e-004 1.0818e-004 2.9460e-005 k = 6 1.1272e-003 3.2831e-004 2.0079e-004 1.7628e-004 2.3828e-004 3.8232e-004 2.1590e-004 1.2863e-004 1.1993e-004 1.8947e-004 2.3580e-007 1.1839e-006 5.7642e-008 8.6352e-008 Example  Let us consider the following fractional nonlinear integro-differential equa(x – t) y(t) dt + g(x),  ≤ x < , where g(x) =    (/) (  x/ – x/) –  x + x – x , and the equation is subject to the initial absolute errors obtained by Euler wavelets and SCW [], respectively. Table  shows the approximate norm- of absolute errors obtained by the Euler wavelet and SCW methods. From Table , we find that the absolute errors become smaller and smaller with k increasing. Table  shows that the Euler wavelet method can reach a higher degree of accuracy than the SCW method. Example  Consider the nonlinear fractional-order Volterra integro-differential equation  ≤ x < , exact solution of this equation is y(x) = x. Table  shows the approximate norm- of absolute errors obtained by the Euler wavelet and SCW methods. The comparisons between approximate and exact solutions for varresults become more accurate and we infer that the approximate solutions converge to the exact solution. Example  Consider the nonlinear Volterra integro-differential equation dt = g(x),  ≤ x < , Figure 2 The approximate solution of Example 2 for some k. Table 3 Comparison of approximate norm-2 of absolute errors with reproducing kernel and CAS k = 3, M = 3 2.5664e-006 2.3107e-006 k = 3, M = 1 5.6933e-004 9.0145e-004 where g(x) = (/) (  x/ – x/) +  x +  x –  x +  x +  x, and the equation is subject to the initial conditions y() = y () = . The exact solution of this equation is y(x) = x – x. Table  shows the approximate solution obtained by our method ( mˆ = (k–)M), reproducing the kernel method ( mˆ = k (M + )) [] and CAS wavelet methods ( mˆ = k (M + )) []. To make each method having the same number of wavelet bases, we select M =  for the Euler wavelet. Under the condition of the same error, our method is closer to the exact solution. Example  Consider this equation: (x + t) y(t) dt = g(x),  ≤ x < , where g(x) = ( ) √x – x – x –  , and the supplementary condition y() = y () = . The exact solution is y(x) = x. Table  shows the approximate solution obtained by the Euler wavelet method, reproducing the kernel method and CAS wavelet methods. From Table  we can see our method is closer to the exact solution. Example  In the following we consider the fourth-order equation [] e–t y(t)  dt = , Table 4 Numerical results for Example 5 with comparison to CAS Figure 3 Numerical and exact solution of Example 5 for mˆ = 8. parison with []. Table  shows the Euler wavelet numerical solution to be in excellent agreement with the solution of CAS method in []. It is worth noticing that the method introduced above only can solve equation () for variable x of equation () is defined on the interval [, ], so we should turn in Figure . The numerical solution is in perfect agreement with the exact solutions. Let us consider examples with non-smooth and singular solutions. Example  Consider the following equation: x [y(t)] Dtαy(x) = –  (x – t)/ dt + g(x), where g(x) = – (–/()/x–)/ + π . In this case, there is a singularity at point x = . The solution around this point is not good (see Figure  with k = , M = ). The Euler wavelet method can be combined with the definition of Riemann-Liouville fractional integral operator to deal with the weakly singular integral. As observed, our method provides a reasonable estimate even in this case with singular solution. Figure 4 The approximate solution of Example 6 for k = 4, M = 2. In the examples above, we do not show the computational times of the different methods. In fact, the Euler wavelet method has the faster computing speed, compared with the CAS wavelet method and the second Chebyshev wavelet method. In Example , for . s, . s, and . s, while the computational times of the Euler wavelet are . s, . s, and . s. The same conclusion can be drawn from the other examples. 6 Conclusion In this paper, we construct the Euler wavelet and derive the wavelet operational matrix of the fractional integration, and we use it to solve the fractional integro-differential equations. By solving the nonlinear system, approximate solutions are got. Graphical illustrations and tables of the numerical results with the aid of Euler wavelets indicate that the numerical results are well in agreement with exact solutions and superior to other results. Also the proposed method can be efficiently applied to a large number of similar fractional problems. Of course, the convergence of this algorithm has not been derived, which will be future research work. Competing interests The authors declare that they have no competing interests. Authors’ contributions All authors contributed equally to the writing of this paper. All authors read and approved the manuscript. Acknowledgements This work was supported by the KC Wong Education Foundation, Hong Kong, and the Project of Education of Zhejiang Province (No. Y201533324). The authors gratefully acknowledge the support of the KC Wong Education Foundation, Hong Kong. The authors are grateful to the editor and the anonymous referees for their constructive and helpful comments. 1. Yuzbasi , S: A numerical approximation for Volterra's population growth model with fractional order . Appl. Math. Model . 37 , 3216 - 3227 ( 2013 ) 2. Sadeghian , H, Salarieh , H, Alasty , A, Meghdari, A: On the fractional-order extended Kalman filter and its application to chaotic cryptography in noisy environment . Appl. Math. Model . 38 , 961 - 973 ( 2014 ) 3. Kumar , S, Kumar, A, Argyros, IK: A new analysis for the Keller-Segel model of fractional order . Numer. Algorithms ( 2016 ). doi:10.1007/s11075- 016 - 0202 -z 4. Li , C, Kumar, A, Kumar , S, Yang, X-J: On the approximate solution of nonlineartime-fractional KdV equation via modified-homotopy analysis Laplace transform method . J. Nonlinear Sci. Appl . 9 , 5463 - 5470 ( 2016 ) 5. Kumar , D, S ingh, J, Kumar, S , et al.: Numerical computation of nonlinear shock wave equation of fractional order . Ain Shams Eng. J . 6 , 605 - 611 ( 2014 ) 6. Kumar , D, S ingh, J, Kumar, S , et al.: Numerical computation of Klein-Gordon equations arising in quantum field theory by using homotopy analysis transform method . Alex. Eng. J . 53 , 469 - 474 ( 2014 ) 7. Kumar , S, Kumar, D, Singh, J: Fractional modelling arising in unidirectional propagation of long waves in dispersive media . Adv. Nonlinear Anal . ( 2016 ). doi:10.1515/anona- 2013 -0033 8. Baleanu , D, Srivastava, HM, Yang, XJ : Local fractional variational iteration algorithms for the parabolic Fokker-Planck equation defined on Cantor sets . Prog. Fract. Differ. Appl . 1 , 1 - 10 ( 2015 ) 9. Yang , XJ, Baleanu , D, Srivastava, HM : Local fractional similarity solution for the diffusion equation defined on Cantor sets . Appl. Math. Lett . 47 , 54 - 60 ( 2015 ) 10. Yang , AM, Han , Y, Mang , YZ , et al.: On local fractional Volterra integro-differential equations in fractal steady heat transfer . Therm. Sci . 20 , S789 - S793 ( 2016 ) 11. Nawaz , Y: Variational iteration method and homotopy perturbation method for fourth-order fractional integro-differential equations. Comput. Math. Appl. 61 , 2330 - 2341 ( 2011 ) 12. Odibat , ZM: A study on the convergence of variational iteration method . Math. Comput. Model . 51 , 1181 - 1192 ( 2010 ) 13. Sayevand , K: Analytical treatment of Volterra integro-differential equations of fractional order . Appl. Math. Model . 39 , 4330 - 4336 ( 2015 ) 14. Momani , S, Noor, MA: Numerical methods for fourth-order fractional integro-differential equations . Appl. Math. Comput . 182 , 754 - 760 ( 2006 ) 15. Arikoglu , A, Ozkol, I: Solution of fractional integro-differential equations by using fractional differential transform method . Chaos Solitons Fractals 34 , 1473 - 1481 ( 2007 ) 16. Wei , J, Tian, T: Numerical solution of nonlinear Volterra integro-differential equations of fractional order by the reproducing kernel method . Appl. Math. Model . 39 , 4871 - 4876 ( 2015 ) 17. Rawashdeh , E: Numerical solution of fractional integro-differential equations by collocation method . Appl. Math. Comput . 176 , 1 - 6 ( 2006 ) 18. Ma, XH, Huang, CM: Numerical solution of fractional integro-differential equations by a hybrid collocation method . Appl. Math. Comput . 219 , 6750 - 6760 ( 2013 ) 19. Saeedi , H, Mohseni Moghadam , M, Mollahasani, N, Chuev, GN: A CAS wavelet method for solving nonlinear Fredholm integro-differential equations of fractional order . Commun. Nonlinear Sci. Numer. Simul . 16 , 1154 - 1163 ( 2011 ) 20. Saeedi , H, Mohseni Moghadam , M: Numerical solution of nonlinear Volterra integro-differential equations of arbitrary order by CAS wavelets . Commun. Nonlinear Sci. Numer. Simul . 16 , 1216 - 1226 ( 2011 ) 21. Zhu , L, Fan, QB : Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet . Commun. Nonlinear Sci. Numer. Simul . 17 , 2333 - 2341 ( 2012 ) 22. Zhu , L, Fan, QB : Numerical solution of nonlinear fractional-order Volterra integro-differential equations by SCW . Commun. Nonlinear Sci. Numer. Simul . 18 , 1203 - 1213 ( 2013 ) 23. Meng , Z, Wang , L, Li , H, Zhang , W: Legendre wavelets method for solving fractional integro-differential equations . Int. J. Comput. Math . 92 , 1275 - 1291 ( 2015 ) 24. Wang , YX, Zhu , L: SCW method for solving the fractional integro-differential equations with a weakly singular kernel . Appl. Math. Comput . 275 , 72 - 80 ( 2016 ) 25. Beylkin , G, Coifman , R, Rokhlin , V: Fast wavelet transforms and numerical algorithms I. Commun. Pure Appl. Math. 44 , 141 - 183 ( 1991 ) 26. Podlubny , I: Fractional Differential Equations: An Introduction to Fractional Derivatives , Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications . Academic Press, New York ( 1999 ) 27. Wang , YX, Fan, QB : The second kind Chebyshev wavelet method for solving fractional differential equation . Appl. Math. Comput . 218 , 8592 - 8601 ( 2012 ) 28. Zhu , L, Wang, YX : Numerical solutions of Volterra integral equation with weakly singular kernel using SCW method . Appl. Math. Comput . 260 , 63 - 70 ( 2015 ) 29. Yuan , H: Some new results on products of Apostol-Bernoulli and Apostol-Euler polynomials . J. Math. Anal. Appl . 431 , 34 - 46 ( 2015 ) 30. Maleknejad , K, Hashemizadeh , E, B asirat, B : Computational method based on Bernstein operational matrices for nonlinear Volterra-Fredholm-Hammerstein integral equations . Commun. Nonlinear Sci. Numer. Simul . 17 , 52 - 61 ( 2012 ) 31. Kilicman , A, Al Zhour, ZAA: Kronecker operational matrices for fractional calculus and some applications . Appl. Math. Comput. 187 , 250 - 265 ( 2007 )


This is a preview of a remote PDF: http://www.advancesindifferenceequations.com/content/pdf/s13662-017-1085-6.pdf

Yanxin Wang, Li Zhu. Solving nonlinear Volterra integro-differential equations of fractional order by using Euler wavelet method, Advances in Difference Equations, 2017, 27, DOI: 10.1186/s13662-017-1085-6