Lecture 6. Numerical Solution of Differential Equations B21/B1. Lecture 6. Lecture 6. Initial value problem. Initial value problem
|
|
- Aubrey Cummings
- 5 years ago
- Views:
Transcription
1 Lecture 6 Numerical Solution of Differential Equations B2/B Initial value problem y = f(x,y), x [x,x M ] y(x ) = y Professor Endre Süli Numerical Solution of Differential Equations p./45 Numerical Solution of Differential Equations p.2/45 Lecture 6 Lecture 6 Initial value problem Initial value problem y = f(x,y), x [x,x M ] y(x ) = y y = f(x,y), x [x,x M ] y(x ) = y General form of a linear multistep method: General form of a linear multistep method: α j y n+j = h β j f n+j f n+j = f(x n+j,y n+j ), α k, α + β Numerical Solution of Differential Equations p.2/45 Numerical Solution of Differential Equations p.2/45
2 Stiff problems: Van der Pol s oscillator Stiff problems: Van der Pol s oscillator y + µ(y 2 )y + y =, y() = 2, y () = with µ =,, 2, 5,, 2, 5. Solve on the interval [, 6]. y + µ(y 2 )y + y =, y() = 2, y () = with µ =,, 2, 5,, 2, 5. Solve on the interval [, 6]. y = y 2, y () = 2 y 2 = µ( y 2 )y 2 y, y 2 () = Numerical Solution of Differential Equations p.3/45 Numerical Solution of Differential Equations p.3/45 Result for µ = : y = y Result for µ = : y = y Numerical Solution of Differential Equations p.4/ Numerical Solution of Differential Equations p.5/45
3 Result for µ = 2: y = y Result for µ = 5: y = y Result for µ = : y = y Numerical Solution of Differential Equations p.6/ Result for µ = 2: y = y Numerical Solution of Differential Equations p.7/ Numerical Solution of Differential Equations p.8/ Numerical Solution of Differential Equations p.9/45
4 Result for µ = 5: y = y A stiff linear system 2.5 Let λ < and consider 2.5 y (λ )y λy =, y() =, y () = λ 2.5 Write y = y, y 2 = y, y = (y,y 2 ) T, y = Ay, y() = y..5 A = ( λ λ ) and y = ( λ 2 ) A stiff linear system Numerical Solution of Differential Equations p./45 A stiff linear system: λ = 45 Numerical Solution of Differential Equations p./45 det(a zi) = z 2 (λ )z λ 45 z =, y (x) = 2e x e λx, z 2 = λ y 2 (x) = 2e x λe λx y y Numerical Solution of Differential Equations p.2/45 5 Numerical Solution of Differential Equations p.3/45
5 Matlab code (for λ = 45) The simplest example linear.m function dy = lineq(x,y) lambda=-45; dy = zeros(2,); dy() = y(2); dy(2) = lambda * y()+(lambda-)*y(2); launchlin.m lambda=-45; [x,y] = ode3( linear,[ 2],[ -lambda-2]); plot(x,y(:,), kd-, x, y(:,2), kp- ) where Exact solution y = λy, y() = Re(λ) <. y(x) = e λx = e xre(λ) e ıxim(λ) Numerical Solution of Differential Equations p.4/45 Numerical Solution of Differential Equations p.5/45 The simplest example: y = λy, y() =, for λ = 45 Van der Pol s equation by Euler s method: h = u Numerical Solution of Differential Equations p.6/ t Numerical Solution of Differential Equations p.7/45
6 Van der Pol s equation by Euler s method: h =.84 Matlab code for van der Pol s eq: y + µ(y 2 )y + y = u x t Numerical Solution of Differential Equations p.8/45 What has gone wrong?... Absolute stability h = input( h? ); mu = 2; Xmax = 5; nmax = round(xmax/h); y = [.5 ]; yvec = [y; zeros(nmax,2)]; xvec = zeros(nmax+,); for i = :nmax yprime = [y(2) -mu*(y().ˆ2-)*y(2)-y()]; y = y + h*yprime; yvec(i+,:) = y; xvec(i+) = i*h; end plot(xvec,yvec(:,), linewidth,), grid Absolute stability xlabel( x ), ylabel( y ) Numerical Solution of Differential Equations p.9/45 Now suppose we apply α j y n+j = h β j f(x n+j,y n+j ) to this model problem (f(x,y) = λy with Re(λ) < ). Then, General solution (α j hβ j )y n+j = y n = p s (n)zs n s (α j λhβ j )y n+j = Let us write h = λh. Clearly, Re( h) <. where z s is a zero of the stability polynomial π(z,h) ρ(z) hσ(z) = (α j hβ j )z j Numerical Solution of Differential Equations p.2/45 Numerical Solution of Differential Equations p.2/45
7 Absolute stability Absolute stability Now, Therefore, we want For this to hold, we need lim y(x) = x lim y n = n A linear multistep method is called absolutely stable in an open set R A of the complex plane, if for all h R A all roots z s = z s ( h), s =,...,k, of π(z, h) satisfy z s <. The set R A is called the region of absolute stability. z s < for all s =,...,k. Numerical Solution of Differential Equations p.22/45 Numerical Solution of Differential Equations p.23/45 Exercise Find the region of absolute stability of: Euler s method Exercise Find the region of absolute stability of: Euler s method The implicit Euler method Numerical Solution of Differential Equations p.24/45 Numerical Solution of Differential Equations p.24/45
8 Region of absolute stability of Euler s method Region of absolute stability of Euler s method y n+ y n = hf(x n,y n ) Apply to y = λy, y() =, where Re(λ) <. y n+ y n = hf(x n,y n ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = Numerical Solution of Differential Equations p.25/45 Numerical Solution of Differential Equations p.25/45 Region of absolute stability of Euler s method Region of absolute stability of Euler s method y n+ y n = hf(x n,y n ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = y n+ y n = hf(x n,y n ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = π(z, h) = z h π(z, h) = z h z = + h Numerical Solution of Differential Equations p.25/45 Numerical Solution of Differential Equations p.25/45
9 Region of absolute stability of Euler s method Region of absolute stability: interior of circle y n+ y n = hf(x n,y n ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = π(z, h) = z h z = + h - z < iff + h < R A = { h C : + h < } Numerical Solution of Differential Equations p.25/45 Numerical Solution of Differential Equations p.26/45 y = λy, y() =, λ =, h =. y = λy, y() =, λ =, h = y.5 y t Numerical Solution of Differential Equations p.27/ t Numerical Solution of Differential Equations p.28/45
10 y = λy, y() =, λ =, h =. y = λy, y() =, λ =, h = y.5 y t Numerical Solution of Differential Equations p.29/ t Numerical Solution of Differential Equations p.3/45 y = λy, y() =, λ =, h =.5 y = λy, y() =, λ =, h = y y t Numerical Solution of Differential Equations p.3/ t Numerical Solution of Differential Equations p.32/45
11 Matlab code for y = λy, y() =, by Euler s method Absolute stability of the implicit Euler method h=input( h? ); lambda = -; Xmax = 3; nmax = round(xmax/h); y = ; yvec = [y; zeros(nmax,)]; xvec = zeros(nmax+,); for i = :nmax yprime = lambda*y(); y = y + h*yprime; yvec(i+,:) = y; xvec(i+) = i*h; end plot(xvec,yvec(:,), k-, linewidth,2), grid xlabel( x ), Absolute stability ylabel( y ) of the implicit Euler method Numerical Solution of Differential Equations p.33/45 y n+ y n = hf(x n+,y n+ ) Apply to y = λy, y() =, where Re(λ) <. Absolute stability of the implicit Euler method Numerical Solution of Differential Equations p.34/45 y n+ y n = hf(x n+,y n+ ) Apply to y = λy, y() =, where Re(λ) <. y n+ y n = hf(x n+,y n+ ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = z ρ(z) = z, σ(z) = z π(z, h) = z hz Numerical Solution of Differential Equations p.34/45 Numerical Solution of Differential Equations p.34/45
12 Absolute stability of the implicit Euler method Absolute stability of the implicit Euler method y n+ y n = hf(x n+,y n+ ) Apply to y = λy, y() =, where Re(λ) <. y n+ y n = hf(x n+,y n+ ) Apply to y = λy, y() =, where Re(λ) <. ρ(z) = z, σ(z) = z ρ(z) = z, σ(z) = z π(z, h) = z hz z = h π(z, h) = z hz z = h R A = { h C : h > } Region of absolute stability: exterior of circle Numerical Solution of Differential Equations p.34/45 A stability Numerical Solution of Differential Equations p.34/45 A linear multistep method is said to be A stable if its region of absolute stability, R A, contains the whole of the open left-hand complex half plane, Re( h) <. Example: The implicit Euler method is A stable. Numerical Solution of Differential Equations p.35/45 Numerical Solution of Differential Equations p.36/45
13 The Dahlquist Barrier Theorem The Dahlquist Barrier Theorem No explicit linear multistep method is A stable. No explicit linear multistep method is A stable. The order of an A stable implicit linear multistep method is 2. Numerical Solution of Differential Equations p.37/45 The Dahlquist Barrier Theorem Numerical Solution of Differential Equations p.37/45 Stepsize control: The Milne device No explicit linear multistep method is A stable. The order of an A stable implicit linear multistep method is 2. The second order A stable linear multistep method with smallest error constant is the trapezium rule method. Suppose we are using the pth order linear k-step method α j y n+j = h β j f(x n+j,y n+j ) How to choose h appropriately? Numerical Solution of Differential Equations p.37/45 Numerical Solution of Differential Equations p.38/45
14 Stepsize control: The Milne device The Milne device Suppose we are using the pth order linear k-step method α j y n+j = h β j f(x n+j,y n+j ) How to choose h appropriately? If y n,y n+,...,y n+k are error free, then y(x n+k ) y n+k = ch p+ y (p+) (x n+k ) + O(h p+2 ) h c is called the local error constant. [Proof at the back.] o Consider another pth order method: α j ỹ n+j = h β j f(x n+j,ỹ n+j ) Numerical Solution of Differential Equations p.38/45 Numerical Solution of Differential Equations p.39/45 The Milne device The Milne device If y n,y n+,...,y n+k are error free, then y(x n+k ) y n+k = ch p+ y (p+) (x n+k ) + O(h p+2 ) h c is called the local error constant. [Proof at the back.] o Similarly, Then, Hence, ỹ n+k y n+k (c c)h p+ y (p+) (x n+k ) y(x n+k ) y n+k c c c (ỹ n+k y n+k ) y(x n+k ) ỹ n+k = ch p+ y (p+) (x n+k ) + O(h p+2 ) h Suppose c c. Numerical Solution of Differential Equations p.39/45 Numerical Solution of Differential Equations p.4/45
15 Then, Hence, The Milne device ỹ n+k y n+k (c c)h p+ y (p+) (x n+k ) y(x n+k ) y n+k c c c (ỹ n+k y n+k ) The Milne device: local refinement criteria Error control per step Given a fixed tolerance δ, require that κ δ Define local refinement indicator: κ = c c c ỹ n+k y n+k Numerical Solution of Differential Equations p.4/45 Numerical Solution of Differential Equations p.4/45 The Milne device: local refinement criteria The Milne device: the adaptive algorithm Error control per step Given a fixed tolerance δ, require that. Set h κ δ Error control unit per step Given a fixed tolerance δ, require that κ hδ Numerical Solution of Differential Equations p.4/45
16 The Milne device: the adaptive algorithm The Milne device: the adaptive algorithm. Set h. Set h 3. Is κ hδ? The Milne device: the adaptive algorithm The Milne device: the adaptive algorithm. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4
17 The Milne device: the adaptive algorithm The Milne device: the adaptive algorithm. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M?. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M? YES: STOP The Milne device: the adaptive algorithm The Milne device: the adaptive algorithm. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M? YES: STOP NO: GOTO 5. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M? YES: STOP NO: GOTO 5 5. Is κ hδ?
18 The Milne device: the adaptive algorithm The Milne device: the adaptive algorithm. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M? YES: STOP NO: GOTO 5 5. Is κ hδ? YES: double h, advance x, remesh, GOTO 2. Set h 3. Is κ hδ? NO: halve h, remesh and GOTO 2 YES: GOTO 4 4. Is x X M? YES: STOP NO: GOTO 5 5. Is κ hδ? YES: double h, advance x, remesh, GOTO 2 NO: advance x, GOTO 2 Proof Proof α k y n+k hβ k f(x n+k,y n+k ) k + k α j y n+j h β j f(x n+j,y n+j ) = By the definition of the truncation error T n, and writing B = k β j, α k y(x n+k ) hβ k f(x n+k,y(x n+k )) k + k α j y(x n+j ) h β j f(x n+j,y(x n+j )) = B ht n Since y n,y n+,...,y n+k are error free, α k [y(x n+k ) y n+k ] hβ k [f(x n+k,y(x n+k )) f(x n+k,y n+k )] = B ht n α k [y(x n+k ) y n+k ] hβ k f y (x n+k,η n+k ) [y(x n+k ) y n+k ] = B ht n [ ] f α k hβ k y (x n+k,η n+k ) (y(x n+k ) y n+k ) = B ht n Numerical Solution of Differential Equations p.43/45 Numerical Solution of Differential Equations p.44/45
19 Proof As h, Therefore, Recalling that α k hβ k f y (x n+k,η n+k ) α k y(x n+k ) y n+k = Const. ht n T n = Const.h p y (p+) (x n+k ) + O(h p+ ) we deduce the desired result. o Numerical Solution of Differential Equations p.45/45
Second Order ODEs. CSCC51H- Numerical Approx, Int and ODEs p.130/177
Second Order ODEs Often physical or biological systems are best described by second or higher-order ODEs. That is, second or higher order derivatives appear in the mathematical model of the system. For
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 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 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 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 informationTHE TRAPEZOIDAL QUADRATURE RULE
y 1 Computing area under y=1/(1+x^2) 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Trapezoidal rule x THE TRAPEZOIDAL QUADRATURE RULE From Chapter 5, we have the quadrature
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 informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
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 informationTS Method Summary. T k (x,y j 1 ) f(x j 1,y j 1 )+ 2 f (x j 1,y j 1 ) + k 1
TS Method Summary Let T k (x,y j 1 ) denote the first k +1 terms of the Taylor series expanded about the discrete approximation, (x j 1,y j 1 ), and ẑ k,j (x) be the polynomial approximation (to y(x))
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationInitial Value Problems
Numerical Analysis, lecture 13: Initial Value Problems (textbook sections 10.1-4, 10.7) differential equations standard form existence & uniqueness y 0 y 2 solution methods x 0 x 1 h h x 2 y1 Euler, Heun,
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 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 information1 Implicit Differentiation
MATH 1010E University Mathematics Lecture Notes (week 5) Martin Li 1 Implicit Differentiation Sometimes a function is defined implicitly by an equation of the form f(x, y) = 0, which we think of as a relationship
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 informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationLinear Multistep Methods
Linear Multistep Methods Linear Multistep Methods (LMM) A LMM has the form α j x i+j = h β j f i+j, α k = 1 i 0 for the approximate solution of the IVP x = f (t, x), x(a) = x a. We approximate x(t) on
More informationMath 10C - Fall Final Exam
Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient
More informationAdditional exercises with Numerieke Analyse
Additional exercises with Numerieke Analyse March 10, 017 1. (a) Given different points x 0, x 1, x [a, b] and scalars y 0, y 1, y, z 1, show that there exists at most one polynomial p P 3 with p(x i )
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
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 informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 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 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 informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS ( WI3097 TU AESB0 ) Thursday April 6
More informationHigher order derivative
2 Î 3 á Higher order derivative Î 1 Å Iterated partial derivative Iterated partial derivative Suppose f has f/ x, f/ y and ( f ) = 2 f x x x 2 ( f ) = 2 f x y x y ( f ) = 2 f y x y x ( f ) = 2 f y y y
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationOrdinary differential equation II
Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
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 informationThe Dynamics of Runge Kutta Methods
The Dynamics of Runge Kutta Methods Julyan H. E. Cartwright & Oreste Piro School of Mathematical Sciences Queen Mary and Westfield College University of London Mile End Road London E1 4NS U.K. Int. J.
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 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 informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More information(6.2) -The Shooting Method for Nonlinear Problems
(6.2) -The Shooting Method for Nonlinear Problems Consider the boundary value problems (BVPs) for the second order differential equation of the form (*) y f x,y,y, a x b, y a and y b. When f x,y,y is linear
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 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 informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
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 informationOrdinary Differential Equations: Initial Value problems (IVP)
Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer
More informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
More informationHIGHER ORDER METHODS. There are two principal means to derive higher order methods. b j f(x n j,y n j )
HIGHER ORDER METHODS There are two principal means to derive higher order methods y n+1 = p j=0 a j y n j + h p j= 1 b j f(x n j,y n j ) (a) Method of Undetermined Coefficients (b) Numerical Integration
More informationCHAPTER 7. An Introduction to Numerical Methods for. Linear and Nonlinear ODE s
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 1 A COLLECTION OF HANDOUTS ON FIRST ORDER ORDINARY DIFFERENTIAL
More informationResearch Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations
Research Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations Ramses van Zon SciNet HPC Consortium November 25, 2014 Ramses van Zon (SciNet HPC Consortium)Research
More informationNumerical Algorithms for ODEs/DAEs (Transient Analysis)
Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no
More informationNumerical Methods for Differential Equations
CHAPTER 5 Numerical Methods for Differential Equations In this chapter we will discuss a few of the many numerical methods which can be used to solve initial value problems and one-dimensional boundary
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
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 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 informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
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 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 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 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 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 informationTHREE- AND FOUR-DERIVATIVE HERMITE BIRKHOFF OBRECHKOFF SOLVERS FOR STIFF ODE
THREE- AND FOUR-DERIVATIVE HERMITE BIRKHOFF OBRECHKOFF SOLVERS FOR STIFF ODE By Njwd Albishi February 2016 A Thesis submitted to the Faculty of Graduate and Postdoctoral Studies in partial fulfillment
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0
More informationVariational Methods & Optimal Control
Variational Methods & Optimal Control lecture 12 Matthew Roughan Discipline of Applied Mathematics School of Mathematical Sciences University of Adelaide April 14, 216
More informationA New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations
A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations Rei-Wei Song and Ming-Gong Lee* d09440@chu.edu.tw, mglee@chu.edu.tw * Department of Applied Mathematics/
More informationFIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS II: Graphical and Numerical Methods David Levermore Department of Mathematics University of Maryland
FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS II: Graphical and Numerical Methods David Levermore Department of Mathematics University of Maryland 9 January 0 Because the presentation of this material in
More informationExercises, module A (ODEs, numerical integration etc)
FYTN HT18 Dept. of Astronomy and Theoretical Physics Lund University, Sweden Exercises, module A (ODEs, numerical integration etc) 1. Use Euler s method to solve y (x) = y(x), y() = 1. (a) Determine y
More informationRunge Kutta Collocation Method for the Solution of First Order Ordinary Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 1, 17-26 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.5823 Runge Kutta Collocation Method for the Solution of First
More information16.7 Multistep, Multivalue, and Predictor-Corrector Methods
740 Chapter 16. Integration of Ordinary Differential Equations 16.7 Multistep, Multivalue, and Predictor-Corrector Methods The terms multistepand multivaluedescribe two different ways of implementing essentially
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 informationSolving ODEs and PDEs in MATLAB. Sören Boettcher
16.02.2009 Introduction Quick introduction to syntax ODE in the form of Initial Value Problems (IVP) what equations can handle how to code into how to choose the right solver how to get the solver to do
More informationMultivariable Calculus and Matrix Algebra-Summer 2017
Multivariable Calculus and Matrix Algebra-Summer 017 Homework 4 Solutions Note that the solutions below are for the latest version of the problems posted. For those of you who worked on an earlier version
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationDefinition of differential equations and their classification. Methods of solution of first-order differential equations
Introduction to differential equations: overview Definition of differential equations and their classification Solutions of differential equations Initial value problems Existence and uniqueness Mathematical
More informationEfficient numerical methods for the solution of stiff initial-value problems and differential algebraic equations
10.1098/rspa.2003.1130 R EVIEW PAPER Efficient numerical methods for the solution of stiff initial-value problems and differential algebraic equations By J. R. Cash Department of Mathematics, Imperial
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 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 informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationLet f(x) = x, but the domain of f is the interval 0 x 1. Note
I.g Maximum and Minimum. Lagrange Multipliers Recall: Suppose we are given y = f(x). We recall that the maximum/minimum points occur at the following points: (1) where f = 0; (2) where f does not exist;
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 informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
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 informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
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 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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
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 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 informationFinal Exam Advanced Mathematics for Economics and Finance
Final Exam Advanced Mathematics for Economics and Finance Dr. Stefanie Flotho Winter Term /5 March 5 General Remarks: ˆ There are four questions in total. ˆ All problems are equally weighed. ˆ This is
More informationDEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS. MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 2018
DEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 208 Version A refers to the regular exam and Version B to the make-up. Version A. A particle
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 informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationPractice Midterm 2 Math 2153
Practice Midterm 2 Math 23. Decide if the following statements are TRUE or FALSE and circle your answer. You do NOT need to justify your answers. (a) ( point) If both partial derivatives f x and f y exist
More informationOn weighted Adams-Bashforth rules
MATHEMATICAL COMMUNICATIONS 127 Math. Commun. 23(218), 127 144 On weighted Adams-Bashforth rules Mohammad Masjed-Jamei 1, Gradimir V. Milovanović 2,3, and Amir Hossein Salehi Shayegan 1 1 Department of
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationLecture 28 The Main Sources of Error
Lecture 28 The Main Sources of Error Truncation Error Truncation error is defined as the error caused directly by an approximation method For instance, all numerical integration methods are approximations
More information