Numerical Methods for Ordinary Differential Equations

Size: px
Start display at page:

Download "Numerical Methods for Ordinary Differential Equations"

Transcription

1 Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics

2 Ordiary Differetial Equatios d d d = F y, y ( t ), y ( t ),..., y ( t dt dt dt ) 0 Ordiary: oly oe idepedet variable Differetial: ukow fuctios eter ito the equatio through its derivatives Order: highest derivative i F Degree: expoet of the highest derivative 3 d y t y t = Example: ( ) ( ) 0 dt

3 What is Solutio of ODE? y = y( t) A problem ivolvig ODE is ot completely specified by its equatio ODE has to be supplemeted with boudary coditios: y Iitial value problem: is give at some startig value, i ad it is desired to fid y at some fial poits tf or at some discrete list of poits (for example, at tabulated itervals). Two poit bodary value problem: Boudary coditios are specified at more tha oe t ; typically some of the coditios will be specified at t i ad some at t. f t

4 What is Numerical Solutio to the Iitial Value Problem? dy( t) dt (, ( )); ( ) 0 0 = f t y t y t = y A umerical solutio to this problem geerates sequece of values for the idepedet variable t, t,, t ad a correspodig sequece of values of the depedet variable y y, y,, y so that each approximates solutio at : y( t ) y, = 0,, t

5 Euler Metod All fiite differece methods start from the same coceptual idea: Add small icremets to your fuctio correspodig to derivatives (right-had side of the equatios) multiplied by the stepsize. Euler method is a implemetatio of this idea i the simplest ad most direct form. y t y E u ler y tru e t Sigle-Step Forward Propagatio

6 Euler Algorithm for First-Order ODE dy f ( t, y) dt = y = f ( t, y) t iitialize t, y y( t ) do while i y = y + f ( t, y) t i + i i i+ ed do i t = t + t

7 Step Size Effects i Radioactive Decay dn N Aalytics: U = U N = N ( t = 0) e U U dt τ Numerics(Euler): dn U NU ( t) = NU (0) + t + O ( t) dt Ni Ni + Ni t τ ( ) Numerical solutio will deped o the step size t t τ

8 Stability of Euler Algorithm Step size if ofte limited by the stability criterio: dy at = ay y(0) =, y = e dt After Euler steps of size t: y = y ay t y = ( a t) + Approximate solutio will decay mootoically oly if is small eough: t t For a sigle decayig expoetial-like solutio (i.e. if there is oly oe max first order equatio) the existece of a stability criterio is ot a problem because has to be small for the reasos of accuracy. t a t

9 Accuracy: Discretizatio ad Roudoff Errors Local: C p+ du Number of steps for roudoff error to be comparable with the discretizatio error: N L = f ( u, t) Lε dt LE = y+ u+ ( t+ ) N u( t) = y t N f = f ( t) y( t) = f ( τ ) dτ h f ( t ) t0 = 0 Global: t+ LE = h f ( t ) f ( τ ) dτ GE = y y( t ) Method is of order iff: LE = O( h ) LE Ch h = t t t Itegrate over iterval: Lε = + h p L t f t0 Full Error: Ch N N GE = h f ( t ) f ( τ ) dτ = 0 GE = N = 0 t t LE t 0

10 Global Discretizatio Error Example Suppose we wat to fid the solutio over the iterval Divide the iterval ito equal steps so that t = T at T y( T) = e, y = a 3 ( at ) ( at) y( T) = at + +! 3! 3 ( ) ( at ) ( )( ) ( at) y = at + + 3! 3! ( at) 3 ( at) a t y( T) y = + + O ate! 3! at [ 0,T ] This is a measure of the global trucatio error, i.e., the error over a fixed rage i t. It is proportioal to the first power of the step size ad hece the Euler method is a first order method (do ot cofuse this with the fact that we are applyig it i this case to a first order equatio).

11 Reducig Higher Order ODE to System of First Order ODE Solve higher order ODEs by splittig them ito sets of first order equatios: d y dy + p( t) + q( t) y = g ( t) dt dt dz = g ( t) p( t) z q( t) y dy dt z = dt dy = z dt There is o uique way to do this: z dz dp( t) g ( t) q( t) y dy = + dt dt = + p( t) y dt dy = z p( t) y dt

