Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
|
|
- Bennett Lawson
- 6 years ago
- Views:
Transcription
1 Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t))
2 First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous on D def = {(t, y) a t b, c y d.} Let (t 0, y 0 ), (t, y) D, t = t t 0, y = y y 0. Ten f (t, y) = P 1 (t, y) + R 1 (t, y), were
3 First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous on D def = {(t, y) a t b, c y d.} Let (t 0, y 0 ), (t, y) D, t = t t 0, y = y y 0. Ten f (t, y) = P 1 (t, y) + R 1 (t, y), were P 1 (t, y) = f (t 0, y 0 ) + t f t (t 0, y 0 ) + y f y (t 0, y 0 ), R 1 (t, y) = 2 t 2 for some point (ξ, µ) D. 2 f t 2 (ξ, µ) + t y 2 f t y (ξ, µ) + 2 y 2 2 f (ξ, µ), y 2
4 Second order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous on D def = {(t, y) a t b, c y d.} Let (t 0, y 0 ), (t, y) D, t = t t 0, y = y y 0. Ten f (t, y) = P 2 (t, y) + R 2 (t, y), were
5 Second order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous on D def = {(t, y) a t b, c y d.} Let (t 0, y 0 ), (t, y) D, t = t t 0, y = y y 0. Ten ( P 2 (t, y) = f (t 0, y 0 ) + + f (t, y) = P 2 (t, y) + R 2 (t, y), were ) t f t (t 0, y 0 ) + y f y (t 0, y 0 ) ( 2 t 2 f 2 t 2 (t 0, y 0 ) + t y 2 f t y (t 0, y 0 ) + 2 y 2 f 2 y 2 (t 0, y 0 ) ),
6 Second order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous on D def = {(t, y) a t b, c y d.} Let (t 0, y 0 ), (t, y) D, t = t t 0, y = y y 0. Ten ( P 2 (t, y) = f (t 0, y 0 ) + + ( 2 t 2 R 2 (t, y) = 1 3 3! f (t, y) = P 2 (t, y) + R 2 (t, y), j=0 were ) t f t (t 0, y 0 ) + y f y (t 0, y 0 ) 2 f t 2 (t 0, y 0 ) + t y ( 3 j ) 3 j t j y 2 f t y (t 0, y 0 ) + 2 y 2 3 f t 3 j y (ξ, µ) j ) 2 f y 2 (t 0, y 0 ),
7 Runge-Kutta Metod of Order Two (III) Midpoint Metod w 0 = α, w j+1 = w j + f ( t j + 2, w j + ) 2 f (t j, w j ), j = 0, 1,, N 1. Two function evaluations for eac j, Second order accuracy. No need for derivative calculations
8 General 2nd order Runge-Kutta Metods w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + (a 1 f (t j, w j ) + a 2 f (t j + α 2, w j + δ 2 f (t j, w j ))).
9 General 2nd order Runge-Kutta Metods w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + (a 1 f (t j, w j ) + a 2 f (t j + α 2, w j + δ 2 f (t j, w j ))). Two function evaluations for eac j, Want to coose a 1, a 2, α 2, δ 2 for igest possible order of accuracy.
10 General 2nd order Runge-Kutta Metods w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + (a 1 f (t j, w j ) + a 2 f (t j + α 2, w j + δ 2 f (t j, w j ))). Two function evaluations for eac j, Want to coose a 1, a 2, α 2, δ 2 for igest possible order of accuracy. local truncation error τ j+1 () = y(t j+1) y(t j ) (a 1 f (t j, y(t j ))+a 2 f (t j +α 2, y(t j )+δ 2 f (t j, y(t j )))) = y (t j ) + 2 y (t j ) + O( 2 ) ( f (a 1 + a 2 ) f (t j, y(t j )) + a 2 α 2 t (t j, y(t j )) + a 2 δ 2 f (t j, y(t j )) f ) y (t j, y(t j )) + O( 2 ).
11 local truncation error τ j+1 () = y (t j ) + 2 y (t j ) + O( 2 ) ( f (a 1 + a 2 ) f (t j, y(t j )) + a 2 α 2 t (t j, y(t j )) + a 2 δ 2 f (t j, y(t j )) f ) y (t j, y(t j )) + O( 2 ). were y (t j ) = f (t j, y(t j )),
12 local truncation error τ j+1 () = y (t j ) + 2 y (t j ) + O( 2 ) ( f (a 1 + a 2 ) f (t j, y(t j )) + a 2 α 2 t (t j, y(t j )) + a 2 δ 2 f (t j, y(t j )) f ) y (t j, y(t j )) + O( 2 ). were y (t j ) = f (t j, y(t j )), y (t j ) = d f d t (t j, y(t j )) = f t (t j, y(t j )) + f (t j, y(t j )) f y (t j, y(t j )). For any coice wit we ave a second order metod a 1 + a 2 = 1, a 2 α 2 = a 2 δ 2 = 2, τ j+1 () = O( 2 ).
13 local truncation error τ j+1 () = y (t j ) + 2 y (t j ) + O( 2 ) ( f (a 1 + a 2 ) f (t j, y(t j )) + a 2 α 2 t (t j, y(t j )) + a 2 δ 2 f (t j, y(t j )) f ) y (t j, y(t j )) + O( 2 ). were y (t j ) = f (t j, y(t j )), y (t j ) = d f d t (t j, y(t j )) = f t (t j, y(t j )) + f (t j, y(t j )) f y (t j, y(t j )). For any coice wit we ave a second order metod a 1 + a 2 = 1, a 2 α 2 = a 2 δ 2 = 2, τ j+1 () = O( 2 ). Four parameters, tree equations
14 General 2nd order Runge-Kutta Metods w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + (a 1 f (t j, w j ) + a 2 f (t j + α 2, w j + δ 2 f (t j, w j ))). a 1 + a 2 = 1, a 2 α 2 = a 2 δ 2 = 2, Midpoint metod: a 1 = 0, a 2 = 1, α 2 = δ 2 = 2, w j+1 = w j + f ( t j + 2, w j + ) 2 f (t j, w j ). Modified Euler metod: a 1 = a 2 = 1 2, α 2 = δ 2 =, w j+1 = w j + 2 (f (t j, w j ) + f (t j+1, w j + f (t j, w j ))).
15 3rd order Runge-Kutta Metod (rarely used in practice) w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + ( ( f (t j, w j ) + 3f t j , w j f (t j + 3, w j + )) 3 f (t j, w j )) def = w j + φ(t j, w j ).
16 3rd order Runge-Kutta Metod (rarely used in practice) w 0 = α; for j = 0, 1,, N 1, w j+1 = w j + ( ( f (t j, w j ) + 3f t j , w j f (t j + 3, w j + )) 3 f (t j, w j )) def = w j + φ(t j, w j ). local truncation error τ j+1 () = y(t j+1)y(t j ) φ(t j, y(t j )) = O( 3 ).
17 4t order Runge-Kutta Metod w 0 = α; for j = 0, 1,, N 1, k 1 = f (t j, w j ), ( k 2 = f t j + 2, w j + 1 ) 2 k 1, ( k 3 = f t j + 2, w j + 1 ) 2 k 2, k 4 = f (t j+1, w j + k 3 ), w j+1 = w j (k 1 + 2k 2 + 2k 3 + k 4 ). 4 function evaluations per step
18 Example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
19 Example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
20 Example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
21 Example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
22 Adaptive Error Control (I) dy = f (t, y(t)), a t b, y(a) = α. dt Consider a variable-step metod wit a well-cosen function φ(t, w, ): w 0 = α. for j = 0, 1,, coose step-size j = t j+1 t j, set wj+1 = w j + j φ (t j, w j, j ).
23 Adaptive Error Control (I) dy = f (t, y(t)), a t b, y(a) = α. dt Consider a variable-step metod wit a well-cosen function φ(t, w, ): w 0 = α. for j = 0, 1,, coose step-size j = t j+1 t j, set wj+1 = w j + j φ (t j, w j, j ). Adaptively coose step-size to satisfy given tolerance
24 Adaptive Error Control (II) Given an order-n metod w0 = α. for j = 0, 1,, w j+1 = w j + φ (t j, w j, ), = t j+1 t j, local truncation error (LTE) τ j+1 () = y(t j+1) y(t j ) φ (t j, y(t j ), ) = O ( n ). Given tolerance τ > 0, we would like to estimate largest step-size for wic τ j+1 () τ.
25 Adaptive Error Control (II) Given an order-n metod w0 = α. for j = 0, 1,, w j+1 = w j + φ (t j, w j, ), = t j+1 t j, local truncation error (LTE) τ j+1 () = y(t j+1) y(t j ) φ (t j, y(t j ), ) = O ( n ). Given tolerance τ > 0, we would like to estimate largest step-size for wic τ j+1 () τ. Approac: Estimate τ j+1 () wit order-(n + 1) metod w j+1 = w j + φ (t j, w j, ), for j 0.
26 LTE Estimation (I) Assume w j y(t j ), w j y(t j ) (only estimating LTE).
27 LTE Estimation (I) Assume w j y(t j ), w j y(t j ) (only estimating LTE). φ(t, w, ) is order-n metod τ j+1 () def = w j y(t j ) = y(t j+1 ) y(t j ) φ (t j, y(t j ), ) y(t j+1 ) (w j + φ (t j, w j, )) y(t j+1 ) w j+1 = O ( n ).
28 LTE Estimation (I) Assume w j y(t j ), w j y(t j ) (only estimating LTE). φ(t, w, ) is order-n metod τ j+1 () def = w j y(t j ) = φ(t, w, ) is order-(n + 1) metod, y(t j+1 ) y(t j ) φ (t j, y(t j ), ) y(t j+1 ) (w j + φ (t j, w j, )) y(t j+1 ) w j+1 = O ( n ). τ j+1 () def = w j y(t j ) = y(t j+1 ) y(t j ) φ (t j, y(t j ), ) ( y(t j+1 ) w j + φ ) (t j, w j, ) y(t j+1 ) w j+1 = O ( n+1).
29 LTE Estimation (III) Assume w j y(t j ), w j y(t j ) (only estimating LTE).
30 LTE Estimation (III) Assume w j y(t j ), w j y(t j ) (only estimating LTE). φ(t, w, ) is order-n metod τ j+1 () y(t j+1) w j+1 = O ( n ).
31 LTE Estimation (III) Assume w j y(t j ), w j y(t j ) (only estimating LTE). φ(t, w, ) is order-n metod τ j+1 () y(t j+1) w j+1 φ(t, w, ) is order-(n + 1) metod, = O ( n ). τ j+1 () y(t j+1) w j+1 = O ( n+1).
32 LTE Estimation (III) Assume w j y(t j ), w j y(t j ) (only estimating LTE). φ(t, w, ) is order-n metod τ j+1 () y(t j+1) w j+1 φ(t, w, ) is order-(n + 1) metod, = O ( n ). τ j+1 () y(t j+1) w j+1 = O ( n+1). terefore τ j+1 () y(t j+1) w j+1 + w j+1 w j+1 = O ( n+1) + w j+1 w j+1 LTE estimate: τ j+1 () w j+1 w j+1 = O ( n ).
33 step-size selection (I) LTE estimate: τ j+1 () w j+1 w j+1
34 step-size selection (I) LTE estimate: τ j+1 () w j+1 w j+1 Since τ j+1 () = O ( n ), assume τ j+1 () K n were K is independent of.
35 step-size selection (I) LTE estimate: τ j+1 () w j+1 w j+1 Since τ j+1 () = O ( n ), assume τ j+1 () K n were K is independent of. K sould satisfy K n w j+1 w j+1. Assume LTE for new step-size q satisfies given tolerance ɛ: τ j+1 (q ) ɛ, need to estimate q.
36 step-size selection (II) LTE estimate: K n τ j+1 () w j+1 w j+1
37 step-size selection (II) LTE estimate: K n τ j+1 () w j+1 w j+1 Assume LTE for q satisfies given tolerance ɛ: τ j+1 (q ) K (q ) n = q n K n q n w j+1 w j+1 ɛ.
38 step-size selection (II) LTE estimate: K n τ j+1 () w j+1 w j+1 Assume LTE for q satisfies given tolerance ɛ: τ j+1 (q ) K (q ) n = q n K n q n w j+1 w j+1 ɛ. 1 new step-size estimate: q ɛ n w j+1 w j+1
39 Summary Given order-n metod w j+1 = w j + φ (t j, w j, ), = t j+1 t j, j 0, and given order-(n + 1) metod w j+1 = w j + φ (t j, w j, ), for j 0. for eac j, compute w j+1 = w j + φ (t j, w j, ), w j+1 = w j + φ (t j, w j, ), new step-size q sould satisfy 1 1 q ɛ n ɛ n w j+1 w j+1 =. φ (t j, w j, ) φ (t j, w j, )
40 Runge-Kutta-Felberg: 4 t order metod, 5 t order estimate w j+1 = w j k k k k 5, w j+1 = w j k k k k k 6, were
41 Runge-Kutta-Felberg: 4 t order metod, 5 t order estimate w j+1 = w j k k k k 5, w j+1 = w j k k k k k 6, were k 1 = f (t j, w j ),
42 Runge-Kutta-Felberg: 4 t order metod, 5 t order estimate w j+1 = w j k k k k 5, w j+1 = w j k k k k k 6, were k 1 = f (t j, w j ), ( k 2 = f t j + 4, w j + 1 ) 4 k 1, ( k 3 = f t j + 3 8, w j k ) 32 k 2, ( k 4 = f t j , w j k k ) 2197 k 3, ( k 5 = f t j +, w j k 1 8 k k ) 4104 k 4, ( k 6 = f t j + 2, w j 8 27 k k k k 4 11 ) 40 k 5.
43 Runge-Kutta-Felberg: step-size selection procedure compute a conservative value for q: 1 q = ɛ 4 2 ( w j+1 w j+1 ).
44 Runge-Kutta-Felberg: step-size selection procedure compute a conservative value for q: 1 q = ɛ 4 2 ( w j+1 w j+1 ). make restricted step-size cange: { 0.1, if q 0.1, = 4, if q 4.
45 Runge-Kutta-Felberg: step-size selection procedure compute a conservative value for q: 1 q = ɛ 4 2 ( w j+1 w j+1 ). make restricted step-size cange: { 0.1, if q 0.1, = 4, if q 4. step-size can t be too big: = min (, max ).
46 Runge-Kutta-Felberg: step-size selection procedure compute a conservative value for q: 1 q = ɛ 4 2 ( w j+1 w j+1 ). make restricted step-size cange: { 0.1, if q 0.1, = 4, if q 4. step-size can t be too big: = min (, max ). step-size can t be too small: if < min ten declare failure.
47 Runge-Kutta-Felberg: example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
48 Runge-Kutta-Felberg: example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
49 Runge-Kutta-Felberg: example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t.
50 Runge-Kutta-Felberg: example dy Initial Value ODE = y t 2 + 1, 0 t 2, y(0) = 0.5, dt exact solution y(t) = (1 + t) e t. t j j y(t j ) w j w j
51 Runge-Kutta-Felberg: solution plots Initial Value ODE dy dt = y t 2 + 1, 0 t 2, y(0) = 0.5.
52 Runge-Kutta-Felberg: solution plots Initial Value ODE dy dt = y t 2 + 1, 0 t 2, y(0) = 0.5.
53 Runge-Kutta-Felberg: solution errors Initial Value ODE dy dt = y t 2 + 1, 0 t 2, y(0) = 0.5.
54 Runge-Kutta-Felberg: solution errors Initial Value ODE dy dt = y t 2 + 1, 0 t 2, y(0) = 0.5.
55 Runge-Kutta-Felberg: truncation errors Initial Value ODE actual def = dy dt y(t j ) w j j = y t 2 + 1, 0 t 2, y(0) = 0.5., def estimate = w j w j j.
56 Runge-Kutta-Felberg: truncation errors Initial Value ODE actual def = dy dt y(t j ) w j j = y t 2 + 1, 0 t 2, y(0) = 0.5., def estimate = w j w j j.
AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
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 informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationCHAPTER 2 MODELING OF THREE-TANK SYSTEM
7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More information(5.5) Multistep Methods
(5.5) Mulstep Metods Consider te inial-value problem for te ordinary differenal equaon: y t f t, y, a t b, y a. Let y t be te unique soluon. In Secons 5., 5. and 5.4, one-step numerical metods: Euler Metod,
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationChapter XI. Solution of Ordinary Differential Equations
dy d 7 Capter XI Solution of Ordinary Differential Equations Several approaces are described for te solution of ODE (ordinary differential equations). Tese include: Te Taylor series Metod Te Euler Metod
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationChapter 1. Density Estimation
Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationApplications of the van Trees inequality to non-parametric estimation.
Brno-06, Lecture 2, 16.05.06 D/Stat/Brno-06/2.tex www.mast.queensu.ca/ blevit/ Applications of te van Trees inequality to non-parametric estimation. Regular non-parametric problems. As an example of suc
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationqwertyuiopasdfghjklzxcvbnmqwerty uiopasdfghjklzxcvbnmqwertyuiopasd fghjklzxcvbnmqwertyuiopasdfghjklzx cvbnmqwertyuiopasdfghjklzxcvbnmq
qwertyuiopasdfgjklzxcbnmqwerty uiopasdfgjklzxcbnmqwertyuiopasd fgjklzxcbnmqwertyuiopasdfgjklzx cbnmqwertyuiopasdfgjklzxcbnmq Projectile Motion Quick concepts regarding Projectile Motion wertyuiopasdfgjklzxcbnmqwertyui
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationFunction Composition and Chain Rules
Function Composition an Cain Rules James K. Peterson Department of Biological Sciences an Department of Matematical Sciences Clemson University November 2, 2018 Outline Function Composition an Continuity
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More information5.1 introduction problem : Given a function f(x), find a polynomial approximation p n (x).
capter 5 : polynomial approximation and interpolation 5 introduction problem : Given a function f(x), find a polynomial approximation p n (x) Z b Z application : f(x)dx b p n(x)dx, a a one solution : Te
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
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 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 informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationLesson 6: The Derivative
Lesson 6: Te Derivative Def. A difference quotient for a function as te form f(x + ) f(x) (x + ) x f(x + x) f(x) (x + x) x f(a + ) f(a) (a + ) a Notice tat a difference quotient always as te form of cange
More informationMA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE
MA64 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS B Dr.T.Kulandaivel Department of Applied Matematics SVCE Numerical ordinar differential equations is te part
More information3.1. COMPLEX DERIVATIVES 89
3.1. COMPLEX DERIVATIVES 89 Teorem 3.1.4 (Cain rule) Let f : D D C and g : D C be olomorpic functions on te domains D and D respectively. Suppose tat bot f and g are C 1 -smoot. 10 Ten (g f) (z) g (f(z))
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationMATH1901 Differential Calculus (Advanced)
MATH1901 Dierential Calculus (Advanced) Capter 3: Functions Deinitions : A B A and B are sets assigns to eac element in A eactl one element in B A is te domain o te unction B is te codomain o te unction
More information1 Introduction to Optimization
Unconstrained Convex Optimization 2 1 Introduction to Optimization Given a general optimization problem of te form min x f(x) (1.1) were f : R n R. Sometimes te problem as constraints (we are only interested
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationChapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1
Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationAn Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method
Int. J. Contemp. Mat. Sciences, Vol. 2, 27, no. 2, 983-989 An Approximation to te Solution of te Brusselator System by Adomian Decomposition Metod and Comparing te Results wit Runge-Kutta Metod J. Biazar
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More information1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -
Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.
More informationA.P. CALCULUS (AB) Outline Chapter 3 (Derivatives)
A.P. CALCULUS (AB) Outline Capter 3 (Derivatives) NAME Date Previously in Capter 2 we determined te slope of a tangent line to a curve at a point as te limit of te slopes of secant lines using tat point
More informationVariable Step Size Differential Equation Solvers
Math55: Differential Equations 1/30 Variable Step Size Differential Equation Solvers Jason Brewer and George Little Introduction The purpose of developing numerical methods is to approximate the solution
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationWalrasian Equilibrium in an exchange economy
Microeconomic Teory -1- Walrasian equilibrium Walrasian Equilibrium in an ecange economy 1. Homotetic preferences 2 2. Walrasian equilibrium in an ecange economy 11 3. Te market value of attributes 18
More informationMATH 155A FALL 13 PRACTICE MIDTERM 1 SOLUTIONS. needs to be non-zero, thus x 1. Also 1 +
MATH 55A FALL 3 PRACTICE MIDTERM SOLUTIONS Question Find te domain of te following functions (a) f(x) = x3 5 x +x 6 (b) g(x) = x+ + x+ (c) f(x) = 5 x + x 0 (a) We need x + x 6 = (x + 3)(x ) 0 Hence Dom(f)
More informationApplied Linear Statistical Models. Simultaneous Inference Topics. Simultaneous Estimation of β 0 and β 1 Issues. Simultaneous Inference. Dr.
Applied Linear Statistical Models Simultaneous Inference Dr. DH Jones Simultaneous Inference Topics Simultaneous estimation of β 0 and β 1 Bonferroni Metod Simultaneous estimation of several mean responses
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More information2.3 Product and Quotient Rules
.3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More informationNotes: DERIVATIVES. Velocity and Other Rates of Change
Notes: DERIVATIVES Velocity and Oter Rates of Cange I. Average Rate of Cange A.) Def.- Te average rate of cange of f(x) on te interval [a, b] is f( b) f( a) b a secant ( ) ( ) m troug a, f ( a ) and b,
More informationOnline Learning: Bandit Setting
Online Learning: Bandit Setting Daniel asabi Summer 04 Last Update: October 0, 06 Introduction [TODO Bandits. Stocastic setting Suppose tere exists unknown distributions ν,..., ν, suc tat te loss at eac
More informationNew Distribution Theory for the Estimation of Structural Break Point in Mean
New Distribution Teory for te Estimation of Structural Break Point in Mean Liang Jiang Singapore Management University Xiaou Wang Te Cinese University of Hong Kong Jun Yu Singapore Management University
More informationApproximation of the Viability Kernel
Approximation of te Viability Kernel Patrick Saint-Pierre CEREMADE, Université Paris-Daupine Place du Marécal de Lattre de Tassigny 75775 Paris cedex 16 26 october 1990 Abstract We study recursive inclusions
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 informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More information