A computationally effective predictor-corrector method for simulating fractional order dynamical control system

Size: px
Start display at page:

Download "A computationally effective predictor-corrector method for simulating fractional order dynamical control system"

Transcription

1 ANZIAM J. 47 (EMA25) pp , A computationally effective predictor-corrector method for simulating fractional order dynamical control system. Yang F. Liu (Received 14 October 25; revised 24 July 26) Abstract Multi-order fractional differential equations are applied to fractional order dynamical controlled systems. The multi-order fractional differential equation is transferred into a system of fractional order differential equations. A new computationally effective fractional predictorcorrector method is proposed for simulating the fractional order systems and controllers. A detailed error analysis is derived. Finally, we give some numerical examples. School of Mathematical Sciences, Xiamen University, Xiamen 3615, hina. mailto:morningsail@163.com School of Mathematical Sciences, Queensland University of Technology, Queensland 41, Australia. mailto:f.liu@qut.edu.au. See for this article, c Austral. Mathematical Soc. 26. Published August 1, 26. ISSN

2 ANZIAM J. 47 (EMA25) pp , ontents 1 Introduction Basic ideas and lemmas A fractional predictor-corrector technique Error analysis for the fractional predictor-corrector method Numerical Results onclusions 182 References Introduction The Fractional alculus is, simultaneously, a new and old research issue. This contradiction comes from the fact that most scientists are not aware of its existence and, often, discover that this new tool provides fruitful perspectives in their studies. For three centuries the theory of fractional derivatives developed mainly as a pure theoretical field of mathematics useful only for mathematicians [1]. However, in the last few decades many have pointed out that derivatives and integrals of non-integer order are very suitable for the description of the properties of various real materials. It has been shown that new fractional order models are more adequate than previously used integer order models [1]. In the past two decades, fractional order control systems attracted the attention of many researchers lately [2]. However, because of the absence of appropriate mathematical methods, fractional order dynamical systems were studied only marginally in theory and practice of control systems.

3 1 Introduction 17 W(S) + E(S) U(s) G c(s) _ G(s) Y(s) Figure 1: Simple unity-feedback control system Let us consider the simple unity-feedback control system [1] shown in Figure 1, where G(s) is the transfer function of the controlled system, G c (s) is the transfer of the controller, W (s) is an input, E(s) is an error, U(s) is the controller s output, and Y (s) is the system s output. ontrary to the traditional model, we consider transfer functions of arbitrary real order. Such systems are called fractional order control systems, which are better described by fractional order mathematical models. Now consider the fractional order transfer function (fotf) [1]: G n (s) = 1 a n s αn + a n 1 s α n a1 s α 1 + a s α, (1) where α n > α n 1 > > α 1 > α, < α n α n 1 < 1, and α k (k =, 1,..., n) are arbitrary real numbers. In the time domain, the fotf (1) corresponds to the (n + 1)-term fractional order differential equation a n D αn y(t) + a n 1 D α n 1 y(t) + + a 1 D α 1 y(t) + a D α y(t) = u(t), (2) where D α y = Dt α y = J m α y (m) (t) is the aputo fractional derivatives (in time) of order α: Dt α y(t) = t 1 Γ(m α) d m y(t) dt, m N, m y (m) (τ) dτ, m 1 < α < m, (t τ) α+1 m where m is the integer defined by the relation m 1 < α < m, and J q is the (3)

4 1 Introduction 171 fractional integral operator J q y(t) = 1 Γ(q) t (t τ) q 1 y(τ) dτ. (4) Numerical solution of a fractional order differential equation is difficult to obtain accurately using standard discretization methods. Hu and Liu [3] considered a four term fractional differential equation corresponding to a fractional order controlled system. The existence and uniqueness of the solution of this model are given, and they proposed three numerical methods for the fractional order control systems. Shen and Liu considered [4] a fractional order Bagley Torvik equation, and proved the existence and uniqueness of solution. The analytical solution of a fractional order Bagley Torvik equation is also derived by the corresponding Green s function. Using the relationship between the Riemann Liouville definition and the Gronwald Letnikov definition, Shen and Liu proposed a computationally effective method for the fractional order Bagley Torvik Equation. Diethelm et al. [5] proposed some techniques for solving fractional order ordinary differential equations and analysed errors. Lin and Liu [6] also proposed a fractional order numerical method for the fractional relaxation equation, and the consistence, convergence and stability of this method is proven. Liu et al. [7] proposed some efficient fractional numerical methods for solving fractional partial differential equation. A special multi-order fractional differential equation was also considered by Miller and Ross [8]. The multi-order fractional differential equations were written as an equivalent system of fractional order differential equations. However, numerical methods for solving fractional order differential systems are limited. Therefore, a new numerical strategy is important in solving these fractional order systems. We propose a new technique for solving fractional order differential equation systems and are apply the technique to simulate a fractional order control system. Some basic ideas and lemmas are presented in Section 2. A computationally effective fractional predictor-corrector method is proposed in Section 3.