12 Example: Realistic Motio of Baseball v v + ω r v ω r iitialize t, y( t ) do while i y = y + f ( t, y ) t i+ i i i t = t + t i+ ed do i F d r a g ω + d r v m = mg B v + S0v ω dt v x = x + v t x i i i x x B x vi+ = vi vvi t m y = y + v t y i+ i i v = v g t y i+ y i z = z + v t z i + i i z z S 0v xω vi+ = vi t m

13 More Realistic Modelig of Air Flow

14 ODE of Liear Harmoic Oscillator d θ g + siθ = 0 dt l for small θ siθ θ d θ g g + θ = 0, Ω = dt l l dθ Etotal = ml + mglθ dt must be coserved!

15 Euler Method for Liear Harmoic Oscillator Switch to dimesioless quatities: d θ + 0 θ = θ = θ 0 si( Ω t + φ ) dt E total dθ = + θ dt Euler discretizatio algorithm: ω+ = ω θ t E = + θ + = θ + ω t t E t t total = E + t + = + ( ω ) + θ + total ( )

16 Euler Fails o θ ( t)

17 Euler Fails o ω( t)

18 Euler Fails o Phase Space Trajectory

19 Ca We Resurrect Euler by Usig Smaller Step Size?

20 Cromer Fixed Euler Method for LHO = t ω ω θ + ω ω θ θ ω = + t = t + t Apparetly trivial trick, but: E + = E + t + O t θ = θ si( t t ), ω = θ cos( t t ) ( ) ( 3 ω θ ) t t 0 ω θ = θ cos( ) = 0 over a period t

21 From Euler to Higher Order Algorithms y = y + f ( t, y ) + t = t + h ( ) ( ) + vs. Mea value theorem exact y t + t = y t + dy dt t t m

22 Midpoit Method: Secod Order Ruge-Kutta s = f ( t, y ) h h s = f ( t +, y + s ) y = y + hs + O( h ) + t = t + h + 3

23 Classical Ruge-Kutta Fourth-order method s = f ( t, y ) h h s = f ( t +, y + s ) h h s3 = f ( t +, y + s ) s = f ( t + h, y + hs ) 4 3 h y = y + ( s + s + s + s ) + O( h 5 ) t = t + h +

24 Classical Ruge-Kutta F90 Subroutie

25 Geeral Sigle-Step Methods Each of the k stages of the algorithm computes slope by evaluatig (, ) for a particular value of ad a value of obtaied by takig liear combiatios of the previous slopes: f t y t y i s = f ( t + α h, y + h β s ), i =,, k i i i, j j j= The proposed step is also a liear combiatio of the slopes: k y = y + h γ s + i i i = Error is estimated from yet aother liear combiatio of the slopes: e = h δ s + i i i = k The parameters are determied by matchig terms i the Taylor series expasio of the slopes the order of the method is the expoet of the smallest power of h that caot be matched. I MATLAB ODE umerical routies are amed as odexx, where idicates the order ad xx is some special feature of the method. s i

26 Example: MATLAB ode3 Fuctio (Bogacki ad Shampie BS3 Algorithm) s = f ( t, y ) h y = y + (s + 3s + 4 s ) h h s = f ( t +, y + s ) 3 3 s3 = f ( t + h, y + hs ) 4 4 h e = ( 5s + 6s + + 8s 9 s ) t = t + h; s = f ( t, y ) 7 3 4

27 Beyod Ruge-Kutta Methods Ruge-Kutta methods propagates a solutio over a iterval by combiig the iformatio from several Euler-style steps (each ivolvig oe evaluatio of the right-had side f s), ad the usig the iformatio obtaied to match Taylor series expasio up to some higher order. Richardso extrapolatio method used the powerful idea of extrapolatig computed result to the value that would have bee obtaied if the stepsize had bee very much smaller tha it actually was. I particular, extrapolatio to zero stepsize is the desired goal implemeted by Burlich-Stoer algorithm. Predictor-corrector methods store the solutio alog the way, ad use those results to extrapolate the solutio oe step advaced; they correct the extrapolatio usig derivative iformatio at the ew poit.

