Topic 6b Finite Difference Approximations

Size: px
Start display at page:

Download "Topic 6b Finite Difference Approximations"

Transcription

1 /8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions? Overview o deriving inite dierence pproximtions Algoritm or clculting inite dierence pproximtions Exmples Clculting inite dierence coeicients using MATLAB Exmples

2 /8/8 Wt re Finite Dierence Approximtions? Wt re Finite Dierence Approximtions? ( o ) Very oten in science nd engineering we must clculte derivtive. Wen we re processing dt rom mesurements or simultions, tere my not be n nlyticl eqution to work wit symboliclly. Typiclly, we only know te unction t discrete points. x x x Slide 4

3 /8/8 Wt re Finite Dierence Approximtions? ( o ) Suppose we wis to numericlly clculte te irst order derivtive t x. Te irst order derivtive is slope. We cn estimte te slope s riserun using inormtion rom surrounding points. rise run x x x run = x rise = x x x Slide 5 Wt re Finite Dierence Approximtions? ( o ) We cn estimte te derivtive t te midpoint between dt points. x x.5 x.5 x Te second order derivtive is te slope o te slope. x.5 x.5 x x x x x x.5.5 x x x Slide 6

4 /8/8 Generl Concept o Finite Dierence Approximtions ( o ) t x, y t x, y t x, y t x, y t x, y t x, y x i i or t xi, y i y t x, y Suppose we wis to estimte te unction (x) or one o its derivtives t loction (x i,y i ). x i or y i It is lwys possible to estimte tis rom liner sum o te unction vlues t surrounding points. Slide 7 Generl Concept o Finite Dierence Approximtions ( o ) Te trick is, ow do we clculte te coeicients n? Tese re unction o te positions o te points. x, y x y x, y,,,, 4, 5, 6, 7 x, y x, y x, y x, y Finite dierence coeicients depend only on te reltive position o te points. Tey do not depend on te bsolute positions. Sme coeicients Dierent coeicients Dierent coeicients Slide 8 4

5 /8/8 Types o Finite Dierences d x x Forwrd Finite Dierence Reces ed to use dt in te orwrd direction. d x.5 d x x x Centrl Finite Dierence Reces symmetriclly to use dt in bot directions or igest ccurcy. Bckwrd Finite Dierence Reces beind to use dt in te bckwrd direction. Slide 9 Continuum o Finite Dierence Approximtions Slide 5

6 /8/8 Two Key Considertions. Te position o te points rom wic te inite dierence pproximtion is clculted. More closely spced points improves ccurcy, but typiclly leds to more computtions.. Te loction o te point were te inite dierence is being evluted. We typiclly wnt to be s centered s possible or best ccurcy. Slide Overview o Deriving Finite Dierence Approximtions 6

7 /8/8 Concept o Using Polynomils We cn it n t order polynomil given + points. x x x x Ater we ve done te curve it, we cn interpolte te unction or ny o its derivtives rom tis polynomil x x x x x x x x x x x x x x x x 4x x Slide Fitting te Polynomil We write our polynomil t ec o our discrete points. Since our t order polynomil contins + coeicients n, we need t lest + discrete points to determine ll o tem. Tis lst point sows tt more points in te pproximtion will give better ccurcy in te inite dierence pproximtion. x xx x x x x x x x x x Slide 4 7

8 /8/8 Esiest Point or Evluting (x) Recll te equtions we will use to evlute (x) or one o its derivtives: 4 x xx x 4x x xxx 44x x x 6 x x x 4 Tese re most esily evluted t x = becuse te bove equtions reduce to 6 Slide 5 How to Mke Any Point Esy ow suppose we wis to evlute (x) or one o its derivtives t te generl point x = x d. To do tis, we sit our x xis by x d beore itting te polynomil. Recll tt te inite dierence coeicients depend only on te reltive position o te points. An oset will not ect teir vlues. ow we write our polynomil t ec sited point. x x x x x x x x x x x x x x x n n d In our sited coordinte system, our inite dierence is being evluted t x. 6 Slide 6 8

9 /8/8 Mtrix Form o Curve Fit Given our sited coordinte system, we write te polynomil t ec discrete point nd put tis lrge set o equtions in mtrix orm. 4 x xx x 4x x xxx 44x x x 6x4x x x x x x x x xx xx x x x x x x x x Slide 7 Algoritm or Clculting Finite Dierence Approximtions 8 9