5 1 Introduction 172 A detailed error analysis is derived in Section 4. Finally, Section 5 gives some numerical results. 2 Basic ideas and lemmas We are concerned with providing good quality methods for the solution of multi-order fractional differential equations of the general form a n D αn y(t) + a n 1 D α n 1 y(t) + + a 1 D α 1 y(t) + a D α y(t) = f(t), y k () = y (k), k =, 1,..., α n (5) where α n > α n 1 > > α 1 > α, < α n α n 1 < 1, α k (k =, 1,..., n) are arbitrary real numbers, a k (k =, 1,..., n) are arbitrary constants. The basic analytical results on existence and uniqueness of solutions to fractional differential equations are given in [1, 8]. Miller and Ross [8] gave a method for calculating the analytic solution to a multi-term fractional differential equation of the form a n D nα y(t) + a n 1 D (n 1)α y(t) + + a 1 D α y(t) + a D y(t) =, (6) where α = 1/M, M N. They interpret the fractional operator D iα as D α applied i times. Under this interpretation they describe the multi-term fractional differential equation as a sequential fractional differential equation. For the sake of simplicity we assume that α i = i/2 = m i N, if i is an even number; (i 1)/2 < α i < (i + 1)/2, if i is an odd number; that is, there is at most one non-integer order derivative between successive integer orders. Lemma 1 Let y k [, T ] for some T > and some k N, and let β / N such that < β < k, then a D β t y() =.

