The goal of this section is to learn how to use a computer to approximate definite integrals, i.e. expressions of the form. Z b
|
|
- Geraldine Hicks
- 5 years ago
- Views:
Transcription
1 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 technique in numericl nlysis) 5. Error nlysis of Trpezoid rule 6. Rtes of convergence 7. A look further field: Newton-Cotes formuls 8. Notes on reding nd further topics. Problem sttement nd motivtion The gol of this section is to lern how to use computer to pproximte definite integrls, i.e. expressions of the form f(x)dx () for finite rel numbers nd b nd smooth rel-vlued functions f. It is worth tking moment to reflect on the question of why we might wnt to use computer to compute definite integrls. After ll, in first yer clculus, we lern how to solve definite integrl using nti-derivtives. Wht is wrong with the nti-derivtive pproch? And in rel life, why would we ctully wnt to solve definite integrls nywy? One nswer to the question of wht s wrong with the nti-derivtive pproch is tht not ll functions hve nti-derivtives. For exmple, the clssic Bell curve hs the functionl form f(x) = p 2 2 e (x µ)2 /2 2 (2) for constnts µ nd. Bell curves cn be used to describe mny kinds of dt, nd show up ll the time in the nturl sciences, but the ntiderivtives of these functions hve no nlyticlly trctble form. Even when integrnds do hve nice nti-derivtes, finding them cn be lborious nd time consuming, wheres numericl methods tend to be esy, quick, nd ccurte. Thus nother nswer to the wht s wrong with the nti-derivtive pproch question is tht it cn ctully be much more cumbersome thn numericl pproch. The question of why we would wnt to solve definite integrls nywy hs myrid nswers. Here re few exmples. Exmple. (Probbility). A density function is positive functions f(x) : R! [0, ) such tht Z f(x)dx =. A rndom vrible X is quntity whose vlue is determined by some rndom process (e.g. X could be the mount of money gined or lost in tomorrow s stock mrket.) A rndom vrible X is sid to be described by density function f if, for ny intervl [, b], the probbility tht X lies in [, b] is given by P ( pple X pple b) = f(x)dx. In other words, probbilities re computed by solving definite integrls. Applictions re ubiquitous, nd include engineering (e.g. the probbility of system filure for geosttionry stellites), biology (e.g. the probbility of genetic survivl), chemistry (the probbility of bond forming), nd mny other fields.
2 .2 First pproches: Riemnn sums Lecture notes for Numericl Anlysis Exmple.2 (Bridge design). Bridges often employ cbles hung in rcs from tll supports. Clculting the length of the cble is importnt for cost estimtes. While function f(x) might be produced tht describes the position of the cble, the length of the cble needs to be ttined by computing n rc-length integrl..2 First pproches: Riemnn sums There re mny wys to evlute definite integrls. Before lerning bout these techniques, it is useful to see wht you cn come up with yourself. Clssroom Activity: Together with prtner, brinstorm bout how you might go bout pproximting definite integrl. Assume tht the tools t your disposl re computer the cpcity to evlute f(x) t ny x 2 [, b]. It s likely tht you stumbled cross the ide of using Riemnn sum to estimte the integrl. Recll tht one wy to pproximte the definite integrl is vi the right hnd endpoint rule. Definition.. The right hnd endpoint pproximtion of R b f(x)dx is f(x)dx h f(x i ), x i = + ih, h =(b )/n. (3) i= For fixed n, the right hnd side is n pproximtion to the left. When n =, the pproximtion hs the form (b )f(b). This pproximtion should get better s n gets lrger, nd indeed converge s n!. A few nturl questions emerge from this nive pproch, however: how lrge should n be in order to chieve given ccurcy? how should tht ccurcy be mesured? how long will the process tke? cn n be mde rbitrrily lrge, or is there threshold beyond which my mchine will blk? is there better pproch? We will ddress these nd other questions in wht follows..3 The Trpezoid Rule One thing to note bout formul (3) is tht its clcultionl requirements re s follows: n function evlutions n dditions multipliction Clssroom Activity: Are there other wys to implement this lgorithm? If so, re they more or less e cient? Discuss. Inccurcies in this formul might rise from severl plces. 2
3 .4 Tylor Series Lecture notes for Numericl Anlysis Clssroom Activity: Where might these inccurcies come from? Discuss. Ides: roundo error in representing numbers on the computer mismtch between f(x) nd the pproximtion f(x i ). One of the problems with the right hnd endpoint rule is tht it tries to cpture the verge vlue of f(x) on the intervl [x i,x i+ ] vi smple tht is chosen to lie t n extreme endpoint of this intervl. A more resonble pproch is to figure out wht f(x) is t the left hnd endpoint nd the right hnd endpoint, nd form the verge. This is the content of the Trpezoid rule. The single trpezoid pproximtion to R b f(x)dx is (b ) f()+f(b), 2 which essentilly tries to pproximte f(x) on [, b] vitheverge of the endpoints. If we divide [, b] into n subintervls nd pply the single trpezoid rule to ech, we get the following: Definition.2. The trpezoid rule pproximtes I = R b f(x)dx vi the formul Z " # b f(x)dx b f()+2 f(x i )+f(b), x i = + ih, h =(b )/n. (4) 2 i= Clssroom Activity: Prove this formul. The trpezoid rule seems more promising thn the right-hnd endpoint rule, in tht we ve mde more intelligent choice for our pproximtion of f(x) in the intervl [x i,x i+ ]. The computtionl requirements re bout the sme: n + function evlutions n + dditions 3 multiplictions In order to figure out if the trpezoid rule relly bets the right hnd endpoint rule, however, we need to do some hrd mthemticl nlysis. Our tool for this nlysis is clled Tylor series..4 Tylor Series The bsic ide behind Tylor series is tht it is possible to pproximte n rbitrry function f(x) with polynomil, s long s f(x) is smooth (see below.) The pproximtion tends to get better s the degree of the pproximting polynomil gets lrger, nd there re forml bounds on how lrge the pointwise error cn be. To stte Tylor s Theorem, let s first define wht we men by smooth function. Definition.3. A function is sid to be of clss C f(x) :! R if f(x) nd ll of its derivtives re defined on the domin. Note: functions in C re relly nice you cn di erentite them s much s you wnt, nd you never run into problems. In pplictions, is generlly tken to be some finite intervl [, b], but it doesn t hve to be. In stting Tylor s theorem (below), I tke = R, just to simplify things. Tylor s theorem cn be extended to other clsses of functions, but this one is good plce to strt. 3
4 .5 Error nlysis of Trpezoid Rule Lecture notes for Numericl Anlysis Theorem.. Suppose f(x) 2CR. Then for ny point 2 R, the vlue of f t the point x cn be expressed s f(x) =f()+f 0 ()(x )+f 00 (x )2 () + f 000 (x )3 () +. 2! 3! If T n (x) represents the first n + terms of the right hnd side nd R n+ (x) represents the reminder, then we hve f(x) =T n (x)+r n+ (x), where the reminder is of the form R n+ (x) =f (n+) ( ) (x )n+ n! for some between nd x. Nottion: Sometimes will write T n (f,,x) nd R n+ (f,,x) if we wish to emphsize the dependence of these terms on f nd. Clssroom Activity: Find T 2 (x) iff(x) =sin(x) nd =. Then find n expression for R 3 (x). Finlly, suppose you use T 2 (x) to estimte sin(5 /4). Clculte bound on the error..5 Error nlysis of Trpezoid Rule At first blush, it s not totlly obvious why Tylor s theorem is so importnt. One wy to think bout it is tht it llows us to use simple polynomil pproximtions to clculte how big errors will be. Those pproximtions will hve their own errors, of course, but the explicit form of the reminder term in the Tylor Theorem llows us to get bounds on these errors. As n exmple of the Tylor Series in ction, we cn clculte bound on the errors involved in the Trpezoid method. Theorem.2. Suppose f(x) 2 C 2 [,b]. Suppose we use the single intervl trpezoid pproximtion to pproximte Then the error of the pproximtion stisfies where Proof. See clss notes for detils. Ĩ =(b I = ) f()+f(b) 2 f(x)dx. I Ĩ pple(b ) 3 M, M = sup f 00 (x) x2[,b] Corollry.. Suppose f(x) 2 C 2 [,b], nd we use the trpezoid pproximtion with n intervls " # n (b ) X Ĩ = f()+2 f(x i )+f(b) n to pproximte I = i f(x)dx. Then the error of the pproximtion stisfies 3 b I Ĩ pple Mn = h 2 M(b n where M is s bove nd h =(b ), )/n is the spcing between smple points. 4
5 .6 Convergence rtes Lecture notes for Numericl Anlysis Proof. Apply Theorem (.2) to ech subintervl, sum the errors, nd simplify. In other words, we hve used Tylor series to show tht if we use the Trpezoid rule on smooth enough function, the error decreses in proportion to the squre of the spce between smple points. This is good news: it mens, mong other things, tht if we decrese tht smple spcing, we decrese the error. Not only tht, the fct tht the error decreses in proportion to the squre of this spcing mens we don t hve to decrese it very much to see big improvement. For exmple, if we double the number of smple points (nd thus decrese h by /2), we will decrese the error by fctor of four (i.e. (/2) 2.) It will turn out tht other lgorithms yield other error formuls tht involve di erent powers of h. To tlk cogently bout how these vrious lgorithms stck up to one nother, we need to define wht we men by convergence rte..6 Convergence rtes Definition.4. Suppose n is sequence converging to, nd if there exists positive k nd n integer N such tht n > 0 is sequence converging to 0. Then n pplek n 8n>N, we sy n converges to with rte of convergence O( n ) (sy big-oh n ), nd write n = + O( n ). If it hppens tht we cn choose sequence K n > 0 such tht K n! 0 nd n pplek n n, we sy n! little-oh n nd write n = + o( n ). Note: The definition is frmed for generl n, but in relity the cse of interest is when n =/n p for some p>0. In generl, we re interested in knowing wht the lrgest p is for which the convergence is O(/n p ). Exmple.3. Let n =(n + )/n 2. Set It follows tht n converges to 0 with O(/n). n =/n. Then n 0 = n + n 2 pple 2n n 2 = 2 n. The sequence definition of convergence rte cn esily be extended to functions: Definition.5. Suppose If there exists K>0 nd lim G(h) =0, lim F (h) =L. h!0 h!0 >0 such tht F (h) L pplek G(h) 8h <, then we write F (h) =L + O(G(h)). Once gin, the cse of interest is when G(h) =h p,p>0. Exmple.4. Let F (h) =( cos h)/h. The second order Tylor pproximtion of cos h ner =0is h 2 /2+R 3 (h), nd by the Tylor reminder theorem, R 3 is O(h 3 ).Itfollowstht cos h h = ( h2 /2+O(h 3 )) h = h + O(h 2 )=h( + O(h)) Since nything O(h) will eventully decrese below s h! 0, we cn set K =2nd tke G(h) =h. Then F (h) pple KG(h) for ll su ciently smll h, sof (h)! 0 O(h). 5
6 .7 A look further field: Newton-Cotes formuls Lecture notes for Numericl Anlysis.7 A look further field: Newton-Cotes formuls As we hve seen, the right hnd endpoint rule tries to pproximte f(x) in n intervl vi horizontl line tht hppens to mtch f(x) t one point in the intervl, while the trpezoid rule tries to mtch f(x) in n intervl with stright line tht mtches f(x) t two points in the intervl. This suggests n ide: if we choose n + points in n intervl, we cn fit n n-degree polynomil to them, nd use tht polynomil s our pproximtion to f(x) over the intervl. Bsed on the Tylor Reminder theorem, we expect the pproximtion to improve s the degree of the pproximting polynomil increses. Formuls derived this wy re clled Newton-Cotes formuls when the n + smple points re uniformly spced. There re two kinds of Newton-Cotes, closed nd open. The closed formuls include the endpoints of the intervl s smple points, the open formuls do not. But in essence, the wy to develop Newton-Cotes formul is s follows: Algorithm How to develop Newton-Cotes formul. Choose n + smple points in [, b]. For closed formuls, these points should be x i, i =0,,n, x i = + ih, h =(b )/n nd for open formuls they should be x i, i =0,,n, x i = +(i + )h, h =(b )/(n + 2) 2. Find the unique n degree polynomil p n (x) tht interpoltes f t the x i, i.e. for which p n (x i )=f(x i ) i =0,,n. 3. Define the (open or closed) n + point Newton-Cotes pproximtion s f(x)dx p n (x)dx It turns out there is clever wy to find the interpolting polynomils p n (x). Definition.6. Given n + distinct smple points x i,theith Lgrnge interpolting polynomil L i,n is the unique nth order polynomil tht stisfies xk = x L i,n (x k )= i, i =0,,n 0 x k 6= x i for ny member x k of the smple points. Theorem.3. Given n + distinct smple points x i nd function f(x), thepolynomilp n (x) tht interpoltes the points (x i,f(x i )) is given by p n (x) = f(x i )L i,n (x). (5) i=0 Proof. Since the right hnd side of (5) is the sum of degree n polynomils, it itself must be degree n polynomil. Moreover, t ech smple point x k, ll terms in the sum dispper except for the kth one, which evlutes to f(x k ). QED. 6
7 .7 A look further field: Newton-Cotes formuls Lecture notes for Numericl Anlysis It turns out tht is ctully relly esy to come up with the Lgrnge polynomils L i,n (x). Here s formul: L i (x) = ny k=0,k6=i (x x k ) (x i x k ). (6) Clssroom Activity: Discuss why this formul works. Clssroom Activity: Find L 0,2 (x), L,2 (x), nd L 2,2 (x) for smple points, 0, nd. Finlly, we cn put ll this together to get ctul, implementble formuls for Newton-Cotes integrtion. Specificlly, substituting formul (6) into (5) nd the result into (3), we get where f(x)dx p n (x)dx c i = = = = f(x i )L i,n (x)dx i=0 f(x i ) i=0 c i f(x i ) i=0 L i,n (x)dx. L i,n (x)dx Note tht the generl Newton-Cotes formul is thus just weighted sum of smples of f(x). To get the weight coe cients, we solve for the L i,n (x) nd integrte them from to b. Clssroom Activity: Using your results from the lst clss ctivity, find the closed Newton-Cotes formul for n = 2 on the intervl [, ]. It is worth noting tht if we get the Newton-Cotes formuls for one intervl, we get them for ll intervls. Clssroom Activity: Use the fct tht the mp (x) x b mps n intervl [, b] into[0, ] to rgue tht if c i re the Newton-Cotes coe cients on [0, ], then (b )c i re the Newton-Cotes coe cients on [, b]. Hint: you ll need u subsitution. Of course, in prctice, we don t spend lot of time solving for Newton Cotes formuls becuse someone else hs lredy done it. Here is list of some common ones: Closed Newton-Cotes formuls with error terms (Burden nd Fires, pg. 99.) n =, Trpezoid: Z x x 0 f(x)dx = h 2 [f(x 0)+f(x )] h 3 2 f 00 ( ). 7
8 .8 Notes on reding nd further topics Lecture notes for Numericl Anlysis n = 2, Simpson: Z x2 x 0 f(x)dx = h 3 [f(x 0)+4f(x )+f(x 2 )] h 5 90 f (4) ( ). n = 3, Simpson s 3/8ths Rule: Z x3 x 0 f(x)dx = 3h 8 [f(x 0)+3f(x )+3f(x 2 )+f(x 3 )] 3h 5 80 f (4) ( ). n = 4: Z x4 x 0 f(x)dx = 2h 45 [7f(x 0) + 32f(x ) + 2f(x 2 ) + 32f(x 3 )+7f(x 4 )] 8h f (6) ( ). Open Newton-Cotes formuls with error terms (Burden nd Fires, pg. 99.) n = 0, Midpoint Rule: n = : n = 3: Z x3 Z x2 x Z x x f(x)dx =2hf(x 0 )+ h3 3 f 00 ( ). f(x)dx = 3h 2 [f(x 0)+f(x )] + 3h3 4 f 00 ( ). x f(x)dx = 4h 3 [2f(x 0) f(x )+2f(x 2 )] + 45h5 45 f (4) ( ). n = 4: Z x4 x f(x)dx = 5h 24 [f(x 0)+f(x )+f(x 2 ) + f(x 3 )] + 95h5 44 f (4) ( )..8 Notes on reding nd further topics The mteril in this section ws tken primrily from three sources:. our textbook 2. Numericl Anlysis, by Burden nd Fires 3. Numericl Anlysis, by Cheney nd Kincid Numericl integrtion is big topic, nd there re lots of innovtions tht we hven t covered. Here is very prtil list of things you might explore for finl project:. Gussin qudrture 2. Adptive qudrture 3. Romber integrtion 4. Extensions to multiple dimensions 5. How to hndle improper integrls 8
The Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More 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 informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
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 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 informationNumerical 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 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 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 informationLecture 1: Introduction to integration theory and bounded variation
Lecture 1: Introduction to integrtion theory nd bounded vrition Wht is this course bout? Integrtion theory. The first question you might hve is why there is nything you need to lern bout integrtion. You
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 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 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 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 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 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 informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More 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 informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
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 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 informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More 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 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 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 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 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 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 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 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 informationIntegrals - Motivation
Integrls - Motivtion When we looked t function s rte of chnge If f(x) is liner, the nswer is esy slope If f(x) is non-liner, we hd to work hrd limits derivtive A relted question is the re under f(x) (but
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 informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
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 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 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 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 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 informationMath 131. Numerical Integration Larson Section 4.6
Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n
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 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 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 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 informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More 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 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 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 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 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 information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
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 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 informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationInterpreting Integrals and the Fundamental Theorem
Interpreting Integrls nd the Fundmentl Theorem Tody, we go further in interpreting the mening of the definite integrl. Using Units to Aid Interprettion We lredy know tht if f(t) is the rte of chnge of
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 informationSYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus
SYDE 112, LECTURES & 4: The Fundmentl Theorem of Clculus So fr we hve introduced two new concepts in this course: ntidifferentition nd Riemnn sums. It turns out tht these quntities re relted, but it is
More informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
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 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 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 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 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 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 informationMA 124 January 18, Derivatives are. Integrals are.
MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationf(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all
3 Definite Integrl 3.1 Introduction In school one comes cross the definition of the integrl of rel vlued function defined on closed nd bounded intervl [, b] between the limits nd b, i.e., f(x)dx s the
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 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 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 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 informationCalculus I-II Review Sheet
Clculus I-II Review Sheet 1 Definitions 1.1 Functions A function is f is incresing on n intervl if x y implies f(x) f(y), nd decresing if x y implies f(x) f(y). It is clled monotonic if it is either incresing
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 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 information7.2 Riemann Integrable Functions
7.2 Riemnn Integrble Functions Theorem 1. If f : [, b] R is step function, then f R[, b]. Theorem 2. If f : [, b] R is continuous on [, b], then f R[, b]. Theorem 3. If f : [, b] R is bounded nd continuous
More information. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. cos(2θ) = sin(2θ) =.
Review of some needed Trig Identities for Integrtion Your nswers should be n ngle in RADIANS rccos( 1 2 ) = rccos( - 1 2 ) = rcsin( 1 2 ) = rcsin( - 1 2 ) = Cn you do similr problems? Review of Bsic Concepts
More informationLECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and
LECTURE 9 Numericl Integrtion Recll from Clculus I tht denite integrl is generlly thought of s representing the re under the grph of fèxè between the points x = nd x = b, even though this is ctully only
More informationThe Evaluation Theorem
These notes closely follow the presenttion of the mteril given in Jmes Stewrt s textook Clculus, Concepts nd Contexts (2nd edition) These notes re intended primrily for in-clss presenttion nd should not
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 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 information38 Riemann sums and existence of the definite integral.
38 Riemnn sums nd existence of the definite integrl. In the clcultion of the re of the region X bounded by the grph of g(x) = x 2, the x-xis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationLecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),
More 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 informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More 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 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 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 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 2-Dimensionl Gussin Qudrture 20 4
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 informationSummer MTH142 College Calculus 2. Section J. Lecture Notes. Yin Su University at Buffalo
Summer 6 MTH4 College Clculus Section J Lecture Notes Yin Su University t Bufflo yinsu@bufflo.edu Contents Bsic techniques of integrtion 3. Antiderivtive nd indefinite integrls..............................................
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More 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 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 informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationBig idea in Calculus: approximation
Big ide in Clculus: pproximtion Derivtive: f (x) = df dx f f(x +h) f(x) =, x h rte of chnge is pproximtely the rtio of chnges in the function vlue nd in the vrible in very short time Liner pproximtion:
More informationMath 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8
Mth 3 Fll 0 The scope of the finl exm will include: Finl Exm Review. Integrls Chpter 5 including sections 5. 5.7, 5.0. Applictions of Integrtion Chpter 6 including sections 6. 6.5 nd section 6.8 3. Infinite
More information. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. sin 2 (θ) =
Review of some needed Trig. Identities for Integrtion. Your nswers should be n ngle in RADIANS. rccos( 1 ) = π rccos( - 1 ) = 2π 2 3 2 3 rcsin( 1 ) = π rcsin( - 1 ) = -π 2 6 2 6 Cn you do similr problems?
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More 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 informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
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 informationAN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
Applied Mthemtics E-Notes, 5(005), 53-60 c ISSN 1607-510 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
More information