10 /8/8 Step Coose x Coordintes Identiy te x coordintes o te points rom wic you wis to pproximte derivte. T x x x x 9 Step Sit x Axis Sit te unction cross te x xis until x = corresponds to te point were you wis to pproximte te derivtive. d x x d x x x x x x x x x x T

11 /8/8 Step Build [X] Mtrix Use te column vector x to build [X]. X x x x x x x x x x x x x x x x x x x x xx Insert s insted o x. Step 4 Invert X Clculte te inverse o X. y y y, y y y, Y X y, y, y,

12 /8/8 Step 5 Extrct Polynomil Coeicients X Y y y y y, y y y y, y y y y, y, y, y, y, y y y y, y y y y, y y y y, y y y y,,,, Step 6 Write Finite Dierence Approximtion x d x d x y y y y, y y y y, y y y y, y y y y,,,, 4

13 /8/8 Exmples 5 Exmple # Derive irst order nd second order inite dierence pproximtions tt spn cross tree points. Te pproximtions sould be evluted t te midpoint. x X x x x Y X x d d x d d xd 6

14 /8/8 Exmple # Derive irst order nd second order inite dierence pproximtions tt spn cross tree points. Te pproximtions sould be evluted t te irst point. x X x x x Y X.5.5 x d d x.5.5 d d xd 7 Exmple # Higer Order Accurcy ( o ) Let s evlute some derivtives t te midpoint o our discrete points Y X x X x x x x 8 4

15 /8/8 Exmple # Higer Order Accurcy ( o ) Te coeicients re ten Y x.5 6 d x x d x.5 4 x Clculting Finite Dierence Approximtions Using MATLAB 5

16 /8/8 Determining Finite Dierence Approximtions wit MATLAB ( o 6) We ve so r derived inite dierence pproximtions symboliclly. Wt i we wnt 6 t order ccurte inite dierences? Tis is unresonble to do symboliclly. Recll our mtrix eqution representing te polynomil written t ec discrete point. It lwys d te ollowing orm were te w s re just constnts. w w w w w w w, w, w, Slide Determining Finite Dierence Approximtions wit MATLAB ( o 6) ow, we re ble to seprte te w terms rom te terms. w w w w w w w, w, w, We were ble to put numbers to ll o tese coeicients. Tis is ully numericl mtrix. It does not contin ny symbolic vribles. Hint: W X wen X Tese re our symbolic vribles. So we build [W] by building nd pretending =. Slide 6

17 /8/8 Determining Finite Dierence Approximtions wit MATLAB ( o 6) Solving our mtrix eqution or [] gives w w w w w w w, w, w, V v v v v v v v v v, v, v, v, W Slide Determining Finite Dierence Approximtions wit MATLAB (4 o 6) ext we multiply our mtrices. v v v v v v v v v, v, v, v, v v v v v v v v v, v, v, v, Slide 4 7

18 /8/8 Determining Finite Dierence Approximtions wit MATLAB (5 o 6) ext, we red o te polynomil coeicients. v v v v v v v v v, v, v, v, vv v v v v v, v, v, Slide 5 Determining Finite Dierence Approximtions wit MATLAB (6 o 6) Lst, we write our inite dierence pproximtions rom te polynomil coeicients. vv v d v v v v v v d String t tese equtions long enoug, we relize tt te v ij coeicients cn be determined completely numericlly. We just ve to remember to divide by nd perps multiply te initedierence expression by constnt. Slide 6 8

19 /8/8 MATLAB Exmples 7 Exmple #4 6 t Order Accurte Finite Dierences ( o ). Here we need seven points to clculte seven polynomil coeicients. x T. To build te [W] mtrix, set = or now. x ˆ T W x ˆ xˆ xˆ xˆ xˆ xˆ xˆ Slide 8 9

20 /8/8 Exmple #4 6 t Order Accurte Finite Dierences ( o ). Invert [W]. VW Write te inite dierence pproximtions, remembering to incorporte te symbolic s bck in x x Slide 9

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),

More information

Math Week 5 concepts and homework, due Friday February 10

Math Week 5 concepts and homework, due Friday February 10 Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

Chapter 3 Solving Nonlinear Equations

Chapter 3 Solving Nonlinear Equations Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,

More information

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature COSC 336 Numericl Anlysis I Numericl Integrtion nd Dierentition III - Guss Qudrture nd Adptive Qudrture Edgr Griel Fll 5 COSC 336 Numericl Anlysis I Edgr Griel Summry o the lst lecture I For pproximting

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

ENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions

ENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions ENGI 44 Engineering Mthemtics Five Tutoril Exmples o Prtil Frctions 1. Express x in prtil rctions: x 4 x 4 x 4 b x x x x Both denomintors re liner non-repeted ctors. The cover-up rule my be used: 4 4 4

More information

Discrete Least-squares Approximations

Discrete Least-squares Approximations Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve

More information

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy. O SPECTRA ACCURACY OF QUADRATURE FORMUAE BASED O PIECEWISE POYOMIA ITERPOATIO A KURGAOV AD S TSYKOV Abstrct It is well-known tt te trpezoidl rule, wile being only second-order ccurte in generl, improves

More information

Math 20C Multivariable Calculus Lecture 5 1. Lines and planes. Equations of lines (Vector, parametric, and symmetric eqs.). Equations of lines

Math 20C Multivariable Calculus Lecture 5 1. Lines and planes. Equations of lines (Vector, parametric, and symmetric eqs.). Equations of lines Mt 2C Multivrible Clculus Lecture 5 1 Lines nd plnes Slide 1 Equtions of lines (Vector, prmetric, nd symmetric eqs.). Equtions of plnes. Distnce from point to plne. Equtions of lines Slide 2 Definition

More information

Orthogonal Polynomials and Least-Squares Approximations to Functions

Orthogonal Polynomials and Least-Squares Approximations to Functions Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

Numerical Integration Problems

Numerical Integration Problems Integrtion COS 33 Nuericl Integrtion Proles Bsic D nuericl integrtion Given ility to evlute or ny, ind Gol: est ccurcy wit ewest sples Clssic prole even nlytic unctions not necessrily integrle in closed

More information

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Best Approximation. Chapter The General Case

Best Approximation. Chapter The General Case Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given

More information

Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS

Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS S Cpter Numericl Integrtion lso clled qudrture Te gol of numericl integrtion is to pproximte numericlly. f(x)dx Tis is useful for difficult integrls like sin(x) ; sin(x ); x + x 4 Or worse still for multiple-dimensionl

More information

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but... Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Chapter 3 Polynomials

Chapter 3 Polynomials Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling

More information

Best Approximation in the 2-norm

Best Approximation in the 2-norm Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion

More information

Numerical Integration - (4.3)

Numerical Integration - (4.3) Numericl Itegrtio - (.). Te Degree of Accurcy of Qudrture Formul: Te degree of ccurcy of qudrture formul Qf is te lrgest positive iteger suc tt x k dx Qx k, k,,,...,. Exmple fxdx 9 f f,,. Fid te degree

More information

Chapter 6 Notes, Larson/Hostetler 3e

Chapter 6 Notes, Larson/Hostetler 3e Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn

More information

Quadratic Forms. Quadratic Forms

Quadratic Forms. Quadratic Forms Qudrtic Forms Recll the Simon & Blume excerpt from n erlier lecture which sid tht the min tsk of clculus is to pproximte nonliner functions with liner functions. It s ctully more ccurte to sy tht we pproximte

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

1 Linear Least Squares

1 Linear Least Squares Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving

More information

Section 6.1 Definite Integral

Section 6.1 Definite Integral Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined

More information

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2. Mth 43 Section 6. Section 6.: Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot

More information

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230 Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

Chapter Direct Method of Interpolation More Examples Civil Engineering

Chapter Direct Method of Interpolation More Examples Civil Engineering Chpter 5. Direct Method of Interpoltion More Exmples Civil Engineering Exmple o mximie ctch of bss in lke, it is suggested to throw the line to the depth of the thermocline. he chrcteristic feture of this

More information

LECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and

LECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and LECTURE 9 Numericl Integrtion Recll from Clculus I tht denite integrl is generlly thought of s representing the re under the grph of fèxè between the points x = nd x = b, even though this is ctully only

More information

University of Houston, Department of Mathematics Numerical Analysis II

University of Houston, Department of Mathematics Numerical Analysis II University of Houston, Deprtment of Mtemtics Numericl Anlysis II 6 Glerkin metod, finite differences nd colloction 6.1 Glerkin metod Consider sclr 2nd order ordinry differentil eqution in selfdjoint form

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ), 1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on

More information

Lecture 20: Numerical Integration III

Lecture 20: Numerical Integration III cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed

More information

Solutions to Assignment #8