6 2 Basic ideas and lemmas 173 Proof: See [9]. Lemma 2 The differentiation operators a Dt α f(t) and a Dt m f(t) satisfy the interchange rule: ( Dt α Dt m f(t) ) ( = Dt m Dt α f(t) ) = Dt α+m f(t), (7) f (s) () =, s = n, n + 1,..., m, m =, 1,..., n 1 < α < n. Proof: See [1]. Let y αn [, T ] for some T >. Using Lemma 1, we have D α 2i α 2i 1 t D α 2i 1 t y(t) = D α 2i t y, D α 2i+1 α 2i t D α 2i t y(t) = D α 2i+1 t y. (8) Now we rewrite the multi-order fractional differential equations (5) in the form of a system of fractional order differential equations: D β 1 t x 1 (t) = D α 1 t x 1 (t) = x 2 (t),. D β n 1 t D βn t x n 1 (t) = D α n 1 α n 2 t x n 1 (t) = x n (t), x n (t) = D αn α n 1 t x n (t) (9) = 1 a n [f(t) a x 1 a 1 x 2 a n 1 x n ], with initial equations x 1 () = x (1) = y (), x 2 () = x (2) =,..., { x i () = x (i) (k) y =, if i = 2k + 1, i n. (1), if i = 2k, Using Lemma 1 and Lemma 2, we obtain the following theorem:

7 2 Basic ideas and lemmas 174 Theorem 3 The multi-order fractional differential equations (5) is equivalent to the system of equations (9) with the initial conditions (1). 3 A fractional predictor-corrector technique In this section, numerical technique for simulate fractional order control systems (1) are presented. Firstly, the fractional order differential systems is decoupled, which is equivalent to solving D β 1 t x 1 (t) = g 1 (t, x 1 ),. D βn t x n (t) = g n (t, x n ). (11) Secondly, we consider a computationally effective fractional predictorcorrector method for solving the following initial-value problem: where < β i < 1. D β i t x i (t) = g i (t, x i ), x i () = x (i), i = 1, 2,..., n, (12) It is well known that the initial-value problem (12) is equivalent to the Volterra integral equation x i (t) = x (i) + 1 Γ(β i ) t (t τ) β i 1 [g i (τ, x i (τ))] dτ. (13) For the sake of simplicity, we assume that we are working on a uniform grid t j = jτ, j =, 1,..., Mτ = T. The issue of stability is very important when implementing the method on a computer in finite precision arithmetic because we must take into account effects introduced by rounding errors. It is known that the classical Adams Bashforth Moulton method for first order ordinary differential equations is

8 3 A fractional predictor-corrector technique 175 a reasonable and practically useful compromise in the sense that its stability properties allow for a safe application to mildly stiff equations without undue propagation of rounding error, whereas the implementation does not require extremely time consuming elements [1]. Thus, a fractional Adams Bashforth method and a fractional Adams Moulton method are chosen as our predictor and corrector formulas. The predictor x P i,k+1 method [5, 11]: is determined by the fractional Adams Bashforth where x P i,k+1 = x (i) + 1 Γ(α) b β j,k+1 = τ β k b β j j,k+1 g i(t j, x i,j ), (14) β [(k + 1 j)β (k j) β ]. (15) The fractional Adams Moulton method determines the corrector formula [5, 11]: ( k ) x i,k+1 = x (i) + 1 a β i j,k+1 Γ(β i ) g i(t j, x i,j ) + a β i k+1,k+1 g i(t k+1, x P i,k+1), (16) where a β j,k+1 = τ β β(β + 1) k β+1 (k β)(k + 1) β, j =, (k j + 2) β+1 + (k j) β+1 2(k j + 1) β+1, 1 j k, 1, j = k + 1. (17) Therefore, we obtain the following new fractional predictor-corrector method for solving the fractional order differential systems (9).

9 3 A fractional predictor-corrector technique 176 Fractional Predictor formulas: x P i,k+1 = x (i) + 1 k b β i j,k+1 Γ(β i ) i+1,j, i = 1, 2,..., n 1, (18) x P n,k+1 = x (i) + 1 Γ(β n ) k b βn j,k+1 1 [u(t j ) a x 1,j a n 1 x n,j ].(19) a n Fractional orrector formulas: ( k ) x i,k+1 = x (i) + 1 a β i j,k+1 Γ(β i ) x i+1,j + a β k+1,k+1 xp i+1,k+1 i = 1, 2,..., n 1, (2) { k x n,k+1 = x (i) a βn j,k+1 [u(t j ) a x 1,j a n 1 x n,j ] Γ(β n ) a n + a βn k+1,k+1 1 [ ] } u(tk+1 ) a x P 1,k+1 a n 1 x P n,k+1 (21) a n, 4 Error analysis for the fractional predictor-corrector method In this section we present the theorems concerning the error of our fractional predictor-corrector method. Lemma 4 Let z 1 [, T ], then tk+1 (t k+1 t) β 1 z(t) dt k b β j,k+1 z(t j) 1 β z t β k+1 τ.

10 4 Error analysis for the fractional predictor-corrector method 177 Proof: See [5, 11]. Lemma 5 If z 2 [, T ], then there is a constant β depending only on β such that tk+1 k+1 (t k+1 t) β 1 z(t) dt a β j,k+1 z(t j) β z t β k+1 τ 2. (22) Proof: See [5, 11]. Theorem 6 If D β i t x i 2 [, T ], (i = 1, 2,..., n), then max j M 1 i n where q = 1 + min 1 i n β i. x i (t j ) x i,j = O(h q ) (23) Proof: Using given condition D β i t x i 2 [, T ], (i = 1, 2,..., n), Lemma 4 and Lemma 5, we have tk+1 k (t k+1 t) β i 1 D β i t x i (t) dt b β j,k+1 D β i t x i (t j ) 1t β i k+1 τ, (24) tk+1 k+1 (t k+1 t) β i 1 D β i t x i (t) dt a β j,k+1 D β i t x i (t j ) 2t β i k+1 τ 2. (25) We show that, for sufficiently small τ = T/M, max j M 1 i n x i (t j ) x i,j = O(τ q ), (26)

11 4 Error analysis for the fractional predictor-corrector method 178 where q = 1 + min 1 i n β i. The proof will be based on mathematical induction. In view of the given initial condition, the induction basis (j = ) is presupposed. Now assume that (26) is true for j =, 1,..., k (k M 1), that is, x i (t j ) x i,j = O(τ q ). (27) max j M 1 j n We must then prove the inequality also holds for j = k + 1. To do this, we first look at the error of the predictor solutions x P i,k+1, i = 1, 2,..., n. By construction of the predictor formulas, using (24), assumption (27), and k b β i j,k+1 = tk+1 (t k+1 t) β i 1 dt = 1 β i t β i k+1 1 β i T β i, we find that x i (t k+1 ) x P i,k+1 1T β i Γ(β i ) τ + T βi Γ(β i + 1) τ q, i = 1, 2,..., n 1,(28) x n (t k+1 ) x P n,k+1 1T βn Γ(β n ) τ + LT βn Γ(β n + 1) τ q. (29) Now we begin the analysis of the corrector error. For j = k + 1, arguing in a similar way to the above, by construction of the predictor formulas, using (25), assumption (27), and (28) (29), we find that x i (t k+1 ) x i,k+1 ( 2 T β i Γ(β i ) + ) T βi Γ(β i + 1) + 1 T βi Γ(β i )Γ(β i + 2) + T βi Γ(β i + 1)Γ(β i + 2) τ β i τ q τ q, i = 1, 2,..., n 1 ; (3) x n (t k+1 ) x n,k+1

12 4 Error analysis for the fractional predictor-corrector method 179 ( 2 T βn Γ(β n ) + τ q. T βn Γ(β n + 1) + 1 T βn Γ(β n )Γ(β n + 2) + T βn Γ(β n + 1)Γ(β n + 2) τ βn ) τ q (31) This completes the proof. Therefore, we obtain not only an approximation for the solution y(t) of the n term fractional order differential equation, but also approximations for its (aputo type) derivatives of order α i (i = 1, 2,..., n 1). Apart from this useful feature, the method has a rather simple structure that makes it very easy to implement. 5 Numerical Results onsider a fractional order controlled system with transfer function [1, 2] G 4 (s) = 1 a 4 s α 4 + a3 s α 3 + a2 s α 2 + a1 s α 1 + a, (32) where we take < α 1 < 1., α 2 = 1., 1. < α 3 < 2., α 4 = 2., a = 1, a 1 =, a 2 =, a 3 =.5, a 4 = 1. The fractional order transfer function (32) corresponds in the time domain to the five term fractional order differential equation a 4 D α 4 y(t) + a 3 D α 4 + a 2 D α 2 y(t) + a 1 D α 1 y(t) + a y(t) = u(t), (33) with the initial conditions y() = y () =, y () = y (1) = 1. (34) Now rewrite the multi-order fractional differential equations (32) as a

13 5 Numerical Results x 1 u(t)=1 1 Values.5 x x 4 x Time,t Figure 2: Unit-step response of the fractional order system with u(t) = 1. system of fractional order differential equations: D β 1 t x 1 (t) = D α 1 t x 1 (t) = x 2 (t), D β 2 t x 2 (t) = D α 2 α 1 t x 2 (t) = x 3 (t), D β 3 t x 3 (t) = D α 3 α 2 t x 3 (t) = x 4 (t), D β 4 t x 4 (t) = D α 4 α 3 t x 4 (t) = 1 [u(t) a x 1 a 1 x 2 a 2 x 3 a 3 x 4 ], a 4 with initial conditions (35) x 1 () = y () =, x 2 () =, x 3 () = y (1) = 1, x 4 () =. (36) Example 1. We take α 1 =.5, α 2 = 1., α 3 = 1.5, α 4 = 2..

14 5 Numerical Results x 1 u(t)=.5 x 2 Value x x Time,t Figure 3: Unit-step response of the fractional order system with u(t) =. Figures 2 and 3 show the simulation results with u(t) = 1 and u(t) =, respectively. The system exhibits behaviors of the solution and its derivatives of order.5, 1, 1.5. From Figures 2 and 3, see that the computed results are in excellent agreement with those of Wang et al. [12]. Example 2. We take α 1 =.5, α 2 = 1., 1.1 α 3 1.9, α 4 = 2.. Figure 4 displays the effect of 1.1 α (that is,.1 β 3.9) using predictor-corrector method for the fractional order control system with u(t) = 1. From Figure 4 shows that when β 3 increases from.1 to.9, the oscillation of solution also increases. Figures 2 4 show that the fractional predictor-corrector method provides a computationally effective method for simulating the behavior of the solution of a fractional order control system.

15 5 Numerical Results 182 Figure 4: Solution of fractional order system for.1 < β 3 <.9, u(t) = 1. 6 onclusions A fractional predictor-corrector method for simulating the fractional order control systems has been described and demonstrated. This method provides a computationally effective method for simulating the behavior of the solution of a fractional control system. A detailed error analysis is discussed. This numerical technique can be applied to simulate other fractional order differential system. Acknowledgments: This research has been supported by the National Natural Science Foundation of hina grant

16 6 onclusions 183 References [1] Podlubny, I., Fractional differential equations. New York: Academic Press, , 17, 172, 173, 179 [2] Petras, I. and Dorcak, L., Fractional-order control systems: modeling and simulation, fractional calculus and applied analysis, 6, 23, , 179 [3] Hu Y. and Liu F., Numerical methods for a fractional-order control system, Journal of Xiamen University (NATURAL Science), 44, 25, URL 171 [4] Shen S., Liu F., A computational effective method for fractional order Bagley Torvik equation, Journal of Xiamen University (NATURAL Science), 45, 24, URL 171 [5] Diethelm K., Ford N. J., and Freed Alen D., Detailed error analysis for a fractional adams method[m], Numerical Algorithms, Kluwer Academic Publishers, 24, , 175, 177 [6] Lin R., Liu F., Analysis of fractional-order numerical method for the fractional relaxation equation, omputational Mechanics, Tsinghua University and Springer-Verlag, 24, R [7] Liu F., Anh V., Turner I., Numerical solution of the space fractional Fokker Planck equation, J. omp. and Appl. Math., 166, 24,

17 References 184 [8] Miler K. S., Ross B., An introduction to the fractional calculus and fractional differential equations, , John Wiley, , 172 [9] Diethelm K., Ford N. J., Numerical solution of the Bagley Torvik equation, BIT 42, 22, [1] Hairer E., Wanner G. Solving ordinary differential equations II: stiff and differential-algebraic problems. Springer, Berlin URL 175 [11] Yang., Liu F., A fractional predictor-corrector method of the fractional relaxation-oscillation equation, Journal of Xiamen University, 44, 25, URL 175, 177 [12] Wang Z., ao G., Two system modeling methods using fractional calculus, J. System Simulation, 16, 24, URL 181

On the fractional-order logistic equation

On 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 information

POSITIVE SOLUTIONS FOR BOUNDARY VALUE PROBLEM OF SINGULAR FRACTIONAL FUNCTIONAL DIFFERENTIAL EQUATION

POSITIVE SOLUTIONS FOR BOUNDARY VALUE PROBLEM OF SINGULAR FRACTIONAL FUNCTIONAL DIFFERENTIAL EQUATION International Journal of Pure and Applied Mathematics Volume 92 No. 2 24, 69-79 ISSN: 3-88 (printed version); ISSN: 34-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/.2732/ijpam.v92i2.3

More information

A Survey of Numerical Methods in Fractional Calculus

A 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 information

Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients

Multi-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 information

Exact Solution of Some Linear Fractional Differential Equations by Laplace Transform. 1 Introduction. 2 Preliminaries and notations

Exact 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 information

DIfferential equations of fractional order have been the

DIfferential 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 information

multistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):

multistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP): MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem

More information

Time Fractional Wave Equation: Caputo Sense

Time Fractional Wave Equation: Caputo Sense Adv. Studies Theor. Phys., Vol. 6, 2012, no. 2, 95-100 Time Fractional Wave Equation: aputo Sense H. Parsian Department of Physics Islamic Azad University Saveh branch, Saveh, Iran h.parsian@iau-saveh.ac.ir

More information

Comparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations

Comparing 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 information

SMOOTHNESS 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 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 information

Higher order numerical methods for solving fractional differential equations

Higher 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 information

Applied Mathematics Letters

Applied 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 information

A generalized Gronwall inequality and its application to a fractional differential equation

A 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 information

On boundary value problems for fractional integro-differential equations in Banach spaces

On 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 information

We 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. International authors and editors 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

More information

Numerical Detection of the Lowest Efficient Dimensions for Chaotic Fractional Differential Systems

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 information

Numerical solution of the Bagley Torvik equation. Kai Diethelm & Neville J. Ford

Numerical 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 information

NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX

NUMERICAL 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 information

arxiv: v1 [math.na] 8 Jan 2019

arxiv: 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 information

Bifurcations of Fractional-order Diffusionless Lorenz System

Bifurcations 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 information

Picard,Adomian and Predictor-Corrector methods for integral equations of fractional order

