4.4 Improper Integrals
|
|
- Michael Moody
- 5 years ago
- Views:
Transcription
1 4.4 Improper Integrls 135 converges (with prmeters s shown bove) on the very first extrpoltion, fter just 5 clls to trpzd, while qsimp requires 8 clls (8 times s mny evlutions of the integrnd) nd qtrp requires 13 clls (mking 256 times s mny evlutions of the integrnd). CITED REFERENCES AND FURTHER READING: Stoer, J., nd Bulirsch, R. 1980, Introduction to Numericl Anlysis (New York: Springer-Verlg), Dhlquist, G., nd Bjorck, A. 1974, Numericl Methods (Englewood Cliffs, NJ: Prentice-Hll), Rlston, A., nd Rbinowitz, P. 1978, A First Course in Numericl Anlysis, 2nd ed. (New York: McGrw-Hill), Improper Integrls For our present purposes, n integrl will be improper if it hs ny of the following problems: its integrnd goes to finite limiting vlue t finite upper nd lower limits, but cnnot be evluted right on one of those limits (e.g., sin x/x t x =0) its upper limit is, or its lower limit is it hs n integrble singulrity t either limit (e.g., x 1/2 t x =0) it hs n integrble singulrity t known plce between its upper nd lower limits it hs n integrble singulrity t n unknown plce between its upper nd lower limits If n integrl is infinite (e.g., x 1 dx), or does not exist in limiting sense 1 (e.g., cos xdx), we do not cll it improper; we cll it impossible. No mount of clever lgorithmics will return meningful nswer to n ill-posed problem. In this section we will generlize the techniques of the preceding two sections to cover the first four problems on the bove list. A more dvnced discussion of qudrture with integrble singulrities occurs in Chpter 18, notbly The fifth problem, singulrity t unknown loction, cn relly only be hndled by the use of vrible stepsize differentil eqution integrtion routine, s will be given in Chpter 16. We need workhorse like the extended trpezoidl rule (eqution ), but one which is n open formul in the sense of 4.1, i.e., does not require the integrnd to be evluted t the endpoints. Eqution (4.1.19), the extended midpoint rule, is the best choice. The reson is tht (4.1.19) shres with (4.1.11) the deep property of hving n error series tht is entirely even in h. Indeed there is formul, not s well known s it ought to be, clled the Second Euler-Mclurin summtion formul, xn h[f 3/2 + f 5/2 + f 7/2 + +f N 3/2 +f N 1/2 ] x 1 + B 2h 2 4 (f N f 1 )+ (4.4.1) + B 2kh 2k (1 2 2k+1 )(f (2k 1) N f (2k 1) 1 )+ (2k)! Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ).
2 136 Chpter 4. Integrtion of Functions This eqution cn be derived by writing out (4.2.1) with stepsize h, then writing it out gin with stepsize h/2, then subtrcting the first from twice the second. It is not possible to double the number of steps in the extended midpoint rule nd still hve the benefit of previous function evlutions (try it!). However, it is possible to triple the number of steps nd do so. Shll we do this, or double nd ccept the loss? On the verge, tripling does fctor 3 of unnecessry work, since the right number of steps for desired ccurcy criterion my in fct fll nywhere in the logrithmic intervl implied by tripling. For doubling, the fctor is only 2, but we lose n extr fctor of 2 in being unble to use ll the previous evlutions. Since < , it is better to triple. Here is the resulting routine, which is directly comprble to trpzd. SUBROUTINE midpnt(func,,b,s,n) REAL,b,s,func EXTERNAL func This routine computes the nth stge of refinement of n extended midpoint rule. func is input s the nme of the function to be integrted between limits nd b, lso input. When clled with n=1, the routine returns s s the crudest estimte of f(x)dx. Subsequent clls with n=2,3,... (in tht sequentil order) will improve the ccurcy of s by dding (2/3) 3 n-1 dditionl interior points. s should not be modified between sequentil clls. REAL ddel,del,sum,tnm,x s=(b-)*func(0.5*(+b)) else it=3**(n-2) tnm=it del=(b-)/(3.*tnm) ddel=del+del The dded points lternte in spcing between del nd ddel. x=+0.5*del sum=0. do 11 j=1,it x=x+ddel x=x+del enddo 11 s=(s+(b-)*sum/tnm)/3. endif return END The new sum is combined with the old integrl to give refined integrl. The routine midpnt cn exctly replce trpzd in driver routine like qtrp ( 4.2); one simply chnges cll trpzd to cll midpnt, nd perhps lso decreses the prmeter JMAX since 3 JMAX 1 (from step tripling) is much lrger number thn 2 JMAX 1 (step doubling). The open formul implementtion nlogous to Simpson s rule (qsimp in 4.2) substitutes midpnt for trpzd nd decreses JMAX s bove, but now lso chnges the extrpoltion step to be Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ). s=(9.*st-ost)/8.
3 4.4 Improper Integrls 137 since, when the number of steps is tripled, the error decreses to 1/9th its size, not 1/4th s with step doubling. Either the modified qtrp or the modified qsimp will fix the first problem on the list t the beginning of this section. Yet more sophisticted is to generlize Romberg integrtion in like mnner: C SUBROUTINE qromo(func,,b,ss,choose) INTEGER JMAX,JMAXP,K,KM REAL,b,func,ss,EPS EXTERNAL func,choose PARAMETER (EPS=1.e-6, JMAX=14, JMAXP=JMAX+1, K=5, KM=K-1) USES polint Romberg integrtion on n open intervl. Returns s ss the integrl of the function func from to b, using ny specified integrting subroutine choose nd Romberg s method. Normlly choose will be n open formul, not evluting the function t the endpoints. It is ssumed tht choose triples the number of steps on ech cll, nd tht its error series contins only even powers of the number of steps. The routines midpnt, midinf, midsql, midsqu, re possible choices for choose. The prmeters hve the sme mening s in qromb. INTEGER j REAL dss,h(jmaxp),s(jmaxp) h(1)=1. do 11 j=1,jmax cll choose(func,,b,s(j),j) if (j.ge.k) then cll polint(h(j-km),s(j-km),k,0.,ss,dss) if (bs(dss).le.eps*bs(ss)) return endif s(j+1)=s(j) h(j+1)=h(j)/9. This is where the ssumption of step tripling nd n even enddo 11 error series is used. puse too mny steps in qromo END The differences between qromo nd qromb ( 4.3) re so slighttht it is perhps grtuitous to list qromo in full. It, however, is n excellent driver routine for solving ll the other problems of improper integrls in our first list (except the intrctble fifth), s we shll now see. The bsic trick for improper integrls is to mke chnge of vribles to eliminte the singulrity, or to mp n infinite rnge of integrtion to finite one. For exmple, the identity 1/ 1/b ( ) 1 1 t 2 f dt b > 0 (4.4.2) t cn be used with either b nd positive, or with nd b negtive, nd works for ny function which decreses towrds infinity fster thn 1/x 2. You cn mke the chnge of vrible implied by (4.4.2) either nlyticlly nd then use (e.g.) qromo nd midpnt to do the numericl evlution, or you cn let the numericl lgorithm mke the chnge of vrible for you. We prefer the ltter method s being more trnsprent to the user. To implement eqution (4.4.2) we simply write modified version of midpnt, clled midinf, which llows b to be infinite (or, more precisely, very lrge number on your prticulr mchine, such s ), or to be negtive nd infinite. Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ).
4 138 Chpter 4. Integrtion of Functions SUBROUTINE midinf(funk,,bb,s,n) REAL,bb,s,funk This routine is n exct replcement for midpnt, i.e., returns s s the nth stge of refinement of the integrl of funk from to bb, except tht the function is evluted t evenly spced points in 1/x rther thn in x. This llows the upper limit bb to be s lrge nd positive s the computer llows, or the lower limit to be s lrge nd negtive, but not both. nd bb must hve the sme sign. REAL,b,ddel,del,sum,tnm,func,x func(x)=funk(1./x)/x**2 This sttement function effects the chnge of vrible. b=1./ These two sttements chnge the limits of integrtion ccordingly. =1./bb From this point on, the routine is exctly identicl to midpnt. s=(b-)*func(0.5*(+b)) else it=3**(n-2) tnm=it del=(b-)/(3.*tnm) ddel=del+del x=+0.5*del sum=0. do 11 j=1,it x=x+ddel x=x+del enddo 11 s=(s+(b-)*sum/tnm)/3. endif return END If you need to integrte from negtive lower limit to positive infinity, you do this by breking the integrl into two pieces t some positive vlue, for exmple, cll qromo(funk,-5.,2.,s1,midpnt) cll qromo(funk,2.,1.e30,s2,midinf) nswer=s1+s2 Where should you choose the brekpoint? At sufficiently lrge positive vlue so tht the function funk is t lest beginning to pproch its symptotic decrese to zero vlue t infinity. The polynomil extrpoltion implicit in the second cll to qromo dels with polynomil in 1/x, not in x. To del with n integrl tht hs n integrble power-lw singulrity t its lower limit, one lso mkes chnge of vrible. If the integrnd diverges s (x ) γ, 0 γ<1, ner x =, use the identity 1 (b ) 1 γ t γ 1 1 γ f(t 1 γ + )dt (b >) (4.4.3) 1 γ 0 Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ). If the singulrity is t the upper limit, use the identity 1 (b ) 1 γ t γ 1 1 γ f(b t 1 γ )dt (b >) (4.4.4) 1 γ 0
5 4.4 Improper Integrls 139 If there is singulrity t both limits, divide the integrl t n interior brekpoint s in the exmple bove. Equtions (4.4.3) nd (4.4.4) re prticulrly simple in the cse of inverse squre-root singulrities, cse tht occurs frequently in prctice: for singulrity t, nd b 0 b 0 2tf( + t 2 )dt (b >) (4.4.5) 2tf(b t 2 )dt (b >) (4.4.6) for singulrity t b. Once gin, we cn implement these chnges of vrible trnsprently to the user by defining substitute routines for midpnt which mke the chnge of vrible utomticlly: SUBROUTINE midsql(funk,,bb,s,n) REAL,bb,s,funk This routine is n exct replcement for midpnt, except tht it llows for n inverse squreroot singulrity in the integrnd t the lower limit. REAL ddel,del,sum,tnm,x,func,,b func(x)=2.*x*funk(+x**2) b=sqrt(bb-) =0. The rest of the routine is exctly like midpnt nd is omitted. Similrly, SUBROUTINE midsqu(funk,,bb,s,n) REAL,bb,s,funk This routine is n exct replcement for midpnt, except tht it llows for n inverse squreroot singulrity in the integrnd t the upper limit bb. REAL ddel,del,sum,tnm,x,func,,b func(x)=2.*x*funk(bb-x**2) b=sqrt(bb-) =0. The rest of the routine is exctly like midpnt nd is omitted. One lst exmple should suffice to show how these formuls re derived in generl. Suppose the upper limit of integrtion is infinite, nd the integrnd flls off exponentilly. Then we wnt chnge of vrible tht mps e x dx into (±)dt (with the sign chosen to keep the upper limit of the new vrible lrger thn the lower limit). Doing the integrtion gives by inspection Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ). t = e x or x = log t (4.4.7)
6 140 Chpter 4. Integrtion of Functions so tht x= x= t=e t=0 f( log t) dt t (4.4.8) The user-trnsprent implementtion would be SUBROUTINE midexp(funk,,bb,s,n) REAL,bb,s,funk This routine is n exct replcement for midpnt, except tht bb is ssumed to be infinite (vlue pssed not ctully used). It is ssumed tht the function funk decreses exponentilly rpidly t infinity. REAL ddel,del,sum,tnm,x,func,,b func(x)=funk(-log(x))/x b=exp(-) =0. The rest of the routine is exctly like midpnt nd is omitted. CITED REFERENCES AND FURTHER READING: Acton, F.S. 1970, Numericl Methods Tht Work; 1990, corrected edition (Wshington: Mthemticl Assocition of Americ), Chpter 4. Dhlquist, G., nd Bjorck, A. 1974, Numericl Methods (Englewood Cliffs, NJ: Prentice-Hll), 7.4.3, p Stoer, J., nd Bulirsch, R. 1980, Introduction to Numericl Anlysis (New York: Springer-Verlg), 3.7, p Gussin Qudrtures nd Orthogonl Polynomils In the formuls of 4.1, the integrl of function ws pproximted by the sum of its functionl vlues t set of eqully spced points, multiplied by certin ptly chosen weighting coefficients. We sw tht s we llowed ourselves more freedom in choosing the coefficients, we could chieve integrtion formuls of higher nd higher order. The ide of Gussin qudrtures is to give ourselves the freedom to choose not only the weighting coefficients, but lso the loction of the bscisss t which the function is to be evluted: They will no longer be eqully spced. Thus, we will hve twice the number of degrees of freedom t our disposl; it will turn out tht we cn chieve Gussin qudrture formuls whose order is, essentilly, twice tht of the Newton-Cotes formul with the sme number of function evlutions. Does this sound too good to be true? Well, in sense it is. The ctch is fmilir one, which cnnot be overemphsized: High order is not the sme s high ccurcy. High order trnsltes to high ccurcy only when the integrnd is very smooth, in the sense of being well-pproximted by polynomil. Smple pge from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN X) visit website or cll (North Americ only),or send emil to trde@cup.cm.c.uk (outside North Americ).
Numerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More 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 informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More 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 informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationConstruction of Gauss Quadrature Rules
Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture
More informationHow can we approximate the area of a region in the plane? What is an interpretation of the area under the graph of a velocity function?
Mth 125 Summry Here re some thoughts I ws hving while considering wht to put on the first midterm. The core of your studying should be the ssigned homework problems: mke sure you relly understnd those
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More information1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More 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 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 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 information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationdifferent methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).
Mth 1A with Professor Stnkov Worksheet, Discussion #41; Wednesdy, 12/6/217 GSI nme: Roy Zho Problems 1. Write the integrl 3 dx s limit of Riemnn sums. Write it using 2 intervls using the 1 x different
More informationMAA 4212 Improper Integrals
Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More 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 informationSuppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.
Mth 43 Section 6. Section 6.: Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationSections 5.2: The Definite Integral
Sections 5.2: The Definite Integrl In this section we shll formlize the ides from the lst section to functions in generl. We strt with forml definition.. The Definite Integrl Definition.. Suppose f(x)
More informationLecture 23: Interpolatory Quadrature
Lecture 3: Interpoltory Qudrture. Qudrture. The computtion of continuous lest squres pproximtions to f C[, b] required evlutions of the inner product f, φ j = fxφ jx dx, where φ j is polynomil bsis function
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
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 information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More 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 informationp(t) dt + i 1 re it ireit dt =
Note: This mteril is contined in Kreyszig, Chpter 13. Complex integrtion We will define integrls of complex functions long curves in C. (This is bit similr to [relvlued] line integrls P dx + Q dy in R2.)
More informationCOT4501 Spring Homework VII
COT451 Spring 1 Homework VII The ssignment is due in clss on Thursdy, April 19, 1. There re five regulr problems nd one computer problem (using MATLAB). For written problems, you need to show your work
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
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 information4181H Problem Set 11 Selected Solutions. Chapter 19. n(log x) n 1 1 x x dx,
48H Problem Set Selected Solutions Chpter 9 # () Tke f(x) = x n, g (x) = e x, nd use integrtion by prts; this gives reduction formul: x n e x dx = x n e x n x n e x dx. (b) Tke f(x) = (log x) n, g (x)
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationThe goal of this section is to learn how to use a computer to approximate definite integrals, i.e. expressions of the form. Z b
Lecture notes for Numericl Anlysis Integrtion Topics:. Problem sttement nd motivtion 2. First pproches: Riemnn sums 3. A slightly more dvnced pproch: the Trpezoid rule 4. Tylor series (the most importnt
More informationUNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3
UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,
More informationDefinite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +
Definite Integrls --5 The re under curve cn e pproximted y dding up the res of rectngles. Exmple. Approximte the re under y = from x = to x = using equl suintervls nd + x evluting the function t the left-hnd
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More 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 informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationNumerical Integration
Numericl Integrtion Wouter J. Den Hn London School of Economics c 2011 by Wouter J. Den Hn June 3, 2011 Qudrture techniques I = f (x)dx n n w i f (x i ) = w i f i i=1 i=1 Nodes: x i Weights: w i Qudrture
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationTaylor Polynomial Inequalities
Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil
More informationINTRODUCTION TO INTEGRATION
INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.
Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More 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 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 informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationx = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b
CHAPTER 5. INTEGRALS 61 where nd x = b n x i = 1 (x i 1 + x i ) = midpoint of [x i 1, x i ]. Problem 168 (Exercise 1, pge 377). Use the Midpoint Rule with the n = 4 to pproximte 5 1 x e x dx. Some quick
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More information6.5 Numerical Approximations of Definite Integrals
Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 6.5 Numericl Approximtions of Definite Integrls Sometimes the integrl of function cnnot be expressed with elementry functions, i.e., polynomil,
More informationNumerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b
Numericl Integrtion Newton Cotes Formuls Given function f : R R nd two rel numbers, b R, < b, we clculte (pproximtely) the integrl I(f,, b) = f (x) dx K. Frischmuth (IfM UR) Numerics for CSE 08/09 8 /
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More informationRiemann is the Mann! (But Lebesgue may besgue to differ.)
Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >
More informationModule 11.4: nag quad util Numerical Integration Utilities. Contents
Qudrture Module Contents Module 11.4: ng qud util Numericl Integrtion Utilities ng qud util provides utility procedures for computtion involving integrtion of functions, e.g., the computtion of the weights
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationNumerical integration. Quentin Louveaux (ULiège - Institut Montefiore) Numerical analysis / 10
Numericl integrtion Quentin Louveux (ULiège Institut Montefiore) Numericl nlysis 2018 1 / 10 Numericl integrtion We consider definite integrls Z b f (x)dx better to it use if known! A We do not ssume tht
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationChapter 5. Numerical Integration
Chpter 5. Numericl Integrtion These re just summries of the lecture notes, nd few detils re included. Most of wht we include here is to be found in more detil in Anton. 5. Remrk. There re two topics with
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More 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 informationTrapezoidal Rule, n = 1, x 0 = a, x 1 = b, h = b a. f (x)dx = h 2 (f (x 0) + f (x 1 )) h3
Trpezoidl Rule, n = 1, x 0 =, x 1 = b, h = b f (x)dx = h 2 (f (x 0) + f (x 1 )) h3 12 f (ξ). Simpson s Rule: n = 3, x 0 =, x 1 = +b 2, x 2 = b, h = b 2. Qudrture Rule, double node t x 1 P 3 (x)dx = f (x
More information1 The Riemann Integral
The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re
More informationAP Calculus Multiple Choice: BC Edition Solutions
AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationName Solutions to Test 3 November 8, 2017
Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
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 informationLecture 12: Numerical Quadrature
Lecture 12: Numericl Qudrture J.K. Ryn@tudelft.nl WI3097TU Delft Institute of Applied Mthemtics Delft University of Technology 5 December 2012 () Numericl Qudrture 5 December 2012 1 / 46 Outline 1 Review
More informationf(a+h) f(a) x a h 0. This is the rate at which
M408S Concept Inventory smple nswers These questions re open-ended, nd re intended to cover the min topics tht we lerned in M408S. These re not crnk-out-n-nswer problems! (There re plenty of those in the
More informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
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 informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More information