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

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

2 Tody: Objectives identify the most widely used qudrture method is it chep? is it effective? how does it compre to Newton-Cotes (Trpezoid, Simpson, etc)? Mteril Section 5.4 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 / 9

3 Qudrture up until now, our qudrture methods were of the form f (x) dx n w j f (x j ) j=0 where x j re eqully spced nodes Trpezoid: Simpson: f (x) dx b 6 f (x) dx b f () + b f (b) (b ) f () + f ( + b Similr for higher order polynomil Newton-Cotes rules ) + b 6 f (b) Dvid Semerro (NCSA) CS 57 Mrch 0, 014 / 9

4 Qudrture with Freedom! These qudrture rules hve one thing in common: they re restrictive e.g. Simpson: f (x) dx b 6 (b ) f () + f ( + b ) + b 6 f (b) Trpezoid, Simpson, etc (Newton-Cotes) re bsed on eqully spced nodes We know one thing lredy from interpoltion: eqully spced nodes result in wiggle. Wht other choice do we hve? (...recll how we fixed wiggle in interpoltion: by moving the loction of the nodes) Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

5 Gussin Qudrture free ourselves from eqully spced nodes combine selection of the nodes nd selection of the weights into one qudrture rule Gussin Qudrture Choose the nodes nd coefficients optimlly to mximize the degree of precision of the qudrture rule: f (x) dx n w j f (x j ) j=0 Gol Seek w j nd x j so tht the qudrture rule is exct for relly high polynomils Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

6 Gussin Qudrture f (x) dx n w j f (x j ) j=0 we hve n + 1 points x j [, b], x 0 < x 1 < < x n < x n b. we hve n + 1 rel coefficients w j so there re n + totl unknowns to tke cre of there were only unknowns in the cse of trpezoid ( weights) there were only unknowns in the cse of Simpson ( weights) there were only n + 1 unknowns in the cse of generl Newton-Cotes (n + 1 weights) Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

7 Gussin Qudrture f (x) dx n w j f (x j ) j=0 we hve n + 1 points x j [, b], x 0 < x 1 < < x n < x n b. we hve n + 1 rel coefficients w j so there re n + totl unknowns to tke cre of there were only unknowns in the cse of trpezoid ( weights) there were only unknowns in the cse of Simpson ( weights) there were only n + 1 unknowns in the cse of generl Newton-Cotes (n + 1 weights) n + unknowns (using n + 1 nodes) cn be used to exctly interpolte nd integrte polynomils of degree up to n + 1 Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

8 Better Nodes Exmple The first thing we do is SIMPLIFY consider the cse of n = 1 (-point) consider [, b] = [, 1] for simplicity we know how the trpezoid rule works Question: cn we possibly do better using only function evlutions? Gol: Find w 0, w 1, x 0, x 1 so tht is s ccurte s possible... f (x) dx w 0 f (x 0 ) + w 1 f (x 1 ) Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

9 Grphicl View Consider x + 1 dx = Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

10 Derive... Agin, we re considering [, b] = [, 1] for simplicity: f (x) dx w 0 f (x 0 ) + w 1 f (x 1 ) Gol: find w 0, w 1, x 0, x 1 so tht the pproximtion is exct up to cubics. So try ny cubic: This implies tht: f (x) = c 0 + c 1 x + c x + c x f (x) dx = ( c0 + c 1 x + c x + c x ) dx ( = w 0 c0 + c 1 x 0 + c x 0 + c x0) + ( ) w 1 c0 + c 1 x 1 + c x 1 + c x 1 Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