Picard,Adomian and Predictor-Corrector methods for integral equations of fractional order Picard,Adomian and Predictor-Corrector methods for integral equations of fractional order WKZahra 1, MAShehata 2 1 Faculty of Engineering, Tanta University, Tanta, Egypt 2 Faculty of Engineering, Delta

More information

An Efficient Numerical Method for Solving. the Fractional Diffusion Equation

An 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 information

Existence of Ulam Stability for Iterative Fractional Differential Equations Based on Fractional Entropy

Existence of Ulam Stability for Iterative Fractional Differential Equations Based on Fractional Entropy Entropy 215, 17, 3172-3181; doi:1.339/e1753172 OPEN ACCESS entropy ISSN 199-43 www.mdpi.com/journal/entropy Article Existence of Ulam Stability for Iterative Fractional Differential Equations Based on

More information

A Numerical Scheme for Generalized Fractional Optimal Control Problems

A Numerical Scheme for Generalized Fractional Optimal Control Problems Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 11, Issue 2 (December 216), pp 798 814 Applications and Applied Mathematics: An International Journal (AAM) A Numerical Scheme for Generalized

More information

Nontrivial solutions for fractional q-difference boundary value problems

Nontrivial solutions for fractional q-difference boundary value problems Electronic Journal of Qualitative Theory of Differential Equations 21, No. 7, 1-1; http://www.math.u-szeged.hu/ejqtde/ Nontrivial solutions for fractional q-difference boundary value problems Rui A. C.

More information

MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS

MULTISTAGE 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 information

A COLLOCATION METHOD FOR SOLVING FRACTIONAL ORDER LINEAR SYSTEM

A COLLOCATION METHOD FOR SOLVING FRACTIONAL ORDER LINEAR SYSTEM J Indones Math Soc Vol 23, No (27), pp 27 42 A COLLOCATION METHOD FOR SOLVING FRACTIONAL ORDER LINEAR SYSTEM M Mashoof, AH Refahi Sheikhani 2, and H Saberi Najafi 3 Department of Applied Mathematics, Faculty

More information

Solving nonlinear fractional differential equation using a multi-step Laplace Adomian decomposition method

Solving 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 information

Homotopy Analysis Method for Nonlinear Differential Equations with Fractional Orders

Homotopy 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 information

Stability Ordinates of Adams Predictor-Corrector Methods

Stability Ordinates of Adams Predictor-Corrector Methods BIT manuscript No. will be inserted by the editor) Stability Ordinates of Adams Predictor-Corrector Methods Michelle L. Ghrist Bengt Fornberg Jonah A. Reeger Received: date / Accepted: date Abstract How