28 Stiff Systems of Differetial Equatios Stiffess arises i systems of ODE where there are two or more very differet scales of the idepedet variable: du dv = 998u v, = 999u 999v u = y z u = e e dt dt v = y + t z u(0) =, v(0) = 0 v = e + e Follow the variatio i the solutio o the shortest legth scale to maitai stability of the itegratio eve though accuracy requiremets allow for a much larger step size use implicit methods: explicit y = cy, c > 0 y = y + ty = ( c t ) y + t > c y as implicit = + = = + c t y cy y y ty y y 000 t t 000t

29 Solutios to Stiffess Beyod Implicit Euler To improve higher-order (tha Euler, which is first-order) methods use: Geeralizatios of Ruge-Kutta methods Rosebrock methods ad Kaps-Retrop methods. Burlich-Stoer algorithm geeralized to Bader-Deuflhard semi-implicit extrapolatio method. Predictor-corrector methods geeralized to Gear backward differetiatio method.

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

Chapter 9: Numerical Differentiation

Chapter 9: Numerical Differentiation 178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at

More information

Math 113, Calculus II Winter 2007 Final Exam Solutions

Math 113, Calculus II Winter 2007 Final Exam Solutions Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

Chapter 2: Numerical Methods

Chapter 2: Numerical Methods Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,

More information

ENGI 9420 Engineering Analysis Assignment 3 Solutions

ENGI 9420 Engineering Analysis Assignment 3 Solutions ENGI 9 Egieerig Aalysis Assigmet Solutios Fall [Series solutio of ODEs, matri algebra; umerical methods; Chapters, ad ]. Fid a power series solutio about =, as far as the term i 7, to the ordiary differetial

More information

Fall 2013 MTH431/531 Real analysis Section Notes

Fall 2013 MTH431/531 Real analysis Section Notes Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters

More information

An efficient time integration method for extra-large eddy simulations

An efficient time integration method for extra-large eddy simulations A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler Departmet of Mathematics Master s Thesis A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler

More information

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said

More information

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations Geeraliig the DTFT The Trasform M. J. Roberts - All Rights Reserved. Edited by Dr. Robert Akl 1 The forward DTFT is defied by X e jω = x e jω i which = Ω is discrete-time radia frequecy, a real variable.

More information

ENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01

ENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01 ENGI 940 Lecture Notes 3 - Numerical Methods Page 3.0 3. Numerical Methods The majority of equatios of iterest i actual practice do ot admit ay aalytic solutio. Eve equatios as simple as = e ad I = e d

More information

ChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS

ChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS I a exothermic reactio the temperature will cotiue to rise as oe moves alog a plug flow reactor util all of the limitig reactat is exhausted. Schematically

More information

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig

More information

Mathematical Methods for Physics and Engineering

Mathematical Methods for Physics and Engineering Mathematical Methods for Physics ad Egieerig Lecture otes Sergei V. Shabaov Departmet of Mathematics, Uiversity of Florida, Gaiesville, FL 326 USA CHAPTER The theory of covergece. Numerical sequeces..

More information

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y

More information

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:

More information

ENGI Series Page 6-01

ENGI Series Page 6-01 ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,

More information

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014 UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)

(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b) Chapter 0 Review 597. E; a ( + )( + ) + + S S + S + + + + + + S lim + l. D; a diverges by the Itegral l k Test sice d lim [(l ) ], so k l ( ) does ot coverge absolutely. But it coverges by the Alteratig

More information

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim Math 3, Sectio 2. (25 poits) Why we defie f(x) dx as we do. (a) Show that the improper itegral diverges. Hece the improper itegral x 2 + x 2 + b also diverges. Solutio: We compute x 2 + = lim b x 2 + =

More information

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5! aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (

More information

Section A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics

Section A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.

More information

Generalized Semi- Markov Processes (GSMP)

Generalized Semi- Markov Processes (GSMP) Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual

More information

Taylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH

Taylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )

More information

Five Steps Block Predictor-Block Corrector Method for the Solution of ( )

Five Steps Block Predictor-Block Corrector Method for the Solution of ( ) Applied Mathematics, 4,, -66 Published Olie May 4 i SciRes. http://www.scirp.org/oural/am http://dx.doi.org/.46/am.4.87 Five Steps Block Predictor-Block Corrector y = f x, y, y Method for the Solutio of

More information

Reading. Required: Witkin, Particle System Dynamics, SIGGRAPH 97 course notes on Physically Based Modeling.

