chap09_8up - Ordinary Differential Equations Numerical...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Scientic 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 Copyright c 2002. Reproduction permitted for noncommercial, educational use only. Michael T. Heath Scientic Computing 1 / 84 Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Outline 1 Ordinary Differential Equations 2 Numerical Solution of ODEs 3 Additional Numerical Methods Michael T. Heath Scientic Computing 2 / 84 Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Differential Equations Initial Value Problems Stability Differential Equations Differential equations involve derivatives of unknown solution function Ordinary differential equation (ODE): all derivatives are with respect to single independent variable, often representing time Solution of differential equation is function in innite-dimensional space of functions Numerical solution of differential equations is based on nite-dimensional approximation Differential equation is replaced by algebraic equation whose solution approximates that of given differential equation Michael T. Heath Scientic Computing 3 / 84 Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Differential Equations Initial Value Problems Stability Order of ODE Order of ODE is determined by highest-order derivative of solution function appearing in ODE ODE with higher-order derivatives can be transformed into equivalent rst-order system We will discuss numerical solution methods only for rst-order ODEs Most ODE software is designed to solve only rst-order equations Michael T. Heath Scientic Computing 4 / 84 Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Differential Equations Initial Value Problems Stability Higher-Order ODEs, continued For k-th order ODE y ( k ) ( t ) = f ( t,y,y ,...,y ( k-1) ) dene k new unknown functions u 1 ( t ) = y ( t ) , u 2 ( t ) = y ( t ) , ..., u k ( t ) = y ( k-1) ( t ) Then original ODE is equivalent to rst-order system u 1 ( t ) u 2 ( t ) . . . u k-1 ( t ) u k ( t ) = u 2 ( t ) u 3 ( t ) . . . u k ( t ) f ( t,u 1 ,u 2 ,...,u k ) Michael T. Heath Scientic Computing 5 / 84 Ordinary Differential Equations Numerical Solution of ODEs Additional Numerical Methods Differential Equations Initial Value Problems Stability Example: Newtons Second Law Newtons Second Law of Motion, F = ma , is second-order ODE, since acceleration a is second derivative of position coordinate, which we denote by y Thus, ODE has form y 00 = F/m where F and m are force and mass, respectively Dening u 1 = y and u 2 = y...
View Full Document

Page1 / 11

chap09_8up - Ordinary Differential Equations Numerical...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online