More information

A Wavelet Method for Solving Bagley-Torvik Equation

A Wavelet Method for Solving Bagley-Torvik Equation Copyright 2014 Tech Science Press CMES, vol.102, no.2, pp.169-182, 2014 A Wavelet Method for Solving Bagley-Torvik Equation Xiaomin Wang 1,2 Abstract: In this paper, an efficient and robust wavelet Laplace

More information

Math 128A Spring 2003 Week 12 Solutions

Math 128A Spring 2003 Week 12 Solutions Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems

More information

Exact Solutions of Fractional-Order Biological Population Model

Exact 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 information

9.6 Predictor-Corrector Methods

9.6 Predictor-Corrector Methods SEC. 9.6 PREDICTOR-CORRECTOR METHODS 505 Adams-Bashforth-Moulton Method 9.6 Predictor-Corrector Methods The methods of Euler, Heun, Taylor, and Runge-Kutta are called single-step methods because they use

More information

An Analytical Scheme for Multi-order Fractional Differential Equations

An Analytical Scheme for Multi-order Fractional Differential Equations Tamsui Oxford Journal of Mathematical Sciences 26(3) (2010) 305-320 Aletheia University An Analytical Scheme for Multi-order Fractional Differential Equations H. M. Jaradat Al Al Bayt University, Jordan

More information

Equilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models

Equilibrium 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 information

The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers

The 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 information

Computational Non-Polynomial Spline Function for Solving Fractional Bagely-Torvik Equation

Computational Non-Polynomial Spline Function for Solving Fractional Bagely-Torvik Equation Math. Sci. Lett. 6, No. 1, 83-87 (2017) 83 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.18576/msl/060113 Computational Non-Polynomial Spline Function for Solving Fractional

More information

Local Polynomial Smoother for Solving Bagley-Torvik Fractional Differential Equations

Local 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 information

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR FOURTH-ORDER BOUNDARY-VALUE PROBLEMS IN BANACH SPACES

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR FOURTH-ORDER BOUNDARY-VALUE PROBLEMS IN BANACH SPACES Electronic Journal of Differential Equations, Vol. 9(9), No. 33, pp. 1 8. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE AND UNIQUENESS

More information

A STUDY OF GENERALIZED ADAMS-MOULTON METHOD FOR THE SATELLITE ORBIT DETERMINATION PROBLEM

A STUDY OF GENERALIZED ADAMS-MOULTON METHOD FOR THE SATELLITE ORBIT DETERMINATION PROBLEM Korean J Math 2 (23), No 3, pp 27 283 http://dxdoiorg/568/kjm232327 A STUDY OF GENERALIZED ADAMS-MOULTON METHOD FOR THE SATELLITE ORBIT DETERMINATION PROBLEM Bum Il Hong and Nahmwoo Hahm Abstract In this

More information

SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD

SOLUTION 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 information

ON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD

ON 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 information

Numerical 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. 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 information

Differential equations with fractional derivative and universal map with memory

Differential equations with fractional derivative and universal map with memory IOP PUBLISHING JOURNAL OF PHYSIS A: MATHEMATIAL AND THEORETIAL J. Phys. A: Math. Theor. 42 (29) 46512 (13pp) doi:1.188/1751-8113/42/46/46512 Differential equations with fractional derivative and universal

More information

Fractional differential equations with integral boundary conditions

Fractional differential equations with integral boundary conditions Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 (215), 39 314 Research Article Fractional differential equations with integral boundary conditions Xuhuan Wang a,, Liping Wang a, Qinghong Zeng

More information

The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers

The 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 information

arxiv: v2 [math.ca] 13 Oct 2015

arxiv: v2 [math.ca] 13 Oct 2015 arxiv:1510.03285v2 [math.ca] 13 Oct 2015 MONOTONICITY OF FUNCTIONS AND SIGN CHANGES OF THEIR CAPUTO DERIVATIVES Kai Diethelm 1,2 Abstract It is well known that a continuously differentiable function is

More information

IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM

IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 8 (1), 17-30 (2002) IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM MAŁGORZATA JANKOWSKA 1, ANDRZEJ MARCINIAK 1,2 1 Poznań University

More information

Oscillatory Solutions of Nonlinear Fractional Difference Equations

Oscillatory Solutions of Nonlinear Fractional Difference Equations International Journal of Difference Equations ISSN 0973-6069, Volume 3, Number, pp. 9 3 208 http://campus.mst.edu/ijde Oscillaty Solutions of Nonlinear Fractional Difference Equations G. E. Chatzarakis

More information

A 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 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 information

The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation

The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM

More information

APPLICATION OF HAAR WAVELETS IN SOLVING NONLINEAR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS

APPLICATION 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 information

Finite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients

Finite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients International Journal of Contemporary Mathematical Sciences Vol. 9, 014, no. 16, 767-776 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.411118 Finite Difference Method of Fractional Parabolic