11 Derive... f (x) dx = ( c0 + c 1 x + c x + c x ) dx ( = w 0 c0 + c 1 x 0 + c x 0 + c x0) + ( ) w 1 c0 + c 1 x 1 + c x 1 + c x 1 Rerrnge into constnt, liner, qudrtic, nd cubic terms: ) 1 c 0 (w 0 + w 1 dx c (w 0 x 0 + w 1 x 1 x dx ) ) 1 + c 1 (w 0 x 0 + w 1 x 1 x dx + ( ) 1 + c w 0 x 0 + w 1 x 1 x dx = 0 Since c 0, c 1, c nd c re rbitrry, then their coefficients must ll be zero. Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

12 Derive... This implies: w 0 + w 1 = w 0 x 0 + w 1 x 1 = Some lgebr leds to: dx = w 0 x 0 + w 1 x 1 = x dx = x dx = 0 w 0 x 0 + w 1 x 1 = x dx = 0 Therefore: w 0 = 1 w 1 = 1 x 0 = x 1 = ( ) ( ) f (x) dx f + f Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

13 Over nother intervl? ( ) ( ) f (x) dx f + f f (x) dx? integrting over [, b] insted of [, 1] needs trnsformtion: chnge of vribles wnt t = c 1 x + c 0 with t = t x = nd t = 1 t x = b let t = b x b+ b (verify) let x = b t + b+ then dx = b dt Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

14 Over nother intervl? f (x) dx? let x = b t + b+ then dx = b dt f (x) dx = now use the qudrture formul over [, 1] ( ) (b )t + b + b f dt note: using two points, n = 1, gve us exct integrtion for polynomils of degree less *1+1 = nd less. Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

15 Over nother intervl? Previous exmple... 1 x + 1 dx = f (x)dx = 1 [ 1 f ( ) t + f dt ( ) + + f ( )] where x = +t Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

16 Over nother intervl? Evluting... 1 [ x + 1 dx 1 f [ ( 1 f 1 x + 1 dx ( ) ) + f + f ( ( 9 6 )] )] 1 [f ( ) + f (1.115)] 1 [ ] = where f (x) = x + 1 Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

17 Extending Guss Qudrture we need more to mke this work for more thn two points A sensible qudrture rule for the intervl [, 1] bsed on 1 node would use the node x = 0. This is root of φ(x) = x Notice: ± 1 re the roots of φ(x) = x 1 generl φ(x)? Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

18 Guss Qudrture Theorem Krl Friedrich Guss proved the following result: Let q(x) be nontrivil polynomil of degree n + 1 such tht x k q(x)dx = 0 (0 k n) nd let x 0, x 1,..., x n be the zeros of q(x). Then f (x)dx n A i f (x i ), A i = i=0 l i (x)dx will be exct for ll polynomils of degree t most n + 1. Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

19 Sketch of Proof Let f (x) be polynomil of degree n + 1. Then we cn write f (x) = p(x)q(x) + r(x), where p(x) nd r(x) re of degree t most n (This is bsiclly dividing f by q with reminder r). Then by the hypothesis, p(x)q(x)dx = 0. Further, f (x i ) = p(x i )q(x i ) + r(x i ) = r(x i ). Thus, f (x)dx = r(x)dx n f (x i ) i=0 l i (x)dx But this is exct becuse r(x) is (t most) degree n polynomil. Thus, we need to find the polynomils q(x). Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

20 Orthogonl Polynomils Orthogonlity of Functions Two functions g(x) nd h(x) re orthogonl on [, b] if g(x)h(x) dx = 0 so the nodes we re using re roots of orthogonl polynomils these re the Legendre Polynomils Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

21 Legendre Polynomils given on the exm φ 0 = 1 φ 1 = x φ = x 1 φ = 5x x. In generl: φ n (x) = n 1 n xφ n (x) n 1 n φ n (x) Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

22 Notes on Legendre Roots The Legendre Polynomils re orthogonl (nice!) The Legendre Polynomils increse in polynomils order (like monomils) The Legendre Polynomils don t suffer from poor conditioning (unlike monomils) The Legendre Polynomils don t hve closed form expression (recursion reltion is needed) The roots of the Legendre Polynomils re the nodes for Gussin Qudrture (GL nodes) Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

23 Qudrture Nodes (see) Often listed in tbles Weights determined by extension of bove Roots re symmetric in [, 1] Exmple: 1 if(n==0) x = 0; w = ; if(n==1) 4 x(1) = -1/sqrt(); x() = -x(1); 5 w(1) = 1; w() = w(1); 6 if(n==) 7 x(1) = -sqrt(/5); x() = 0; x() = -x(1); 8 w(1) = 5/9; w() = 8/9; w() = w(1); 9 if(n==) 10 x(1) = ; x(4) = -x(1); 11 x() = ; x() = -x(); 1 w(1) = ; w(4) = w(1); 1 w() = ; w() = w(); 14 if(n==4) 15 x(1) = ; x(5) = -x(1); 16 x() = ; x(4) = -x(); 17 x() = 0; 18 w(1) = ; w(5) = w(1); 19 w() = ; w(4) = w(); 0 w() = ; 1 if(n==5) x(1) = ; x(6) = -x(1); x() = ; x(5) = -x(); 4 x() = ; x(4) = -x(); 5 w(1) = ; w(6) = w(1); 6 w() = ; w(5) = w(); 7 w() = ; w(4) = w(); Dvid Semerro (NCSA) CS 57 Mrch 0, 014 / 9

24 View of Nodes Order Order Order Order Order Dvid Semerro (NCSA) CS 57 Mrch 0, 014 / 9

25 Theory The connection between the roots of the Legendre polynomils nd exct integrtion of polynomils is estblished by the following theorem. Theorem Suppose tht x 0, x 1,..., x n re roots of the nth Legendre polynomil P n (x) nd tht for ech i = 0, 1,..., n the numbers w i re defined by w i = n j = 0 j i x x 1 j dx = l i (x) dx x i x j Then f (x)dx = n w i f (x i ), where f (x) is ny polynomil of degree less or equl to n + 1. i=0 Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

26 Do not!!!! When evluting qudrture rule f (x)dx = n w i f (x i ). i=0 do not generte the nodes nd weights ech time. Use lookup tble... Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

27 Exmple.5 Approximte x ln x dx using Gussin qudrture with n = 1. 1 ( 1 ) ( ) Solution As derived erlier we wnt to use f (x) dx f + f From erlier we know tht we re interested in.5 ( ) (1.5 1)t + ( ) f (x) dx = f dt 1 Therefore, we re looking for the integrl of 1 1 ( ) x + 5 f dx = 1 ( ) ( ) x + 5 x + 5 ln dx Using Gussin qudrture, our numericl integrtion becomes: ( ) ( ) ( ) ( ) ln + ln = Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

28 Exmple Approximte 0 x e x dx using Gussin qudrture with n = 1. Solution We gin wnt to convert our limits of integrtion to -1 to 1. Using the sme process s the erlier exmple, we get: 0 x e x dx = 1 ( ) t + 1 e (t+1)/ dt. Using the Gussin roots we get: ( 1 x e x dx ) ( e ( +1)/ ) e ( +1)/ = Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

29 Numericl Question How does n point Guss qudrture compre with n point Newton-Cotes... Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

30 Exmples with Python... Exmple int guss.py: bse routine for Guss qudrture Exmple int guss test.py: integrte 5 0 xe x dx with 1 1 subintervl, incresing number of nodes nodes, increses number of intervls (pnels) Result: fewer totl evlutions in GL qudrture with 1 subintervl nd mny nodes versus nodes nd mny subpnels. Also more ccurte. Dvid Semerro (NCSA) CS 57 Mrch 0, / 9

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

### 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

### 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

### Interpolation. Gaussian Quadrature. September 25, 2011

Gussin Qudrture September 25, 2011 Approximtion of integrls Approximtion of integrls by qudrture Mny definite integrls cnnot be computed in closed form, nd must be pproximted numericlly. Bsic building

### 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

### 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

### 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

### 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

### 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

### 20 MATHEMATICS POLYNOMIALS

0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of

### Lecture Notes: Orthogonal Polynomials, Gaussian Quadrature, and Integral Equations

18330 Lecture Notes: Orthogonl Polynomils, Gussin Qudrture, nd Integrl Equtions Homer Reid My 1, 2014 In the previous set of notes we rrived t the definition of Chebyshev polynomils T n (x) vi the following

### 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

### 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

### Continuous Random Variables

STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

### Math 8 Winter 2015 Applications of Integration

Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

### 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

### The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

Mth 104: finl informtion The finl exm will tke plce on Fridy My 11th from 8m 11m in Evns room 60. The exm will cover ll prts of the course with equl weighting. It will cover Chpters 1 5, 7 15, 17 21, 23

### 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 =

### Introduction to Numerical Analysis

Introduction to Numericl Anlysis Doron Levy Deprtment of Mthemtics nd Center for Scientific Computtion nd Mthemticl Modeling (CSCAMM) University of Mrylnd June 14, 2012 D. Levy CONTENTS Contents 1 Introduction

### Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions

### 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

### MTH 122 Fall 2008 Essex County College Division of Mathematics Handout Version 10 1 October 14, 2008

MTH 22 Fll 28 Essex County College Division of Mthemtics Hndout Version October 4, 28 Arc Length Everyone should be fmilir with the distnce formul tht ws introduced in elementry lgebr. It is bsic formul

### We know that if f is a continuous nonnegative function on the interval [a, b], then b

1 Ares Between Curves c 22 Donld Kreider nd Dwight Lhr We know tht if f is continuous nonnegtive function on the intervl [, b], then f(x) dx is the re under the grph of f nd bove the intervl. We re going

### Line and Surface Integrals: An Intuitive Understanding

Line nd Surfce Integrls: An Intuitive Understnding Joseph Breen Introduction Multivrible clculus is ll bout bstrcting the ides of differentition nd integrtion from the fmilir single vrible cse to tht of

### Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

IOSR Journl of Mthemtics (IOSR-JM) e-issn: - p-issn: 9-X. Volume Issue Ver. IV (My. - Jun. 0) PP - www.iosrjournls.org Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture T. Rmchndrn D.Udykumr nd

### Numerical Methods I. Olof Widlund Transcribed by Ian Tobasco

Numericl Methods I Olof Widlund Trnscribed by In Tobsco Abstrct. This is prt one of two semester course on numericl methods. The course ws offered in Fll 011 t the Cournt Institute for Mthemticl Sciences,

### 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

### (0.0)(0.1)+(0.3)(0.1)+(0.6)(0.1)+ +(2.7)(0.1) = 1.35

7 Integrtion º½ ÌÛÓ Ü ÑÔÐ Up to now we hve been concerned with extrcting informtion bout how function chnges from the function itself. Given knowledge bout n object s position, for exmple, we wnt to know

### 7. Numerical evaluation of definite integrals

7. Numericl evlution of definite integrls Tento učení text yl podpořen z Operčního progrmu Prh - Adptilit Hn Hldíková Numericl pproximtion of definite integrl is clled numericl qudrture, the formuls re

### Math 113 Exam 2 Practice

Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number

### Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

### Math 360: A primitive integral and elementary functions

Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:

### The Wave Equation I. MA 436 Kurt Bryan

1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

### 7.6 The Use of Definite Integrals in Physics and Engineering

Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 7.6 The Use of Definite Integrls in Physics nd Engineering It hs been shown how clculus cn be pplied to find solutions to geometric problems

### Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

Improper Integrls The First Fundmentl Theorem of Clculus, s we ve discussed in clss, goes s follows: If f is continuous on the intervl [, ] nd F is function for which F t = ft, then ftdt = F F. An integrl

### 7 - Continuous random variables

7-1 Continuous rndom vribles S. Lll, Stnford 2011.01.25.01 7 - Continuous rndom vribles Continuous rndom vribles The cumultive distribution function The uniform rndom vrible Gussin rndom vribles The Gussin

### 10. AREAS BETWEEN CURVES

. AREAS BETWEEN CURVES.. Ares etween curves So res ove the x-xis re positive nd res elow re negtive, right? Wrong! We lied! Well, when you first lern out integrtion it s convenient fiction tht s true in

### Lecture 7 notes Nodal Analysis

Lecture 7 notes Nodl Anlysis Generl Network Anlysis In mny cses you hve multiple unknowns in circuit, sy the voltges cross multiple resistors. Network nlysis is systemtic wy to generte multiple equtions

### x = a To determine the volume of the solid, we use a definite integral to sum the volumes of the slices as we let!x " 0 :

Clculus II MAT 146 Integrtion Applictions: Volumes of 3D Solids Our gol is to determine volumes of vrious shpes. Some of the shpes re the result of rotting curve out n xis nd other shpes re simply given

### 31.2. Numerical Integration. Introduction. Prerequisites. Learning Outcomes

Numericl Integrtion 3. Introduction In this Section we will present some methods tht cn be used to pproximte integrls. Attention will be pid to how we ensure tht such pproximtions cn be gurnteed to be

### Notes on Calculus II Integral Calculus. Miguel A. Lerma

Notes on Clculus II Integrl Clculus Miguel A. Lerm November 22, 22 Contents Introduction 5 Chpter. Integrls 6.. Ares nd Distnces. The Definite Integrl 6.2. The Evlution Theorem.3. The Fundmentl Theorem

### 5: The Definite Integral

5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce

### PDE Notes. Paul Carnig. January ODE s vs PDE s 1

PDE Notes Pul Crnig Jnury 2014 Contents 1 ODE s vs PDE s 1 2 Section 1.2 Het diffusion Eqution 1 2.1 Fourier s w of Het Conduction............................. 2 2.2 Energy Conservtion.....................................

### Lecture 1. Functional series. Pointwise and uniform convergence.

1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

### Test , 8.2, 8.4 (density only), 8.5 (work only), 9.1, 9.2 and 9.3 related test 1 material and material from prior classes

Test 2 8., 8.2, 8.4 (density only), 8.5 (work only), 9., 9.2 nd 9.3 relted test mteril nd mteril from prior clsses Locl to Globl Perspectives Anlyze smll pieces to understnd the big picture. Exmples: numericl

### The Definite Integral

CHAPTER 3 The Definite Integrl Key Words nd Concepts: Definite Integrl Questions to Consider: How do we use slicing to turn problem sttement into definite integrl? How re definite nd indefinite integrls

### Anonymous Math 361: Homework 5. x i = 1 (1 u i )

Anonymous Mth 36: Homewor 5 Rudin. Let I be the set of ll u (u,..., u ) R with u i for ll i; let Q be the set of ll x (x,..., x ) R with x i, x i. (I is the unit cube; Q is the stndrd simplex in R ). Define

### Not for reproduction

AREA OF A SURFACE OF REVOLUTION cut h FIGURE FIGURE πr r r l h FIGURE A surfce of revolution is formed when curve is rotted bout line. Such surfce is the lterl boundry of solid of revolution of the type

### MA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.

MA123, Chpter 1: Formuls for integrls: integrls, ntiderivtives, nd the Fundmentl Theorem of Clculus (pp. 27-233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.

### Week 10: Riemann integral and its properties

Clculus nd Liner Algebr for Biomedicl Engineering Week 10: Riemnn integrl nd its properties H. Führ, Lehrstuhl A für Mthemtik, RWTH Achen, WS 07 Motivtion: Computing flow from flow rtes 1 We observe the

### Problem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:

(x + y ) = y + (x + y ) = x + Problem Set 9 Discussion: Nov., Nov. 8, Nov. (on probbility nd binomil coefficients) The nme fter the problem is the designted writer of the solution of tht problem. (No one

### 2

1 Notes for Numericl Anlysis Mth 5466 by S. Adjerid Virgini Polytechnic Institute nd Stte University (A Rough Drft) 2 Contents 1 Differentition nd Integrtion 5 1.1 Introduction............................

### Distance And Velocity

Unit #8 - The Integrl Some problems nd solutions selected or dpted from Hughes-Hllett Clculus. Distnce And Velocity. The grph below shows the velocity, v, of n object (in meters/sec). Estimte the totl

### Math 211A Homework. Edward Burkard. = tan (2x + z)

Mth A Homework Ewr Burkr Eercises 5-C Eercise 8 Show tht the utonomous system: 5 Plne Autonomous Systems = e sin 3y + sin cos + e z, y = sin ( + 3y, z = tn ( + z hs n unstble criticl point t = y = z =

### Resistors. Consider a uniform cylinder of material with mediocre to poor to pathetic conductivity ( )

10/25/2005 Resistors.doc 1/7 Resistors Consider uniform cylinder of mteril with mediocre to poor to r. pthetic conductivity ( ) ˆ This cylinder is centered on the -xis, nd hs length. The surfce re of the

### Numerical Methods 2007

Physics Mster Course: Numericl Methods 2007 Hns Mssen Rdboud Universiteit Nijmegen Onderwijsinstituut Wiskunde, Ntuur- en Sterrenkunde Toernooiveld 1 6525 ED Nijmegen September 2007 1 Introduction In these

### MATH 222 Second Semester Calculus. Fall 2015

MATH Second Semester Clculus Fll 5 Typeset:August, 5 Mth nd Semester Clculus Lecture notes version. (Fll 5) This is self contined set of lecture notes for Mth. The notes were written by Sigurd Angenent,

### Math 554 Integration

Mth 554 Integrtion Hndout #9 4/12/96 Defn. A collection of n + 1 distinct points of the intervl [, b] P := {x 0 = < x 1 < < x i 1 < x i < < b =: x n } is clled prtition of the intervl. In this cse, we

### Numerical Methods in Economics

Numericl Methods in Economics Alexis Akir Tod November 5, 07 Polynomil interpoltion Polynomils re gret becuse we cn differentite nd integrte them very esily. Since degree n polynomil is determined by n

### CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES

Filomt 27:4 (2013) 649 658 DOI 10.2298/FIL1304649M Published by Fculty of Sciences nd Mthemtics University of Niš Serbi Avilble t: http://www.pmf.ni.c.rs/filomt CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED

### OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS

IJRET: Interntionl Journl of Reserch in Engineering nd Technology eissn: 9-6 pissn: -78 OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS T. Rmchndrn R.Priml

### MATH 174A: PROBLEM SET 5. Suggested Solution

MATH 174A: PROBLEM SET 5 Suggested Solution Problem 1. Suppose tht I [, b] is n intervl. Let f 1 b f() d for f C(I; R) (i.e. f is continuous rel-vlued function on I), nd let L 1 (I) denote the completion

### Topic 1 Notes Jeremy Orloff

Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble

### Calculus 2: Integration. Differentiation. Integration

Clculus 2: Integrtion The reverse process to differentition is known s integrtion. Differentition f() f () Integrtion As it is the opposite of finding the derivtive, the function obtined b integrtion is

### l 2 p2 n 4n 2, the total surface area of the

Week 6 Lectures Sections 7.5, 7.6 Section 7.5: Surfce re of Revolution Surfce re of Cone: Let C be circle of rdius r. Let P n be n n-sided regulr polygon of perimeter p n with vertices on C. Form cone

### Practice final exam solutions

University of Pennsylvni Deprtment of Mthemtics Mth 26 Honors Clculus II Spring Semester 29 Prof. Grssi, T.A. Asher Auel Prctice finl exm solutions 1. Let F : 2 2 be defined by F (x, y (x + y, x y. If

### Problem. Statement. variable Y. Method: Step 1: Step 2: y d dy. Find F ( Step 3: Find f = Y. Solution: Assume

Functions of Rndom Vrible Problem Sttement We know the pdf ( or cdf ) of rndom r vrible. Define new rndom vrible Y = g. Find the pdf of Y. Method: Step : Step : Step 3: Plot Y = g( ). Find F ( y) by mpping

### Precalculus Spring 2017

Preclculus Spring 2017 Exm 3 Summry (Section 4.1 through 5.2, nd 9.4) Section P.5 Find domins of lgebric expressions Simplify rtionl expressions Add, subtrct, multiply, & divide rtionl expressions Simplify

### Linear Systems with Constant Coefficients

Liner Systems with Constnt Coefficients 4-3-05 Here is system of n differentil equtions in n unknowns: x x + + n x n, x x + + n x n, x n n x + + nn x n This is constnt coefficient liner homogeneous system

### Math 0230 Calculus 2 Lectures

Mth Clculus Lectures Chpter 7 Applictions of Integrtion Numertion of sections corresponds to the text Jmes Stewrt, Essentil Clculus, Erly Trnscendentls, Second edition. Section 7. Ares Between Curves Two

### Topic 6b Finite Difference Approximations

/8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?

### MATH 409 Advanced Calculus I Lecture 19: Riemann sums. Properties of integrals.

MATH 409 Advnced Clculus I Lecture 19: Riemnn sums. Properties of integrls. Drboux sums Let P = {x 0,x 1,...,x n } be prtition of n intervl [,b], where x 0 = < x 1 < < x n = b. Let f : [,b] R be bounded

### Mapping the delta function and other Radon measures

Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support

### ECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance

Generl structure ECO 37 Economics of Uncertinty Fll Term 007 Notes for lectures 4. Stochstic Dominnce Here we suppose tht the consequences re welth mounts denoted by W, which cn tke on ny vlue between

### Homework Assignment 3 Solution Set

Homework Assignment 3 Solution Set PHYCS 44 6 Ferury, 4 Prolem 1 (Griffiths.5(c The potentil due to ny continuous chrge distriution is the sum of the contriutions from ech infinitesiml chrge in the distriution.

### FINALTERM EXAMINATION 2011 Calculus &. Analytical Geometry-I

FINALTERM EXAMINATION 011 Clculus &. Anlyticl Geometry-I Question No: 1 { Mrks: 1 ) - Plese choose one If f is twice differentible function t sttionry point x 0 x 0 nd f ''(x 0 ) > 0 then f hs reltive...

### Keywords : Generalized Ostrowski s inequality, generalized midpoint inequality, Taylor s formula.

Generliztions of the Ostrowski s inequlity K. S. Anstsiou Aristides I. Kechriniotis B. A. Kotsos Technologicl Eductionl Institute T.E.I.) of Lmi 3rd Km. O.N.R. Lmi-Athens Lmi 3500 Greece Abstrct Using

### 1 i n x i x i 1. Note that kqk kp k. In addition, if P and Q are partition of [a, b], P Q is finer than both P and Q.

Chpter 6 Integrtion In this chpter we define the integrl. Intuitively, it should be the re under curve. Not surprisingly, fter mny exmples, counter exmples, exceptions, generliztions, the concept of the

### Theoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4

WiSe 1 8.1.1 Prof. Dr. A.-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Mtthis Sb m Lehrstuhl für Theoretische Physik I Deprtment für Physik Friedrich-Alexnder-Universität Erlngen-Nürnberg Theoretische

### Convex Sets and Functions

B Convex Sets nd Functions Definition B1 Let L, +, ) be rel liner spce nd let C be subset of L The set C is convex if, for ll x,y C nd ll [, 1], we hve 1 )x+y C In other words, every point on the line

### Unit #10 De+inite Integration & The Fundamental Theorem Of Calculus

Unit # De+inite Integrtion & The Fundmentl Theorem Of Clculus. Find the re of the shded region ove nd explin the mening of your nswer. (squres re y units) ) The grph to the right is f(x) = -x + 8x )Use

### Applied Physics Introduction to Vibrations and Waves (with a focus on elastic waves) Course Outline

Applied Physics Introduction to Vibrtions nd Wves (with focus on elstic wves) Course Outline Simple Hrmonic Motion && + ω 0 ω k /m k elstic property of the oscilltor Elstic properties of terils Stretching,

### UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE

UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE 1. Pointwise Convergence of Sequence Let E be set nd Y be metric spce. Consider functions f n : E Y for n = 1, 2,.... We sy tht the sequence

### STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

### Normal Distribution. Lecture 6: More Binomial Distribution. Properties of the Unit Normal Distribution. Unit Normal Distribution

Norml Distribution Lecture 6: More Binomil Distribution If X is rndom vrible with norml distribution with men µ nd vrince σ 2, X N (µ, σ 2, then P(X = x = f (x = 1 e 1 (x µ 2 2 σ 2 σ Sttistics 104 Colin

### CSCI FOUNDATIONS OF COMPUTER SCIENCE

1 CSCI- 2200 FOUNDATIONS OF COMPUTER SCIENCE Spring 2015 My 7, 2015 2 Announcements Homework 9 is due now. Some finl exm review problems will be posted on the web site tody. These re prcqce problems not

### METHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS

Journl of Young Scientist Volume III 5 ISSN 44-8; ISSN CD-ROM 44-9; ISSN Online 44-5; ISSN-L 44 8 METHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS An ALEXANDRU Scientific Coordintor: Assist

### Week 7 Riemann Stieltjes Integration: Lectures 19-21

Week 7 Riemnn Stieltjes Integrtion: Lectures 19-21 Lecture 19 Throughout this section α will denote monotoniclly incresing function on n intervl [, b]. Let f be bounded function on [, b]. Let P = { = 0

### 4.1. Probability Density Functions

STT 1 4.1-4. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile - vers - discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of

### Ordinary Differential Equations- Boundary Value Problem

Ordinry Differentil Equtions- Boundry Vlue Problem Shooting method Runge Kutt method Computer-bsed solutions o BVPFD subroutine (Fortrn IMSL subroutine tht Solves (prmeterized) system of differentil equtions

### Non-Linear & Logistic Regression

Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find

### Section 7.1 Area of a Region Between Two Curves

Section 7.1 Are of Region Between Two Curves White Bord Chllenge The circle elow is inscried into squre: Clcultor 0 cm Wht is the shded re? 400 100 85.841cm White Bord Chllenge Find the re of the region

### Fundamental Methods of Numerical Extrapolation With Applications

Fundmentl Methods of Extrpoltion 1 Fundmentl Methods of Numericl Extrpoltion With Applictions Eric Hung-Lin Liu ehliu@mit.edu Keywords: numericl nlysis, extrpoltion, richrdson, romberg, numericl differentition,

### Cu 3 (PO 4 ) 2 (s) 3 Cu 2+ (aq) + 2 PO 4 3- (aq) circle answer: pure water or Na 3 PO 4 solution This is the common-ion effect.

CHEM 1122011 NAME: ANSWER KEY Vining, Exm # SHORT ANSWER QUESTIONS ============= 4 points ech ============ All work must be shown. 1. Wht re [H O ] nd [OH ] for solution tht hs ph of 9.0? Choose one. ()

### Efficient Computation of a Class of Singular Oscillatory Integrals by Steepest Descent Method

Applied Mthemticl Sciences, Vol. 8, 214, no. 31, 1535-1542 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/1.12988/ms.214.43166 Efficient Computtion of Clss of Singulr Oscilltory Integrls by Steepest Descent

### Improper Integrals. Introduction. Type 1: Improper Integrals on Infinite Intervals. When we defined the definite integral.

Improper Integrls Introduction When we defined the definite integrl f d we ssumed tht f ws continuous on [, ] where [, ] ws finite, closed intervl There re t lest two wys this definition cn fil to e stisfied:

### 1 Nondeterministic Finite Automata

1 Nondeterministic Finite Automt Suppose in life, whenever you hd choice, you could try oth possiilities nd live your life. At the end, you would go ck nd choose the one tht worked out the est. Then you