Almost Runge-Kutta Methods Of Orders Up To Five
|
|
- Emil Palmer
- 6 years ago
- Views:
Transcription
1 Almost Runge-Kutta Methods Of Orders Up To Five ABRAHAM OCHOCHE Department of Information Technology, Federal University of Technology, Minna, Niger State PMB 65, Nigeria PETER NDAJAH Graduate School of Science and Technology, Niigata University, 85 Ikarashi -no-cho, Nishi-ku Niigata, 95-8, Japan Abstract: In this paper, we have sought to investigate the viability of a type of general linear methods called Almost Runge-Kutta (ARK) methods, as a means of obtaining acceptable numerical approximations of the solution of problems in continuous mathematics We have outlined the derivation and implementation of this class of methods up to order five Extensive numerical experiments were carried out and the results clearly show that ARK methods are indeed a viable alternative to existing traditional methods Key Words: Almost, Order, Alternative, Euler, Runge Kutta, General Linear Methods Introduction Ordinary differential equations (ODEs) can be used to model many different types of physical behaviors, form chemical reactions to the notion of planets around the sun they are at the very heart of our understanding of the physical world today[see, 5] Unfortunately, as the systems being modeled get more complex, so do the solutions This can mean that we are left in the position of being unable to solve the problem analytically This is where the need for numerical methods become obvious Even though they do not give exactly solutions, they do enable us to find good approximations to the solution at any point where it is required General Linear Methods The name General linear methods applies to a large family of numerical methods for ordinary differential equations They were introduced in [9] to provide a unifying theory of the basic questions of consistency, stability and convergence Later they were used as a framework for studying accuracy questions and later the phenomena associated with nonlinear convergence They combine the multi-stage nature of Runge-Kutta methods with the notion of passing more than one piece of information between steps that is used in linear multi-step methods Their discovery opened the possibility of obtaining essentially new methods that were neither Runge-Kutta nor linear multi-step methods would exist which are practical and have advantages over the traditional methods These extremely broad class of methods, besides containing Runge-Kutta and linear multi-step methods as special cases also contain hybrid methods, cyclic composite linear multi-step methods and pseudo Runge-Kutta methods Examples of general linear methods As noted above, this class of methods is a large one It includes the traditional methods such as Runge Kutta methods and linear multi-step methods, along with methods that have been developed within the general linear methods framework, such as DIMSIMs, IRKS methods (with DESIRE and ESIRK as special cases) and Almost Runge-Kutta (ARK) methods 3 Almost Runge-Kutta Methods Almost Runge-Kutta (ARK) methods are a very special class of general linear methods introduced by Butcher [] The basic idea of these methods, is to retain the multi-stage nature of Runge-Kutta methods, while allowing more than one value to be passed from step to step Hence, they have a multi-value nature These methods have advantages over traditional methods, which are to be found in low-cost local estimation and dense output These latter features will be a consequence of the higher stage orders that are possible because of the multi-value nature of the new methods This multi-value nature brings its own difficulties The number of values passed between steps varies among general linear methods; the number of val- ISSN: Issue 5, Volume, May
2 ues is three for ARK methods Of the three input and output values in ARK methods, one approximates the solution value ie y(x n ) the second approximates the scaled first derivative (hy (x n )) while the third approximates second, derivatives respectively (h y (x n )) To simplify the starting procedure, the second derivative is required to be accurate only to within O(h 3 ), where h as usual is the In order to ensure that this low order does not adversely affect the solution value, the method has unbuilt annihilation conditions As a result of these extra input values, we are able to obtain stage order two, unlike explicit Runge-Kutta methods that are only able to obtain stage order one The advantage of this higher order is that we are able to interpolate and obtain an estimate at little extra cost [] 3 Motivation for the Methods When ARK methods were first proposed, it was hoped to develop a method that retained the spirit of R- K methods, but had advantages over the traditional methods The method was required to retain the simple stability functions of R-K methods as well as pass through very little information between steps The information that was passed from step to step was required to be easy to adjust for changes The main advantages that could be seen were a cheap estimator and dense output as well as the higher stage order that could be achieved [] 3 Structure of ARK Methods The general form of ARK methods is given by Y Y Y s 3 hf (Y ) hf (Y ) hf (Y s ) y [n ] y [n ] y [n ] r () where s is the number of internal stages A is an s s strictly lower triangular matrix, such that e T s A b T which implies that the last row of the A matrix is the same as the b vector in the first row of the B matrix s is the number of internal stages The vector e is a vector of ones, and is also of length s e s is a vector of length s consisting entirely of zeros except the sth component which is For an order p method the three output values are given by the equations below y(x n ) + O(h p+ ), hy (x n ) + O(h p+ ), () h y (x n ) + O(h 3 ) In choosing the coefficients of the method, we are careful to ensure that the simple stability properties of Runge-Kutta methods are retained In this introductory paper, we would focus on the case where the matrix A is strictly lower triangular by which we mean the methods are explicit Explicit ARK methods have the general form: Y Y Y s 3 a a 3 a 3 e c Ae c Ac a s, a s, a s,3 b b b 3 b s b b b 3 b s b β β β 3 β s β s β hf (Y ) hf (Y ) hf (Y s ) y [n ] y [n ] y [n ] r (3) As in traditional Runge-Kutta theory, b is a vector of length s representing the weights and c is a vector of length s representing the positions at which the function f is evaluated The vector e is a vector of ones, and is also of length s The form of the U matrix is given by a Taylor series expansion of the internal stages The internal stages of an ARK method is given by y(x + hc i ) u i y + u i hy + u i3 h y ISSN: Issue 5, Volume, May
3 a ij y (x + hc j ) (4) i +h j By carrying out the Taylor series expansion on both sides of equation (4) and equate the coefficients of y, y and y, we find u i, i,,, s u i c i j a ij u i3 c i j a ijc j (5) We must note some very important features of ARK methods: The first row of the B matrix is the same as the last row of the A matrix, we denote this as b T and the first row of the V matrix is the same as the last row of the U matrix, because we wish the final internal stage to give us the same quantity that is to be exported as the first outgoing approximation It also means that c s We also wish the second outgoing approximations to be h times the derivative of the final stage This means that the second row of the B and V matrices consists of zeros, with the exception of a in the (, s) position of B that is B is e T s The last row of B isβ T There are three types of conditions to consider when deriving Almost Runge-Kutta methods: Order conditions Annihilation conditions 3 Runge-Kutta stability conditions Detailed explanation of how these conditions are implemented, can be found in [], [] Theorem 3: An ARK method of order p with p stages has RK- stability if and only if β T (I + β s A) β s e T s (6) ( + β sc )b T A s c s! c exp s( β s ) β s exp s ( β s ) (7) (8) Proof: A proof of this theorem can be found in a number of literature including [], [3] 33 Methods with s p This section, considers the derivation procedure for methods which have the same number of stages as the order of the method Methods with this property are considered so as to minimize computational costs, and because, it is not possible to satisfy all the order conditions for s < p Also, the stability function is always a polynomial of degree s but the stability function is an approximation to exp(z) with an of (o{z p+ })as z This means that the stability function must have degree at least p and therefore s p 33 Derivation of Methods With s p 4 The general form of a fourth order, four stage ARK method is [ ] c c a c a c a c a 3 a 3 c 3 a 3 a 3 c 3 a 3c a 3 c b b b 3 b b b b 3 b β β β 3 β 4 β (9) The order conditions for a fourth order method are: b + b T e () b T c b T c 3 b T c 3 4 () () (3) b T Ac 6 b T Ac (4) (5) β T e + β (6) β T (I + β 4 A) β 4 e T 4 (7) c exp 4( β 4 ) β 4 exp 3 ( β 4 ) (8) ISSN: Issue 5, Volume, May
4 ( + β 4c )b T A c 4! (9) and (5) Provided c c, These methods have three free parameters, and we will take them to be c, c 3 and β 4 Once we assign values for these parameters, the method can then be uniquely defined using equations ()-(9) c can be calculated from equation (8), the b T vector can be found from the quadrature conditions (), () and (3), from where we get b 3 4c 4c 3 + 6c c 3 c (c c )(c c 3 ) b 3 4c 4c 3 + 6c c 3 c (c c )(c c 3 ) b 3 3 4c 4c + 6c c c 3 (c 3 c )(c 3 c ) () () () The next step is to obtain b from equation (), and a 3 from the linear combination of equations (4) a 3 c b 3 c (c c ) From conditions (4) and (9), a a 3 4b 3 a 3 c ( + β 4c ) (3) (4) 6 b 3a 3 c b a c b 3 c (5) Finally, β T can be found from equation (7) Based on the complete classification of fourth order methods with four stages given in [4], we present two methods based on Case 4 (ark4a) and Case 5 (ark4b) c T [,,, ], b ct [,,, ], β 4 [ ] , (6) 33 Derivation of Methods with s p 5 A fifth order, five stage ARK method takes the form c c a c a c a c a 3 a 3 c 3 a 3 a 3 c 3 a 3c a 3 c a 4 a 4 a 43 c 4 a 4 a 4 a 43 c 4 a 4c a 4 c a 43 c 3 b b b 3 b 4 b b b b 3 b 4 b β β β 3 β 4 β 5 β (7) with c [c, c, c 3 c 4, ] T The order conditions for a fifth order five stage method are: b + b T e (8) b T c (9) b T c 4 5 (3) c 4 (33) b T A b T (I C) (34) b T c 3 (3) b T (I C)Ac 4 (35) b T c 3 4 (3) b T (I C)Ac 6 (36) ISSN: Issue 5, Volume, May
5 and β T e + β (37) β T (I + β 5 A) β 5 e T 5 (38) c exp 5( β 5 ) β 5 exp 4 ( β 5 ) (39) ( + β 5c )b T A 3 c 5! (4) 333 Some Derived Methods Not much is known about the best choice of parameters for order five methods; however we will attempt to present some methods ARK5a and ARK5a respectively c T [ 5 65, 4,,, ], β c T [ 5 65,, 4 5,, ], β (4) (4) 34 Methods with s p+ From traditional Runge-Kutta theory, we know that we can achieve enhanced performance if there are more stages than are required for a particular order; the same applies to ARK methods The focus is on the case where there exists one more stage than is required to attain the required order 34 Methods with s 5, p 4 A fourth order method with five stages takes the form [ with stability function ] a a 3 a 3 e c Ae a 4 a 4 a 43 c Ae b b b 3 b 4 b b b 3 b 4 b β β β 3 β 4 β 5 β R(z) + z + z + z3 3! + z4 4! + Kz5 (43) ISSN: Issue 5, Volume, May
6 34 Order Conditions The following conditions must be satisfied by an ARK method with five stages to have order four b b T e (44) b T c b T c 3 (45) (46) b T c 3 4 (47) c 4 (48) b T Ac 6 b T Ac (49) (5) β β T e (5) β T c (5) β T Ac β 5(θ + ) φ θβ 5 (53) β 5 e T 5 (I + θa) β T (I + φa + β 5 θa ) (54) K( β 5c θα 4 α 5 ) ( + β 5c )(+α + α + α 3 3! + α 4 4! ) (55) where the values of α i are determined by expanding and + (φ β 5 )z + φz + β 5 θz α i z i (56) i b T A c 4 θ(bt A 3 c K) (57) β 5 (b T A c K) (β 5 φ)(b T A 3 c K) (58) 343 Some Examples Two methods are presented here; ark45a and ark45b respectively We have chosen L 5 and K for both methods, ensuring zero for both the bushy tree and the tall tree c T [ 4,, 3 4,, ], φ 4, µ 8, L 5, K c T [ 4, 3, 3 4,, ], φ 4, µ 8, L 5, K Implementation of ARK Methods (59) (6) The use of an ARK method is very similar to that of an RK method, with the main difference being that three pieces of information is now passed between steps The first two starting values are y(x )and hf(y(x ))respectively and the third starting value is obtained by taking a single Euler step forward and then taking the difference between the derivatives at these two points Therefore, the starting vector is ISSN: Issue 5, Volume, May
7 given by [y(x ), hf(y(x )), hf (y(x ) + hf(y(x ))) hf(y(x ))] (6) This choice of starting method was chosen for its simplicity, but it is adequate The method for computing the three starting approximations can be written in the form of the generalized Runge Kutta tableau where the zero in the first column of the last two rows indicates the fact that the term y n is absent from the output approximation This can be interpreted in the same way as a Runge Kutta method, but with three output approximations To change the we simply scale the vector in the same way we would scale a Nordsieck vector [5] Setting r h j /h j means the y vector needs to be scaled by [, r, r ] Once the starting vector is obtained from equation (6), the next thing is to calculate the internal stages; Y [n] i s j (hf (Y n i )a ij ) + Uy [n ] i (6) as well as the output approximations; i s j (hf (Y n i )B ij ) + V y [n ] i (63) 5 Numerical Experiments The problems that have been chosen are part of the DETest set of problems [6]; a set of standard test problems were suggested for testing ODE solvers Solutions have all been done with variable A: Logistics Curve: Single Equations y y 4 ( y ); y() Exact Solution: y E (t) e t 4 9+e t 4 B: The Radioactive Decay Chain: Moderate Systems y y y 9 y 9 9 with y() [,,,, ] T y y y 9 y C: A Non-linear Chemical Reaction Problem y y, y (), y y y, y (), y 3 y, y 3() ; (64) (65) D: A Problem Derived From Van der Pol s Equation y y, y (), y ( y )y y, y () 6 Discussions of Results (66) Figures to (see, next pages) show the results when the ARK methods we have derived in this paper, are used to solve problems A D From the s given by these methods, they may be considered as very good approximations of the exact results and we can say that they performed very well 7 Conclusion In this paper an effort has been made to present this unique class of methods and show that they are quite viable and reliable for solving not just single ODEs but large system as well For a more detailed analysis of Ark methods see [7] ISSN: Issue 5, Volume, May
8 Results for The Logistics Curve Problem ARK4a Results for The x Radioactive Decay Chain Problem ARK4a Figure : Solving Problem A using ark4a and ark4b Results for The Logistics Curve Problem ARK45a 6 3 Figure 4: Solving Problem B using ARK4a and 3 Results for The x Radioactive Decay Chain Problem ARK45a Figure : Solving Problem A using ark45a and ark45b 6 3 Results for The Logistics Curve Problem ARK5a Figure 5: Solving Problem B using ARK45a and 3 Results for The x Radioactive Decay Chain Problem ARK5a Figure 3: Solving Problem A using ARK5a and Figure 6: Solving Problem B using ARK5a and ISSN: Issue 5, Volume, May
9 6 8 Result for The Non linear Chemical Reaction Problem ARK4a Result for The Problem Derived From Van der Pol Equation ARK4a Figure 7: Solving Problem C using ARK4a and Figure : Solving Problem D using ARK4a and 7 8 Result for The Non linear Chemical Reaction Problem ARK45a 4 Result for The Problem Derived From Van der Pol Equation ARK45a Figure 8: Solving Problem C using ARK45a and Figure : Solving Problem D using ARK45a and 8 9 Result for The Non linear Chemical Reaction Problem ARK5a 4 Result for The Problem Derived From Van der Pol Equation ARK5a Figure 9: Solving Problem C using ARK5a and Figure : Solving Problem D using ARK5a and ISSN: Issue 5, Volume, May
10 Special Acknowledgement The authors wish to thank the ISAT fund of New Zealand for funding the recent visit of Abraham Ochoche, to the Department of Mathematics, The University of Auckland, Auckland, New Zealand where he had a one month stimulating exposure under Prof J C Butcher, who facilitated the visit Also, to acknowledge the immense contributions and cooperation received from Dr Nicolette Rattenbury, of the Manchester Metropolitan University, UK and Dr Helmut Podhaisky of the Institute of Mathematics, Martin-Luther-University, Halle- Wittenberg, Germany We also want to acknowledge Dr Allison Heard of the Department of Mathematics, University of Auckland, New Zealand for her immense contributions that led to the success of the Auckland visit References: [] Abraham, O, Improving the Improved Modified Euler Method for Better Performance on Autonomous Initial Value Problems, Leonardo J of Sciences, Issue, January June 8, [] Euler H, Institutiones Calculi Integralis Volumen Primum (768), Opera Omnia, Vol XI, B G Teubneri Lipsiae et Berolini MCMXIII [3] Butcher, J C, General Linear Methods: A Survey, Appl Numer Math (985) [4] Runge, C, Uber die Numerische Auflosung von Differentialgleichungen, Math Ann,46 (895), [5] Heun, K Neue Methoden zur Approximativen Integration Der Differentialgleichungen Einer UnabhAngigen VerAnderlichen, Z Math Phys, 45 (9), 3 38 [6] Kutta, W, Beitrag fur NAherungsweisen Integration Totaler Differentialgleichungen, Z Math Phys, 46 (9), [7] Burrage, K & Butcher, J C Non-Linear Stability for a General Class of Differential Equation Methods, BIT, (98), 85-3 [8] Bashforth, F and Adams, J C, An attempt to Test the Theories of Capillary Action by Comparing the Theoretical and Measured Forms of Drops of Fluid, with an Explanation of the Method of Integration Employed in Constructing the Tables which Give the Theoretical Forms of Such Drops, Cambridge University Press, Cambridge (883) [9] Butcher,J C, On the Convergence of Numerical Solutions of Ordinary Differential Equations, Math Comp, (966), [] Butcher, J C, An Introduction to Almost Runge Kutta Methods, Appl Numer Math, 4 (997),33 34 [] Rattenbury, N, Almost Runge-Kutta Methods for Stiff and Non-Stiff Problems, PhD Thesis, The University of Auckland, Department of Mathematics, 5 [] Rattenbury, N, Almost Runge-Kutta Methods: A Theoretical and Practical Investigation, MSc Thesis, The University of Auckland, department of Mathematics, (998) [3] Butcher, J C, Numerical Methods for Ordinary Differential Equations, John Wiley & Sons, Second Edition (8) [4] Butcher, J C, ARK Method Up to Order Five,Numer Algorithm, 7 (998), 93- [5] Nordseick, A, On Numerical Integration of Ordinary Differential Equations, Math Comp, 6 (96), 49 [6] Hull, T E, Enright, W H, Fellen, B M and Sedgwick, A E, Comparing Numerical Methods for Ordinary Differential Equations, SIAM J Numer Anal, 9 (97), [7] Abraham, O, The Dynamics of Almost Runge- Kutta Methods, Abacus: J of the Math Assoc of Nig, Vol 33, Num A (6) 6 7 ISSN: Issue 5, Volume, May
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3 J.H. Verner August 3, 2004 Abstract. In [5], Jackiewicz and Verner derived formulas for, and tested the implementation of two-step
More informationJournal of Quality Measurement and Analysis JQMA 4(1) 2008, 1-9 Jurnal Pengukuran Kualiti dan Analisis NUMERICAL ANALYSIS JOHN BUTCHER
Journal of Quality Measurement and Analysis JQMA 4(1) 2008, 1-9 Jurnal Pengukuran Kualiti dan Analisis NUMERICAL ANALYSIS JOHN BUTCHER ABSTRACT Mathematics has applications in virtually every scientific
More informationStarting Methods for Two-Step Runge Kutta Methods of Stage-Order 3 and Order 6
Cambridge International Science Publishing Cambridge CB1 6AZ Great Britain Journal of Computational Methods in Sciences and Engineering vol. 2, no. 3, 2, pp. 1 3 ISSN 1472 7978 Starting Methods for Two-Step
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 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 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 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 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 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 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 informationA NOTE ON EXPLICIT THREE-DERIVATIVE RUNGE-KUTTA METHODS (ThDRK)
BULLETIN OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN 303-4874 (p), ISSN (o) 303-4955 www.imvibl.org / JOURNALS / BULLETIN Vol. 5(015), 65-7 Former BULLETIN OF THE SOCIETY OF MATHEMATICIANS
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 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 informationA Class of an Implicit Stage-two Rational Runge-Kutta Method for Solution of Ordinary Differential Equations
Journal of Applied Mathematics & Bioinformatics, vol.2, no.3, 2012, 17-31 ISSN: 1792-6602 (print), 1792-6939 (online) Scienpress Ltd, 2012 A Class of an Implicit Stage-two Rational Runge-Kutta Method for
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 informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
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 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 informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
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 informationExplicit General Linear Methods with quadratic stability
Explicit with quadratic stability Angelamaria Cardone a, Zdzislaw Jackiewicz b ancardone@unisa.it, jackiewi@math.la.asu.edu a Dipartimento di Matematica e Informatica, Università degli Studi di Salerno,
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 - 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 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 informationGeneralized RK Integrators for Solving Ordinary Differential Equations: A Survey & Comparison Study
Global Journal of Pure and Applied Mathematics. ISSN 7-78 Volume, Number 7 (7), pp. Research India Publications http://www.ripublication.com/gjpam.htm Generalized RK Integrators for Solving Ordinary Differential
More informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
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 informationAug Vol. 5. No. 03 International Journal of Engineering and Applied Sciences EAAS & ARF. All rights reserved
FORMULATION OF PREDICTOR-CORRECTOR METHODS FROM 2-STEP HYBRID ADAMS METHODS FOR THE SOLUTION OF INITIAL VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS 1 ABUBAKAR M. BAKOJI, 2 ALI M. BUKAR, 3 MUKTAR
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 informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
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 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 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 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 informationA SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new
More informationNumerical solution of stiff ODEs using second derivative general linear methods
Numerical solution of stiff ODEs using second derivative general linear methods A. Abdi and G. Hojjati University of Tabriz, Tabriz, Iran SciCADE 2011, Toronto 0 5 10 15 20 25 30 35 40 1 Contents 1 Introduction
More informationEnergy-Preserving Runge-Kutta methods
Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,
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 informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
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 informationExplicit Runge Kutta Pairs with Lower Stage-Order *
myjournal manuscript No. (will be inserted by the editor) Explicit Runge Kutta Pairs with Lower Stage-Order * J.H. Verner Dedicated to John C. Butcher in celebration of his eightieth birthday October,
More informationNew Third Order Runge Kutta Based on Contraharmonic Mean for Stiff Problems
Applied Mathematical Sciences, Vol., 2009, no. 8, 65-76 New Third Order Runge Kutta Based on Contraharmonic Mean for Stiff Problems Osama Yusuf Ababneh 1 and Rokiah Rozita School of Mathematical Sciences
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 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 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 informationParallel general linear methods for stiff ordinary differential and differential algebraic equations
~-~H-~-~ ~ APPLIED ~ NUMERICAL MATHEMATICS ELSEVIER Applied Numerical Mathematics 17 (1995) 213-222 Parallel general linear methods for stiff ordinary differential and differential algebraic equations
More information(again assuming the integration goes left to right). Standard initial value techniques are not directly applicable to delay problems since evaluation
Stepsize Control for Delay Differential Equations Using Continuously Imbedded Runge-Kutta Methods of Sarafyan Skip Thompson Radford University Radford, Virginia Abstract. The use of continuously imbedded
More informationSeyyed Ahmad Hosseini
Seyyed Ahmad Hosseini Curriculum Vitae Assistant Professor of Applied Mathematics-Numerical Analysis at Golestan University Gender Personal Details Male Date of birth 21 st March, 1982 Place of birth Present
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
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 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 informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationSANGRADO PAGINA 17CMX24CM. PhD Thesis. Splitting methods for autonomous and non-autonomous perturbed equations LOMO A AJUSTAR (AHORA 4CM)
SANGRADO A3 PAGINA 17CMX24CM LOMO A AJUSTAR (AHORA 4CM) We have considered the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients.
More informationNumerical Solution of Linear Fredholm Integro-Differential Equations by Non-standard Finite Difference Method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISS: 1932-9466 Vol. 10, Issue 2 (December 2015), pp. 1019-1026 Applications and Applied Mathematics: An International Journal (AAM) umerical Solution
More information9.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 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 informationDownloaded 12/15/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 40, No. 4, pp. 56 537 c 2002 Society for Industrial and Applied Mathematics PREDICTOR-CORRECTOR METHODS OF RUNGE KUTTA TYPE FOR STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE
More informationAdebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction
Malaya Journal of Matematik Vol. No. 73-73 8 https://doi.org/.37/mjm/ An accurate five-step trigonometrically-fitted numerical scheme for approximating solutions of second order ordinary differential equations
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationOn the reliability and stability of direct explicit Runge-Kutta integrators
Global Journal of Pure and Applied Mathematics. ISSN 973-78 Volume, Number 4 (), pp. 3959-3975 Research India Publications http://www.ripublication.com/gjpam.htm On the reliability and stability of direct
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 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 informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität Halle Wittenberg Institut für Mathematik Implicit peer methods for large stiff ODE systems St. Beck, R. Weiner, H. Podhaisky and B. A. Schmitt Report No. 07 (2010) Editors: Professors
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 informationA New Numerical Approach for Solving Initial Value Problems of Ordinary Differential Equations
Annals of Pure and Applied Mathematics Vol. 17, No. 2, 2018, 157-162 ISSN: 2279-087X (P), 2279-0888(online) Published on 28 May 2018 www.researchmathsci.org DOI: http://dx.doi.org/10.22457/apam.v17n2a2
More informationExponentially Fitted Error Correction Methods for Solving Initial Value Problems
KYUNGPOOK Math. J. 52(2012), 167-177 http://dx.doi.org/10.5666/kmj.2012.52.2.167 Exponentially Fitted Error Correction Methods for Solving Initial Value Problems Sangdong Kim and Philsu Kim Department
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 informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
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 informationInitial-Value Problems for ODEs. Introduction to Linear Multistep Methods
Initial-Value Problems for ODEs Introduction to Linear Multistep Methods Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University
More informationYou may not use your books, notes; calculators are highly recommended.
Math 301 Winter 2013-14 Midterm 1 02/06/2014 Time Limit: 60 Minutes Name (Print): Instructor This exam contains 8 pages (including this cover page) and 6 problems. Check to see if any pages are missing.
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Numerical Methods Gaëtan Kerschen Space Structures & Systems Lab (S3L) Why Different Propagators? Analytic propagation: Better understanding of the perturbing forces. Useful
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationExplicit One-Step Methods
Chapter 1 Explicit One-Step Methods Remark 11 Contents This class presents methods for the numerical solution of explicit systems of initial value problems for ordinary differential equations of first
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
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 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 informationA Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital Problems
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal A Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
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 informationIntegration of Ordinary Differential Equations
Integration of Ordinary Differential Equations Com S 477/577 Nov 7, 00 1 Introduction The solution of differential equations is an important problem that arises in a host of areas. Many differential equations
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 informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
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 informationSolving Delay Differential Equations (DDEs) using Nakashima s 2 Stages 4 th Order Pseudo-Runge-Kutta Method
World Applied Sciences Journal (Special Issue of Applied Math): 8-86, 3 ISSN 88-495; IDOSI Publications, 3 DOI:.589/idosi.wasj.3..am.43 Solving Delay Differential Equations (DDEs) using Naashima s Stages
More information-Stable Second Derivative Block Multistep Formula for Stiff Initial Value Problems
IAENG International Journal of Applied Mathematics, :3, IJAM 3_7 -Stable Second Derivative Bloc Multistep Formula for Stiff Initial Value Problems (Advance online publication: 3 August ) IAENG International
More informationa 32 γ 0 1 a 31 a 42 a 43 γ 1 a 41 v n = u n + haf(v n ) (3) u n+1 = u n + hbf(v n ) (4)
A Generalized Runge Kutta Method of Order three Per G. Thomsen April 22, 2002 Abstract We present a numerical method for the solution of stiff systems of ODE's and index one DAE's. The type of method is
More informationStability and Convergence of Numerical Computations
Stability and Convergence of Numerical Computations Pavla Sehnalová Department of Intelligent Systems Faculty of Information Technology Brno University of Technology Božetěchova 2, 621 66 Brno, Czech Republic
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 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 informationExtended Runge Kutta-like formulae
Applied Numerical Mathematics 56 006 584 605 www.elsevier.com/locate/apnum Extended Runge Kutta-like formulae Xinyuan Wu a Jianlin Xia b a State Key Laboratory for Novel Software Technology Department
More informationThe variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
Cent. Eur. J. Eng. 4 24 64-7 DOI:.2478/s353-3-4-6 Central European Journal of Engineering The variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
More informationQuadratic SDIRK pair for treating chemical reaction problems.
Quadratic SDIRK pair for treating chemical reaction problems. Ch. Tsitouras TEI of Chalkis, Dept. of Applied Sciences, GR 34400 Psahna, Greece. I. Th. Famelis TEI of Athens, Dept. of Mathematics, GR 12210
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationThe Derivation of Interpolants for Nonlinear Two-Point Boundary Value Problems
European Society of Computational Methods in Sciences and Engineering ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics JNAIAM) vol. 1, no. 1, 2006, pp. 49-58 ISSN 1790 8140 The
More informationc 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