More information

EXISTENCE OF SOLUTIONS TO FRACTIONAL ORDER ORDINARY AND DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS

EXISTENCE 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 information

Second Order Step by Step Sliding mode Observer for Fault Estimation in a Class of Nonlinear Fractional Order Systems

Second Order Step by Step Sliding mode Observer for Fault Estimation in a Class of Nonlinear Fractional Order Systems Second Order Step by Step Sliding mode Observer for Fault Estimation in a Class of Nonlinear Fractional Order Systems Seyed Mohammad Moein Mousavi Student Electrical and Computer Engineering Department

More information

Fourth Order RK-Method

Fourth Order RK-Method Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)

More information

Chaos Control and Synchronization of a Fractional-order Autonomous System

Chaos Control and Synchronization of a Fractional-order Autonomous System Chaos Control and Snchronization of a Fractional-order Autonomous Sstem WANG HONGWU Tianjin Universit, School of Management Weijin Street 9, 37 Tianjin Tianjin Universit of Science and Technolog College

More information

The geometric and physical interpretation of fractional order derivatives of polynomial functions

The geometric and physical interpretation of fractional order derivatives of polynomial functions The geometric and physical interpretation of fractional order derivatives of polynomial functions M.H. Tavassoli, A. Tavassoli, M.R. Ostad Rahimi Abstract. In this paper, after a brief mention of the definitions

More information

Projective synchronization of a complex network with different fractional order chaos nodes

Projective synchronization of a complex network with different fractional order chaos nodes Projective synchronization of a complex network with different fractional order chaos nodes Wang Ming-Jun( ) a)b), Wang Xing-Yuan( ) a), and Niu Yu-Jun( ) a) a) School of Electronic and Information Engineering,

More information

This 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. 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 information

Boundary value problems for fractional differential equations with three-point fractional integral boundary conditions

Boundary 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 information

Solution of Conformable Fractional Ordinary Differential Equations via Differential Transform Method

Solution of Conformable Fractional Ordinary Differential Equations via Differential Transform Method Solution of Conformable Fractional Ordinary Differential Equations via Differential Transform Method Emrah Ünal a, Ahmet Gödoğan b a Department of Elementary Mathematics Education, Artvin Çoruh University,

More information

Positive solutions for a class of fractional boundary value problems

Positive 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 information

Fractional Order Heat Equation in Higher Space-Time Dimensions

Fractional Order Heat Equation in Higher Space-Time Dimensions Fractional Order Heat Equation in Higher Space-Time Dimensions Dimple Singh a,, Bhupendra Nath Tiwari b,, Nunu Yadav c, 3 a, b, c Amity School of Applied Sciences, Amity University Haryana Gurgaon, India

More information

Mathematics for chemical engineers. Numerical solution of ordinary differential equations

Mathematics for chemical engineers. Numerical solution of ordinary differential equations Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler

More information

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

More information

THE objective of this paper is to introduce a comparative

THE 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 information

THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands

THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay

More information

Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig

Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig A CAS Wavelet Method for Solving Nonlinear Fredholm Integro- Differential Equations of Fractional Order by Habib Allah Saeedi, Nasibeh

More information

Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials

Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department

More information

Nonlocal problems for the generalized Bagley-Torvik fractional differential equation

Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)

More information

SECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

SECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or

More information

c 2012 Society for Industrial and Applied Mathematics

c 2012 Society for Industrial and Applied Mathematics SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG

More information

Lecture 4: Numerical solution of ordinary differential equations

Lecture 4: Numerical solution of ordinary differential equations Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor

More information

Fractional Trigonometric Functions in Complexvalued Space: Applications of Complex Number to Local Fractional Calculus of Complex Function

Fractional Trigonometric Functions in Complexvalued Space: Applications of Complex Number to Local Fractional Calculus of Complex Function From the SelectedWorks of Xiao-Jun Yang June 4, 2 Fractional Trigonometric Functions in omplevalued Space: Applications of omple Number to Local Fractional alculus of omple Function Yang Xiao-Jun Available

More information

ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL

ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS

More information

FRACTIONAL BOUNDARY VALUE PROBLEMS ON THE HALF LINE. Assia Frioui, Assia Guezane-Lakoud, and Rabah Khaldi

FRACTIONAL BOUNDARY VALUE PROBLEMS ON THE HALF LINE. Assia Frioui, Assia Guezane-Lakoud, and Rabah Khaldi Opuscula Math. 37, no. 2 27), 265 28 http://dx.doi.org/.7494/opmath.27.37.2.265 Opuscula Mathematica FRACTIONAL BOUNDARY VALUE PROBLEMS ON THE HALF LINE Assia Frioui, Assia Guezane-Lakoud, and Rabah Khaldi

More information

Complete Identification of a Dynamic Fractional Order System Under Non-ideal Conditions Using Fractional Differintegral Definitions

