MULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS. Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University
|
|
- Earl Stephens
- 5 years ago
- Views:
Transcription
1 MULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University HU-1117 Budapest, Pázmány Péter sétány 1/c. Abstract Here we present a numerical method to solve MBVP s of ODE s and its computer adaptation developed as software packages written in Turbo Pascal and C languages. 1. Introduction Boundary Value Problems of Ordinary Differential Equations appear in a large domain of sciences and many practically important problems lead to Multipoint Boundary Value Problems (MBVP s). Such is the case, e.g., for modelling and analysing problems arising from electric power networks, electric railway systems, telecommunication lines and also in chemistry, analysing kinetical reaction problems. Applications for difference frequency domain problems, e.g. harmonic penetration in unbalanced power networks, railway circuits with auto- or booster transformers, telecommunication circuits have been studied at the Technical University in Budapest and at the Swedish Transmission Research Institute in Sweden and problems have been solved ([1],[2]) using the algorithm and computer programs presented in this article. 1
2 The class of problems studied here has the following form: (1) x (t) = f(t, x(t)), t [t 1, t m ] R, x : [t 1, t m ] R n, f : [t 1, t m ] R n R n, n, m N +, m > 1, subject to multipoint boundary conditions (2) r(x(t 1 ), x(t 2 +),..., x(t m 1 +), x(t m )) = 0, r : (R n ) m R n, t 1 < t 2 <... < t m, (3) g j (x(t j ), x(t j +)) = 0, g j : (R n ) 2 R n, (j =2, 3,..., m 1), where x(t j ) and x(t j +) mean the left and right limits of x at point t j, respectively (j =2, 3,..., m 1). We assume that f is continuous and twice continuously differentiable with respect to its second variable on the regions [t j, t j+1 ] R n (j = 1, 2,..., m 1), and the functions r and g j are continuously differentiable. These requierements can be weakened, but for simplicity and since the applications generally satisfy these assumptions, we suppose them to be fulfilled. 2. Numerical algorithm To solve the MBVP (1) (2) (3) we apply an iteration method, based on a modified shooting method ([3],[4]). In each iteration step, we solve the initial value problems (4) x (t) = f(t, x(t)), t [t j, t j+1 ], (5) x(t j ) = s j, (j =1, 2,..., m 1), where the initial values s j R n are given in the first iteration step using one of the strategies described below (at paragraph 3.), and determined in the following iteration steps by the formula (12). 2
3 Let us denote the solutions of (4) (5) by ξ j ( ; s j ), j =1, 2,..., m 1. We have to determine the elements s j in such a way that the function ξ : [t 1, t m ] R n, defined by (6) ξ(t) := ξ j ( t ; s j ), t [ t j, t j+1 ), ξ(t m ) := ξ m 1 ( t m ; s m 1 ), satisfies the boundary conditions (2) (3), that is (7) r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ) = 0, (8) g j ( ξ j 1 ( t j ; s j 1 ), s j ) = 0 (j =2, 3,..., m 1), hence ξ is a solution of the problem (1) (2) (3). Thus we have to find a solution of the equation F (s 1, s 2,..., s m 1 ) = 0, where F is defined by F : (R n ) m 1 (R n ) m 1, (9) F (s 1, s 2,..., s m 1 ) := r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ), g 2 ( ξ 1 ( t 2 ; s 1 ), s 2 ),.. g m 1 ( ξ m 2 ( t m 1 ; s m 2 ), s m 1 ) To solve this equation we use the general Newton s iteration (10) (s 1, s 2,..., s m 1 ) (i+1) := (s 1, s 2,..., s m 1 ) (i) (DF ((s 1, s 2,..., s m 1 ) (i) )) 1 F ((s 1, s 2,..., s m 1 ) (i) ). Therefore, we have to compute F (s 1, s 2,..., s m 1 ), DF (s 1, s 2,..., s m 1 ) in each iteration step, and to solve the system of linear equations (11) DF (s 1, s 2,..., s m 1 ) ( s 1, s 2,..., s m 1 ) = = F (s 1, s 2,..., s m 1 ) and finally, using the solution of the above system we get the initial values s j R n for the next iteration step by 3
4 (12) s (i+1) j := s (i) j + s j (j =1, 2,..., m 1). To compute DF (s 1, s 2,..., s m 1 ) we determine the partial derivatives of F with respect to the k th components of its j th variable s j, (j = 1, 2,..., m 1, k = 1, 2,..., n). Since the functions ξ j ( ; s j ) do not depend on s l if j l and 2k ξ ji ( t j ; s j ) = δ ki, the [(j 1) n + k] th column of the Jacobian matrix can be written (13) for j =1 as 1k F (s 1, s 2,..., s m 1 ) = 1k r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ) n 1i g 2 ( ξ 1 ( t 2 ; s 1 ), s 2 ) 2k ξ 1i (t 2 ; s 1 ) i=1 0 n }.. (m 3) 0 n (14) for j =2,..., m 2 as jk F (s 1, s 2,..., s m 1 ) = jk r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ) 0 n }.. (j 2) 0 n 2k g j ( ξ j 1 ( t j ; s j 1 ), s j ) n 1i g j+1 ( ξ j ( t j+1 ; s j ), s j+1 ) 2k ξ ji (t j+1 ; s j ) 0 n }.. (m 2 j) 0 i=1 4 n
5 (15) and for j =m 1 as (m 1)k F (s 1, s 2,..., s m 1 ) = (m 1)k r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ) + + n mi r( s 1, s 2,..., s m 1, ξ m 1 ( t m ; s m 1 ) ) i=1 2k ξ (m 1)i (t m ; s m 1 ) 0 n }. (m 3) 0 n 2k g m 1 ( ξ m 2 ( t m 1 ; s m 2 ), s m 1 ) where [ 0 ] is the zero column of n Rn ; ji r, 1i g j and 2i g j are the partial derivatives of r and g j with respect to the i th component of their j th, first and second variables, respectively. To compute the partial derivatives 2k ξ ji (t j+1 ; s j ) for j =1, 2,..., m 1, k = 1, 2,..., n, i = 1, 2,..., n, we solve the following linear initial value problems (16) y (t) = 2 f( t, ξ j ( t ; s j ) ) y(t), t [ t j, t j+1 ], (17) y(t j ) = e k, (j =1, 2,..., m 1), (k =1, 2,..., n), where e k is the k th unit vector in R n. Denoting the solutions by η j,k we have 2k ξ ji ( ; s j ) = (η j,k ) i ( j = 1, 2,..., m 1, k = 1, 2,..., n, i=1, 2,..., n). Let us summarize the main steps of the algorithm: We have to prepare the first iteration by choosing the starting elements s j, (the initial values of the Cauchy problems (4) (5)), then 5
6 in each iteration step we have to (i) solve the Cauchy problems (4) (5) and (16) (17) simultaneously, (ii) compute F (s 1, s 2,..., s m 1 ) using (9) and the solutions of the Cauchy problems (4) (5), (iii) compute DF (s 1, s 2,..., s m 1 ) using (13), (14), (15) and the solutions of the Cauchy problems (4) (5) and (16) (17), (iv) determine the elements s 1, s 2,..., s m 1 by solving the system of linear equations (11), (v) compute the new elements s j by the formula (12), (vi) test if these new elements s j can be taken as good approximations for the solution of the problem (1) (2) (3) or initial values (5) for the next iteration step. Since the Newton method is only locally convergent, we replace (12) by the correspondent formula of the modified Newton method ([4]), that is, in the step (v) we compute the new elements s j by (18) s (i+1) j := s (i) j + λ s j (j =1, 2,..., m 1), where λ is determined so that the method converge to the solution of the problem (1) (2) (3), provided it exists ([4],[5]). 3. Computer adaptation The main steps of crucial importance in the algorithm which have to be worked out very carefully in the computer programmation are as follows: (3.1) the choice of the values s j (j =1, 2,..., m 1) in the first iteration, (3.2) the numerical integration to solve the Cauchy problems (4) (5) and (16) (17), (3.3) the numerical method to solve the system of linear equations (11), (3.4) the numerical method to determine λ in the formula (18) of the modified Newton method. 6
7 The choice of the values s j (j =1, 2,..., m 1) in the first iteration step, and the computation of these values in the following iteration steps, have crucial importance not only for the convergence of the method. In some practical problems the partial derivative 2 f is unbounded on the regions [t j, t j+1 ] R n, and in such cases the solutions of the initial value problems (4) (5) and (16) (17) for certain values of s j can only be defined in a neighborhood of t j, and not on the whole intervals [t j, t j+1 ]. So the method would break down in these cases. Moreover, even if the initial value problems (4) (5) and (16) (17) can be solved in principle, in practice if the solutions depend very sensitively on the initial values, there can be considerable inaccuracies at the end of the intervals [t j, t j+1 ]. The computation even to full machine accuracy and by quasi-exact initial value problem solvers (codes) does not guarantee that the computed solutions can be determined accurately. These problems can be demonstrated by the well-known estimate (3.5) ξ j (t; s 1 j) ξ j (t; s 2 j) s 1 j s 2 j e L j t t j, where ξ j (t; s 1 j ) and ξ j(t; s 2 j ) are the values at t of the solutions of the initial value problem (4) (5) with s j := s 1 j and s2 j, respectively; and L j is a Lipschitz constant for the function f with respect to its second variable on the region [t j, t j+1 ] R n. The estimate (3.5) shows that the influence of inaccurate initial values s j can be reduced arbitrarily small if we define sufficiently close points in the intervals [t j, t j+1 ], t j =: t j1 < t j2 <... < t jl := t j+1, such that the values e L j t jk+1 t jk are sufficiently small. We can consider these intermediate points as if were given points (together with the endpoints t j and t j+1 ) for the MBVP problem (1) (2) (3). Simple continuity conditions have to be given for these newly defined points, such as (3.6) x(t jk +) x(t jk ) = 0, thus the boundary conditions (3) have to be completed by the mapping (3.7) g jk : (R n ) 2 R n, (u, v) u v. for the newly generated points t jk. 7
8 The function r in the boundary condition (2) has to be also modified. Its domain has to be completed corresponding to the newly defined intermediate points, but the mapping really does not change. Since we do not know the Lipschitz constants L j, in each iteration step we generate new t jk points if the magnitude of the solution of the initial value problems (4) (5) or (16) (17) at some point t becomes greater than a user prescribed value. While integrating the differential equations (4) and (16) on the intervals [t j, t j+1 ], we examine whether the values of ξ j (t; s j ) and η j,k (t) become too large to s j and 1, respectively. If (3.8) ξ j ( t; s j ) > γ ( s j + 1 ) or η j,k (t) > γ, then we terminate the integration and generate the point t as a new intermediate point t jk. The parameter γ has to be be given by the user. To continue the integrations on the interval [t jk, t j+1 ] the initial values in (5) are computed by linear interpolation from the output trajectory of the previous iteration or from the starting trajectory in the first iteration. It is important to mention that (i) (ii) (iii) for numerical reasons, the program uses greater γ in the second and in the following iterations than in the first iteration, sometimes, especially when large derivatives are expected for the trajectories, it is recommended to put intermediate points as given points of the multipoint boundary conditions (2) (3), the given points t j together with the program generated intermediate points t jk will be called shooting points in the following. The first estimations for the solution at the given points, i.e. the values of s j (j = 1, 2,..., m 1), have to be given by the user. This can easily be done when the qualitative behavior of the solution is known by the phisycal meaning of the problem, so at least a rough approximation can be found. Howewer, in complicated cases, the first approximation must be treated very carefully. In many cases the so-called continuation 8
9 or homotopy method is the most effective method, where we solve the problem gradually. Starting with a neighboring simpler problem, we go toward the solution of the actually posed problem, step by step, using the output result of the previous step as input data for the actually executed step. The homotopy procedure can easily be programmed by the following possibilities to give the first approximations. (3.9) The first approximative trajectory for the differential equation (1) is the constant zero function, thus, all components of s j equal zero in the first iteration (j = 1, 2,..., m 1). In this case the user does not have to prepare data before entering the program. (3.10) The first estimations of s j for the initial value problems (4) (5) are computed by interpolation from a data file which has to be given by the user. Data must be written in matrix form, the rows should contain points τ k (at which approximations are given for the trajectory) and the estimated components of ξ(τ k ). At least two points, t b and t e have to be given (together with the corresponding estimations of the trajectory) in such a way that the interval [t b, t e ] includes all given points t j. (3.11) The first estimations of s j for the Cauchy problems (4) (5) are interpolated from the results at the shooting points of the previous run, so this can be chosen only when the program run previously. The user does not have to give estimations, since the program automatically creates the necessary input data file from the results of the previous run. (3.12) The first approximations of the trajectory have to be found in a data file like in case (3.10). The set of given points will be enlarged so that the shooting points of the previous run will be embedded in the set of the actually given points. The estimations at these points are interpolated from the input data file, which is recommended to be the output results of the previous run. 9
10 (3.13) The program creates the given points t j together with the approximations s j automatically so that they will be exactly the same as they were in the last iteration of the previous run. The user does not have to give approximations, since the program utilizes the output of the previous run as input for the actual run. We emphasize that in the cases (3.9), (3.11) and (3.13) the program automatically creates the necessary input data file (by the constant zero function in (3.9) or by the output results at the shooting points of the previous run in (3.11) and (3.13) ), while in the cases (3.10) and (3.12) the input data file has to be created by the user. When the user applies the homotopy method, we propose to use the case (3.13). However, if the solutions of the differential equation (4) depend very sensitively on the initial values s j, the case (3.12) can be used the most efficiently (by more detailed approximations in the data file). For the integration of the differential equations (4) and (16) different one-step methods are implemented in the program package, such as (i) third, fourth, fifth-sixth and seventh-eighth order Runge-Kutta methods ([4],[6],[7]), (ii) fourth order Runge-Kutta-Fehlberg methods ([4]), (iii) fourth, fifth and sixth order ROW methods for stiff ODE s ([8],[9]). Each code is developed to work with automatic stepsize control (two algorithms are implemented for each code), and forward or backward direction can be chosed for the integrations. 10
11 Reference: [ 1 ] G.Varjú,K.Károlyi: Calculating screening effect of a metal cable sheath with nonlinearity, Int.Symp.EMC. Wroclaw,1990. ( ) [ 2 ] F.Jonas,G.Varjú: Gen.Model & Num.Method for Multicond.Systems Frequ.Dom., IEEE/KTH Pow.Tech.Conf.Stockholm,1995. ( ) [ 3 ] K.Károlyi: An interactive code to solve MBVP s, Int.Conference on Diff.Equations, Barcelona,1991. ( ) [ 4 ] J.Stoer,R.Bulirsch: Introduction to Numerical Analysis, Springer-Verlag [ 5 ] R.E.Bank,D.J.Rose: Glob.Appr.Newton Meth. Num.Math ( ) [ 6 ] P.J.Prince,J.R.Dormand: High order embedded Runge-Kutta formulae, J.Comp.Appl.Math.vol.7.no (67-75) [ 7 ] J.H.Verner: Explicit Runge-Kutta Methods with Estimates of the Local Error, Report 92.Univ.Auckland, New Zealand,1976. [ 8 ] P.Kaps,P.Rentrop: Generalized Runge-Kutta Meth.of Ord.Four with Steps Contr. for Stiff ODE s, Numer.Math.vol.33,1979 (55-68) [ 9 ] P.Kaps,G.Wanner: A Study of Rosenbrock-Type Methods of High Order, Numer.Math.vol.38,1981 ( ) 11
Initial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More 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 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 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 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 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 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 Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
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 2 Dianne P. O Leary c 2008 The Plan
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 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 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 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 informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More 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 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 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 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 informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationDelay Differential Equations with Constant Lags
Delay Differential Equations with Constant Lags L.F. Shampine Mathematics Department Southern Methodist University Dallas, TX 75275 shampine@smu.edu S. Thompson Department of Mathematics & Statistics Radford
More 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 informationPh 22.1 Return of the ODEs: higher-order methods
Ph 22.1 Return of the ODEs: higher-order methods -v20130111- Introduction This week we are going to build on the experience that you gathered in the Ph20, and program more advanced (and accurate!) solvers
More informationThree Body Problem using High-Order Runge-Kutta Interpolation
Three Body Problem using High-Order Runge-Kutta Interpolation Lawrence Mulholland Contents 1 NAG Library Mark 26 Reverse Communication Runge-Kutta Routines 1 1.1 Reverse Communication.....................
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
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 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 informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
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 information5. Ordinary Differential Equations. Indispensable for many technical applications!
Indispensable for many technical applications! Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 30 5.1. Introduction Differential Equations One of the most important fields of application of
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
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 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 informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
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 information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
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 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 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 informationMatrix Reduction Techniques for Ordinary Differential Equations in Chemical Systems
Matrix Reduction Techniques for Ordinary Differential Equations in Chemical Systems Varad Deshmukh University of California, Santa Barbara April 22, 2013 Contents 1 Introduction 3 2 Chemical Models 3 3
More informationEuler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.
Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be
More informationRicharson Extrapolation for Runge-Kutta Methods
Richarson Extrapolation for Runge-Kutta Methods Zahari Zlatevᵃ, Ivan Dimovᵇ and Krassimir Georgievᵇ ᵃ Department of Environmental Science, Aarhus University, Frederiksborgvej 399, P. O. 358, 4000 Roskilde,
More informationCONVERGENCE AND STABILITY CONSTANT OF THE THETA-METHOD
Conference Applications of Mathematics 2013 in honor of the 70th birthday of Karel Segeth. Jan Brandts, Sergey Korotov, Michal Křížek, Jakub Šístek, and Tomáš Vejchodský (Eds.), Institute of Mathematics
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
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 informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.
Logistics Week 12: Monday, Apr 18 HW 6 is due at 11:59 tonight. HW 7 is posted, and will be due in class on 4/25. The prelim is graded. An analysis and rubric are on CMS. Problem du jour For implicit methods
More informationExcel for Scientists and Engineers Numerical Method s. E. Joseph Billo
Excel for Scientists and Engineers Numerical Method s E. Joseph Billo Detailed Table of Contents Preface Acknowledgments About the Author Chapter 1 Introducing Visual Basic for Applications 1 Chapter
More informationEfficient path tracking methods
Efficient path tracking methods Daniel J. Bates Jonathan D. Hauenstein Andrew J. Sommese April 21, 2010 Abstract Path tracking is the fundamental computational tool in homotopy continuation and is therefore
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 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 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 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 informationNumerical Questions in ODE Boundary Value Problems
Numerical Questions in ODE Boundary Value Problems M.R.Osborne Contents 1 Introduction 2 2 ODE stability 6 2.1 Initial value problem stability.................. 6 2.2 Boundary value problem stability................
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 informationThe Chemical Kinetics Time Step a detailed lecture. Andrew Conley ACOM Division
The Chemical Kinetics Time Step a detailed lecture Andrew Conley ACOM Division Simulation Time Step Deep convection Shallow convection Stratiform tend (sedimentation, detrain, cloud fraction, microphysics)
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 informationVariable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems
Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems J. R. Cash and S. Girdlestone Department of Mathematics, Imperial College London South Kensington London SW7 2AZ,
More informationMATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction
MATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction Differential equations are a convenient way to express mathematically a change of a dependent variable (e.g. concentration
More informationSolving Models With Off-The-Shelf Software. Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings
An Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings Ric D. Herbert 1 Peter J. Stemp 2 1 Faculty of Science and Information Technology, The University of Newcastle,
More informationONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 5,53-65 (1999) ONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE ANDRZEJ MARCINIAK 1, BARBARA SZYSZKA 2 1 Institute of Computing Science,
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
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 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 informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationMAT300/500 Programming Project Spring 2019
MAT300/500 Programming Project Spring 2019 Please submit all project parts on the Moodle page for MAT300 or MAT500. Due dates are listed on the syllabus and the Moodle site. You should include all neccessary
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 informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Newton-Euler formulation) Dimitar Dimitrov Örebro University June 8, 2012 Main points covered Newton-Euler formulation forward dynamics inverse dynamics
More informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction 1 Introduction In this module, we develop solution techniques for numerically solving ordinary
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 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 information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationOpleiding Wiskunde & Informatica
Opleiding Wiskunde & Informatica Numerics and continuation for Reaction-Diffusion equations Renzo Baasdam s1524054 Supervisors: Martina Chirilus-Bruckner (MI) & Michael Emmerich (LIACS) BACHELOR THESIS
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 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 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 informationAn Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints
An Active Set Strategy for Solving Optimization Problems with up to 200,000,000 Nonlinear Constraints Klaus Schittkowski Department of Computer Science, University of Bayreuth 95440 Bayreuth, Germany e-mail:
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 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 informationINTRODUCTION TO COMPUTER METHODS FOR O.D.E.
INTRODUCTION TO COMPUTER METHODS FOR O.D.E. 0. Introduction. The goal of this handout is to introduce some of the ideas behind the basic computer algorithms to approximate solutions to differential equations.
More informationThe Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1
The Definition and Numerical Method of Final Value Problem and Arbitrary Value Problem Shixiong Wang 1*, Jianhua He 1, Chen Wang 2, Xitong Li 1 1 School of Electronics and Information, Northwestern Polytechnical
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter VII: Numerical Differentiation and Solution of Ordinary Differential Equations Greg Fasshauer Department of Applied Mathematics Illinois Institute
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationParallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
More informationDevelopment of a New One-Step Scheme for the Solution of Initial Value Problem (IVP) in Ordinary Differential Equations
International Journal of Theoretical and Applied Mathematics 2017; 3(2): 58-63 http://www.sciencepublishinggroup.com/j/ijtam doi: 10.11648/j.ijtam.20170302.12 Development of a New One-Step Scheme for the
More informationA Brief Introduction to Numerical Methods for Differential Equations
A Brief Introduction to Numerical Methods for Differential Equations January 10, 2011 This tutorial introduces some basic numerical computation techniques that are useful for the simulation and analysis
More 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 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 informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
Lab 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 informationF-TRANSFORM FOR NUMERICAL SOLUTION OF TWO-POINT BOUNDARY VALUE PROBLEM
Iranian Journal of Fuzzy Systems Vol. 14, No. 6, (2017) pp. 1-13 1 F-TRANSFORM FOR NUMERICAL SOLUTION OF TWO-POINT BOUNDARY VALUE PROBLEM I. PERFILIEVA, P. ŠTEVULIÁKOVÁ AND R. VALÁŠEK Abstract. We propose
More informationTHREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 6, 41-59 (2000) THREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE KAROL GAJDA 1, ANDRZEJ MARCINIAK 2, BARBARA SZYSZKA 1 1 Institute of
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 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 informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationFIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS III: Numerical and More Analytic Methods David Levermore Department of Mathematics University of Maryland
FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS III: Numerical and More Analytic Methods David Levermore Department of Mathematics University of Maryland 30 September 0 Because the presentation of this material
More informationRemark on the Sensitivity of Simulated Solutions of the Nonlinear Dynamical System to the Used Numerical Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 55, 2749-2754 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.59236 Remark on the Sensitivity of Simulated Solutions of
More informationSTATE VARIABLE (SV) SYSTEMS
Copyright F.L. Lewis 999 All rights reserved Updated:Tuesday, August 05, 008 STATE VARIABLE (SV) SYSTEMS A natural description for dynamical systems is the nonlinear state-space or state variable (SV)
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More information