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


 Scott Jordan
 9 months ago
 Views:
Transcription
1 Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t UrbnChmpign 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 NewtonCotes (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 NewtonCotes 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 (NewtonCotes) 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 NewtonCotes (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 NewtonCotes (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 NewtonCotes... 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 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 informationCOSC 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 informationNumerical 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 informationInterpolation. 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
More informationNUMERICAL 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 informationIII. Lecture on Numerical Integration. File faclib/dattab/lecturenotes/numericalinter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecturenotes/numericalinter03.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 informationTHE 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 2Dimensionl Gussin Qudrture 20 4
More informationOrthogonal Polynomials and LeastSquares Approximations to Functions
Chpter Orthogonl Polynomils nd LestSqures Approximtions to Functions **4/5/3 ET. Discrete LestSqures Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationDiscrete Leastsquares Approximations
Discrete Lestsqures 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 information20 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
More informationLecture 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
More informationPolynomial 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 information1 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 informationContinuous 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 relvlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationMath 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
More informationIf 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 ntiderivtives. We now introduce some methods for finding ntiderivtives: If u = g(x) is differentible
More informationThe 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
More information63. 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 informationIntroduction 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
More informationSection 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
More informationSection 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 informationMTH 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
More informationWe 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
More informationLine 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
More informationHarmonic Mean Derivative  Based Closed Newton Cotes Quadrature
IOSR Journl of Mthemtics (IOSRJM) eissn:  pissn: 9X. Volume Issue Ver. IV (My.  Jun. 0) PP  www.iosrjournls.org Hrmonic Men Derivtive  Bsed Closed Newton Cotes Qudrture T. Rmchndrn D.Udykumr nd
More informationNumerical 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,
More informationNumerical 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(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
More information7. 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
More informationMath 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
More informationChapter 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
More informationMath 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:
More informationThe 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
More information7.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
More informationImproper 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
More information7  Continuous random variables
71 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
More information10. AREAS BETWEEN CURVES
. AREAS BETWEEN CURVES.. Ares etween curves So res ove the xxis 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
More informationLecture 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
More informationx = 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
More information31.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
More informationNotes 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
More information5: 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
More informationPDE 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.....................................
More informationLecture 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
More informationTest , 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
More informationThe 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
More informationAnonymous 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
More informationNot 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
More informationMA123, 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. 27233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.
More informationWeek 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
More informationProblem 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
More information2
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............................
More informationDistance And Velocity
Unit #8  The Integrl Some problems nd solutions selected or dpted from HughesHllett Clculus. Distnce And Velocity. The grph below shows the velocity, v, of n object (in meters/sec). Estimte the totl
More informationMath 211A Homework. Edward Burkard. = tan (2x + z)
Mth A Homework Ewr Burkr Eercises 5C 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 =
More informationResistors. 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
More informationNumerical 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
More informationMATH 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,
More informationMath 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
More informationNumerical 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
More informationCLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTONCOTES 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
More informationOPEN NEWTON  COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS
IJRET: Interntionl Journl of Reserch in Engineering nd Technology eissn: 96 pissn: 78 OPEN NEWTON  COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS T. Rmchndrn R.Priml
More informationMATH 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 relvlued function on I), nd let L 1 (I) denote the completion
More informationTopic 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
More informationCalculus 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
More informationl 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 nsided regulr polygon of perimeter p n with vertices on C. Form cone
More informationPractice 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
More informationProblem. 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
More informationPrecalculus 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
More informationLinear Systems with Constant Coefficients
Liner Systems with Constnt Coefficients 4305 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
More informationMath 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
More informationTopic 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?
More informationMATH 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
More informationMapping 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
More informationECO 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
More informationHomework 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.
More informationFINALTERM EXAMINATION 2011 Calculus &. Analytical GeometryI
FINALTERM EXAMINATION 011 Clculus &. Anlyticl GeometryI 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...
More informationKeywords : 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. LmiAthens Lmi 3500 Greece Abstrct Using
More information1 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
More informationTheoretische 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 FriedrichAlexnderUniversität ErlngenNürnberg Theoretische
More informationConvex 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
More informationUnit #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
More informationApplied 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,
More informationUNIFORM 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
More informationSTEP 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
More informationNormal 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
More informationCSCI 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
More informationMETHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS
Journl of Young Scientist Volume III 5 ISSN 448; ISSN CDROM 449; ISSN Online 445; ISSNL 44 8 METHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS An ALEXANDRU Scientific Coordintor: Assist
More informationWeek 7 Riemann Stieltjes Integration: Lectures 1921
Week 7 Riemnn Stieltjes Integrtion: Lectures 1921 Lecture 19 Throughout this section α will denote monotoniclly incresing function on n intervl [, b]. Let f be bounded function on [, b]. Let P = { = 0
More information4.1. Probability Density Functions
STT 1 4.14. 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
More informationOrdinary Differential Equations Boundary Value Problem
Ordinry Differentil Equtions Boundry Vlue Problem Shooting method Runge Kutt method Computerbsed solutions o BVPFD subroutine (Fortrn IMSL subroutine tht Solves (prmeterized) system of differentil equtions
More informationNonLinear & Logistic Regression
NonLiner & 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
More informationSection 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
More informationFundamental Methods of Numerical Extrapolation With Applications
Fundmentl Methods of Extrpoltion 1 Fundmentl Methods of Numericl Extrpoltion With Applictions Eric HungLin Liu ehliu@mit.edu Keywords: numericl nlysis, extrpoltion, richrdson, romberg, numericl differentition,
More informationCu 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 commonion 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. ()
More informationEfficient Computation of a Class of Singular Oscillatory Integrals by Steepest Descent Method
Applied Mthemticl Sciences, Vol. 8, 214, no. 31, 15351542 HIKARI Ltd, www.mhikri.com http://dx.doi.org/1.12988/ms.214.43166 Efficient Computtion of Clss of Singulr Oscilltory Integrls by Steepest Descent
More informationImproper 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:
More information1 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
More informationSynoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?
Synoptic Meteorology I: Finite Differences 1618 September 2014 Prtil Derivtives (or, Why Do We Cre About Finite Differences?) With the exception of the idel gs lw, the equtions tht govern the evolution
More informationAP Calculus BC Review Applications of Integration (Chapter 6) noting that one common instance of a force is weight
AP Clculus BC Review Applictions of Integrtion (Chpter Things to Know n Be Able to Do Fin the re between two curves by integrting with respect to x or y Fin volumes by pproximtions with cross sections:
More information