Numerical Analysis II. Problem Sheet 12
|
|
- Aleesha Dorsey
- 5 years ago
- Views:
Transcription
1 P. Grohs S. Hosseini Ž. Kereta Spring Term 2015 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 12 Problem 12.1 The Exponential Runge-Kutta Method Consider the exponential Runge Kutta single-step method k i := φ(cha) ( i 1 ) f(u i ) + ha c ij k j, i = 1,..., s, j=1 i 1 u i := y 0 + h a ij k j i = 1,..., s, (12.1.1) Ψ h y 0 := y 0 + h j=1 s b i k i. i=1 for the autonomous differential equation ẏ = f(y), where A := Df(y 0 ) and (12.1a) φ(z) = exp(z) 1. z Show that the two-step exponential Runge Kutta method, with the parameter values c = 1 2, a 21 = α, c 21 = 3 4 α2 α, b 1 = 1 1 3α, b 2 2 = 1 3α, (12.1.2) 2 and α free to be chosen, solves the linear inhomogeneous differential equation ẏ = Ay + g, A R d d, g R d exactly. Solution: We must show that the method applied to the linear inhomogeneous differential equation ẏ = Ay + g, A R d d, g R d gives us the exact solution y(t) = exp(at)y 0 + t 0 exp(a(t τ))g dτ = exp(at)y 0 + (e ta 1)A 1 g. The procedure of the 2 step Runge Kutta method is (for the linear case, Df(y 0 ) = A) u 1 = y 0 k 1 = φ(cha)f(y 0 ) = φ(cha)(ay 0 + g) u 2 = y 0 + ha 21 k 1 k 2 = φ(cha)(f(y 0 + ha 21 k 1 ) + hac 21 k 1 ) = φ(cha)(ay 0 + h (a 21 + c 21 ) }{{} = 3 4 α2 Ak 1 + g), Problem Sheet 12 Page 1 Problem 12.1
2 and 2 Ψ h y 0 = y 0 + h b i k i. i=1 With the parameter values we get c = 1 2, a 21 = α, c 21 = 3 4 α2 α, b 1 = 1 1 3α 2, b 2 = 1 3α 2, ( Ψ h y 0 = y 0 + hφ(cha) (1 1 3α )(Ay g)+ ) 1 3α (Ay hα2 Ak 1 + g) ( = y 0 + hφ(cha) Ay 0 + g + h ) 4 Ak 1 ( = y 0 + hφ(cha) Ay 0 + g + h ) 4 Aφ(chA)(Ay 0 + g) = y 0 + 2(exp (ha/2) I)A 1 (Ay 0 + g + ) h 2 AA 1 (exp (ha/2) I)(Ay 0 + g) = y 0 + 2(exp (ha/2) I)y 0 + 2(exp (ha/2) I)A 1 g+ (exp (ha/2) I)A 1 A(exp ha/2 I)(y 0 + A 1 g) = exp(ah)y 0 + (e ha 1)A 1 g. (12.1b) Implement a MATLAB function function [t,y] = exponentialrk2(f,df,y0,t,h,alpha) which solves the autonomous differential equation ẏ = f(y), y(0) = y 0 with the 2-step Runge Kutta method (12.1.1) and parameter values (12.1.2). HINT: Even though the method (12.1.1) solves autonomous differential equations, write your implementation exponentialrk2 for the general system ẏ = f(t, y) so it is compatible with MATLAB ODE-solvers. Solution: See Listing 12.1 Listing 12.1: Implementation of exponentialrk2.m from subproblem (12.1b) 1 f u n c t i o n [t,y] = ExponentialRK2(f,df,y0,T,h,alpha) 2 % this code solves the differential equation y =F(y) and Jacobi DF using a 3 % 2-stage Exponential Runge-Kutta method. Problem Sheet 12 Page 2 Problem 12.1
3 4 % 5 6 i f nargin < 6; 7 alpha = 1; 8 end 9 10 % Initialize variables 11 N = c e i l ((T(2)-T(1))/h); % number of steps 12 h = (T(2)-T(1))/N; % adjust step size to fit interval length 13 d = s i z e (y0,1); % dimension of solution 14 t = l i n s p a c e (T(1),T(2),N+1); % time grid 15 y = nan(d,n+1); % solution 16 y(:,1) = y0; % set initial value % Compute solution 19 f o r j=1:n 20 y(:,j+1) = ExponentialRK2Step(f,df,y(:,j),t,h,alpha); 21 end end f u n c t i o n y = ExponentialRK2Step(f,df,y0,t,h,alpha) c = 0.5; 28 a21 = alpha; 29 c21 = 0.75*alphaˆ2 - alpha; b = nan( s i z e (y0)); 32 b(2) = 1/(3*alphaˆ2); 33 b(1) = 1 - b(2); DF = df(t,y0); phi = ( expm(c*h*df)-eye( s i z e (DF)) )/(c*h*df); u = nan( s i z e (y0)); 40 k = nan( s i z e (y0)); u(:,1) = y0; 43 k(:,1) = phi*f(t,u(:,1)); u(:,2) = y0 + h*a21*k(:,1); 46 k(:,2) = phi*( f(t,u(:,2)) + h*df*c21*k(:,1) ); y = y0 + h*(b(1)*k(:,1) + b(2)*k(:,2)); end Problem Sheet 12 Page 3 Problem 12.1
4 (12.1c) Write a MATLAB function which determines the convergence rate of the method numerically. For your numerical experiment, solve Zeeman s heartbeat model l = l 3 + α heart l p ṗ = βl (12.1.3) with the parameter values α heart = 3, β = 0.01, initial values (p(0), l(0)) = (0, 1) up to the end point T = 100. Use the step sizes h = 1, k = 2,..., 6, and as a reference solution, use the 2 k results given by ode15s with very high accuracy, f.e. AbsTol=1e-12 and RelTol=1e-12. Solution: See Listing 12.2 for the implementation. See Figure 12.1 and Figure 12.2 for the solution graph and the convergence loglog plot respectively. We get an estimated convergence rate of f u n c t i o n ExponentialRK2Conv Listing 12.2: Implementation of subproblem (12.1c) 2 % this code calculates the order of exponentialrk2 using a reference 3 % solution for Zeeman s Heartbeat model obtained using ode23s. 4 5 a l p h a = 0.75; % free-parameter in ExponentialRK beta = 0.01; 8 heartalpha = 3; 9 10 % Zeeman s Heartbeat model: Example 11, Section f [ -y(1)ˆ3 + heartalpha*y(1) - y(2); 12 beta*y(1)]; df [ -3*y(1)ˆ2 + heartalpha, -1; 15 beta, 0 ]; y0 = [1; 0]; 18 T = [0 100]; 19 h = 2.ˆ-[2:6]; options = odeset( RelTol,1e-12, AbsTol,1e-12, Jacobian,df); 22 [tgood,ygood] = ode23s(f,t,y0,options); disp( ode23s done... ); ygood = ygood. ; % matlab compatibility err = nan( s i z e (h)); f o r j=1: l e n g t h (h) Problem Sheet 12 Page 4 Problem 12.1
5 31 [t,y] = exponentialrk2(f,df,y0,t,h(j), a l p h a); 32 err(j) = norm( y(1,end) - ygood(1,end) ); 33 end f i g u r e ; 36 p l o t (t,y(1,:), ro ); 37 hold on; 38 p l o t (tgood, ygood(1,:), g-, LineWidth,1.1); 39 t i t l e ( Solution to Zeeman s Heartbeat Model using Exponential RK2 ); 40 x l a b e l ( t ); 41 y l a b e l ( y(t) ); 42 legend([ Exponential RK2, h= num2str(h(end))], ode23s ); f i g u r e ; 45 l o g l o g (h,err, bx: ) 46 grid on; 47 t i t l e ( Log-Log Plot of Error against Step-Size for Exponential RK2 ); 48 x l a b e l ( h ); 49 y l a b e l ( error ); fit = p o l y f i t ( l o g (h), l o g (err),1); 52 f p r i n t f ( calculated numerical convergence orders: %f\n,fit(1)); Problem 12.2 Robustness of L-Stable Method In this exercise we investigate the robustness of L-stable 1-step methods for the scalar linear model problem [NUMODE, Eq. (3.1.1)]. Robustness denotes the special property of a numerical integrator, namely that the integration error can be reasonably estimated independently of a parameter in the differential equation. We apply the two step Radau-2-method ( [NUMODE, Eq. (3.6.4)]) to the model problem ẏ = λy, y(0) = 1. (12.2.1) Hence we construct an equidistant mesh with step size h = 1 N on the time interval [0, 1]. (12.2a) Determine an analytic expression for e N (λ) := y N y(1), (12.2.2) by expressing the numerical solution using the stability function ( [NUMODE, Thm ]) of the method. Solution: By [NUMODE, Eq. (3.1.16)] we have y N = (Ψ h ) N y(0) = S(hλ) N y(0) = S( λ N )N, Problem Sheet 12 Page 5 Problem 12.2
6 Solution to Zeeman s Heartbeat Model using Exponential RK2 Exponential RK2, h= ode23s y(t) t Figure 12.1: subproblem (12.1c), Solution graph 10 4 Log Log Plot of Error against Step Size for Exponential RK error h Figure 12.2: subproblem (12.1c), convergence rate plot Problem Sheet 12 Page 6 Problem 12.2
7 where S( ) denotes the stability function of the Radau-2-method. This gives e N (λ) = S( λ N )N e λ. For the stability function [NUMODE, Thm ] implies S(z) = z z z2. (12.2b) Write a MATLAB function which determines a lower bound for by sampling on the interval [0, N]. Solution: See Listing f u n c t i o n e = sup_e(n) e = sup e(n), sup e N (λ) 0 λ N Listing 12.3: subproblem (12.2b) 2 % determine lower bound for the sup over e_n(lambda) for 3 % 0 <= lambda <= N 4 % 5 % Input: 6 % N: number of steps, see above 7 % 8 % Output: 9 % e: lower bound for the sup 10 % 11 % call: 12 % e = sup_e(10); % stability function of the Radau-2- method 15 S (1 + z/3)./ (1 + (z/2-2).* z/3); % generate net for the sampling 18 lambda = l i n s p a c e (0, N, 5000); % find max 21 e = max( abs( S(lambda/N).ˆN - exp(lambda) ) ); (12.2c) Give a constant C, independent of h and λ, such that y N y(1) Ch. (12.2.3) Solution: No C with those properties exists because sup e N (λ) e N (N) = S(1) N e N for N = 1. h 0 λ N Problem Sheet 12 Page 7 Problem 12.2
8 (12.2d) State a 1-step method which when applied to (12.2.1) with uniform time step size h has the property (12.2.3) with C independent of h and λ. Solution: The exponential Euler method has the desired property since its stability function is given by S(z) = exp(z) and hence Problem 12.3 e N (λ) = S( λ N )N exp(λ) = exp( λ N )N exp(λ) = 0. Exact Quadrature on P 2s 2 Implies Positive Weights It is often important to know whether the coefficients b i in the Butcher-scheme of a RK single step method are positive. Let the quadrature formula 1 0 f(x) dx s b i f(c i ) i=1 c i c j be exact for polynomials f P 2s 2. Show that the weights b i are necessarily positive. HINT: Define to the s different supporting points suitable polynomials p i (x) P 2s 2 with p i (x) 0. Solution: We consider the polynomials p k (x) = s i=1,i k (x c i ) 2 0 x. Then we have 0 < 1 0 p k (x) dx = s b i p k (c i ) = b k p k (c k ). i=1 So that b k > 0, because p k (c k ) > 0. Published on 12 May To be submitted by 19 May References [NUMODE] Lecture Slides for the course Numerical Methods for Ordinary Differential Equations, SVN revision # Last modified on May 15, 2015 Problem Sheet 12 Page 8 Problem 12.3
Numerical Analysis II. Problem Sheet 9
H. Ammari W. Wu S. Yu Spring Term 08 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 9 Problem 9. Extrapolation of the Implicit Mid-Point Rule. Taking the implicit mid-point rule as a basis method,
More informationNumerical Analysis II. Problem Sheet 8
P. Grohs S. Hosseini Ž. Kereta Spring Term 15 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 8 Problem 8.1 The dierential equation Autonomisation and Strang-Splitting ẏ = A(t)y, y(t ) = y, A(t)
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
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 informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 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
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 informationNumerical Methods for CSE. Examination - Solutions
R. Hiptmair A. Moiola Fall term 2010 Numerical Methods for CSE ETH Zürich D-MATH Examination - Solutions February 10th, 2011 Total points: 65 = 8+14+16+16+11 (These are not master solutions but just a
More informationTest you method and verify the results with an appropriate function, e.g. f(x) = cos 1 + sin
Advanced methods for ODEs - Lab exercises Verona, 20-30 May 2015 1. FFT and IFFT in MATLAB: Use the Fast Fourier Transform to compute compute the k-th derivative of an appropriate function. In MATLAB the
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationLast name Department Computer Name Legi No. Date Total. Always keep your Legi visible on the table.
First name Last name Department Computer Name Legi No. Date 02.02.2016 1 2 3 4 5 Total Grade Fill in this cover sheet first. Always keep your Legi visible on the table. Keep your phones, tablets and computers
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationVector Fields and Solutions to Ordinary Differential Equations using Octave
Vector Fields and Solutions to Ordinary Differential Equations using Andreas Stahel 6th December 29 Contents Vector fields. Vector field for the logistic equation...............................2 Solutions
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
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 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 informationDifferential Equations FMNN10 Graded Project #1 c G Söderlind 2017
Differential Equations FMNN10 Graded Project #1 c G Söderlind 2017 Published 2017-10-30. Instruction in computer lab 2017-11-02/08/09. Project report due date: Monday 2017-11-13 at 10:00. Goals. The goal
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
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 informationSymplectic integration with Runge-Kutta methods, AARMS summer school 2015
Symplectic integration with Runge-Kutta methods, AARMS summer school 2015 Elena Celledoni July 13, 2015 1 Hamiltonian systems and their properties We consider a Hamiltonian system in the form ẏ = J H(y)
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 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 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 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 information24, B = 59 24, A = 55
Math 128a - Homework 8 - Due May 2 1) Problem 8.4.4 (Page 555) Solution: As discussed in the text, the fourth-order Adams-Bashforth formula is a formula of the type x n+1 = x n + h[af n + Bf n 1 + Cf n
More informationMAT 275 Laboratory 6 Forced Equations and Resonance
MAT 275 Laboratory 6 Forced Equations and Resonance In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate on equations with a periodic harmonic forcing
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 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 / 33 Almost Done! Last
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 informationNumerical Methods - Boundary Value Problems for ODEs
Numerical Methods - Boundary Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Boundary Value Problems for ODEs 2013 1 / 14 Outline 1 Boundary Value
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
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 informationChecking the Radioactive Decay Euler Algorithm
Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present
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 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 informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
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 informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationLaboratory 10 Forced Equations and Resonance
MATLAB sessions: Laboratory 9 Laboratory Forced Equations and Resonance ( 3.6 of the Edwards/Penney text) In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE201: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE21: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationMath 308 Week 8 Solutions
Math 38 Week 8 Solutions There is a solution manual to Chapter 4 online: www.pearsoncustom.com/tamu math/. This online solutions manual contains solutions to some of the suggested problems. Here are solutions
More informationTheory, Solution Techniques and Applications of Singular Boundary Value Problems
Theory, Solution Techniques and Applications of Singular Boundary Value Problems W. Auzinger O. Koch E. Weinmüller Vienna University of Technology, Austria Problem Class z (t) = M(t) z(t) + f(t, z(t)),
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 informationIntroduction to Functional Analysis
Introduction to Functional Analysis Carnegie Mellon University, 21-640, Spring 2014 Acknowledgements These notes are based on the lecture course given by Irene Fonseca but may differ from the exact lecture
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 informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #5
University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #5 M.G. Lipsett, Updated 00 Integration Methods with Higher-Order Truncation Errors with MATLAB MATLAB is capable
More information2D1240 Numerical Methods II / André Jaun, NADA, KTH NADA
ND 2D24 Numerical Methods II / ndré Jaun, ND, KTH 2. INITIL- / OUNDRY VLUE PROLEMS 2.. Remember: what we saw during the last lesson Euler and Runge-Kutta methods for solving ODEs Predator-prey model using
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 informationModified Defect Correction Algorithms for ODEs. Part II: Stiff Initial Value Problems
Modified Defect Correction Algorithms for ODEs. Part II: Stiff Initial Value Problems Winfried Auzinger Harald Hofstätter Wolfgang Kreuzer Ewa Weinmüller to appear in Numerical Algorithms ANUM Preprint
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 informationHomogeneous Equations with Constant Coefficients
Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
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 informationRunge - Kutta Methods for first and second order models
Runge - Kutta Methods for first and second order models James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 3, 2013 Outline 1 Runge -
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More information13 Numerical Solution of ODE s
13 NUMERICAL SOLUTION OF ODE S 28 13 Numerical Solution of ODE s In simulating dynamical systems, we frequently solve ordinary differential equations. These are of the form dx = f(t, x), dt where the function
More information2 Solving Ordinary Differential Equations Using MATLAB
Penn State Erie, The Behrend College School of Engineering E E 383 Signals and Control Lab Spring 2008 Lab 3 System Responses January 31, 2008 Due: February 7, 2008 Number of Lab Periods: 1 1 Objective
More information4th Preparation Sheet - Solutions
Prof. Dr. Rainer Dahlhaus Probability Theory Summer term 017 4th Preparation Sheet - Solutions Remark: Throughout the exercise sheet we use the two equivalent definitions of separability of a metric space
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationMAT 275 Laboratory 6 Forced Equations and Resonance
MATLAB sessions: Laboratory 6 MAT 275 Laboratory 6 Forced Equations and Resonance In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate on equations with
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 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 informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (35 points) Leslie Leroy Irvin bails out of an airplane at the altitude of 16,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s
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 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 informationRunge-Kutta and Collocation Methods Florian Landis
Runge-Kutta and Collocation Methods Florian Landis Geometrical Numetric Integration p.1 Overview Define Runge-Kutta methods. Introduce collocation methods. Identify collocation methods as Runge-Kutta methods.
More informationA First Course on Kinetics and Reaction Engineering Supplemental Unit S5. Solving Initial Value Differential Equations
Supplemental Unit S5. Solving Initial Value Differential Equations Defining the Problem This supplemental unit describes how to solve a set of initial value ordinary differential equations (ODEs) numerically.
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 informationLecture 17: Ordinary Differential Equation II. First Order (continued)
Lecture 17: Ordinary Differential Equation II. First Order (continued) 1. Key points Maple commands dsolve dsolve[interactive] dsolve(numeric) 2. Linear first order ODE: y' = q(x) - p(x) y In general,
More informationSection 2.1 Differential Equation and Solutions
Section 2.1 Differential Equation and Solutions Key Terms: Ordinary Differential Equation (ODE) Independent Variable Order of a DE Partial Differential Equation (PDE) Normal Form Solution General Solution
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 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 information1 Systems of First Order IVP
cs412: introduction to numerical analysis 12/09/10 Lecture 24: Systems of First Order Differential Equations Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore 1 Systems
More informationSolving systems of ODEs with Matlab
Solving systems of ODEs with Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 20, 2013 Outline 1 Systems of ODEs 2 Setting Up
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
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 informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationInverse Problems D-MATH FS 2013 Prof. R. Hiptmair. Series 1
Inverse Problems D-MATH FS 2013 Prof. R. Hiptmair Series 1 1. Differentiation in L 2 per([0, 1]) In class we investigated the inverse problem related to differentiation of continous periodic functions.
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 informationCS 205b / CME 306. Application Track. Homework 2
CS 205b / CME 306 Application Track Homework 2 1. ALE An Eulerian formulation of conservation of mass uses control volumes that are fixed in space as material flows freely through the control volumes.
More informationRunge - Kutta Methods for first order models
Runge - Kutta Methods for first order models James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 30, 2013 Outline 1 Runge - Kutte Methods
More information2.1 Differential Equations and Solutions. Blerina Xhabli
2.1 Math 3331 Differential Equations 2.1 Differential Equations and Solutions Blerina Xhabli Department of Mathematics, University of Houston blerina@math.uh.edu math.uh.edu/ blerina/teaching.html Blerina
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
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 informationProblem Sheet 0: Numerical integration and Euler s method. If you find any typos/errors in this problem sheet please
Problem Sheet 0: Numerical integration and Euler s method If you find any typos/errors in this problem sheet please email jk08@ic.ac.uk. The material in this problem sheet is not examinable. It is merely
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
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 informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
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 informationLösning: Tenta Numerical Analysis för D, L. FMN011,
Lösning: Tenta Numerical Analysis för D, L. FMN011, 090527 This exam starts at 8:00 and ends at 12:00. To get a passing grade for the course you need 35 points in this exam and an accumulated total (this
More information