Reading. Required: Witkin, Particle System Dynamics, SIGGRAPH 97 course notes on Physically Based Modeling. Particle Systems Readig Required: Witki, Particle System Dyamics, SIGGRAPH 97 course otes o Physically Based Modelig. Optioal Witki ad Baraff, Differetial Equatio Basics, SIGGRAPH 97 course otes o Physically

More information

CS321. Numerical Analysis and Computing

CS321. Numerical Analysis and Computing CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca

More information

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4

More information

Calculus with Analytic Geometry 2

Calculus with Analytic Geometry 2 Calculus with Aalytic Geometry Fial Eam Study Guide ad Sample Problems Solutios The date for the fial eam is December, 7, 4-6:3p.m. BU Note. The fial eam will cosist of eercises, ad some theoretical questios,

More information

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet

More information

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.

More information

Topic 5 [434 marks] (i) Find the range of values of n for which. (ii) Write down the value of x dx in terms of n, when it does exist.

Topic 5 [434 marks] (i) Find the range of values of n for which. (ii) Write down the value of x dx in terms of n, when it does exist. Topic 5 [44 marks] 1a (i) Fid the rage of values of for which eists 1 Write dow the value of i terms of 1, whe it does eist Fid the solutio to the differetial equatio 1b give that y = 1 whe = π (cos si

More information

Math 257: Finite difference methods

Math 257: Finite difference methods Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...

More information

Numerical Methods for Partial Differential Equations

Numerical Methods for Partial Differential Equations Numerical Methods for Partial Differetial Equatios CAAM 45 Sprig 005 Lecture 4 -step time-steppig methods: stability, accuracy Ruge-Kutta Methods, Istructor: Tim Warburto Today Recall AB stability regios

More information

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m? MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle

More information

SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)

SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE) QUESTION BANK 8 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayaavaam Road 5758 QUESTION BANK (DESCRIPTIVE Subject with Code : (6HS6 Course & Brach: B.Tech AG Year & Sem: II-B.Tech& I-Sem

More information

MATH 10550, EXAM 3 SOLUTIONS

MATH 10550, EXAM 3 SOLUTIONS MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,

More information

MATH2007* Partial Answers to Review Exercises Fall 2004

MATH2007* Partial Answers to Review Exercises Fall 2004 MATH27* Partial Aswers to Review Eercises Fall 24 Evaluate each of the followig itegrals:. Let u cos. The du si ad Hece si ( cos 2 )(si ) (u 2 ) du. si u 2 cos 7 u 7 du Please fiish this. 2. We use itegratio

More information

In algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:

In algebra one spends much time finding common denominators and thus simplifying rational expressions. For example: 74 The Method of Partial Fractios I algebra oe speds much time fidig commo deomiators ad thus simplifyig ratioal epressios For eample: + + + 6 5 + = + = = + + + + + ( )( ) 5 It may the seem odd to be watig

More information

CS537. Numerical Analysis and Computing

CS537. Numerical Analysis and Computing CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be

More information

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES Diola Bagayoko (0 PARTAL DFFERENTAL EQUATONS SEPARATON OF ARABLES. troductio As discussed i previous lectures, partial differetial equatios arise whe the depedet variale, i.e., the fuctio, varies with

More information

1 Generating functions for balls in boxes

1 Generating functions for balls in boxes Math 566 Fall 05 Some otes o geeratig fuctios Give a sequece a 0, a, a,..., a,..., a geeratig fuctio some way of represetig the sequece as a fuctio. There are may ways to do this, with the most commo ways

More information

Assignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1

Assignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1 Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

CHAPTER 1 SEQUENCES AND INFINITE SERIES

CHAPTER 1 SEQUENCES AND INFINITE SERIES CHAPTER SEQUENCES AND INFINITE SERIES SEQUENCES AND INFINITE SERIES (0 meetigs) Sequeces ad limit of a sequece Mootoic ad bouded sequece Ifiite series of costat terms Ifiite series of positive terms Alteratig

More information

Time-Domain Representations of LTI Systems

Time-Domain Representations of LTI Systems 2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable

More information

Assignment 5: Solutions

Assignment 5: Solutions McGill Uiversity Departmet of Mathematics ad Statistics MATH 54 Aalysis, Fall 05 Assigmet 5: Solutios. Let y be a ubouded sequece of positive umbers satisfyig y + > y for all N. Let x be aother sequece

More information

Chapter 4. Fourier Series

Chapter 4. Fourier Series Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,

More information

Z ß cos x + si x R du We start with the substitutio u = si(x), so du = cos(x). The itegral becomes but +u we should chage the limits to go with the ew

Z ß cos x + si x R du We start with the substitutio u = si(x), so du = cos(x). The itegral becomes but +u we should chage the limits to go with the ew Problem ( poits) Evaluate the itegrals Z p x 9 x We ca draw a right triagle labeled this way x p x 9 From this we ca read off x = sec, so = sec ta, ad p x 9 = R ta. Puttig those pieces ito the itegralrwe

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Reading. Particle Systems. What are particle systems? Overview

Reading. Particle Systems. What are particle systems? Overview Readig Particle Systems Required: Witki, Particle System Dyamics, SIGGRAPH 97 course otes o Physically Based Modelig. Optioal Witki ad Baraff, Differetial Equatio Basics, SIGGRAPH 97 course otes o Physically

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

FINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,

More information

Lecture 19: Convergence

Lecture 19: Convergence Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may

More information

sin(n) + 2 cos(2n) n 3/2 3 sin(n) 2cos(2n) n 3/2 a n =

sin(n) + 2 cos(2n) n 3/2 3 sin(n) 2cos(2n) n 3/2 a n = 60. Ratio ad root tests 60.1. Absolutely coverget series. Defiitio 13. (Absolute covergece) A series a is called absolutely coverget if the series of absolute values a is coverget. The absolute covergece

More information

Numerical Methods in Fourier Series Applications

Numerical Methods in Fourier Series Applications Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets

More information

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations Chapter The Solutio of Numerical Algebraic ad Trascedetal Equatios Itroductio I this chapter we shall discuss some umerical methods for solvig algebraic ad trascedetal equatios. The equatio f( is said

More information

Math 312 Lecture Notes One Dimensional Maps

Math 312 Lecture Notes One Dimensional Maps Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,

More information

Lecture 2: Finite Difference Methods in Heat Transfer

Lecture 2: Finite Difference Methods in Heat Transfer Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (

More information

ME NUMERICAL METHODS Fall 2007

ME NUMERICAL METHODS Fall 2007 ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course

More information

MATH 31B: MIDTERM 2 REVIEW

MATH 31B: MIDTERM 2 REVIEW MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate x (x ) (x 3).. Partial Fractios Solutio: The umerator has degree less tha the deomiator, so we ca use partial fractios. Write x (x ) (x 3) = A x + A (x ) +

More information

Seunghee Ye Ma 8: Week 5 Oct 28

Seunghee Ye Ma 8: Week 5 Oct 28 Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value

More information

Question 1: The magnetic case

Question 1: The magnetic case September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013 Large Deviatios for i.i.d. Radom Variables Cotet. Cheroff boud usig expoetial momet geeratig fuctios. Properties of a momet

More information

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0, 4. Series Solutios of Differetial Equatios:Special Fuctios 4.. Illustrative examples.. 5. Obtai the geeral solutio of each of the followig differetial equatios i terms of Maclauri series: d y (a dx = xy,

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

Phys 6303 Final Exam Solutions December 19, 2012

Phys 6303 Final Exam Solutions December 19, 2012 Phys 633 Fial Exam s December 19, 212 You may NOT use ay book or otes other tha supplied with this test. You will have 3 hours to fiish. DO YOUR OWN WORK. Express your aswers clearly ad cocisely so that

More information

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5 CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio

More information

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0. THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of

More information

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics: Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets

More information

Chapter 6 Infinite Series

Chapter 6 Infinite Series Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat

More information

Honors Calculus Homework 13 Solutions, due 12/8/5

Honors Calculus Homework 13 Solutions, due 12/8/5 Hoors Calculus Homework Solutios, due /8/5 Questio Let a regio R i the plae be bouded by the curves y = 5 ad = 5y y. Sketch the regio R. The two curves meet where both equatios hold at oce, so where: y

More information

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0, Math Activity 9( Due with Fial Eam) Usig first ad secod Taylor polyomials with remaider, show that for, 8 Usig a secod Taylor polyomial with remaider, fid the best costat C so that for, C 9 The th Derivative

More information

Review Problems for the Final

Review Problems for the Final Review Problems for the Fial Math - 3 7 These problems are provided to help you study The presece of a problem o this hadout does ot imply that there will be a similar problem o the test Ad the absece

More information

Solutions to quizzes Math Spring 2007

Solutions to quizzes Math Spring 2007 to quizzes Math 4- Sprig 7 Name: Sectio:. Quiz a) x + x dx b) l x dx a) x + dx x x / + x / dx (/3)x 3/ + x / + c. b) Set u l x, dv dx. The du /x ad v x. By Itegratio by Parts, x(/x)dx x l x x + c. l x

More information

Chapter 7: Numerical Series

Chapter 7: Numerical Series Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals

More information

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

A 2nTH ORDER LINEAR DIFFERENCE EQUATION A 2TH ORDER LINEAR DIFFERENCE EQUATION Doug Aderso Departmet of Mathematics ad Computer Sciece, Cocordia College Moorhead, MN 56562, USA ABSTRACT: We give a formulatio of geeralized zeros ad (, )-discojugacy

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will

More information

Practice Problems: Taylor and Maclaurin Series

Practice Problems: Taylor and Maclaurin Series Practice Problems: Taylor ad Maclauri Series Aswers. a) Start by takig derivatives util a patter develops that lets you to write a geeral formula for the -th derivative. Do t simplify as you go, because

More information

Løsningsførslag i 4M

Løsningsførslag i 4M Norges tekisk aturviteskapelige uiversitet Istitutt for matematiske fag Side 1 av 6 Løsigsførslag i 4M Oppgave 1 a) A sketch of the graph of the give f o the iterval [ 3, 3) is as follows: The Fourier