Solutions to Assignment #8 Mth 1 Numericl Anlysis (Bueler) December 9, 29 Solutions to Assignment #8 Problems 64, exercise 14: The nswer turns out to be yes, which mens tht I hve to be orgnized in writing it up There re lot of fcts

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor 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 information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

CHAPTER 6b. NUMERICAL INTERPOLATION

CHAPTER 6b. NUMERICAL INTERPOLATION CHAPTER 6. NUMERICAL INTERPOLATION A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering y Dr. Irhim A. Asskk Spring ENCE - Computtion s in Civil Engineering II Deprtment o Civil

More information

On Second Derivative-Free Zero Finding Methods

On Second Derivative-Free Zero Finding Methods 010 Americn Control Conerence Mrriott Wterront, Bltimore, MD, USA June 30-July 0, 010 FrC07.4 On Second Derivtive-Free Zero Finding Methods Mohmmed A. Hsn Deprtment o Electricl & Computer Engineering University

More information

Midpoint Approximation

Midpoint Approximation Midpoint Approximtion Sometimes, we need to pproximte n integrl of the form R b f (x)dx nd we cnnot find n ntiderivtive in order to evlute the integrl. Also we my need to evlute R b f (x)dx where we do

More information

Derivatives, Finding Roots, Interpolating and Integration

Derivatives, Finding Roots, Interpolating and Integration Derivtives, Finding Roots, Interpolting nd Integrtion Doing Numericl Derivtives Altoug tking derivtives is very esy opertion for ll of you, tere re situtions were it must e done numericlly. For exmple,

More information

Chapter 2 Differentiation

Chapter 2 Differentiation Cpter Differentition. Introduction In its initil stges differentition is lrgely mtter of finding limiting vlues, wen te vribles ( δ ) pproces zero, nd to begin tis cpter few emples will be tken. Emple..:

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Lecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014

Lecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014 Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method

More information

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique? XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4

More information

1 Error Analysis of Simple Rules for Numerical Integration

1 Error Analysis of Simple Rules for Numerical Integration cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion

More information

1 The Riemann Integral

1 The Riemann Integral The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re

More information

Numerical Integration

Numerical Integration Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl

More information

Numerical Methods I Orthogonal Polynomials

Numerical Methods I Orthogonal Polynomials Numericl Methods I Orthogonl Polynomils Aleksndr Donev Cournt Institute, NYU 1 donev@cournt.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fll 2014 Nov 6th, 2014 A. Donev (Cournt Institute) Lecture IX

More information

Chapter 18 Two-Port Circuits

Chapter 18 Two-Port Circuits Cpter 8 Two-Port Circuits 8. Te Terminl Equtions 8. Te Two-Port Prmeters 8.3 Anlysis of te Terminted Two-Port Circuit 8.4 nterconnected Two-Port Circuits Motivtion Tévenin nd Norton equivlent circuits

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur Module Anlysis of Stticlly Indeterminte Structures by the Mtrix Force Method Version CE IIT, Khrgpur esson 8 The Force Method of Anlysis: Bems Version CE IIT, Khrgpur Instructionl Objectives After reding

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

CHAPTER 4a. ROOTS OF EQUATIONS

CHAPTER 4a. ROOTS OF EQUATIONS CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Sections 5.2: The Definite Integral

Sections 5.2: The Definite Integral Sections 5.2: The Definite Integrl In this section we shll formlize the ides from the lst section to functions in generl. We strt with forml definition.. The Definite Integrl Definition.. Suppose f(x)

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

1 Review: Volumes of Solids (Stewart )

1 Review: Volumes of Solids (Stewart ) Lecture : Some Bic Appliction of Te Integrl (Stewrt 6.,6.,.,.) ul Krin eview: Volume of Solid (Stewrt 6.-6.) ecll: we d provided two metod for determining te volume of olid of revolution. Te rt w by dic

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Lecture Outline. Dispersion Relation Electromagnetic Wave Polarization 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3c

Lecture Outline. Dispersion Relation Electromagnetic Wave Polarization 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3c Course Instructor Dr. Rymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mil: rcrumpf@utep.edu EE 4347 Applied Electromgnetics Topic 3c Wve Dispersion & Polriztion Wve Dispersion These notes & Polriztion

More information

MATH1901 Differential Calculus (Advanced)

MATH1901 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 information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial Derivatives. Limits. For a single variable function f (x), the limit lim Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Average Rate of Change (AROC) The average rate of change of y over an interval is equal to change in