Complete Identification of a Dynamic Fractional Order System Under Non-ideal Conditions Using Fractional Differintegral Definitions Complete Identification of a Dynamic Fractional Order System Under Non-ideal Conditions Using Fractional Differintegral Definitions Deepyaman Maiti # Ayan Acharya #2 R. Janarthanan *3 Amit Konar #4 # Department

More information

EXACT TRAVELING WAVE SOLUTIONS FOR NONLINEAR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS USING THE IMPROVED (G /G) EXPANSION METHOD

EXACT 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 information

Fractional-order linear time invariant swarm systems: asymptotic swarm stability and time response analysis

Fractional-order linear time invariant swarm systems: asymptotic swarm stability and time response analysis Cent Eur J Phys 11(6) 013 845-854 DOI: 10478/s11534-013-074-5 Central European Journal of Physics Fractional-order linear time invariant swarm systems: asymptotic swarm stability and time response analysis

More information

The Fractional-order SIR and SIRS Epidemic Models with Variable Population Size

The Fractional-order SIR and SIRS Epidemic Models with Variable Population Size Math. Sci. Lett. 2, No. 3, 195-200 (2013) 195 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.12785/msl/020308 The Fractional-order SIR and SIRS Epidemic Models with Variable

More information

1 Error Analysis for Solving IVP

1 Error Analysis for Solving IVP cs412: introduction to numerical analysis 12/9/10 Lecture 25: Numerical Solution of Differential Equations Error Analysis Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore

More information

YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL

YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions

More information

DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION

DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL

More information

Initial value problems for ordinary differential equations

Initial value problems for ordinary differential equations AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential

More information

A note on the uniform perturbation index 1

A note on the uniform perturbation index 1 Rostock. Math. Kolloq. 52, 33 46 (998) Subject Classification (AMS) 65L5 M. Arnold A note on the uniform perturbation index ABSTRACT. For a given differential-algebraic equation (DAE) the perturbation

More information

Jim Lambers MAT 772 Fall Semester Lecture 21 Notes

Jim 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 information

A Gauss Lobatto quadrature method for solving optimal control problems

A Gauss Lobatto quadrature method for solving optimal control problems ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes

More information

Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions

Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Mouffak Benchohra a,b 1 and Jamal E. Lazreg a, a Laboratory of Mathematics, University

More information

Mahmoud M. El-Borai a, Abou-Zaid H. El-Banna b, Walid H. Ahmed c a Department of Mathematics, faculty of science, Alexandria university, Alexandria.

Mahmoud 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 information

SOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER

SOME 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 information

Numerical Differential Equations: IVP

Numerical 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 information

On The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions

On The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions On The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions Xiong Wang Center of Chaos and Complex Network, Department of Electronic Engineering, City University of

More information

Solution of System of Linear Fractional Differential Equations with Modified Derivative of Jumarie Type

Solution of System of Linear Fractional Differential Equations with Modified Derivative of Jumarie Type American ournal of Mathematical Analysis, 015, Vol 3, No 3, 7-84 Available online at http://pubssciepubcom/ajma/3/3/3 Science and Education Publishing DOI:101691/ajma-3-3-3 Solution of System of Linear

More information

Numerical study of time-fractional hyperbolic partial differential equations

Numerical study of time-fractional hyperbolic partial differential equations Available online at wwwisr-publicationscom/jmcs J Math Computer Sci, 7 7, 53 65 Research Article Journal Homepage: wwwtjmcscom - wwwisr-publicationscom/jmcs Numerical study of time-fractional hyperbolic

More information

ERROR ANALYSIS, STABILITY, AND NUMERICAL SOLUTIONS OF FRACTIONAL-ORDER DIFFERENTIAL EQUATIONS

ERROR ANALYSIS, STABILITY, AND NUMERICAL SOLUTIONS OF FRACTIONAL-ORDER DIFFERENTIAL EQUATIONS International Journal of Pure and Applied Mathematics Volume 76 No. 5 2012, 647-659 ISSN: 1311-8080 (printed version) url: http://www.ijpam.eu PA ijpam.eu ERROR ANALYSIS, STABILITY, AN NUMERICAL SOLUTIONS

More information

The dynamical rigid body with memory

The dynamical rigid body with memory The dynamical rigid body with memory Ion Doru Albu, Mihaela Neamţu and Dumitru Opriş Abstract. In the present paper we describe the dynamics of the revised rigid body, the dynamics of the rigid body with

More information

Economic Interpretation of Fractional Derivatives

Economic Interpretation of Fractional Derivatives Progr. Fract. Differ. Appl. 3, No. 1, 1-6 (217) 1 Progress in Fractional Differentiation and Applications An International Journal http://dx.doi.org/1.18576/pfda/311 Economic Interpretation of Fractional

More information

EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS

EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL

More information