Lecture 14: Quadrature


 Easter Mason
 9 months ago
 Views:
Transcription
1 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 relvlues nd smooth The pproximtion of n integrl by numericl method is commonly referred to s qudrture In Clculus one discusses mny techniques for determining n ntiderivtive of fx) The integrl is then computed by evluting the ntiderivtive t the endpoints of the intervl However, mny integrnds of interest in science nd engineering do not hve known ntiderivtive Moreover, in some pplictions the integrnd is only known t few points in the intervl We would like to be ble to determine ccurte pproximtions of integrls lso in these situtions We will pproximte integrls 1) by sums fx j )w j ) These sums re referred to s qudrture rules The x j re the nodes nd the w j the weights of the qudrture rule We re interested in determining weights so tht the qudrture rules gives ccurte pproximtions of the integrl 1 for lrge clsses of integrnds fx) The nodes often cnnot be chosen freely, becuse the integrnd my only be known t certin points For instnce, the function fx) might not be explicitly known, only mesured vlues fx j ) t the nodes x j my be vilble We remrk tht you lredy encountered the pproximtion of integrls by sums in Clculus There, however, one ws less concerned with how smll the error E N f) = fx)dx fx j )w j 3) is for smll vlues of N; it ws sufficient tht the error E N f) converged to zero s N incresed to infinity In this lecture, we would like the error E N f) to be smll lredy for smll to modest number of terms N in the qudrture rule ) Method of undetermined coefficients Let the nodes x j, 1 j N, be given Throughout this lecture the nodes will be ordered so tht x 1 < x < < x N b 4) We would like to determine the weights w j so tht qudrture rule is exct for polynomils of s high degree s possible, ie, we would like the weights be such tht px)dx = for ll polynomils px) of s high degree s possible px j )w j 5) 1
2 The powers {1, x, x, x 3, } form bsis for polynomils Our objective therefore cn be expressed s follows: We would like the equlity 5) to hold for s mny powers px) = x j, j =, 1,,, s possible This requirement gives rise to the liner system of equtions for the weights, px) = 1 : w 1 + w + + w N = px) = x : x 1 w 1 + x w + + x N w N = px) = x : x 1w 1 + x w + + x N w N = px) = x N 1 : x N 1 1 w 1 + x N 1 w + + w N x N 1 N This system conveniently cn be expressed in the form x 1 x x N x N 1 1 x N 1 x N 1 N w 1 w w N = dx = b, xdx = 1 b ), x 1 dx = 3 b3 3 ), = xn 1 dx = 1 N bn N ) b 1 b ) 1 N bn N ) 6) 7) We recognize the mtrix s the trnspose of Vndermonde mtrix Vndermonde mtrices were encountered in the fll semester in connection with polynomil interpoltion, when we showed tht Vndermonde mtrices determined by distinct nodes re nonsingulr Determinnts re invrint under trnsposition This follows from the fct tht determinnts cn be expnded either by rows or by columns We conclude tht the mtrix in 7) is nonsingulr when the nodes x j re distinct Hence the bove liner system of equtions hs unique solution Thus, given n distinct nodes, we cn determine n unique weights w j, such tht the qudrture rule ) integrtes ll polynomils of degree strictly less thn n exctly Exmple 1 Consider the midpoint qudrture rule fx)dx fx 1 )w 1, x 1 = 1 + b) We only hve one weight to determine The liner system of equtions 6) reduces to px) = 1 : w 1 = dx = b This determines the weight w 1 Exmple Apply the midpoint rule to integrte the function fx) = exp 1 x/) ) 8) on the intervl [, 1] The midpoint rule pproximtes the integrnd by the constnt function f1/) = exp 7/8) on the intervl [, b] nd integrtes the ltter function exctly The blue grph of Figure 1 shows fx) nd the vlue of the integrl is the re below this grph The dshed red line displyes the constnt function tht is integrted exctly by the midpoint rule The vlue determined by the midpoint rule is the re below the dshed red line Sometimes Nnode qudrture rules determined by solving the liner system of equtions 7) lso integrte powers x k for k N exctly The following exmple illustrtes this
3 Figure 1: Integrnd of Exmple nd pproximtion used by the midpoint rule Exmple 3 Consider the midpoint qudrture rule of Exmple 1 Appliction of this rule to fx) = x yields x 1 w 1 = 1 + b)b ) = 1 b ) The righthnd side equls xdx, ie, the midpoint rule integrtes not only constnts but lso liner polynomils exctly Exercise 1 Determine the weights of the node qudrture rule fx)dx fx 1 )w 1 + fx )w, x 1 =, x = b This rule is known s the trpezoidl rule By construction, it integrtes liner polynomils exctly Does it integrte polynomils of higher degree exctly? Justify your nswer Exercise Determine the weights of the 3node qudrture rule 1 fx)dx fx 1 )w 1 + fx )w + fx 3 )w 3, x 1 =, x = 1/, x 3 = 1 This rule is known s Simpson s rule It integrtes qudrtic polynomils exctly by construction Does it integrte polynomils of higher degree exctly? Justify your nswer Wht is the nlogous qudrture rule for the intervl [1, ]? Hint: Do not recompute the qudrture rule, just mke chnge of vribles Integrting Lgrnge polynomils When discussing polynomil interpoltion lst fll, we considered different polynomil bses, such s monomils, Lgrnge polynomils, nd Newton polynomils We cn lso in the present context use bses different 3
4 from the monomil one For instnce, the nodes x 1, x,,x N determine the Lgrnge polynomils l k x) = N j k x x j x k x j, k = 1,,, N, 9) which form bsis for ll polynomils of degree t most N 1 We therefore my require tht n Nnode qudrture rule integrtes the Lgrnge polynomils 9) exctly Substituting px) = l k x) into 5) yields l k x)dx = l k x j )w j = l k x k )w k = w k, k = 1,,,N, where the simplifictions of the righthnd side follow from the fct tht { 1, k = j, l k x j ) =, k j Thus, the weights cn be obtined by integrting the Lgrnge polynomils However, the evlution of these integrls is tedious unless N is smll Composite qudrture rules nd Tylor expnsion Assume tht the midpoint rule of Exmple does not yield sufficiently ccurte pproximtion of the integrl We my then subdivide the intervl [, 1] into subintervls nd pply the midpoint rule on ech subintervl This defines the composite midpoint rule Figure : Integrnd of Exmple 4 nd pproximtion used by the composite midpoint rule obtined by dividing the intervl [, 1] into two subintervls of equl length 4
5 Exmple 4 Divide the intervl [, 1] into the subintervls [, 1/] nd [1/, 1] nd pply the midpoint rule on ech subintervl to the integrnd fx) defined by 8) This yields the composite midpoint rule 1 fx)dx fx 1 )w 1 + fx )w with x 1 = 1/4, x = 3/4, nd w 1 = w = 1/ The integrl is pproximted by the re below the dshed curve of Figure Compring the grphs of Figures 1 nd suggests tht subdivision of the intervl should increse the ccurcy of the computed pproximtion While the computtion of the weights by solving the liner system of equtions 7) is esily done in MATLAB or Octve, this pproch does not shed ny light on how the error behves when we increse the number of nodes Further insight on the behvior of qudrture rules cn be gined by expnding the integrnd into Tylor series Consider the pproximtion of the integrl fx)dx by the midpoint rule nd use the Tylor expnsion of fx) t x = h/ Here x h/ should be thought of s firly smll, nd we ssume tht fx) is continuously differentible s mny times s required Then ) ) h h fx) = f + f x h ) + f ) h! x h ) + f ) h x h 3! Integrting the lefthnd nd righthnd sides from to h yields ) h ) h fx)dx = f dx + f x h + f ) h x h ) 3 dx + f ) h 3! 4! ) dx + f h! ) 3 + f ) h x h 4 + 4! ) ) x h ) 4 dx + x h ) dx The integrl over odd powers of x h/ vnishes due to symmetry The righthnd side therefore simplifies to ) h fx)dx = f h + f ) h x h ) dx + f ) h x h 4 dx +! 4! ) ) h = f h + f ) h h 3 + f ) h h The midpoint rule pplied to the integrl on the lefthnd side gives the first term in the righthnd side The remining terms express the qudrture error Hence, this error is given by ) h fx)dx f h = f ) h h 3 + f ) h h 5 + 1) 1 8 Since ll derivtives of order nd higher of liner function vnish, the righthnd side vnishes for such functions It follows tht the midpoint rule is exct for liner functions We know this lredy, nd here it is consequence of the expnsion of the integrl in powers of h Exercise 3 Use Tylor expnsions round x = nd x = h to determine how the pproximtion of the integrl fx)dx determined by the trpezoidl rule depends on h Thus, express the qudrture error fx)dx h f) + fh)) 5
6 s function of h in similrly mnner s 1) Consider the pproximtion of the integrl 1) by the Npoint composite midpoint rule with the nodes nd weights x j = + j 1 ) h, h = b N 1, w j = h, j = 1,,,N 11) Thus, the qudrture rule is given by M h f) = h fx j ) Anlogously to 1), we obtin for ech subintervl [x j h/, x j + h/] of length h the expression xj+h/ x j h/ nd summing over j = 1,,,N yields fx) M h f) = fx)dx fx j )h = f x j ) h 3 +, 1) 1 f x j ) h 3 + = 1 1 N f x j ) Nh ) The verge of the second derivtive vlues 1 N N f x j ) is not smller thn the smllest of the vlues f x j ) Similrly, the verge is not lrger thn lrgest of the vlues f x j ) This is expressed by the inequlities min 1 j N f x j ) 1 N f x j ) mx 1 j N f x j ) We ssumed f x) to be continuous Therefore there is vlue ξ in [, b], such tht f ξ) = 1 N f x j ) Substituting this expression into the righthnd side of 13) yields Multiplying the expression for h in 11) by the denomintor gives fx) M h f) = f ξ) Nh ) which when substituted into 14) yields Nh = b + h, fx) M h f) = f ξ) b 1 h + 15) The righthnd side indictes tht we cn expect the error to decrese by fctor 4 when h is hlved nd the number of nodes x j is doubled) 6
7 Generlly, one does not know in dvnce how mny nodes to use in order to chieve desired ccurcy We therefore my be interested in hlving h until consecutive pproximtions M h f), M h/ f), M h/4 f),, of the integrl do not vry much when h is further reduced Exercise 4 In pplictions with complicted functions, the evlution of the function vlues my dominte the rithmetic work required to evlute qudrture rules The computtion of M h f) requires the evlution of N function vlues How mny of these function vlues cn be used gin when evluting M h/ f)? Exercise 5 Consider the Npoint composite trpezoidl rule for the pproximtion of 1), 1 T h f) = h fx 1) + fx ) + fx 3 ) + + fx N 1 ) + 1 ) fx N) with ) Use the representtion x j = + j 1)h, T h f) = N 1 h = b, j = 1,,, N N 1 h fx j) + fx j+1 )) nd the result of Exercise 3 to determine n expnsion of the qudrture error similr to 13) b) Wht is the nlog of formul 15)? How much is the error reduced when h is hlved? c) Assume tht T h f) hs been evluted, nd we would like to compute T h/ f) How mny dditionl function evlutions re required? Singulr integrnds, infinite intervls, nd dptive qudrture rules The composite midpoint nd trpezoidl rules cn be used lso for functions tht re not twice continuously differentible However, the error then my converge to zero slower when h is reduced nd mny nodes might be required to give smll qudrture error We discuss n pproch to remedy this sitution Consider the evlution of the integrl 1 fx) lnx)dx 16) where fx) is ssumed to be smooth function We note tht the trpezoidl rule cnnot be used since the integrnd is infinite t x = Nevertheless, the integrl exists for mny smooth functions fx) We therefore my consider using composite midpoint rule; see Exercise 6 Alterntively, the method of undetermined coefficients cn be pplied to the function fx) only Thus, we would like to determine qudrture rule of the form fx j )w j with given nodes, nd we seek to determine weights w j so tht the reltion px) lnx)dx = px j )w j 17) 7
8 holds for ll polynomils of s high degree s possible This yields the liner system of equtions, which is nlogous to 6), px) = 1 : w 1 + w + + w N = px) = x : x 1 w 1 + x w + + x N w N = px) = x : x 1w 1 + x w + + x N w N = px) = x N 1 : x N 1 1 w 1 + x N 1 w + + w N x N 1 N lnx)dx, b xlnx)dx, b x lnx)dx, = xn 1 lnx)dx The righthnd side expressions cn be evluted by integrtion by prts Note tht the lloction of nodes x j is quite rbitrry; in prticulr, we my put node t the origin Exercise 6 Compute n pproximtion of 16) with fx) = expx ) by the node composite midpoint rule Use the MATLAB function qud to determine the exct vlue Note tht the MATLAB function for the integrnd hs to written to llow vector rguments, eg, f=expx^)*logx); How lrge is the error? The MATLAB function qud pplies composite Simpson rule The MATLAB function qud is n exmple of n dptive qudrture rule Adptive rules re composite rules tht estimte the qudrture error by compring the result obtined with different mesh sizes h New nodes re llocted in subintervls [x j, x j+1 ] in which the qudrture error is deemed to be lrger thn prescribed tolernce Exercise 7 Determine point qudrture rule with the sme nodes s in Exercise 6 by solving liner system of equtions of the form 18) Apply it to the function fx) of Exercise 6 Is this rule more ccurte thn the one in Exmple 6? Exercise 8 Evlute the integrl exp x 3 )dx with t lest 5 correct deciml digits using the MATLAB function qud This function requires the intervl of integrtion to be finite Therefore the integrl hs to be split, exp x 3 )dx = c exp x 3 )dx + c exp x 3 )dx, where the constnt c > is chosen lrge enough so tht the second integrl cn be neglected nd the first integrl is evluted with the function qud For instnce, we my choose c lrge enough so tht c exp x 3 )dx Determine such vlue of c nd justify your choice The error bound for the integrl evluted with qud then should not be lrger thn In this lecture, we hve fixed the nodes nd then determined the weights so tht we integrte polynomils of s high degree s possible A clever choice of nodes s well s weights in n Npoint qudrture rule mkes it possible to integrte polynomils of degree up to N 1 exctly The midpoint rule is n exmple These rules re known s Gussin qudrture rules The nodes nd weights generlly re not known in closed form, but they cn be computed quite efficiently by solving n eigenvlue problem symmetric N N tridigonl mtrix This is often done with the QRlgorithm or modifiction thereof 18) 8
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 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 informationLecture 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 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 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 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 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 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 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 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 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 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 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 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 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 informationa a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.
Section 9 The Lplce Expnsion In the lst section, we defined the determinnt of (3 3) mtrix A 12 to be 22 12 21 22 2231 22 12 21. In this section, we introduce generl formul for computing determinnts. Rewriting
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 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 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 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 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 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 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 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 informationThe Trapezoidal Rule
_.qd // : PM Pge 9 SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion
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 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 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 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 information2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).
AP Clculus BC Review Chpter 8 Prt nd Chpter 9 Things to Know nd Be Ale to Do Know everything from the first prt of Chpter 8 Given n integrnd figure out how to ntidifferentite it using ny of the following
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 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 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 informationMATH 573 FINAL EXAM. May 30, 2007
MATH 573 FINAL EXAM My 30, 007 NAME: Solutions 1. This exm is due Wednesdy, June 6 efore the 1:30 pm. After 1:30 pm I will NOT ccept the exm.. This exm hs 1 pges including this cover. There re 10 prolems.
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 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 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 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 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 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 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 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 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 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 informationMath 4200: Homework Problems
Mth 4200: Homework Problems Gregor Kovčič 1. Prove the following properties of the binomil coefficients ( n ) ( n ) (i) 1 + + + + 1 2 ( n ) (ii) 1 ( n ) ( n ) + 2 + 3 + + n 2 3 ( ) n ( n + = 2 n 1 n) n,
More informationProblem Set 3
14.102 Problem Set 3 Due Tuesdy, October 18, in clss 1. Lecture Notes Exercise 208: Find R b log(t)dt,where0
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 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 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 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 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 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 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 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 informationSummary of Elementary Calculus
Summry of Elementry Clculus Notes by Wlter Noll (1971) 1 The rel numbers The set of rel numbers is denoted by R. The set R is often visulized geometriclly s numberline nd its elements re often referred
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 informationMACsolutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences  Finite Elements  Finite Volumes  Boundry Elements MACsolutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
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 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 informationUSA Mathematical Talent Search Round 1 Solutions Year 21 Academic Year
1/1/21. Fill in the circles in the picture t right with the digits 18, one digit in ech circle with no digit repeted, so tht no two circles tht re connected by line segment contin consecutive digits.
More informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCKKURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When relvlued
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 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 informationA basic logarithmic inequality, and the logarithmic mean
Notes on Number Theory nd Discrete Mthemtics ISSN 30 532 Vol. 2, 205, No., 3 35 A bsic logrithmic inequlity, nd the logrithmic men József Sándor Deprtment of Mthemtics, BbeşBolyi University Str. Koglnicenu
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 informationS. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:
FACTA UNIVERSITATIS NIŠ) Ser Mth Inform 9 00) 6 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Dedicted to Prof G Mstroinni for his 65th birthdy
More informationChapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY
Chpter 3 MTRIX In this chpter: Definition nd terms Specil Mtrices Mtrix Opertion: Trnspose, Equlity, Sum, Difference, Sclr Multipliction, Mtrix Multipliction, Determinnt, Inverse ppliction of Mtrix in
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The twodimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More informationContinuous probability distributions
Chpter 1 Continuous probbility distributions 1.1 Introduction We cll x continuous rndom vrible in x b if x cn tke on ny vlue in this intervl. An exmple of rndom vrible is the height of dult humn mle, selected
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 informationPart III. Calculus of variations
Prt III. Clculus of vritions Lecture notes for MA342H P. Krgeorgis pete@mths.tcd.ie 1/43 Introduction There re severl pplictions tht involve expressions of the form Jy) = Lx,yx),y x))dx. For instnce, Jy)
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 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 informationSCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics
SCHOOL OF ENGINEERING & BUIL ENVIRONMEN Mthemtics An Introduction to Mtrices Definition of Mtri Size of Mtri Rows nd Columns of Mtri Mtri Addition Sclr Multipliction of Mtri Mtri Multipliction 7 rnspose
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 informationFirst variation. (onevariable problem) January 14, 2013
First vrition (onevrible problem) Jnury 14, 2013 Contents 1 Sttionrity of n integrl functionl 2 1.1 Euler eqution (Optimlity conditions)............... 2 1.2 First integrls: Three specil cses.................
More informationFarey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University
U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions
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 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 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 informationMatrix Eigenvalues and Eigenvectors September 13, 2017
Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues
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 informationThe problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.
ADVANCED CALCULUS PRACTICE PROBLEMS JAMES KEESLING The problems tht follow illustrte the methods covered in clss. They re typicl of the types of problems tht will be on the tests. 1. Riemnn Integrtion
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 informationD i k B(n;D)= X LX ;d= l= A (il) d n d d+ (A ) (kl)( + )B(n;D); where c B(:::;n c;:::) = B(:::;n c ;:::), in prticulr c n = n c. Using the reltions [n
Explicit solutions of the multi{loop integrl recurrence reltions nd its ppliction? P. A. BAKOV ; nstitute of ucler Physics, Moscow Stte University, Moscow 9899, Russi The pproch to the constructing explicit
More informationQUADRATIC EQUATIONS OBJECTIVE PROBLEMS
QUADRATIC EQUATIONS OBJECTIVE PROBLEMS +. The solution of the eqution will e (), () 0,, 5, 5. The roots of the given eqution ( p q) ( q r) ( r p) 0 + + re p q r p (), r p p q, q r p q (), (d), q r p q.
More informationCalculus MATH 172Fall 2017 Lecture Notes
Clculus MATH 172Fll 2017 Lecture Notes These notes re concise summry of wht hs been covered so fr during the lectures. All the definitions must be memorized nd understood. Sttements of importnt theorems
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 informationMatrix & Vector Basic Linear Algebra & Calculus
Mtrix & Vector Bsic Liner lgebr & lculus Wht is mtrix? rectngulr rry of numbers (we will concentrte on rel numbers). nxm mtrix hs n rows n m columns M x4 M M M M M M M M M M M M 4 4 4 First row Secon row
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 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 informationUniversitaireWiskundeCompetitie. Problem 2005/4A We have k=1. Show that for every q Q satisfying 0 < q < 1, there exists a finite subset K N so that
Problemen/UWC NAW 5/7 nr juni 006 47 Problemen/UWC UniversitireWiskundeCompetitie Edition 005/4 For Session 005/4 we received submissions from Peter Vndendriessche, Vldislv Frnk, Arne Smeets, Jn vn de
More information1 Online Learning and Regret Minimization
2.997 DecisionMking in LrgeScle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
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 informationA. Limits  L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.
A. Limits  L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c
More informationC1M14. Integrals as Area Accumulators
CM Integrls s Are Accumultors Most tetbooks do good job of developing the integrl nd this is not the plce to provide tht development. We will show how Mple presents Riemnn Sums nd the ccompnying digrms
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 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 informationThe Trapezoidal Rule
SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion Approimte
More information