Average Rate of Change (AROC) The average rate of change of y over an interval is equal to change in Averge Rte o Cnge AROC Te verge rte o cnge o y over n intervl is equl to b b y y cngein y cnge in. Emple: Find te verge rte o cnge o te unction wit rule 5 s cnges rom to 5. 4 4 6 5 4 0 0 5 5 5 5 & 4 5

More information

Lecture 4: Piecewise Cubic Interpolation

Lecture 4: Piecewise Cubic Interpolation Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion

More information

Section 14.3 Arc Length and Curvature

Section 14.3 Arc Length and Curvature Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in

More information

Chapters 4 & 5 Integrals & Applications

Chapters 4 & 5 Integrals & Applications Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

Construction of Gauss Quadrature Rules

Construction of Gauss Quadrature Rules Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture

More information

Chapter 14. Matrix Representations of Linear Transformations

Chapter 14. Matrix Representations of Linear Transformations Chpter 4 Mtrix Representtions of Liner Trnsformtions When considering the Het Stte Evolution, we found tht we could describe this process using multipliction by mtrix. This ws nice becuse computers cn

More information

If u = g(x) is a differentiable function whose range is an interval I and f is continuous on I, then f(g(x))g (x) dx = f(u) du

If u = g(x) is a differentiable function whose range is an interval I and f is continuous on I, then f(g(x))g (x) dx = f(u) du Integrtion by Substitution: The Fundmentl Theorem of Clculus demonstrted the importnce of being ble to find nti-derivtives. We now introduce some methods for finding ntiderivtives: If u = g(x) is differentible

More information

Advanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.

Advanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration. Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

III. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9

III. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9 III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Improper Integrals, and Differential Equations

Improper Integrals, and Differential Equations Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted

More information

Numerical quadrature based on interpolating functions: A MATLAB implementation

Numerical quadrature based on interpolating functions: A MATLAB implementation SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in

More information

( ) as a fraction. Determine location of the highest

( ) as a fraction. Determine location of the highest AB Clculus Exm Review Sheet - Solutions A. Preclculus Type prolems A1 A2 A3 A4 A5 A6 A7 This is wht you think of doing Find the zeros of f ( x). Set function equl to 0. Fctor or use qudrtic eqution if

More information

Disclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.

Disclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know. Disclimer: This is ment to help you strt studying. It is not necessrily complete list of everything you need to know. The MTH 33 finl exm minly consists of stndrd response questions where students must

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS

THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS CARLOS SUERO, MAURICIO ALMANZAR CONTENTS 1 Introduction 1 2 Proof of Gussin Qudrture 6 3 Iterted 2-Dimensionl Gussin Qudrture 20 4

More information

AB Calculus Review Sheet

AB Calculus Review Sheet AB Clculus Review Sheet Legend: A Preclculus, B Limits, C Differentil Clculus, D Applictions of Differentil Clculus, E Integrl Clculus, F Applictions of Integrl Clculus, G Prticle Motion nd Rtes This is

More information

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30 Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function

More information

Logarithms and Exponential Functions. Gerda de Vries & John S. Macnab. match as necessary, or to work these results into other lessons.

Logarithms and Exponential Functions. Gerda de Vries & John S. Macnab. match as necessary, or to work these results into other lessons. Logritms nd Eponentil Functions Gerd de Vries & Jon S. Mcn It is epected tt students re lred fmilir wit tis mteril. We include it ere for completeness. Te tree lessons given ere re ver sort. Te tecer is

More information

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x). AB Clculus Exm Review Sheet A. Preclculus Type prolems A1 Find the zeros of f ( x). This is wht you think of doing A2 A3 Find the intersection of f ( x) nd g( x). Show tht f ( x) is even. A4 Show tht f

More information

INTRODUCTION TO INTEGRATION

INTRODUCTION TO INTEGRATION INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide

More information

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of

More information

3.4 Numerical integration

3.4 Numerical integration 3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,

More information

UNIT #6 EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS

UNIT #6 EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS Answer Key Name: Date: UNIT # EXPONENTS, EXPONENTS, AND MORE EXPONENTS REVIEW QUESTIONS Part I Questions. Te epression 0 can be simpliied to () () 0 0. Wic o te ollowing is equivalent to () () 8 8? 8.

More information

NUMERICAL INTEGRATION

NUMERICAL INTEGRATION NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls

More information