Solving Delay Differential Equations (DDEs) using Nakashima s 2 Stages 4 th Order Pseudo-Runge-Kutta Method
|
|
- Lizbeth Fitzgerald
- 5 years ago
- Views:
Transcription
1 World Applied Sciences Journal (Special Issue of Applied Math): 8-86, 3 ISSN ; IDOSI Publications, 3 DOI:.589/idosi.wasj.3..am.43 Solving Delay Differential Equations (DDEs) using Naashima s Stages 4 th Order Pseudo-Runge-Kutta Method Nur Ain Ayunni Sabri and Mustafa bin Mamat Department of Business and Finance, Faculty of Entrepreneurship and Business, Universiti Malaysia Kelantan (UMK), 6 Kota Bharu, Kelantan Department of Mathematics, Faculty of Science and Technology, Universiti Malaysia Terengganu, 3 K. Terengganu, Malaysia Abstract: This paper is investigate numerically the problem on solving Delay Differential Equations (DDEs) using Naashima s Stages 4th Order Pseudo-Runge-Kutta Method and also developed the algorithm of Naashima s Stages 4th Order Pseudo-Runge-Kutta Method incorporates with Hermite Interpolation. Numerical result is applied to a real life problem which is Food-Limited Model. Mathematica 7 software and Microsoft Excel are used to develop the calculation. Key words: Delay Differential Equations (DDEs) Pseudo-Runge Kutta method Hermite interpolation INTRODUCTION Delay differential equations (DDEs) are used to describe many phenomena in physics, biology and chemistry [, 7]. Pseudo-Runge-Kutta (PRK) Method is also one of the Runge-Kutta (RK) Methods family. In this article, Stages 4th Order Pseudo-Runge-Kutta Method, which is a type of Pseudo Runge-Kutta (PRK) Methods will be discussed. According to F.Constabile (99), Pseudo-Runge Kutta (PRK) Methods require fewer evaluations than Runge-Kutta (RK) Method of the same order. Fourth Order Runge-Kutta (RK4) Method will be used to solve Initial Value Problems (IVPs) in order to get the first term. Many real life problem phenomena in physics, engineering, biology, medicine, economics etc can be modelled by an Initial Value Problem (IVP) or Cauchy problem for Ordinary Differential Equations (ODEs ) of the type ( ) y y' t gt,y t,t t yt (.) where the function y(t) called the state variable represents some physical quantity that evolves over the time []. ODEs are usually solved using the step-sie methods that. The method that commonly been used to solve ODEs are RK Method, Euler Method and any other step-sie methods. RK Method is the most In mathematics, DDEs are differential equations in which the derivative of the unnown function at certain time is given in terms of the values of the function at previous times. DDEs differ from ODEs in that the derivatives at any time depend on the solution (and in the case of neutral equations on the derivatives) at prior times. The simplest constant delay equation has the form ( ( τ) ( τ) ( τ) ) ' y t f t,y t,y t,yt,,yt (.) where the time delay (lags) τ j are positive constant. More generally, state dependent delays may depends on the solution, that is ττ i i( t,yt ). In more general models, the derivative y (t) may depend on y and y itself at some past value-τ. In this case, Equation (.) changes into the form ' ' φ y t f t,y t,y t τ,y t τ, t t yt t, t t popular method that been used to solve ODEs. Corresponding Author: Nur Ain Ayunni Sabri, Department of Business and Finance, Faculty of Entrepreneurship and Business, Universiti Malaysia Kelantan (UMK), 6 Kota Bharu, Kelantan 8 (.3) where the function φ(t) is assumed to be at least continuous. Equation (.3) is called a Delay Differential Equation of neutral type (NDDE). The methods that used to solve DDEs in this paper are RK Method and one of the step-sie methods, which is Hermite Interpolation. For the recent years, there are many papers have been written regarding on solving DDEs [5-5]. In this
2 World Appl. Sci. J., (Special Issue of Applied Math): 8-86, 3 research, we solve DDEs using Naashima.s Stages 4th Order Pseudo-Runge-Kutta Method and we used this method incorporates with Hermite Interpolation. PSEUDO-RUNGE-KUTTA AND HERMITE INTERPOLATION Pseudo-Runge-Kutta formulae that we are considering is Naashima s Stages 4 th Order Pseudo- Runge-Kutta [3]: 7 5 yn+ yn + h ( ) f x,y n n f x,y n n x +.7h,.56y +.56y n n n f +.833h +.3h (.) Interpolation is to obtain an approximation solution for a specific value of x, it is possible to shorten the final step, if necessary, to complete the step exactly at the right place. It is usually more convenient to rely on step sie control mechanism that is independent of output requirements and to produce require output results by interpolation as the opportunity arises. The use of interpolation maes it also possible to produce output at multiple and arbitrary points [4]. Hermite Interpolation is a method of interpolating data points as a polynomial function, in the field of numerical analysis. The generated Hermite Polynomial is closely related to the Newton Polynomial in that both are derived from the calculation of Divided Differences. Hermite Interpolation matches data points in both value and first derivative. This means that n derivatives values ( x,y ),,( x,y ) must be given in addition n n to n data points ( x,y ),,( x,y ). n n NAKASHIMA S STAGES 4TH ORDER PSEUDO-RUNGE-KUTTA METHOD INCORPORATE WITH HERMITE INTERPOLATION AND THE ALGORITHM y n y + + n ( 3 4) forj,,,,m (3.) where hf t,y n n h hf t n+,yn + h 3 hf t n +,yn + hf t + h,y + 4 n n 3. From the Classical RK4 Method which is Equation (3.), the equations becomes y y n n 3 4 ( ) f t,y,y t τ n n n h h h f t n +,yn +,y tn + τ h h h f t,y,y t τ n n n ( ) f t + h,y + h,yt + h τ 4 n n 3 n (3.) where y is the delay argument which depends on the question and τ given in the question. 3. From the first step, we will get the delay argument that is not fix point. Therefore, here we need do some approximation using Hermite Interpolation. The mesh point,,, n that we going to choose must: a) Depends on the delay argument that we want to approximate and the value of h. For each point that we choose, that point is multiple becomes two points. b) Contain all the argument in the interval of mesh points. 4. Next, substitute the value of mesh point that already choose to find new argument of From the meaning of pseudo and RK Method, PRK Method is the method that is not actually RK Method but having the appearance of RK Method. It has been proof that PRK Method is one of the family members of RK Method. Algorithm. Runge-Kutta Method use the formula t + t +h and 8 [ ] y ( ) [ ] y ( ) y f f f y f n 5. Then, using divided difference to find the approximation. The finite difference formula that need to find are: n
3 World Appl. Sci. J., (Special Issue of Applied Math): 8-86, 3 [ ] y ( ) f[ ] f [ ] f, f[ ] f,,f, y 3 f, n n f f, y n n n n n n n 6. Continue the calculation for, f,, f,, 3 f,, n n n f, f, f, f, 3 3 f, f, n n n n n n Next, Until f,,, 3 f,,, 3 4 f,,, n 3 n n n [ ] f,,,, n f,, f,, 3 3 f,, f,, f,, f,, n n n n 3 n n n n 3 [ ] [ ] f,,, f,,,, n n n 7. Finally, substitute into Hermite Interpolation formulae. ( ) [ ] [ ]( n ) [ ]( n )( n ) [ 3]( n )( n )( n ) [ ] f t,y,y t τ n n n f + f, t τ + f,, t τ t τ + f,,, t τ t τ t τ + + f,,,, t τ t τ t τ t τ n n n n n n (3.3) h h h f t n +,yn +,y tn + τ h h h f[ ] + f [,] tn + τ + f,, [ ] tn + τ tn + τ h h h + f,,, [ 3] tn + τ tn + τ tn + τ + + h h h h + f,,, [,n] tn + τ tn + τ tn + τ tn + τ n h h h 3 f t n +,yn +,y tn + τ h h h f[ ] + f [,] tn + τ + f,, [ ] tn + τ tn + τ h h h + f,,, [ 3] tn + τ tn + τ tn + τ + + h h h h + f,,, [,n] tn + τ tn + τ tn + τ tn + τ n ( ) [ ] [ ]( n ) [ ]( n )( n ) f,,, [ 3]( tn h )( tn h )( tn h ) [ ] f t + h,y + h,yt + h τ 4 n n 3 n f + f, t + h τ + f,, t + h τ t + h τ + + τ + τ + τ f,,,, t + h τ t + h τ t + h τ t + h τ n n n n n n (3.4) (3.5) (3.6) 8. Substitute all the value of,, 3 and 4 into equation 83
4 World Appl. Sci. J., (Special Issue of Applied Math): 8-86, 3 h y n y + + n ( 3 4) to get y. 9. To calculate y, y 3,, we will use Naashima s Stages 4 th Order Pseudo-Runge-Kutta Method incorporate with Hermite Interpolation. From the formula of Naashima s Stages 4 th Order Pseudo-Runge-Kutta Method and Equation (.) becomes, 7 5 yn+ yn + h (3.7) ( ( )) ( ) ( ) f t,y,yt τ n n n f t,y,y t τ n n n f t +.7h,.56y +.56y +.833h +.3h,y t +.7h τ n n n n. From this step, we will get the delay argument that is not fix point. Therefore, here we need do some approximation using Hermite Interpolation. The mesh point,,, n that we going to choose must : a) pends on the delay argument that we want to approximate and the value of h. For each point that we choose, that point is multiple becomes two points. b) Contain all the argument in the interval of mesh point. Repeat steps 4, 5 and 6 for calculations of finite difference.. Next, substitute the values that we get previously in Hermite Interpolation formula. ( ( )) [ ] [ ]( n ) [ ]( n )( n ) f,,, [ 3]( tn )( tn )( tn ) f,,,, [, ]( t )( t )( t ) ( t ) f t,y,yt τ n n n f + f, t τ + f,, t τ t τ + τ τ τ τ τ τ τ 3 n n n n n n ( ) [ ] [ ]( n ) [ ]( n )( n ) f,,, [ 3]( tn )( tn )( tn ) [ ] f t,y,y t τ n n n f + f, t τ + f,, t τ t τ + τ τ τ f,,,, 3,n tn τ tn τ tn τ tn τ n ( ) [ ] [ ]( n ) [ ]( n )( n ) f,,, [ 3]( tn.7h )( tn.7h )( tn.7h ) f,,,, [, ]( t.7h )( t + τ )( + τ ) ( + τ ) f t +.7h,.56y +.56y +.833h +.3h,y t +.7h τ n n n n f + f, t +.7h τ + f,, t +.7h τ t +.7h τ + + τ + τ + τ τ 3 n n n.7h t.7h t.7h n n n (3.8) (3.9) (3.). Finally, substitute all the values of, and into to find the next y. 7 5 yn+ yn + h SOLVING DDE S USING NAKASHIMA S STAGES 4TH ORDER PSEUDO-RUNGE-KUTTA METHOD We have chosen a model in order to apply Naashima s Stages 4 th Order Pseudo-Runge-Kutta Methods. The model that we choose was Food Limited Model. We also mae a comparison with 4 th Order Runge-Kutta Method or also nown as Classical Runge-Kutta Method in order to see the accuracy between these two methods. 84
5 World Appl. Sci. J., (Special Issue of Applied Math): 8-86, 3 Table 4.: Results for application food limited model using RK4 method and Naashima s stages 4 th order Pseudo-Runge- Kutta method incorporate with Hermite interpolation Naashima s stages 4 th order t RK4 method Pseudo-Runge-Kutta method Food limited model: Food-limited model is a model for population growth which is proposed by Smith in 963. This model is an alternative to the logistic equation for food-limited population dx K xt r( x ),t dt K + rcn t (4.) Here N, r and K are the mass of the population, the rate of increase with unlimited food and value of N at saturation respectively. The constant /c is the rate replacement of mass in the population at saturation (this includes both the replacement of metabolic loss and of dead organism) [6]. The delay Food-Limited is where ( τ) dx K x t rx( t) dt K+ rcx t τ r.5,k.,c and τ 8 for t [,] and x t e for t<. 85 h.5. The initial function.5t We see that the solution from both methods is almost same although have small error between both method, this shows that Naashima.s Stages 4th Order Pseudo-Runnge-Kutta Mehod also can be applied to real life problem of DDEs. By comparing both methods, we can see the validity of the new method in solving the real life problem of DDEs. We conclude that Naashima.s Stages 4th Order Pseudo-Runge- Kutta Method can be used to solve DDEs. We used to compare both methods because we want to find other method to solve real life problems of DDEs and we found that Naashima s Stages 4 th Order Pseudo-Runge-Kutta Method is one of the method that we can use to solve DDEs. CONCLUSIONS In this study, the method that been used to solve Delay Differential Equations (DDEs) is Naashima s Stages 4 th Order Pseudo-Runge-Kutta method incorporate with Hermite Interpolation. From numerical results, we have observed that PRK method perform a good result in solving DDEs of single delay because the answer is accurate as Runge- Kutta method. PRK method is accurate as RK method and it is also one step method. From here we can mae a conclusion that Naashima s Stages 4 th Order Pseudo-Runge-Kutta method is another method to solve DDEs problem. We also conclude that, the solution using Naashima s Stages 4 th Order Pseudo-Runge- Kutta method need less calculation than RK4 method. This is because, using Naashima s Stages 4 th Order Pseudo-Runge-Kutta method, we only need to calculate until, but using RK4 method we need to calculate until 4. Therefore, we note that Naashima s Stages 4 th Order Pseudo-Runge-Kutta method is also can be used to solve DDEs and this is another method to solve DDEs. REFERENCES. Sharma, J.N., 4. Numerical Methods for Engineers and Scientist. Deemed University: Alpha Science International Ltd. Pangbourne England.. Bellen, A. and M. Zennaro, 3. Numerical Methods for Delay Differential Equations. United States: Oxford Science Publications. 3. Masaharu Naashima, 98. On Pseudo-Runge- Kutta Methods with and 3 Stages. RIMS Publications, 8 (3): Butcher, J.C., 8. Numerical Methods for Ordinary Differential Equation. nd Edn. John Willey and Sons Ltd.
6 World Appl. Sci. J., (Special Issue of Applied Math): 8-86, 3 5. Al-Mutib, A.N., 984. One-Step Implicit Methods for Solving Delay Differential Equations. International Journal of Computer Mathematics, 6: Isha, F., M.B. Suleiman and Z. Omar, 8. Two-Point Predictor-Corrector Bloc Methods for Solving Delay Differential Equations. Matematia, 4 (): Lwin, A.S., 4. Solving Delay Differential Equations Using Explicit Runge-Kutta Methods. Master Thesis. Universiti Putra Malaysia. 8. Wong, W.S., Y. Zhang and S. Li, 9. Stability of Continuous Runge-Kutta Methods for Non Linear Neutral Delay Differential Equations. Applied Mathematical Modeling, 33: Ismail, F. and M. Suleiman,. Solving Delay Differential Equations Using Intervalwise Partitioning by Runge-Kutta Method. Applied Mathematics and Computation, : Ismail, F., R.A. Al-Khasawneh, S.L. Aung and M. Suleiman,. Numerical Treatment of Delay Differential Equations by Runge-Kutta Method using Hermite Interpolation. Matematia, 8 (): Huang, C., S. Li, H. Fu and G. Chen,. D-Convergence of General Linear Methods for Stiff Delay Differential Equations. Computers and Mathematics with Applications, 4: Verheyden, K., T. Luyanina and D. Roose, 8. Efficient Computation of Characteristic Roots of Delay Differential Equations Using Linear Multistep Method. Journal of Computational and Applied Mathematics, 4: Zhang, C. and G. Sun, 4. Nonlinear Stability of Runge-Kutta Methods Applied to Infinite Delay Differential Equations. Mathematical and Computer Modeling, 39: Hout, K.J., 996. On The Stability of Adaptations of Runge-Kutta Methods to System of Delay Differential Equations. Applied Numerical Mathematics, : Bartosewsi, Z. and Z. Jaciewic,. Stability Analysis of Two-Step Runge-Kutta Method for Delay Differential Equations. Computers and Mathematics with Application, 44: Arino, O., M.L. Hbid and E.A. Dads, 6. Delay Differential Equations and Applications. Berlin:Springer. pp: Hafshejani, M.S., S.K. Vanani and J.S. Hafshejani,. Numerical Solution of Delay Differential equations Using Legendre Wavelet Method. World Applied Sciences Journal, 3 (Special Issue of Applied Maths):
Modified Milne Simpson Method for Solving Differential Equations
Modified Milne Simpson Method for Solving Differential Equations R. K. Devkate 1, R. M. Dhaigude 2 Research Student, Department of Mathematics, Dr. Babasaheb Ambedkar Marathwada University, Aurangabad,
More informationFourth 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 informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationmultistep 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 informationResearch Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations
International Mathematics and Mathematical Sciences Volume 212, Article ID 767328, 8 pages doi:1.1155/212/767328 Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
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 informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationORBIT 14 The propagator. A. Milani et al. Dipmat, Pisa, Italy
ORBIT 14 The propagator A. Milani et al. Dipmat, Pisa, Italy Orbit 14 Propagator manual 1 Contents 1 Propagator - Interpolator 2 1.1 Propagation and interpolation............................ 2 1.1.1 Introduction..................................
More informationApplied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations
Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin
More informationFifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations
Australian Journal of Basic and Applied Sciences, 6(3): 9-5, 22 ISSN 99-88 Fifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations Faranak Rabiei, Fudziah Ismail Department
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More information5.6 Multistep Methods
5.6 Multistep Methods 1 Motivation: Consider IVP: yy = ff(tt, yy), aa tt bb, yy(aa) = αα. To compute solution at tt ii+1, approximate solutions at mesh points tt 0, tt 1, tt 2, tt ii are already obtained.
More informationSolving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method
World Applied Sciences Journal 21 (6): 939-944, 2013 ISSN 11-4952 IDOSI Publications, 2013 DOI: 10.529/idosi.wasj.2013.21.6.2915 Solving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method 1 1 3
More informationSolution of First Order Initial Value Problem by Sixth Order Predictor Corrector Method
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 2277 2290 Research India Publications http://www.ripublication.com/gjpam.htm Solution of First Order Initial
More informationCOSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods
COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:
More informationOne-Step Hybrid Block Method with One Generalized Off-Step Points for Direct Solution of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 10, 2016, no. 29, 142-142 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.6127 One-Step Hybrid Block Method with One Generalized Off-Step Points for
More informationMath 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 informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationLecture 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 informationOn interval predictor-corrector methods
DOI 10.1007/s11075-016-0220-x ORIGINAL PAPER On interval predictor-corrector methods Andrzej Marcinia 1,2 Malgorzata A. Janowsa 3 Tomasz Hoffmann 4 Received: 26 March 2016 / Accepted: 3 October 2016 The
More informationA CLASS OF CONTINUOUS HYBRID LINEAR MULTISTEP METHODS FOR STIFF IVPs IN ODEs
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LVIII, 0, f. A CLASS OF CONTINUOUS HYBRID LINEAR MULTISTEP METHODS FOR STIFF IVPs IN ODEs BY R.I. OKUONGHAE Abstract.
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationInitial 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 informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationDelay Differential Equations with Constant Lags
Delay Differential Equations with Constant Lags L.F. Shampine Mathematics Department Southern Methodist University Dallas, TX 75275 shampine@smu.edu S. Thompson Department of Mathematics & Statistics Radford
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationMath 225 Differential Equations Notes Chapter 1
Math 225 Differential Equations Notes Chapter 1 Michael Muscedere September 9, 2004 1 Introduction 1.1 Background In science and engineering models are used to describe physical phenomena. Often these
More informationarxiv: v1 [math.na] 31 Oct 2016
RKFD Methods - a short review Maciej Jaromin November, 206 arxiv:60.09739v [math.na] 3 Oct 206 Abstract In this paper, a recently published method [Hussain, Ismail, Senua, Solving directly special fourthorder
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 information1 Ordinary Differential Equations
Ordinary Differential Equations.0 Mathematical Background.0. Smoothness Definition. A function f defined on [a, b] is continuous at ξ [a, b] if lim x ξ f(x) = f(ξ). Remark Note that this implies existence
More informationIntroduction to standard and non-standard Numerical Methods
Introduction to standard and non-standard Numerical Methods Dr. Mountaga LAM AMS : African Mathematic School 2018 May 23, 2018 One-step methods Runge-Kutta Methods Nonstandard Finite Difference Scheme
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
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 informationNumerical solution of delay differential equations via operational matrices of hybrid of block-pulse functions and Bernstein polynomials
Computational Methods for Differential Equations http://cmdetabrizuacir Vol, No, 3, pp 78-95 Numerical solution of delay differential equations via operational matrices of hybrid of bloc-pulse functions
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationSolving Delay Differential Equations By Elzaki Transform Method
SCITECH Volume 3, Issue 1 RESEARCH ORGANISATION February 22, 2017 Boson Journal of Modern Physics www.scitecresearch.com Solving Delay Differential Equations By Elzaki Transform Method Ebimene, J. Mamadu
More informationResearch Article Solving Delay Differential Equations of Small and Vanishing Lag Using Multistep Block Method
Applied Mathematics, Article ID 348912, 10 pages http://dx.doi.org/10.1155/2014/348912 Research Article Solving Delay Differential Equations of Small and Vanishing Lag Using Multistep Bloc Method Nurul
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationModule 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008
Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.
Logistics Week 12: Monday, Apr 18 HW 6 is due at 11:59 tonight. HW 7 is posted, and will be due in class on 4/25. The prelim is graded. An analysis and rubric are on CMS. Problem du jour For implicit methods
More informationSolving scalar IVP s : Runge-Kutta Methods
Solving scalar IVP s : Runge-Kutta Methods Josh Engwer Texas Tech University March 7, NOTATION: h step size x n xt) t n+ t + h x n+ xt n+ ) xt + h) dx = ft, x) SCALAR IVP ASSUMED THROUGHOUT: dt xt ) =
More informationA Zero-Stable Block Method for the Solution of Third Order Ordinary Differential Equations.
A Zero-Stable Block Method for the Solution of Third Order Ordinary Differential Equations. K. Rauf, Ph.D. 1 ; S.A. Aniki (Ph.D. in view) 2* ; S. Ibrahim (Ph.D. in view) 2 ; and J.O. Omolehin, Ph.D. 3
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 informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationNotes for Numerical Analysis Math 5466 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) Contents Numerical Methods for ODEs 5. Introduction............................
More information8.1 Introduction. Consider the initial value problem (IVP):
8.1 Introduction Consider the initial value problem (IVP): y dy dt = f(t, y), y(t 0)=y 0, t 0 t T. Geometrically: solutions are a one parameter family of curves y = y(t) in(t, y)-plane. Assume solution
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationOn linear and non-linear equations.(sect. 2.4).
On linear and non-linear equations.sect. 2.4). Review: Linear differential equations. Non-linear differential equations. Properties of solutions to non-linear ODE. The Bernoulli equation. Review: Linear
More informationNUMERICAL ANALYSIS 2 - FINAL EXAM Summer Term 2006 Matrikelnummer:
Prof. Dr. O. Junge, T. März Scientific Computing - M3 Center for Mathematical Sciences Munich University of Technology Name: NUMERICAL ANALYSIS 2 - FINAL EXAM Summer Term 2006 Matrikelnummer: I agree to
More informationThe 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 informationProject on Runge Kutta method
Project on Runge Kutta method Nguyen Quan Ba Hong Doan Tran Nguyen Tung Nguyen An Thinh 3 Students at Faculty of Math and Computer Science Ho Chi Minh University of Science, Vietnam email. email. nguyenquanbahong@gmail.com
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationMathematics 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 informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationZ. Omar. Department of Mathematics School of Quantitative Sciences College of Art and Sciences Univeristi Utara Malaysia, Malaysia. Ra ft.
International Journal of Mathematical Analysis Vol. 9, 015, no. 46, 57-7 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.015.57181 Developing a Single Step Hybrid Block Method with Generalized
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationOrdinary Differential Equations
Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationResearch Article An Accurate Block Hybrid Collocation Method for Third Order Ordinary Differential Equations
Applied Mathematics, Article ID 549597, 9 pages http://dx.doi.org/1.1155/14/549597 Research Article An Accurate Bloc Hybrid Collocation Method for Third Order Ordinary Differential Equations Lee Ken Yap,
More informationMath 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.
Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.8: 1b, 3b, 4 Burden & Faires 5.6. Multistep Methods 1. Use all the Adams-Bashforth
More informationUniform Order Legendre Approach for Continuous Hybrid Block Methods for the Solution of First Order Ordinary Differential Equations
IOSR Journal of Mathematics (IOSR-JM) e-issn: 8-58, p-issn: 319-65X. Volume 11, Issue 1 Ver. IV (Jan - Feb. 015), PP 09-14 www.iosrjournals.org Uniform Order Legendre Approach for Continuous Hybrid Block
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More information2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l.
Numerical Methods II UNIT.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS.1.1 Runge-Kutta Method of Fourth Order 1. Let = f x,y,z, = gx,y,z be the simultaneous first order
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 informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationOn Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
Nonlinear Analysis and Differential Equations, Vol. 6, 08, no., 5-64 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.08.8 On Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
More informationA Brief Introduction to Numerical Methods for Differential Equations
A Brief Introduction to Numerical Methods for Differential Equations January 10, 2011 This tutorial introduces some basic numerical computation techniques that are useful for the simulation and analysis
More informationDelay Differential Equations Part II: Time and State dependent Lags
Delay Differential Equations Part II: Time and State dependent Lags L.F. Shampine Department of Mathematics Southern Methodist University Dallas, Texas 75275 shampine@smu.edu www.faculty.smu.edu/shampine
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
More informationImproved Extended Runge-Kutta-like Method for Solving First Order IVPs
Indian Journal of Science and Technology, Vol 9(48), DOI: 0.7485/ijst/06/v9i48/0938, December 06 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 Improved Extended Runge-Kutta-like Method for Solving
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationMTH 452/552 Homework 3
MTH 452/552 Homework 3 Do either 1 or 2. 1. (40 points) [Use of ode113 and ode45] This problem can be solved by a modifying the m-files odesample.m and odesampletest.m available from the author s webpage.
More informationCONVERGENCE PROPERTIES OF THE MODIFIED NUMEROV RELATED BLOCK METHODS
Vol 4 No 1 January 213 CONVERGENCE PROPERTIES OF THE MODIFIED NUMEROV RELATED BLOCK METHODS YS, Awari 1 and AA, Abada 2 1 Department of Mathematics & 2 Department of Statistics, BinghamUniversity, Karu,
More informationContinuous Block Hybrid-Predictor-Corrector method for the solution of y = f (x, y, y )
International Journal of Mathematics and Soft Computing Vol., No. 0), 5-4. ISSN 49-8 Continuous Block Hybrid-Predictor-Corrector method for the solution of y = f x, y, y ) A.O. Adesanya, M.R. Odekunle
More informationMethod of Successive Approximations for Solving the Multi-Pantograph Delay Equations
Gen. Math. Notes, Vol. 8, No., January, pp.3-8 ISSN 9-784; Copyright c ICSRS Publication, www.i-csrs.org Available free online at http://www.geman.in Method of Successive Approximations for Solving the
More informationChapter 10 Skill biased technological change
Chapter 10 Skill biased technological change O. Afonso, P. B. Vasconcelos Computational Economics: a concise introduction O. Afonso, P. B. Vasconcelos Computational Economics 1 / 29 Overview 1 Introduction
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationNumerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
Available online at http://ijim.srbiau.ac.ir Int. J. Industrial Mathematics Vol. 1, No. 2 (2009)147-161 Numerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
More informationA modification of Kaldor-Kalecki model and its analysis
A modification of Kaldor-Kalecki model and its analysis 1 Introduction Jan Kodera 1, Jarmila Radová 2, Tran Van Quang 1 Abstract. This paper studies the investment cycle as an endogenous phenomenon using
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2
More information