We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
|
|
- Clementine Roberts
- 5 years ago
- Views:
Transcription
1 We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4, 6, 2M Open access books available International authors and editors Downloads Our authors are among the 54 Countries delivered to TOP % most cited scientists 2.2% Contributors from top 5 universities Selection of our books indexed in the Book Citation Index in Web of Science Core Collection (BKCI Interested in publishing with us? Contact book.department@intechopen.com Numbers displayed above are based on latest data collected. For more information visit
2 6 Numerical Schemes for Fractional Ordinary Differential Equations Weihua Deng and Can Li School of Mathematics and Statistics, Lanzhou University, Lanzhou People s Republic of China. Introduction Fractional calculus, which has almost the same history as classic calculus, did not attract enough attention for a long time. However, in recent decades, fractional calculus and fractional differential equations become more and more popular because of its powerful potential applications. A large number of new differential equations (models that involve fractional calculus are developed. These models have been applied successfully, e.g., in mechanics (theory of viscoelasticity, biology (modelling of polymers and proteins, chemistry (modelling the anomalous diffusion behavior of Brownian particles, electrical engineering (electromagnetic waves etc (Bouchaud & Georges, 99; Hilfer, 2; Kirchner et al., 2; Metzler & Klafter, 2; Zaslavsky, 22; Mainrdi, 28. Meanwhile, some rich fractional dynamical behavior which reflect the inherenature of realistic physical systems are observed. In short, fractional calculus and fractional differential equations have played more and more important role in almost all the scientific fields. One of the most important fractional models is the following initial value problem D α y(t f(t, y(t, (x, k,,, α, ( where α (,, can be any real numbers, and D α denotes the fractional derivative in the Caputo sense, defined by D α y(tj n α D n y(t t Γ(n α y (n (τ dτ, (2 (t τ α n here n : α is the smallest integer not less than α, y (n (x is the classical nth-order derivative of y(x and for µ >,J µ is the µ-order Riemann-Liouville integral operator expressed as follows J µ y(t t Γ(µ y(τ dτ. (t τ µ The use of Caputo derivative in above equation is partly because of the convenience to specify the initial conditions. Since the initial conditions are expressed in terms of values of the unknown function and its integer-order derivatives which have clear physical meaning (Podlubny, 999; 22. However, from pure mathematical viewpoint, the Riemann-Liouville derivative is more welcome and many earlier research papers use it instead of Caputo derivative (Podlubny, 999. In general, specifying some additional conditions is necessary
3 356 Numerical Modelling 2 Will-be-set-by-IN-TECH to make sure the discussed equation has a unique solution. These additional conditions, in many situations, describe some properties of the solution at the initial time (Heymans & Podlubny, 25, the fractional derivative does not have convenient used physical meaning (there are already some progress in the geometric and the physical interpretation of the fractional calculus (Podlubny, 22 and the physical interpretation of the initial conditions in terms of the Riemann-Liouville fractional derivatives of the unknown function has also been discussed in (Podlubny, 22. Just like the classic calculus and differential equations, the theories of fractional differentials, integrals and differential equations have been developing. With the development of the theories of fractional calculus, many research monographs are published, e.g., (Oldham & Spanier, 974; Podlubny, 999; Samko et al., 993. In the literatures, several analytical methodologies, such as, Laplace transform, Mellin transform, Fourier transform, are restored to obtain the analytical solutions of the fractional equations by many authors (Metzler & Klafter, 2; Podlubny, 999; Samko et al., 993; Zaslavsky, 22; Mainrdi, 28, however, similar to treating classical differential equations, they can mainly deal with linear fractional differential equation with constant coefficients. Usually, for nonlinear systems these techniques do not work. So in many cases the more reasonable option is to find its numerical solution. As is well known, the difficulty of solving fractional differential equations is essentially because fractional calculus are non-local operators. This non-local property means that the next state of a system not only depends on its current state but also on its historical states starting from the initial time. This property is closer to reality and is the main reason why fractional calculus has become more and more useful and popular. In other words, this non-local property is good for modeling reality, but a challenge for numerical computations. Much effort has been devoted during the recent years to the numerical investigations of fractional calculus and fractional dynamics of ( (Lubich, 985; 986; Podlubny, 999. More recently, Diethelm et al successfully presented the numerical approximation of ( using Adams-type predictor-corrector approach (Diethelm & Ford, 22a and the detailed error analysis of this method was given in (Diethelm et al., 24. The convergent order of Diethelm s predictor-corrector approach was proved to be min(2, α (Diethelm et al., 24. Because of the non-local property of the fractional derivatives, the arithmetic complexity of their algorithm with step size h is O(h 2, whereas a comparable algorithm for a classical initial value problem only gives rise to O(h. To improve the accuracy and reduce the arithmetic complexity, some techniques such as the Richardson extrapolation, short memory principle and corresponding mixed numerical schemes are developed. In (Deng, 27a, we present an improved version of the predictor-corrector algorithm with the accuracy increased to min(2, 2α and half of the computational cost is reduced comparing to the original one in (Diethelm et al., 24. Furthermore, we apprehend the short memory principle from a new viewpoint (Deng, 27b; after using the nested meshes presented in (Ford & Simpson, 2 and combining the short memory principle and the predictor-corrector approach, we minimize the computational complexity to O(h log(h at preserving the order of accuracy 2. This chapter briefly reviews the recent development of the predictor-corrector approach for fractional dynamic systems. The plan of this chapter is as follows. In Section 2, we briefly discuss the short memory principle and the nested meshes. In Section 3, the predictor-corrector schemes and its improved versions are presented, meanwhile the convergent order and arithmetic complexity are also proposed. In Section 4, we provide two
4 } Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations numerical examples to illustrate the performance of our numerical schemes. Conclusions are given in the last section. 2. Short memory principle We can see that the fractional derivative (2 is an operator depending on the past states of the process y(t (see Fig. However, for t and < α <, the behavior of y(t far a D* y( t the " past " of yt ( t Fig.. The fractional derivative operating on the "past" of y(t (the red part. from upper terminal suggests that the original fractional derivative can possibly be replaced by the fractional derivative with a moving lower terminal (Podlubny, 999. This means that the history of the process y(t can be approached only over a fixed period of recent history. Inspired by this kind of idea, Podlubny in (Podlubny, 999 introduces a fixed integral length memory principle for Riemann-Liouville derivative (see Fig 2. He shows that the truncation a a Dt y( t t- LDty( t } t- L t Fig. 2. The "short-memory" principle with "fixed memory length" L(< t, where D α t denote the Riemann-Liouville derivative operator. error gives E < ML α /Γ( α with a fixed integral length L. To accelerate the computation without significant loss of accuracy, in (Ford & Simpson, 2 Ford and Simpson present a short-memory principle for Caputo derivative. In (Deng, 27b we apprehend the short memory principle from a new viewpoint, and then it is closer to reality and extends the effective range of short memory principle from α (, to α (, 2. In view of the scaling property of fractional integral, the nested meshes are possibly used. And the nested meshes can be produced by splitting the integral interval [, ] into [. ][, p m T ] [ p m T, p m T ] [ p 2 T, pt ] [ pt, ], (3 where T ωh, h R, m, ω, p N and p m T < p m T. Denote M h {hn, n N} and l, l 2 N, l > l 2, then M l2 h M l h. 3. Predictor-corrector schemes The problem of determining y(t by means of the information of initial value y, is called initial value problem. As the classic theory of ordinary differential equation, if a function y(t satisfies the initial value problem ( pointwisely, then y(t is called the solution of the fractional differential equation (. For the existence and uniqueness of the solutions of fractional differential equations, one can see (Podlubny, 999. But the explicit formula of the solution y(t can t usually be obtained in spite of we can prove the existence of solution. The most
5 358 Numerical Modelling 4 Will-be-set-by-IN-TECH reasonable way is to use numerical methods, and the obtained solution is called the numerical approximation of the solution of the differential equations ( and denoted by y h in the following sections. 3. Numerical schemes In this section, we show the predictor-corrector schemes of (. Using the Laplace transform formula for the Caputo fractional derivative (Podlubny, 999 From (, we have or L { D y(t α } s α n Y(s s α n Y(s k k s α k (, n < α n. (4 s α k (F(s, Y(s, (5 Y(ss α n F(s, Y(s k s k (, (6 where F(s, Y(s L( f(t, y(t. Applying the inverse Laplace transform gives where the fact and y(t α k L { J µ y(t } L { Γ(µ t k t (t τ α f(τ, y(τdτ, (7 Γ(α t y(τ (t τ µ dτ} L { t µ Γ(µ y(t} s µ Y(s, L { t µ } s µ Γ(µ. are used. The approximation is based on the equivalent form of the Volterra integral equation (7. A fractional Adams-predictor-corrector approach was firstly developed by Diethelm et al (Diethelm et al., 22b to numerically solve the problem (7. Using the standard quadrature techniques for the integral in (7, denote g(τ f(τ, y(τ, the integral is replaced by the trapezoidal quadrature formula at point ( τ α tn g(τdτ ( τ α g n (τdτ, (8 where g n is the piecewise linear interpolation of g with nodes and knots chosen at t j, j,, 2,, n. After some elementary calculations, the right hand side of (8 gives ( τ α g n (τdτ n j a j,n g(t j, (9 where the uniform mesh is used and h is the stepsize. And if we use the product rectangle rule, the right hand of (8 can be written as ( τ α g n (τdτ n b j,n g(t j, ( j
6 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations where n α (n α(n α, if j, a j,n (n j2 α 2(n j α (n j α, if j n,, if j n and b j,n hα α [(n jα (n j α ], if j n. Then the predictor and corrector formulae for solving (7 are given, respectively, by and y h ( α k y P h ( t k n α k t k n n Γ(α b j,n f(t j, y h (t j ( j Γ(2α f(, y P h ( Γ(2α n a j,n f(t j, y h (t j. j (2 The approximation accuracy of the scheme (-(2 is O(h min{2,α} (Diethelm et al., 24. Now we make some improvements for the scheme (-(2. We modify the approximation of (8 as (Deng, 27a ( τ α g(τdτ ( τ α g n (τdτ ( τ α g( dτ, (3 where g n is the piecewise linear interpolation for g with nodes and knots chosen at t j, j,, 2,, n. Then using the standard quadrature technique, the right hand of (3 can be recast as where ( τ α g n (τdτ bj,n ( τ α g( dτ { a j,n, if j n, 2 α if n >,, if j n. b, α, if n. Hence, this algorithm for the predictor step can be improved as (Deng, 27b y P h ( α k t k n Γ(2α n bj,n g(t j (4 j n bj,n f(t j, y h (t j, (5 j The new predictor-corrector approach (5 and (2 has the numerical accuracy O(h min{2,2α} (the detailed analysis is left in Section 3.2. Obviously half of the
7 36 Numerical Modelling 6 Will-be-set-by-IN-TECH computational cost can be reduced, for < α, if we reformulate (5 and (2 as h y α Γ(α f(y h(t, t, if n, yh P(t h n y α Γ(α2 (2α f(y h (, hα Γ(α2 n a j,n f(y h (t j, t j, j if n (6 and y h ( y ( f(y P Γ(α2 h (t, t α f(y h (t, t, if n, y Γ(α2 hα Γ(α2 ( f(yh P(, (2 α 2 f(y h (, n a j,n f(y h (t j, t j, j if n. (7 The arithmetic complexity of the above two predictor-corrector schemes ((-(2, and (6-(7 is O(h 2, where h is step size. For further reducing the computational cost, we understand the short memory principle from a new viewpoint, and then go to design the predictor-corrector scheme (Deng, 27b. For α (, and α (,, we rewrite (7 as, respectively, and y( y( (t Γ(α n τ α f(τ, y(τdτ ( (tn τ α ( τ α f(τ, y(τdτ, α (, (8 Γ(α y( α k (t k n tk ny( (t Γ(α n τ α f(τ, y(τdτ ( (tn τ α ( τ α f(τ, y(τdτ, α (,. (9 Γ(α By observation of (8 and (9, we can see that the non-local property of D α induces the term tn ( (tn τ α ( τ α f(τ, y(τdτ. (2 Γ(α
8 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations 7 36 In fact, if α (, 2 the integration kernel of (2 fades faster when the time history becomes longer, more concretely, tn ( (tn τ α ( τ α f(τ, y(τdτ Γ(α ( τ Γ(α(α Γ(α(α Γ(α(α 2 τ ( τ τ ( τ τ ( τ z α 2 dz f(τ, y(τdτ z α 2 dz f(τ, y(τdτ z α 2 dz f(τ, y(τdτ t2 z α 2 dz f(τ, y(τdτ Γ(α(α t τ t ( τ z α 2 dz f(τ, y(τdτ (2 Γ(α(α Γ(α(α t τ (z (τα 2 f(τ, y(τdτ t2 (z n Γ(α(α (τα 2 f(τ, y(τdτ Γ(α(α t t t (z n(τ α 2 f(τ, y(τdτ, tn (z2 Γ(α(α (τα 2 f(τ, y(τdτ 2 where z m(τ (t m, t m, m, 2,, n. Obviously, we can see that the integration (2 s kernel ( τ α ( τ α decays (algebraically by the order 2 α when α (, 2, but in (Ford & Simpson, 2 the integral kernel ( τ α decays with the order α when α (,. This is the main reason why we can extend the range of the short memory principle of fractional differential equations from α (, to α (, 2. For the nested mashes defined by (3, we can take the step length h in the integral[ pt, ] and in the subsequent intervals[ p 2 T, pt ],[ p 3 T, p 2 T ],,[ p m T, p m T ],[t, p m T ], step lengths ph, p 2 h,, p m h, p m h are used, respectively. Noting that ( p m T can t be divided by p m h, so the integral in the interval[, l](l ( p m T ( p m T /(p m h (p m h is ignored, it does not destroy the computational accuracy in general. We will pay our attention to the numerical approximation by using short memory principle in the integration (2 in the following part. With the similar approximation of (3, the product trapezoidal quadrature formula is applied to replace the integral of (9, where nodes t j, j n, n are taken with respect to the weighted function ( α for the first integral and nodes t j, j,,, n are used w.r.t. the weighted function ( α ( α for the second integral, and it yields ( τ α f(τ, y(τdτ ( τ α fn (τ, y(τdτ ( α f(tn, y( f(t n, y(, α (, (22
9 362 Numerical Modelling 8 Will-be-set-by-IN-TECH and ( (tn τ α ( τ α f(τ, y(τdτ ( (tn τ α ( τ α fn (τ, y(τdτ n ã j,n f(t j, y(t j (23 j where (n α (α nn α (2n α (n α, if j, ã j,n (n j2 α 3(n j α 3(n j α (n j α, if j n, 2 α α 3, if j n. (24 For the first integral of (8 or (9, the product rectangle formula is used ( τ α f(τ, y(τdτ ( τ α f(, y( dτ hα α f(, y(. (25 Combing (22-(25, for α (, 2, the new predictor-corrector approach gives y P h (y ( hy h ( Γ(α f(, y h ( hα n Γ(α2 ã j,n f(t j, y h (t j, α (, 2 (26 j and y h ( y ( hy h ( f(, y P h ( ( α f(tn, y Γ(α2 h ( n Γ(α2 ã j,n f(t j, y h (t j, α (, 2. (27 j For α (2,, the predictor-corrector approach is y P h ( α k hα Γ(α f(, y h ( (t k n tk ny h ( Γ(α2 n ã j,n f(t j, y h (t j, α (2, (28 j and y h ( α k f(, y P h ( (t k n tk ny h ( ( α f(tn, y Γ(α2 h ( n Γ(α2 ã j,n f(t j, y h (t j, α (2,. (29 j
10 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations The predictor and corrector formulae based on the analytical formula (9 is fully described by (28 and (29 with the weighted ã j,n defined by (24. It is apparent that the same term n j ãj,n f(t j, y h (t j exists in both predictor (28 and corrector formulae (29, we juseed to compute one times at each predictor-corrector iteration step. In order to reduce the computational cost, in conjunction with the nested meshes memory principle (3, decompose the integral and still use the product trapezoidal quadrature formula at each subinterval but with different step lengths. ( (tn τ α ( τ α f(τ, y(τdτ ( pωh m i p i ωh p i ωh p m ωh ((tn τ α ( τ α fn (τ, y(τdτ n b j,p,n f(t j, y(t j (3 jn pω m i (p i h α ( (p ω j (pm h α n/p m ω j b j,p i,n f( p i (ωjh, y( p i (ωjh b j,p m,n f( p m (ωjh, y( p m (ωjh, (3 where (pω α (α pω( pω α (2pω α (pω α, if j n pω, (n j2 α 3(n j α 3(n j α b j,p,n (n j α, if n pω j n, 2 α α 3, if j n and for i, 2,, m, (/p i ω α ω α (/p i ω α (ω α ((/p i ω α ω α (α, if j, b j,p i,n (ωj /p i α (ωj α 2(ωj/p i α, 2(ωj α (ωj/p i α (ωj α, if j (p ω, (p i ω /p i α (p i ω α (p i ω /p i α (p i ω α ((p i ω /p i α (p i ω α (α, if j (p ω. Employing above analysis we obtain the following predictor-corrector algorithm y P h (y ( hy h ( m i (p i h α ( (p ω Γ(α f(, y h ( j (pm h α n/p m ω j n b j,p,n f(t j, y(t j jn pω b j,p i,n f( p i (ωjh, y( p i (ωjh b j,p m,n f( p m (ωjh, y( p m (ωjh (32
11 364 Numerical Modelling Will-be-set-by-IN-TECH and y h ( y ( hy h ( f(, y P h ( m i (p i h α ( (p ω j (pm h α n/p m ω j ( α f(tn, y Γ(α2 h ( n b j,p,n f(t j, y(t j jn pω b j,p i,n f( p i (ωjh, y( p i (ωjh b j,p m,n f( p m (ωjh, y( p m (ωjh. (33 In this case, the nested meshes predictor-corrector algorithm has the computational cost of O(h log(h for α (, Error analysis and convergent order In this section, we make the local truncation error and convergent order analysis for the improved predictor-corrector approaches (6-(7, (26-(29, and (32-(33. First we present several lemmas which will be used later. Lemma. (Diethelm et al., 24 Suppose g C 2 [, T], ( τ α tn g(τdτ ( τ α g n (τdτ C αh 2, (34 where C α only depends on α. Proof. g 2 ( τ α ( g(τ g n (τ dτ g 2 2 g 2 2 n tj ( τ α (t j τ(τ t j dτ j t j n( (n j2 α (n j α 2 ( (n j2 α2 (n j α2 α2 j n j ( (j α j α 2 ( (n α2 α2 g 2 ( n 2 τ α dτ 2 g 2 24 n j jα, if α <, g 2 24 n j (jα, if α, n( (j α j α j { g 2 n 24 g 2 24 τ α dτ, if α <, n2 2 τ α dτ, if α, C α h 2. (35
12 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations 365 With the similar method, we can obtain the following lemma. Lemma 2. Suppose g C 2 [, T], ( τ α g(τdτ where C α only depends on α. The error of the product rectangle rule is given as ( τ α g n (τdτ C αh 2, (36 Lemma 3. Suppose that f(τ, y(τ/ t C[, t, for some suitable t, then we have tn (t Γ(α n τ α ( f(τ, y(τ f(, y( dτ C α, (37 where C α only depends on α. Proof. Γ(α ( τ α ( f(τ, y(τ f(, y( dτ f(τ, y(τ/ τ Γ(α f(τ, y(τ/ τ Γ(α ( τ α (τ dτ hα C α where C α f(τ, y(τ/ τ. (38 Γ(α2 Lemma 4. Suppose that 2 f(τ, y(τ/ 2 τ C[, t, for some suitable t, then we have tn (t Γ(α n τ α ( f(τ, y(τ f n (τ, y(τ dτ C α2, (39 where C α only depends on α. Proof. According to the property of linear interpolation polynomials, f(τ, y(τ f n (τ, y(τ f[τ,, ](τ (τ, (4 where f[τ,, ] is second divided differences. And using the fact tj ( τ α (τ t j (τ t j dτ t j tj t j tj t j ( τ α [ (t j τ(t j τ ] dτ [ tn j j ( τ α ( j j ( τ α ( τ α] dτ (tα n j j t α n j j 2 (α2 (tα2 n j tα2 n j (α(α2 (tα2 n j tα2 n j h (tα n j tα n j for all j, (4
13 366 Numerical Modelling 2 Will-be-set-by-IN-TECH we have tn (t Γ(α n τ α ( f(τ, y(τ f n (τ, y(τ dτ tn (t Γ(α n τ α f[τ,, ](τ (τ dτ f[ξ,, ] Γ(α ( τ α (τ (τ dτ f (η, y(η 2Γ(α ( τ α (τ (τ dτ f (η, y(η h 2Γ(α α2 (taking j n in (4 (α(α2 C α 2 where C α f (η, y(η 2Γ(α(α(α2. (42 Here ξ, η [, ] and in the above equalities the second integral mean value theorem and the properties of second divided differences are used. Lemma 5. Suppose that 2 f(τ, y(τ/ 2 τ C[, t, for some suitable t, then we have tn ( (tn τ α ( τ α ( f(τ, y(τ Γ(α f n (τ, y(τ dτ C αh min{α2,3}, (43 where C α only depends on α. Proof. The idea of this lemma s proof is similar to the above lemmas, namely tn ( (tn τ α ( τ α ( f(τ, y(τ Γ(α f n (τ, y(τ dτ f (η, y(η n tj ( 2Γ(α (tn τ α ( τ α (τ t j (τ t j dτ j t j f (η, y(η n tj 2Γ(α ( τ α (τ t j (τ t j dτ j t j tj ( τ α (τ t j (τ t j dτ t j f (η, y(η n [ 2 2Γ(α (α2 (tα2 n j tα2 n j h (tα n j tα n j j 2 (α2 (tα2 f (η, y(η 2Γ(α f (η, y(η 2Γ(α n j tα2 n j n { 2(2t α2 n j tα2 j { 2(t α2 h (tα n j tα n j ] (using(4 n j tα2 n j (α2 n t α2 (α2 n tα2 } h (tα n j tα n j h (tα n tα n t } α
14 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations 3 f (η, y(η 2Γ(α f (η, y(η 2Γ(α f (η, y(η 2Γ(α f (η, y(η 2Γ(α f (η, y(η 2Γ(α Lemma 6. (Diethelm et al., 24 Assume that the solution of the initial value problem satisfies and ( τ α D α y(τdτ n k ( τ α D α n y(τdτ k bj,n D y(t α j C t γ n hδ (45 a j,n D y(t α j C 2t γ 2 n hδ 2, (46 with some γ, γ 2 and δ, δ 2. Then, for some suitably chosen T, we have max y(t j y h (t j O(h q, (47 j N with q min{δ α, δ 2 } and N T/h. Theorem. For the fractional initial problem (, if D α y(t C 2 [, T] for some suitable T, then the convergent order of our algorithm with the predictor and corrector formulae (6 and (7 gives max j,,,n y(t j y h (t j O(h min{2,2α}. Proof. Noting { 2(t α2 n t α2 n (α2 h (tα n h tα n (α(α2} α2 { h(t α n tα n 2(z α h α2 } (mean value theorem (α2 (α(α2 { h((t α n (z α ((z α t α n h α2 ]} (α(α2 { h 2 ((z α ( z α h α2 } (mean value theorem α (α(α2 { h 3 ((z α h α2 } (mean value theorem (α(α2 C α h min{α2,3}, (44 where z [, ], z [z, ] [, ], z [, z ] [, ], z [ z, z ] [, ] and C α f (η, y(η 2Γ(α (z α (α(α2. ( τ α ( g(τ g( dτ h g α, (48 and applying lemmas 2, 4, 6, we have the above result. Theorem 2. When α >,if 2 f(τ, y(τ/ 2 τ C[, t for some suitable t, then the local truncation error of our algorithm with the predictor and corrector formulae (26-(27 (α (, 2 and (28-(29 (α (2, iso(h 3, and the convergent order is 2, i.e., max j,,,n y(t j y h (t j O(h
15 368 Numerical Modelling 4 Will-be-set-by-IN-TECH Proof. This proof will be used based on mathematical induction. In view of the given initial condition, the induction basis j is presupposed, it has convergent order 2. Now, assuming that the convergent order is 2 for j,,, k, k n, we have the local truncation error y( { α k (t k n tk ny( hα ( α f(tn, y( f(t Γ(α2 n, yh P( { α k Γ(α2 n ã j,n f(t j, y h (t j } j (t k n tk ny( (t Γ(α n τ α f(τ, y(τdτ ( (tn τ α ( τ α } { α f(τ, y(τdτ Γ(α k hα ( α f(tn, y( f(t Γ(α2 n, yh P( { Γ(α ( τ α f(τ, y(τdτ hα ( f(tn, y(t Γ(α2 n f(, yh P( Γ(α2 { tn ( (tn τ α ( τ α f(τ, y(τdτ Γ(α { Γ(α2 Then we have y( n ã j,n f(t j, y(t j j { α k Γ(α2 (t k n tk ny( hα ( α f(tn, y( f(t Γ(α2 n, yh P( C 2 Ch 3, αl Γ(α2 hαmin{α,3} (t k n tk ny( n ã j,n f(t j, y h (t j j ( α f(tn, y( f(t Γ(α2 n, y( } Γ(α2 n ã j,n f(t j, y h (t j }. j Γ(α2 αl Γ(α2 hαmin{α2,3} } n } ã j,n f(t j, y(t j j n ã j,n f(t j, y h (t j } j ( 2 hα (z α h Lh 2 where z (,, Lemmas 4 and 5 in the above proof are used, and we also utilize the result y( y P h ( O(h min{α,3} which can be proved by using Lemmas 4 and 5
16 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations and the similar idea to above proof, its sketch proof is given as y( { α k hα Γ(α f(, y( { α k (t k n tk ny( Γ(α2 n ã j,n f(t j, y h (t j } j (t k n tk ny( (t Γ(α n τ α f(τ, y(τdτ ( (tn τ α ( τ α f(τ, y(τdτ Γ(α n } ã j,n f(t j, y h (t j j hα Γ(α f(t h n, y( α Γ(α2 { tn (t Γ(α n τ α f(τ, y(τdτ { Γ(α Γ(α Γ(α Ch min{α,3}. } { α k } Γ(α2 f(, y( ( (tn τ α ( τ α f(τ, y(τdτ ( τ α ( f(τ, y(τ f(, y( dτ ( (tn τ α ( τ α f(τ, y(τdτ (t k n tk ny( n } Γ(α2 ã j,n f(t j, y h (t j j Γ(α2 n j ã j,n f(t j, y h (t j We have proved that the local truncation error of our algorithm (26-(27 and (28-(29 is O(h 3 when α >, so the convergent order is 2. Lemma 7. (Ford Simpson [6,Theorem ] The nested mesh scheme preserves the order of the underlying quadrature rule on which it is based. Because of Theorem 2, Lemma 7, and the analysis in above section, we have Theorem 3. When α >,if 2 f(τ, y(τ/ 2 τ C[, t for some suitable t, then the local truncation error of our algorithm with the predictor and corrector formulae (32-(33 (α (, 2 iso(h 3 and the convergent order is 2, i.e., max j,,,n y(t j y h (t j O(h 2. The comparison of the local truncation error, convergent order and arithmetic complexity for different predictor-corrector schemes are presented in the following table. So far, our convergence results are obtained under the smoothness assumptions of D y(t α and f, which obviously depend on the smoothness properties of the solution y(t (Deng, Numerical results In this section we present two numerical examples to illustrate the performance of our proposed predictor-corrector schemes. We show the order of convergence in the the absolute
17 37 Numerical Modelling 6 Will-be-set-by-IN-TECH Schemes Local truncation error Convergent order Arithmetic complexity (,(2 O(h 3 O(h min{2,α} O(h 2 (6,(7 O(h 3 O(h min{2,2α} O(h 2 (26-(29 O(h 3 O(h 2 O(h 2 (32,(33 O(h 3 O(h 2 O(h log(h Table. Comparison of different predictor-corrector schemes. error. And the convergence order is measured by ( error(h Order log 2 error(h/2 where error(h is the absolute error y(t y h (t with the step h. Example. Consider the following fractional differential equation with the initial conditions or note that the exact solution to this problem is D α y(t Γ(5 Γ(5 α t4 α y(tt 4, α (, 2, (49 y(, α (,, (5 y(, y (, α (, 2, (5 y(tt 4. (52 The numerical results for the predictor-corrector scheme (6-(7 at time t with different steps and different α are reported in Tables 2 and 3. Table shows the numerical errors at time t between the exact solution and the numerical solution of the predictor-corrector schemes (6-(7 for (49 with different α (, for various step sizes. It is seen that the rate of convergence of the numerical results is of order O(τ min{2,2α} for the predictor-corrector schemes (6-(7. Tables 2 and 3 show the ratio of the error reduction as the grid is refined. The last row states the theoretical orders of convergence (abbreviated as TOC, which are the results we theoretically prove in the above theorems. From Tables 2 and 3, it can be noted that the numerical results are in excellent agreement with the theoretical ones given in Theorem. Example 2. Now, consider the following fractional differential equation (Diethelm et al., 24 2 Γ(3 α t2 α Γ(2 α t α y(tt 2 t, if < α <, D y(t α 2 Γ(3 α t2 α y(tt 2 t, if α. with the initial conditions or the exact solution to this problem is (53 y(, α (,, (54 y(, y (, α (, 2, (55 y(tt 2 t. (56
18 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations Order α.5 Order α.9 Order / / / / / / / TOC Table 2. Absolute errors and convergence orders for predictor-corrector schemes (6-(7 at time t with different < α <..25 Order α.5 Order α.85 Order / / / / / / / TOC Table 3. Absolute errors and convergence orders for predictor-corrector schemes (6-(7 at time t with different < α < 2. Table 4 shows the numerical errors at time t between the exact solution and the predictor-corrector schemes (6-(7 for (53 with different α (, for various step sizes. Again we numerically verify that the rate of convergence for predictor-corrector schemes (6-(7 is O(τ min{2,2α}. From Table 4, it can be seen that both the convergence orders. Order α.5 Order α.25 Order / / / / / / TOC Table 4. Absolute errors and convergence orders for predictor-corrector schemes (6-(7 at time t with different α.
19 372 Numerical Modelling 8 Will-be-set-by-IN-TECH p T Computing value Absolute errors Relative error (% Table 5. Error behavior versus the variation of p andt (the definition of p andt are given in (3 at time t 5 with exact (analytical value y(5245, fractional order α.5, step length h /8. and the errors of the improved predictor-corrector scheme (6-(7 are improved significantly compared with the Tables 3 and 4 given in (Diethelm et al., 24. Table 5 shows the computing values, the absolute numerical errors, and the relative numerical errors by using the scheme (32 and (33 to compute Example 2 for different values of p and T defined in (3. According to the numerical results we can see that the computing errors are generally acceptable in engineering when the computational cost is greatly minimized, especially the computing error is not sensitive to the value of p. On the other hand, this numerical example also illuminates that the algorithm is numerically stable. 5. Conclusions We briefly review the numerical techniques for efficiently solving fractional ordinary differential equations. The possible improvements for the fractional predictor-corrector algorithm are presented. Even though the algorithm is designed for scalar fractional ordinary differential equation, it can be easily extended to the systems. In fact, based on this algorithm, we have simulated the dynamics of fractional systems (Deng, 27c;d. In addition, the fractional predictor-corrector algorithm combining with the scheme of generating stochastic variables also works well for the stochastic fractional ordinary differential equations (Deng & Barkai, Acknowledgments This research was partly supported by the Program for New Century Excellent Talents in University under Grant No. NCET-9-438, the National Natural Science Foundation of China under Grant No. 867, and the Fundamental Research Funds for the Central Universities under Grant No. lzujbky-2-63 and No. lzujbky-22-k26.
20 Numerical Schemes for Fractional Ordinary Differential Equations Numerical Schemes for Fractional Ordinary Differential Equations References Bouchaud, J. & Georges, A. (99. Anomalous diffusion in disordered media: Statistical mechanisms models and physical applications. Phys. Rep., Vol. 95, No. 4-5, pp Deng, W.H. (27a. Numerical algorithm for the time fractional Fokker-Planck equation. J. Comput. Phys., Vol. 227, No. 2, pp Deng, W.H. (27b. Short memory principle and a predictor-corrector approach for fractional differential equations. J. Comput. Appl. Math., Vol. 26, No., pp Deng, W.H. (27c. Generating 3-D scroll grid attractors of fractional differential systems via stair function. Int. J. Bifurcation Chaos Appl. Sci. Eng., Vol. 7, No., pp Deng, W.H. (27d. Generalized synchronization in fractional order systems. Phys. Rev. E, Vol. 75, 562. Deng, W.H. & Barkai, E. (29. Ergodic properties of fractional Brownian-Langevin motion. Phys. Rev. E, Vol. 79, 2. Deng, W.H. (2. Smoothness and stability of the solutions for nonlinear fractional differential equations. Nonl. Anal.:TMA, Vol. 72, No. 2, pp Diethelm, K. & Ford, N.J. (22. Analysis of fractional differential equations. Nonl. Anal.:TMA, Vol. 265, No. 2, pp Diethelm, K., Ford, N.J. & Freed, A.D. (22b. A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dynam., Vol. 29, No. -4, pp Diethelm, K., Ford, N.J. & Freed, A.D. (24. Detailed error analysis for a fractional Adams method. Numer. Algorithms, Vol. 36, No., pp Ford, N.J. & Simpson, A.C. (2. The numerical solution of fractional differential equations: Speed versus accuracy. Numer. Algorithms, Vol. 26, No. 4, pp Heymans, N., & Podlubny, I. (25. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol. Acta, Vol. 45, No. 5, pp Hilfer, H. (2. Applications of Fractional Calculus in Physics, World Scientific Press, Singapore. Kirchner, J. W., Feng, X. & Neal, C. (2. Fractal stream chemistry and its implications for contaminant transport in catchments. Nature, Vol. 43, No. 3, pp Lubich, Ch. (985. Fractional linear multistep methods for Abel-Volterra integral equations of the second kind. Math. Comp., Vol. 45, No. 72, pp Lubich, Ch. (986. Discretized fractional calculus, SIAM J. Appl. Math., Vol. 7, No. 3, pp Metzler, R. & Klafter, J. (2. The random walk s guide to anomalous diffusion: a fractional dynamics approach. Phys. Rep., Vol. 339, No., pp Oldham, K.B. & Spanier, J. (974. The fractional calculus, Academic Press, New York. Podlubny, I. (999. Fractional differential equations, Academic Press, San Diego. Podlubny, I. (22. Geometric and Physical interpretation of fractional integration and fractional differentiation. Fract. Calc. Appl. Anal., Vol. 5, No. 4, pp Shlesinger, M. F., Zaslavsky, G.M. & Klafter, J. (993. Strange kinetics, Nature, Vol. 363 No. 6, pp Samko, S.G., Kilbas A.A. & Marichev, O.I. (993. Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach, London.
21 374 Numerical Modelling 2 Will-be-set-by-IN-TECH Zaslavsky, G.M. (22. Chaos, fractional kinetic, and anomalous transport, Phys. Rep., Vol. 37, No. 6, pp Mainrdi, F. (28. Fractional calculus and waves in linear viscoelasticity: an introduction to mathematical models. Imperial College Press, London.
22 Numerical Modelling Edited by Dr. Peep Miidla ISBN Hard cover, 398 pages Publisher InTech Published online 23, March, 22 Published in print edition March, 22 This book demonstrates applications and case studies performed by experts for professionals and students in the field of technology, engineering, materials, decision making management and other industries in which mathematical modelling plays a role. Each chapter discusses an example and these are ranging from wellknown standards to novelty applications. Models are developed and analysed in details, authors carefully consider the procedure for constructing a mathematical replacement of phenomenon under consideration. For most of the cases this leads to the partial differential equations, for the solution of which numerical methods are necessary to use. The term Model is mainly understood as an ensemble of equations which describe the variables and interrelations of a physical system or process. Developments in computer technology and related software have provided numerous tools of increasing power for specialists in mathematical modelling. One finds a variety of these used to obtain the numerical results of the book. How to reference In order to correctly reference this scholarly work, feel free to copy and paste the following: Weihua Deng and Can Li (22. Numerical Schemes for Fractional Ordinary Differential Equations, Numerical Modelling, Dr. Peep Miidla (Ed., ISBN: , InTech, Available from: InTech Europe University Campus STeP Ri Slavka Krautzeka 83/A 5 Rijeka, Croatia Phone: 385 ( Fax: 385 ( InTech China Unit 45, Office Block, Hotel Equatorial Shanghai No.65, Yan An Road (West, Shanghai, 24, China Phone: Fax:
23 22 The Author(s. Licensee IntechOpen. This is an open access article distributed under the terms of the Creative Commons Attribution 3. License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Numerical Detection of the Lowest Efficient Dimensions for Chaotic Fractional Differential Systems
The Open Mathematics Journal, 8, 1, 11-18 11 Open Access Numerical Detection of the Lowest Efficient Dimensions for Chaotic Fractional Differential Systems Tongchun Hu a, b, and Yihong Wang a, c a Department
More informationHigher order numerical methods for solving fractional differential equations
Noname manuscript No. will be inserted by the editor Higher order numerical methods for solving fractional differential equations Yubin Yan Kamal Pal Neville J Ford Received: date / Accepted: date Abstract
More informationNumerical solution of the Bagley Torvik equation. Kai Diethelm & Neville J. Ford
ISSN 1360-1725 UMIST Numerical solution of the Bagley Torvik equation Kai Diethelm & Neville J. Ford Numerical Analysis Report No. 378 A report in association with Chester College Manchester Centre for
More informationA computationally effective predictor-corrector method for simulating fractional order dynamical control system
ANZIAM J. 47 (EMA25) pp.168 184, 26 168 A computationally effective predictor-corrector method for simulating fractional order dynamical control system. Yang F. Liu (Received 14 October 25; revised 24
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 108,000 1.7 M Open access books available International authors and editors Downloads Our
More informationA Survey of Numerical Methods in Fractional Calculus
A Survey of Numerical Methods in Fractional Calculus Fractional Derivatives in Mechanics: State of the Art CNAM Paris, 17 November 2006 Kai Diethelm diethelm@gns-mbh.com, k.diethelm@tu-bs.de Gesellschaft
More informationComparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations
Comparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations Farhad Farokhi, Mohammad Haeri, and Mohammad Saleh Tavazoei Abstract This paper is a result of comparison of some
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4, 6, M Open access books available International authors and editors Downloads Our authors are
More informationHOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION. 1. Introduction
Fractional Differential Calculus Volume 1, Number 1 (211), 117 124 HOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION YANQIN LIU, ZHAOLI LI AND YUEYUN ZHANG Abstract In this paper,
More informationFRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS
FRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS L. Boyadjiev*, B. Al-Saqabi** Department of Mathematics, Faculty of Science, Kuwait University *E-mail: boyadjievl@yahoo.com **E-mail:
More informationA truncation regularization method for a time fractional diffusion equation with an in-homogeneous source
ITM Web of Conferences, 7 18) ICM 18 https://doi.org/1.151/itmconf/187 A truncation regularization method for a time fractional diffusion equation with an in-homogeneous source Luu Vu Cam Hoan 1,,, Ho
More informationEXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1. Yong Zhou. Abstract
EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1 Yong Zhou Abstract In this paper, the initial value problem is discussed for a system of fractional differential
More informationIMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES
Dynamic Systems and Applications ( 383-394 IMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES M ANDRIĆ, J PEČARIĆ, AND I PERIĆ Faculty
More informationSmart Bolometer: Toward Monolithic Bolometer with Smart Functions
Smart Bolometer: Toward Monolithic Bolometer with Smart Functions Matthieu Denoual, Gilles Allègre, Patrick Attia, Olivier De Sagazan To cite this version: Matthieu Denoual, Gilles Allègre, Patrick Attia,
More informationACTA UNIVERSITATIS APULENSIS No 20/2009 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS. Wen-Hua Wang
ACTA UNIVERSITATIS APULENSIS No 2/29 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS Wen-Hua Wang Abstract. In this paper, a modification of variational iteration method is applied
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationMulti-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
International Journal of Difference Equations ISSN 0973-6069, Volume 0, Number, pp. 9 06 205 http://campus.mst.edu/ijde Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4,000 116,000 120M Open access books available International authors and editors Downloads Our
More informationResearch Article New Method for Solving Linear Fractional Differential Equations
International Differential Equations Volume 2011, Article ID 814132, 8 pages doi:10.1155/2011/814132 Research Article New Method for Solving Linear Fractional Differential Equations S. Z. Rida and A. A.
More informationOn boundary value problems for fractional integro-differential equations in Banach spaces
Malaya J. Mat. 3425 54 553 On boundary value problems for fractional integro-differential equations in Banach spaces Sabri T. M. Thabet a, and Machindra B. Dhakne b a,b Department of Mathematics, Dr. Babasaheb
More informationThe Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers
The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers Chengbin Ma and Yoichi Hori Electrical Control System Engineering, Information & System Division, Institute of Industrial
More informationDIfferential equations of fractional order have been the
Multistage Telescoping Decomposition Method for Solving Fractional Differential Equations Abdelkader Bouhassoun Abstract The application of telescoping decomposition method, developed for ordinary differential
More informationHandling the fractional Boussinesq-like equation by fractional variational iteration method
6 ¹ 5 Jun., COMMUN. APPL. MATH. COMPUT. Vol.5 No. Å 6-633()-46-7 Handling the fractional Boussinesq-like equation by fractional variational iteration method GU Jia-lei, XIA Tie-cheng (College of Sciences,
More informationBoundary value problems for fractional differential equations with three-point fractional integral boundary conditions
Sudsutad and Tariboon Advances in Difference Equations 212, 212:93 http://www.advancesindifferenceequations.com/content/212/1/93 R E S E A R C H Open Access Boundary value problems for fractional differential
More informationResearch Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
International Differential Equations Volume 2010, Article ID 764738, 8 pages doi:10.1155/2010/764738 Research Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
More informationA generalized Gronwall inequality and its application to a fractional differential equation
J. Math. Anal. Appl. 328 27) 75 8 www.elsevier.com/locate/jmaa A generalized Gronwall inequality and its application to a fractional differential equation Haiping Ye a,, Jianming Gao a, Yongsheng Ding
More informationSOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER
Dynamic Systems and Applications 2 (2) 7-24 SOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER P. KARTHIKEYAN Department of Mathematics, KSR College of Arts
More informationNumerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes. Kai Diethelm. International Symposium on Fractional PDEs June 3 5, 2013
Numerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes Gesellschaft für numerische Simulation mbh Braunschweig AG Numerik Institut Computational Mathematics Technische Universität
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationA Fractional Spline Collocation Method for the Fractional-order Logistic Equation
A Fractional Spline Collocation Method for the Fractional-order Logistic Equation Francesca Pitolli and Laura Pezza Abstract We construct a collocation method based on the fractional B-splines to solve
More informationAnalytic solution of fractional integro-differential equations
Annals of the University of Craiova, Mathematics and Computer Science Series Volume 38(1), 211, Pages 1 1 ISSN: 1223-6934 Analytic solution of fractional integro-differential equations Fadi Awawdeh, E.A.
More informationThis work has been submitted to ChesterRep the University of Chester s online research repository.
This work has been submitted to ChesterRep the University of Chester s online research repository http://chesterrep.openrepository.com Author(s): Kai Diethelm; Neville J Ford Title: Volterra integral equations
More informationBifurcations of Fractional-order Diffusionless Lorenz System
EJTP 6, No. 22 (2009) 123 134 Electronic Journal of Theoretical Physics Bifurcations of Fractional-order Diffusionless Lorenz System Kehui Sun 1,2 and J. C. Sprott 2 1 School of Physics Science and Technology,
More informationOn the fractional-order logistic equation
Applied Mathematics Letters 20 (2007) 817 823 www.elsevier.com/locate/aml On the fractional-order logistic equation A.M.A. El-Sayed a, A.E.M. El-Mesiry b, H.A.A. El-Saka b, a Faculty of Science, Alexandria
More informationLocal Polynomial Smoother for Solving Bagley-Torvik Fractional Differential Equations
Preprints (wwwpreprintsorg) NOT PEER-REVIEWED Posted: 3 August 216 doi:12944/preprints2168231v1 Article Local Polynomial Smoother for Solving Bagley-Torvik Fractional Differential Equations Tianshun Yan
More informationOn the Finite Caputo and Finite Riesz Derivatives
EJTP 3, No. 1 (006) 81 95 Electronic Journal of Theoretical Physics On the Finite Caputo and Finite Riesz Derivatives A. M. A. El-Sayed 1 and M. Gaber 1 Faculty of Science University of Alexandria, Egypt
More informationConstruction of a New Fractional Chaotic System and Generalized Synchronization
Commun. Theor. Phys. (Beijing, China) 5 (2010) pp. 1105 1110 c Chinese Physical Society and IOP Publishing Ltd Vol. 5, No. 6, June 15, 2010 Construction of a New Fractional Chaotic System and Generalized
More informationA generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives
A generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives Deliang Qian Ziqing Gong Changpin Li Department of Mathematics, Shanghai University,
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
More informationPositive solutions for a class of fractional boundary value problems
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 1, 1 17 ISSN 1392-5113 http://dx.doi.org/1.15388/na.216.1.1 Positive solutions for a class of fractional boundary value problems Jiafa Xu a, Zhongli
More informationMahmoud M. El-Borai a, Abou-Zaid H. El-Banna b, Walid H. Ahmed c a Department of Mathematics, faculty of science, Alexandria university, Alexandria.
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:01 52 On Some Fractional-Integro Partial Differential Equations Mahmoud M. El-Borai a, Abou-Zaid H. El-Banna b, Walid H. Ahmed c
More informationElena Gogovcheva, Lyubomir Boyadjiev 1 Dedicated to Professor H.M. Srivastava, on the occasion of his 65th Birth Anniversary Abstract
FRACTIONAL EXTENSIONS OF JACOBI POLYNOMIALS AND GAUSS HYPERGEOMETRIC FUNCTION Elena Gogovcheva, Lyubomir Boyadjiev 1 Dedicated to Professor H.M. Srivastava, on the occasion of his 65th Birth Anniversary
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationExistence of Solutions for Nonlocal Boundary Value Problems of Nonlinear Fractional Differential Equations
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 7, Number 1, pp. 31 4 (212) http://campus.mst.edu/adsa Existence of Solutions for Nonlocal Boundary Value Problems of Nonlinear Fractional
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationAbstract We paid attention to the methodology of two integral
Comparison of Homotopy Perturbation Sumudu Transform method and Homotopy Decomposition method for solving nonlinear Fractional Partial Differential Equations 1 Rodrigue Batogna Gnitchogna 2 Abdon Atangana
More informationEXACT TRAVELING WAVE SOLUTIONS FOR NONLINEAR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS USING THE IMPROVED (G /G) EXPANSION METHOD
Jan 4. Vol. 4 No. 7-4 EAAS & ARF. All rights reserved ISSN5-869 EXACT TRAVELIN WAVE SOLUTIONS FOR NONLINEAR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS USIN THE IMPROVED ( /) EXPANSION METHOD Elsayed M.
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationEXISTENCE AND UNIQUENESS OF POSITIVE SOLUTIONS TO HIGHER-ORDER NONLINEAR FRACTIONAL DIFFERENTIAL EQUATION WITH INTEGRAL BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 212 (212), No. 234, pp. 1 11. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE AND UNIQUENESS
More informationFrom Unstructured 3D Point Clouds to Structured Knowledge - A Semantics Approach
From Unstructured 3D Point Clouds to Structured Knowledge - A Semantics Approach Christophe Cruz, Helmi Ben Hmida, Frank Boochs, Christophe Nicolle To cite this version: Christophe Cruz, Helmi Ben Hmida,
More informationThe Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers
Nonlinear Dynamics 38: 171 18, 24. C 24 Kluwer Academic Publishers. Printed in the Netherlands. The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers CHENGBIN MA and YOICHI
More informationExact Solutions of Fractional-Order Biological Population Model
Commun. Theor. Phys. (Beijing China) 5 (009) pp. 99 996 c Chinese Physical Society and IOP Publishing Ltd Vol. 5 No. 6 December 15 009 Exact Solutions of Fractional-Order Biological Population Model A.M.A.
More informationStability Analysis and Numerical Solution for. the Fractional Order Biochemical Reaction Model
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 11, 51-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.6531 Stability Analysis and Numerical Solution for the Fractional
More informationA SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION
THERMAL SCIENCE: Year 7, Vol., Suppl., pp. S-S7 S A SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION by Wei XU a, Wen CHEN a*, Ying-Jie LIANG a*, and Jose WEBERSZPIL b a State Key Laboratory
More informationSolving nonlinear fractional differential equation using a multi-step Laplace Adomian decomposition method
Annals of the University of Craiova, Mathematics and Computer Science Series Volume 39(2), 2012, Pages 200 210 ISSN: 1223-6934 Solving nonlinear fractional differential equation using a multi-step Laplace
More informationOn the numerical solutions of some fractional ordinary differential equations by fractional Adams-Bashforth-Moulton method
Open Math. 215; 13: 547 556 Open Mathematics Open Access Research Article Haci Mehmet Baskonus* and Hasan Bulut On the numerical solutions of some fractional ordinary differential equations by fractional
More informationFractional Evolution Integro-Differential Systems with Nonlocal Conditions
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 5, Number 1, pp. 49 6 (21) http://campus.mst.edu/adsa Fractional Evolution Integro-Differential Systems with Nonlocal Conditions Amar
More informationRIEMANN-LIOUVILLE FRACTIONAL COSINE FUNCTIONS. = Au(t), t > 0 u(0) = x,
Electronic Journal of Differential Equations, Vol. 216 (216, No. 249, pp. 1 14. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu RIEMANN-LIOUVILLE FRACTIONAL COSINE FUNCTIONS
More informationResearch Article Solving Fractional-Order Logistic Equation Using a New Iterative Method
International Differential Equations Volume 2012, Article ID 975829, 12 pages doi:10.1155/2012/975829 Research Article Solving Fractional-Order Logistic Equation Using a New Iterative Method Sachin Bhalekar
More informationV. G. Gupta 1, Pramod Kumar 2. (Received 2 April 2012, accepted 10 March 2013)
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9(205 No.2,pp.3-20 Approimate Solutions of Fractional Linear and Nonlinear Differential Equations Using Laplace Homotopy
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationFractional Diffusion Theory and Applications Part II
Fractional Diffusion Theory and Applications Part II p. 1/2 Fractional Diffusion Theory and Applications Part II 22nd Canberra International Physics Summer School 28 Bruce Henry (Trevor Langlands, Peter
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
More informationResearch Article Applying GG-Convex Function to Hermite-Hadamard Inequalities Involving Hadamard Fractional Integrals
International Journal of Mathematics and Mathematical Sciences Volume 4, Article ID 3635, pages http://dx.doi.org/.55/4/3635 Research Article Applying GG-Convex Function to Hermite-Hadamard Inequalities
More informationExistence of Minimizers for Fractional Variational Problems Containing Caputo Derivatives
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 3 12 (2013) http://campus.mst.edu/adsa Existence of Minimizers for Fractional Variational Problems Containing Caputo
More information2 Numerical Methods for Initial Value Problems
Numerical Analysis of Differential Equations 44 2 Numerical Methods for Initial Value Problems Contents 2.1 Some Simple Methods 2.2 One-Step Methods Definition and Properties 2.3 Runge-Kutta-Methods 2.4
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationHomotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders
Homotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders Yin-Ping Liu and Zhi-Bin Li Department of Computer Science, East China Normal University, Shanghai, 200062, China Reprint
More informationCRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 2. Jan. 2012, No. 2, pp. 1-9. ISSN: 2090-5858. http://www.fcaj.webs.com/ CRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION
More informationChapter 1 Introduction
Chapter 1 Introduction During recent years the interest of physicists in nonlocal field theories has been steadily increasing. The main reason for this development is the expectation that the use of these
More informationSIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE. Yuriy Povstenko. Abstract
SIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE Yuriy Povstenko Abstract The time-fractional diffusion-wave equation is considered in a half-plane. The Caputo fractional derivative
More informationHIGHER ORDER METHODS. There are two principal means to derive higher order methods. b j f(x n j,y n j )
HIGHER ORDER METHODS There are two principal means to derive higher order methods y n+1 = p j=0 a j y n j + h p j= 1 b j f(x n j,y n j ) (a) Method of Undetermined Coefficients (b) Numerical Integration
More informationApplied Mathematics Letters. A reproducing kernel method for solving nonlocal fractional boundary value problems
Applied Mathematics Letters 25 (2012) 818 823 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml A reproducing kernel method for
More informationEXISTENCE OF SOLUTIONS TO FRACTIONAL ORDER ORDINARY AND DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS
Electronic Journal of Differential Equations, Vol. 211 (211), No. 9, pp. 1 11. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF SOLUTIONS
More informationExact Solution of Some Linear Fractional Differential Equations by Laplace Transform. 1 Introduction. 2 Preliminaries and notations
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.16(213) No.1,pp.3-11 Exact Solution of Some Linear Fractional Differential Equations by Laplace Transform Saeed
More informationA novel difference schemes for analyzing the fractional Navier- Stokes equations
DOI: 0.55/auom-207-005 An. Şt. Univ. Ovidius Constanţa Vol. 25(),207, 95 206 A novel difference schemes for analyzing the fractional Navier- Stokes equations Khosro Sayevand, Dumitru Baleanu, Fatemeh Sahsavand
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationEquilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models
J. Math. Anal. Appl. 325 (2007) 542 553 www.elsevier.com/locate/jmaa Equilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models E. Ahmed a, A.M.A. El-Sayed
More informationarxiv:cs/ v1 [cs.na] 23 Aug 2004
arxiv:cs/0408053v1 cs.na 23 Aug 2004 WEIGHTED AVERAGE FINITE DIFFERENCE METHODS FOR FRACTIONAL DIFFUSION EQUATIONS Santos B. Yuste,1 Departamento de Física, Universidad de Extremadura, E-0071 Badaoz, Spain
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More informationNEW GENERAL FRACTIONAL-ORDER RHEOLOGICAL MODELS WITH KERNELS OF MITTAG-LEFFLER FUNCTIONS
Romanian Reports in Physics 69, 118 217 NEW GENERAL FRACTIONAL-ORDER RHEOLOGICAL MODELS WITH KERNELS OF MITTAG-LEFFLER FUNCTIONS XIAO-JUN YANG 1,2 1 State Key Laboratory for Geomechanics and Deep Underground
More informationExistence of triple positive solutions for boundary value problem of nonlinear fractional differential equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 5, No. 2, 217, pp. 158-169 Existence of triple positive solutions for boundary value problem of nonlinear fractional differential
More informationarxiv: v1 [math.ca] 28 Feb 2014
Communications in Nonlinear Science and Numerical Simulation. Vol.18. No.11. (213) 2945-2948. arxiv:142.7161v1 [math.ca] 28 Feb 214 No Violation of the Leibniz Rule. No Fractional Derivative. Vasily E.
More informationarxiv: v1 [math.na] 8 Jan 2019
arxiv:190102503v1 [mathna] 8 Jan 2019 A Numerical Approach for Solving of Fractional Emden-Fowler Type Equations Josef Rebenda Zdeněk Šmarda c 2018 AIP Publishing This article may be downloaded for personal
More informationTHE objective of this paper is to introduce a comparative
Proceedings of the World Congress on Engineering 23 Vol I, WCE 23, July 3-5, 23, London, U.K. Numerical Solution of Fractional Integro-differential Equation by Using Cubic B-spline Wavelets Khosrow Maleknejad,
More informationResearch Article A New Fractional Integral Inequality with Singularity and Its Application
Abstract and Applied Analysis Volume 212, Article ID 93798, 12 pages doi:1.1155/212/93798 Research Article A New Fractional Integral Inequality with Singularity and Its Application Qiong-Xiang Kong 1 and
More informationarxiv: v1 [cond-mat.stat-mech] 23 Apr 2010
Fractional Fokker-Planck Equations for Subdiffusion with Space-and-Time-Dependent Forces. B.. Henry Department of Applied Mathematics, School of Mathematics and Statistics, University of New South Wales,
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationOn Local Asymptotic Stability of q-fractional Nonlinear Dynamical Systems
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 11, Issue 1 (June 2016), pp 174-183 Applications and Applied Mathematics: An International Journal (AAM) On Local Asymptotic Stability
More informationA Cauchy Problem for Some Local Fractional Abstract Differential Equation with Fractal Conditions
From the SelectedWorks of Xiao-Jun Yang 2013 A Cauchy Problem for Some Local Fractional Abstract Differential Equation with Fractal Conditions Yang Xiaojun Zhong Weiping Gao Feng Available at: https://works.bepress.com/yang_xiaojun/32/
More informationDynamic Response and Oscillating Behaviour of Fractionally Damped Beam
Copyright 2015 Tech Science Press CMES, vol.104, no.3, pp.211-225, 2015 Dynamic Response and Oscillating Behaviour of Fractionally Damped Beam Diptiranjan Behera 1 and S. Chakraverty 2 Abstract: This paper
More informationNumerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag Leffler Functions
mathematics Article Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag Leffler Functions Marina Popolizio ID Dipartimento di Matematica e Fisica Ennio De Giorgi,
More informationA new Mittag-Leffler function undetermined coefficient method and its applications to fractional homogeneous partial differential equations
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017), 4515 4523 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa A new Mittag-Leffler function
More informationHybrid Functions Approach for the Fractional Riccati Differential Equation
Filomat 30:9 (2016), 2453 2463 DOI 10.2298/FIL1609453M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Hybrid Functions Approach
More informationAPPLICATION OF HAAR WAVELETS IN SOLVING NONLINEAR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS
APPLICATION OF HAAR WAVELETS IN SOLVING NONLINEAR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS H. SAEEDI DEPARTMENT OF MATHEMATICS, SAHID BAHONAR UNIVERSITY OF KERMAN, IRAN, 76169-14111. E-MAILS:
More informationAN EFFICIENT PARALLEL ALGORITHM FOR THE NUMERICAL SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm 1,2
RESEARCH PAPER AN EFFICIENT PARALLEL ALGORITHM FOR THE NUMERICAL SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm 1,2 Abstract The numerical solution of differential equations of fractional order
More informationMULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS
MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS Hossein Jafari & M. A. Firoozjaee Young Researchers club, Islamic Azad University, Jouybar Branch, Jouybar, Iran
More informationA Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional Derivative
Commun. Theor. Phys. 68 2017 301 308 Vol. 68, No. 3, September 1, 2017 A Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional
More information