More information

BESSEL EQUATION and BESSEL FUNCTIONS

BESSEL EQUATION and BESSEL FUNCTIONS BESSEL EQUATION ad BESSEL FUNCTIONS Bessel s Equatio Summary of Bessel Fuctios d y dy y d + d + =. If is a iteger, the two idepedet solutios of Bessel s Equatio are J J, Bessel fuctio of the first kid,

More information

Math 12 Final Exam, May 11, 2011 ANSWER KEY. 2sinh(2x) = lim. 1 x. lim e. x ln. = e. (x+1)(1) x(1) (x+1) 2. (2secθ) 5 2sec2 θ dθ.

Math 12 Final Exam, May 11, 2011 ANSWER KEY. 2sinh(2x) = lim. 1 x. lim e. x ln. = e. (x+1)(1) x(1) (x+1) 2. (2secθ) 5 2sec2 θ dθ. Math Fial Exam, May, ANSWER KEY. [5 Poits] Evaluate each of the followig its. Please justify your aswers. Be clear if the it equals a value, + or, or Does Not Exist. coshx) a) L H x x+l x) sihx) x x L

More information

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1. SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad

More information

The Numerical Solution of Singular Fredholm Integral Equations of the Second Kind

The Numerical Solution of Singular Fredholm Integral Equations of the Second Kind WDS' Proceedigs of Cotributed Papers, Part I, 57 64, 2. ISBN 978-8-7378-39-2 MATFYZPRESS The Numerical Solutio of Sigular Fredholm Itegral Equatios of the Secod Kid J. Rak Charles Uiversity, Faculty of

More information

Chapter 10: Power Series

Chapter 10: Power Series Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because

More information

Sequences and Series of Functions

Sequences and Series of Functions Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges

More information

MH1101 AY1617 Sem 2. Question 1. NOT TESTED THIS TIME

MH1101 AY1617 Sem 2. Question 1. NOT TESTED THIS TIME MH AY67 Sem Questio. NOT TESTED THIS TIME ( marks Let R be the regio bouded by the curve y 4x x 3 ad the x axis i the first quadrat (see figure below. Usig the cylidrical shell method, fid the volume of

More information

Part I: Covers Sequence through Series Comparison Tests

Part I: Covers Sequence through Series Comparison Tests Part I: Covers Sequece through Series Compariso Tests. Give a example of each of the followig: (a) A geometric sequece: (b) A alteratig sequece: (c) A sequece that is bouded, but ot coverget: (d) A sequece

More information

Discrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 22

Discrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 22 CS 70 Discrete Mathematics for CS Sprig 2007 Luca Trevisa Lecture 22 Aother Importat Distributio The Geometric Distributio Questio: A biased coi with Heads probability p is tossed repeatedly util the first

More information