Numerical Integration

Size: px
Start display at page:

Download "Numerical Integration"

Transcription

1 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 pproximtions to the integrl of the function. An importnt use of both types of methods is estimtion of derivtives nd integrls for functions tht re only known t isolted points, s is the cse with for exmple mesurement dt. An importnt difference between differentition nd integrtion is tht for most functions it is not possible to determine the integrl vi symbolic methods, but we cn still compute numericl pproximtions to virtully ny definite integrl. Numericl integrtion methods re therefore more useful thn numericl differentition methods, nd re essentil in mny prcticl situtions. We use the sme generl strtegy for deriving numericl integrtion methods s we did for numericl differentition methods: We find the polynomil tht interpoltes the function t some suitble points, nd use the integrl of the polynomil s n pproximtion to the function. This mens tht the trunction error cn be nlysed in bsiclly the sme wy s for numericl differentition. However, when it comes to round-off error, integrtion behves differently from differentition: Numericl integrtion is very insensitive to round-off errors, so we will ignore round-off in our nlysis. The mthemticl definition of the integrl is bsiclly vi numericl integrtion method, nd we therefore strt by reviewing this definition. We then derive the simplest numericl integrtion method, nd see how its error cn be nlysed. We then derive two other methods tht re more ccurte, but for these we just indicte how the error nlysis cn be done. We emphsise tht the generl procedure for deriving both numericl dif- 79

2 Figure 1.1. The re under the grph of function. ferentition nd integrtion methods with error nlyses is the sme with the exception tht round-off errors re not of must interest for the integrtion methods. 1.1 Generl bckground on integrtion Recll tht if f (x) is function, then the integrl of f from x = to x = b is written f (x)d x. The integrl gives the re under the grph of f, with the re under the positive prt counting s positive re, nd the re under the negtive prt of f counting s negtive re, see figure 1.1. Before we continue, we need to define term which we will use repetedly in our description of integrtion. Definition 1.1 (Prtition). Let nd b be two rel numbers with < b. A prtition of [,b] is finite sequence {x i } n of incresing numbers in [,b] i= with x = nd x n = b, = x < x 1 < x < x n 1 < x n = b. The prtition is sid to be uniform if there is fixed number h, clled the step length, such tht x i x i 1 = h = (b )/n for i = 1,..., n. 8

3 () (b) (c) (d) Figure 1.. The definition of the integrl vi inscribed nd circumsribed step functions. The trditionl definition of the integrl is bsed on numericl pproximtion to the re. We pick prtition {x i } n of [,b], nd in ech subintervl i= [x i 1, x i ] we determine the mximum nd minimum of f (for convenience we ssume tht these vlues exist), m i = min f (x), M i = mx f (x), x [x i 1,x i ] x [x i 1,x i ] for i = 1,,..., n. We cn then compute two obvious pproximtions to the integrl by pproximting f by two different functions which re both ssumed to be constnt on ech intervl [x i 1, x i ]: The first hs the constnt vlue m i nd the other the vlue M i. We then sum up the res under ech of the two step functions n end up with the two pproximtions I = m i (x i x i 1 ), I = M i (x i x i 1 ), (1.1) to the totl re. In generl, the first of these is too smll, the other too lrge. To define the integrl, we consider lrger prtitions (smller step lengths) nd consider the limits of I nd I s the distnce between neighbouring x i s goes 81

4 to zero. If those limits re the sme, we sy tht f is integrble, nd the integrl is given by this limit. Definition 1. (Integrl). Let f be function defined on the intervl [, b], nd let {x i } n i= be prtition of [,b]. Let m i nd M i denote the minimum nd mximum vlues of f over the intervl [x i 1, x i ], respectively, ssuming they exist. Consider the two numbers I nd I defined in (1.1). If sup I nd inf I both exist nd re equl, where the sup nd inf re tken over ll possible prtitions of [,b], the function f is sid to be integrble, nd the integrl of f over [,b] is defined by I = f (x)d x = sup I = inf I. This process is illustrted in figure 1. where we see how the piecewise constnt pproximtions become better when the rectngles become nrrower. The bove definition cn be used s numericl method for computing pproximtions to the integrl. We choose to work with either mxim or minim, select prtition of [,b] s in figure 1., nd dd together the res of the rectngles. The problem with this technique is tht it cn be both difficult nd time consuming to determine the mxim or minim, even on computer. However, it cn be shown tht the integrl hs property tht is very useful when it comes to numericl computtion. Theorem 1.3. Suppose tht f is integrble on the intervl [,b], let {x i } n i= be prtition of [,b], nd let t i be number in [x i 1, x i ] for i = 1,..., n. Then the sum Ĩ = f (t i )(x i x i 1 ) (1.) will converge to the integrl when the distnce between ll neighbouring x i s tends to zero. Theorem 1.3 llows us to construct prcticl, numericl methods for computing the integrl. We pick prtition of [,b], choose t i equl to x i 1 or x i, nd compute the sum (1.). It turns out tht n even better choice is the more symmetric t i = (x i + x i 1 )/ which leds to the pproximtion I f ( (x i + x i 1 )/ ) (x i x i 1 ). (1.3) 8

5 This is the so-clled midpoint rule which we will study in the next section. In generl, we cn derive numericl integrtion methods by splitting the intervl [, b] into smll subintervls, pproximte f by polynomil on ech subintervl, integrte this polynomil rther thn f, nd then dd together the contributions from ech subintervl. This is the strtegy we will follow for deriving more dvnced numericl integrtion methods, nd this works s long s f cn be pproximted well by polynomils on ech subintervl. Exercises 1 In this exercise we re going to study the definition of the integrl for the function f (x) = e x on the intervl [,1]. ) Determine lower nd upper sums for uniform prtition consisting of 1 subintervls. b) Determine the bsolute nd reltive errors of the sums in () compred to the exct vlue e 1 = of the integrl. c) Write progrm for clculting the lower nd upper sums in this exmple. How mny subintervls re needed to chieve n bsolute error less thn 3 1 3? 1. The midpoint rule for numericl integrtion We hve lredy introduced the midpoint rule (1.3) for numericl integrtion. In our stndrd frmework for numericl methods bsed on polynomil pproximtion, we cn consider this s using constnt pproximtion to the function f on ech subintervl. Note tht in the following we will lwys ssume the prtition to be uniform. Algorithm 1.4. Let f be function which is integrble on the intervl [,b], nd let {x i } n be uniform prtition of [,b]. In the midpoint rule, the integrl of f is pproximted i= by f (x)d x I mid (h) = h f (x i 1/ ), (1.4) where x i 1/ = (x i 1 + x i )/ = + (i 1/)h. This my seem like strngely formulted lgorithm, but ll there is to do is to compute the sum on the right in (1.4). The method is illustrted in figure 1.3 in the cses where we hve 1 nd 5 subintervls. 83

6 () (b) Figure 1.3. The midpoint rule with one subintervl () nd five subintervls (b) A detiled lgorithm Algorithm 1.4 describes the midpoint rule, but lcks lot of detil. In this section we give more detiled lgorithm. Whenever we compute quntity numericlly, we should try nd estimte the error, otherwise we hve no ide of the qulity of our computtion. We did this when we discussed lgorithms for finding roots of equtions in chpter 1, nd we cn do exctly the sme here: We compute the integrl for decresing step lengths, nd stop the computtions when the difference between two successive pproximtions is less thn the tolernce. More precisely, we choose n initil step length h nd compute the pproximtions I mid (h ), I mid (h 1 ),..., I mid (h k ),..., where h k = h / k. Suppose I mid (h k ) is our ltest pproximtion. Then we estimte the reltive error by the number I mid (h k ) I mid (h k 1 ), I mid (h k ) nd stop the computtions if this is smller thn ɛ. To void potentil division by zero, we use the test I mid (h k ) I mid (h k 1 ) ɛ I mid (h k ). As lwys, we should lso limit the number of pproximtions tht re computed, so we count the number of times we divide the subintervls, nd stop when we rech predefined limit which we cll M. 84

7 Algorithm 1.5. Suppose the function f, the intervl [,b], the length n of the intitil prtition, positive tolernce ɛ < 1, nd the mximum number of itertions M re given. The following lgorithm will compute sequence of pproximtions to f (x)d x by the midpoint rule, until the estimted reltive error is smller thn ɛ, or the mximum number of computed pproximtions rech M. The finl pproximtion is stored in I. n := n ; h := (b )/n; I := ; x := + h/; for k := 1,,..., n I := I + f (x); x := x + h; j := 1; I := h I ; bser r := I ; while j < M nd bser r > ɛ I j := j + 1; I p := I ; n := n; h := (b )/n; I := ; x := + h/; for k := 1,,..., n I := I + f (x); x := x + h; I := h I ; bser r := I I p ; Note tht we compute the first pproximtion outside the min loop. This is necessry in order to hve meningful estimtes of the reltive error the first two times we rech the while loop (the first time we rech the while loop we will lwys get pst the condition). We store the previous pproximtion in I p nd use this to estimte the error in the next itertion. In the coming sections we will describe two other methods for numericl integrtion. These cn be implemented in lgorithms similr to Algorithm 1.5. In fct, the only difference will be how the ctul pproximtion to the integrl is computed. Exmple 1.6. Let us try the midpoint rule on n exmple. As usul, it is wise to test on n exmple where we know the nswer, so we cn esily check the qulity 85

8 of the method. We choose the integrl cos x d x = sin where the exct nswer is esy to compute by trditionl, symbolic methods. To test the method, we split the intervl into k subintervls, for k = 1,,..., 1, i.e., we hlve the step length ech time. The result is By error, we here men h I mid (h) Error f (x)d x I mid (h). Note tht ech time the step length is hlved, the error seems to be reduced by fctor of The error Algorithm 1.5 determines numericl pproximtion to the integrl, nd even estimtes the error. However, we must remember tht the error tht is computed is not lwys relible, so we should try nd understnd the error better. We do this in two steps. First we nlyse the error in the sitution where we use very simple prtition with only one subintervl, the so-clled locl error. Then we use this result to obtin n estimte of the error in the generl cse this is often referred to s the globl error. Locl error nlysis Suppose we use the midpoint rule with just one subintervl. We wnt to study the error f (x)d x f ( 1/ ) (b ), 1/ = ( + b)/. (1.5) 86

9 Once gin, Tylor polynomil with reminder helps us out. We expnd f (x) bout the midpoint 1/ nd obtin, f (x) = f ( 1/ ) + (x 1/ )f ( 1/ ) + (x 1/) f (ξ), where ξ is number in the intervl ( 1/, x) tht depends on x. Next, we integrte the Tylor expnsion nd obtin f (x)d x = (f ( 1/ ) + (x 1/ )f ( 1/ ) + (x 1/) = f ( 1/ )(b ) + f ( 1/ )[ (x 1/ ) ] b + 1 ) f (ξ) d x (x 1/ ) f (ξ)d x = f ( 1/ )(b ) + 1 (x 1/ ) f (ξ)d x, (1.6) since the middle term is zero. This leds to n expression for the error, f (x)d x f ( 1/ )(b ) = 1 (x 1/ ) f (ξ)d x. (1.7) Let us simplify the right-hnd side of this expression nd explin fterwords. We hve 1 b (x 1/ ) f (ξ)d x 1 (x 1/ ) f (ξ) d x = 1 (x 1/ ) f (ξ) d x M (x 1/ ) d x = M 1 [ (x 1/ ) 3] b 3 = M ( (b 1/ ) 3 ( 1/ ) 3) 6 = M 4 (b )3, (1.8) where M = mx x [,b] f (x). The first inequlity is vlid becuse when we move the bsolute vlue sign inside the integrl sign, the function tht we integrte becomes nonnegtive everywhere. This mens tht in the res where the integrnd in the originl expression is negtive, everything is now positive, nd hence the second integrl is lrger thn the first. Next there is n equlity which is vlid becuse (x 1/ ) is never negtive. The next inequlity follows becuse we replce f (ξ) with its mximum on the 87

10 intervl [,b]. The next step is just the evlution of the integrl of (x 1/ ), nd the lst equlity follows since (b 1/ ) 3 = ( 1/ ) 3 = (b ) 3 /8. This proves the following lemm. Lemm 1.7. Let f be continuous function whose first two derivtives re continuous on the intervl [.b]. The error in the midpoint rule, with only one intervl, is bounded by f (x)d x f ( 1/ ) (b ) M 4 (b )3, where M = mx x [,b] f (x) nd 1/ = ( + b)/. Before we continue, let us sum up the procedure tht led up to lemm 1.7 without focusing on the detils: Strt with the error (1.5) nd replce f (x) by its liner Tylor polynomil with reminder. When we integrte the Tylor polynomil, the liner term becomes zero, nd we re left with (1.7). At this point we use some stndrd techniques tht give us the finl inequlity. The importnce of lemm 1.7 lies in the fctor (b ) 3. This mens tht if we reduce the size of the intervl to hlf its width, the error in the midpoint rule will be reduced by fctor of 8. Globl error nlysis Above, we nlysed the error on one subintervl. Now we wnt to see wht hppens when we dd together the contributions from mny subintervls. We consider the generl cse where we hve prtition tht divides [,b] into n subintervls, ech of width h. On ech subintervl we use the simple midpoint rule tht we just nlysed, I = The totl error is then f (x)d x = I I mid = ( xi xi x i 1 f (x)d x f (x i 1/ )h. x i 1 f (x)d x f (x i 1/ )h We note tht the expression inside the prenthesis is just the locl error on the 88 ).

11 intervl [x i 1, x i ]. We therefore hve ( xi I I mid = f (x)d x f (x i 1/ )h x i 1 xi f (x)d x f (x i 1/ )h x i 1 h 3 4 M i (1.9) where M i is the mximum of f (x) on the intervl [xi 1, x i ]. The first of these inequlities is just the tringle inequlity, while the second inequlity follows from lemm 1.7. To simplify the expression (1.9), we extend the mximum on [x i 1, x i ] to ll of [,b]. This cnnot mke the mximum smller, so for ll i we hve M i = mx f (x) mx f (x) = M. x [x i 1,x i ] Now we cn simplify (1.9) further, h 3 4 M i h 3 x [,b] ) 4 M = h3 nm. (1.1) 4 Here, we need one finl little observtion. Recll tht h = (b )/n, so hn = b. If we insert this in (1.1), we obtin our min error estimte. Theorem 1.8. Suppose tht f nd its first two derivtives re continuous on the intervl [,b], nd tht the integrl of f on [,b] is pproximted by the midpoint rule with n subintervls of equl width, I = Then the error is bounded by where x i 1/ = + (i 1/)h. f (x)d x I mid = f (x i 1/ )h. I I mid (b ) h 4 mx f (x), (1.11) x [,b] This confirms the error behviour tht we sw in exmple 1.6: If h is reduced by fctor of, the error is reduced by fctor of = 4. 89

12 One notble omission in our discussion of the error in the midpoint rule is round-off error, which ws mjor concern in our study of numericl differentition. The good news is tht round-off error is not usully problem in numericl integrtion. The only sitution where round-off my cuse problems is when the vlue of the integrl is. In such sitution we my potentilly dd mny numbers tht sum to, nd this my led to cncelltion effects. However, this is so rre tht we will not discuss it here Estimting the step length The error estimte (1.11) lets us ply stndrd gme: If someone demnds tht we compute n integrl with error smller thn ɛ, we cn find step length h tht gurntees tht we meet this demnd. To mke sure tht the error is smller thn ɛ, we enforce the inequlity which we cn esily solve for h, (b ) h 4 mx f (x) ɛ x [,b] 4ɛ h (b )M, M = mx f (x). x [,b] This is not quite s simple s it my look since we will hve to estimte M, the mximum vlue of the second derivtive, over the whole intervl of integrtion [,b]. This cn be difficult, but in some cses it is certinly possible, see exercise 3. Exercises 1 Clculte n pproximtion to the integrl π/ sin x 1 + x d x = with the midpoint rule. Split the intervl into 6 subintervls. In this exercise you re going to progrm lgorithm 1.5. If you cnnot progrm, use the midpoint lgorithm with 1 subintervls, check the error, nd skip (b). ) Write progrm tht implements the midpoint rule s in lgorithm 1.5 nd test it on the integrl e x d x = e 1. 9

13 x 1 () x 1 4 (b) Figure 1.4. The trpezoidl rule with one subintervl () nd five subintervls (b). b) Determine vlue of h tht gurntees tht the bsolute error is smller thn 1 1. Run your progrm nd check wht the ctul error is for this vlue of h. (You my hve to djust lgorithm 1.5 slightly nd print the bsolute error.) 3 Repet the previous exercise, but compute the integrl 6 ln x d x = ln(11664). 4 Redo the locl error nlysis for the midpoint rule, but replce both f (x) nd f ( 1/ ) by liner Tylor polynomils with reminders bout the left end point. Wht hppens to the error estimte? 1.3 The trpezoidl rule The midpoint rule is bsed on very simple polynomil pproximtion to the function f to be integrted on ech subintervl; we simply use constnt pproximtion tht interpoltes the function vlue t the middle point. We re now going to consider nturl lterntive; we pproximte f on ech subintervl with the secnt tht interpoltes f t both ends of the subintervl. The sitution is shown in figure 1.4. The pproximtion to the integrl is the re of the trpezoidl polygon under the secnt so we hve f () + f (b) f (x)d x (b ). (1.1) To get good ccurcy, we will hve to split [,b] into subintervls with prtition nd use the trpezoidl pproximtion on ech subintervl, s in figure 1.4b. If we hve uniform prtition {x i } n with step length h, we get the pproximtion i= f (x)d x = xi x i 1 f (x)d x 91 f (x i 1 ) + f (x i ) h. (1.13)

14 We should lwys im to mke our computtionl methods s efficient s possible, nd in this cse n improvement is possible. Note tht on the intervl [x i 1, x i ] we use the function vlues f (x i 1 ) nd f (x i ), nd on the next intervl we use the vlues f (x i ) nd f (x i+1 ). All function vlues, except the first nd lst, therefore occur twice in the sum on the right in (1.13). This mens tht if we implement this formul directly we do lot of unnecessry work. From this the following observtion follows. Observtion 1.9 (Trpezoidl rule). Suppose we hve function f defined on n intervl [,b] nd prtition {x i } n of [,b]. If we pproximte f by its i= secnt on ech subintervl nd pproximte the integrl of f by the integrl of the resulting piecewise liner pproximtion, we obtin the pproximtion ( f () + f (b) f (x)d x I trp (h) = h n 1 + ) f (x i ). (1.14) Once we hve the formul (1.14), we cn esily derive n lgorithm similr to lgorithm 1.5. In fct the two lgorithms re identicl except for the prt tht clcultes the pproximtions to the integrl, so we will not discuss this further. Exmple 1.1. We test the trpezoidl rule on the sme exmple s the midpoint rule, cos x d x = sin As in exmple 1.6 we split the intervl into k subintervls, for k = 1,,..., 1. The resulting pproximtions re h I trp (h) Error

15 where the error is defined by f (x)d x I trp (h). We note tht ech time the step length is hlved, the error is reduced by fctor of 4, just s for the midpoint rule. But we lso note tht even though we now use two function vlues in ech subintervl to estimte the integrl, the error is ctully twice s big s it ws for the midpoint rule The error Our next step is to nlyse the error in the trpezoidl rule. We follow the sme strtegy s for the midpoint rule nd use Tylor polynomils. Becuse of the similrities with the midpoint rule, we skip some of the detils. The locl error We first study the error in the pproximtion (1.1) where we only hve one secnt. In this cse the error is given by f () + f (b) f (x)d x (b ), (1.15) nd the first step is to expnd the function vlues f (x), f (), nd f (b) in Tylor series bout the midpoint 1/, f (x) = f ( 1/ ) + (x 1/ )f ( 1/ ) + (x 1/) f (ξ 1 ), f () = f ( 1/ ) + ( 1/ )f ( 1/ ) + ( 1/) f (ξ ), f (b) = f ( 1/ ) + (b 1/ )f ( 1/ ) + (b 1/) f (ξ 3 ), where ξ 1 ( 1/, x), ξ (, 1/ ), nd ξ 3 ( 1/,b). The integrtion of the Tylor series for f (x) we did in (1.6) so we just quote the result here, f (x)d x = f ( 1/ )(b ) + 1 (x 1/ ) f (ξ 1 )d x. (1.16) We note tht 1/ = (b )/ nd b 1/ = (b )/, so the sum of the Tylor series for f () nd f (b) is (b ) f () + f (b) = f ( 1/ ) + f (b ) (ξ ) + f (ξ 3 ). (1.17)

16 If we insert (1.16) nd (1.17) in the expression for the error (1.15), the first two terms cncel, nd we obtin f () + f (b) f (x)d x (b ) = 1 b (x 1/ ) f (b )3 (ξ 1 )d x f (b )3 (ξ ) f (ξ 3 ) b (x 1/ ) f (ξ 1 )d x (b )3 + f (b )3 (ξ ) + f (ξ 3 ) The lst reltion is just n ppliction of the tringle inequlity. The first term we estimted in (1.8), nd in the lst two we use the stndrd trick nd tke mximum vlues of f (x) over ll of [,b]. Then we end up with f () + f (b) f (x)d x (b ) M 4 (b )3 + M 16 (b )3 + M (b )3 16 = M 6 (b )3. Let us sum this up in lemm. Lemm Let f be continuous function whose first two derivtives re continuous on the intervl [, b]. The error in the trpezoidl rule, with only one secnt bsed t nd b, is bounded by where M = mx x [,b] f (x). f () + f (b) f (x)d x (b ) M 6 (b )3, This lemm is completely nlogous to lemm 1.7 which describes the locl error in the midpoint rule. We prticulrly notice tht even though the trpezoidl rule uses two vlues of f, the error estimte is slightly lrger thn the estimte for the midpoint rule. The most importnt feture is the exponent on (b ), which tells us how quickly the error goes to when the intervl width is reduced, nd from this point of view the two methods re the sme. In other words, we hve gined nothing by pproximting f by liner function insted of constnt. This does not men tht the trpezoidl rule is bd, it rther mens tht the midpoint rule is surprisingly good. 94

17 Globl error We cn find n expression for the globl error in the trpezoidl rule in exctly the sme wy s we did for the midpoint rule, so we skip the proof. Theorem 1.1. Suppose tht f nd its first two derivtives re continuous on the intervl [,b], nd tht the integrl of f on [,b] is pproximted by the trpezoidl rule with n subintervls of equl width h, ( f () + f (b) I = f (x)d x I trp = h Then the error is bounded by n 1 + ) f (x i ). I Itrp h (b ) 6 mx f (x). (1.18) x [,b] The error estimte for the trpezoidl rule is not best possible in the sense tht it is possible to derive better error estimte (using other techniques) with the smller constnt 1/1 insted of 1/6. However, the fct remins tht the trpezoidl rule is bit disppointing compred to the midpoint rule, just s we sw in exmple 1.1. Exercises 1 Clculte n pproximtion to the integrl π/ sin x 1 + x d x = with the trpezoidl rule. Split the intervl into 6 subintervls. In this exercise you re going to progrm n lgorithm like lgorithm 1.5 for the trpezoidl rule. If you cnnot progrm, use the trpezoidl rule mnully with 1 subintervls, check the error, nd skip the second prt of (b). ) Write progrm tht implements the midpoint rule s in lgorithm 1.5 nd test it on the integrl e x d x = e 1. b) Determine vlue of h tht gurntees tht the bsolute error is smller thn 1 1. Run your progrm nd check wht the ctul error is for this vlue of h. (You my hve to djust lgorithm 1.5 slightly nd print the bsolute error.) 95

18 3 Fill in the detils in the derivtion of lemm 1.11 from (1.16) nd (1.17). 4 In this exercise we re going to do n lterntive error nlysis for the trpezoidl rule. Use the sme procedure s in section 1.3.1, but expnd both the function vlues f (x) nd f (b) in Tylor series bout. Compre the resulting error estimte with lemm When h is hlved in the trpezoidl rule, some of the function vlues used with step length h/ re the sme s those used for step length h. Derive formul for the trpezoidl rule with step length h/ tht mkes it esy to void recomputing the function vlues tht were computed on the previous level. 1.4 Simpson s rule The finl method for numericl integrtion tht we consider is Simpson s rule. This method is bsed on pproximting f by prbol on ech subintervl, which mkes the derivtion bit more involved. The error nlysis is essentilly the sme s before, but becuse the expressions re more complicted, we omit it here Derivtion of Simpson s rule As for the other methods, we derive Simpson s rule in the simplest cse where we pproximte f by one prbol on the whole intervl [,b]. We find the polynomil p tht interpoltes f t, 1/ = ( + b)/ nd b, nd pproximte the integrl of f by the integrl of p. We could find p vi the Newton form, but in this cse it is esier to use the Lgrnge form. Another simplifiction is to first construct Simpson s rule in the cse where = 1, 1/ =, nd b = 1, nd then generlise fterwrds. Simpson s rule on [ 1,1] The Lgrnge form of the polynomil tht interpoltes f t 1,, nd 1 is given by x(x 1) (x + 1)x p (x) = f ( 1) f ()(x + 1)(x 1) + f (1), nd it is esy to check tht the interpoltion conditions hold. To integrte p, we must integrte ech of the three polynomils in this expression. For the first one we hve 1 1 x(x 1)d x = 1 1 1(x x)d x = 1 [ 1 3 x3 1 x] 1 = Similrly, we find 1 (x + 1)(x 1)d x = 4 3, (x + 1)x d x = 1 3.

19 On the intervl [ 1, 1], Simpson s rule therefore corresponds to the pproximtion f (x)d x 1 ) f ( 1) + 4f () + f (1). (1.19) 3( 1 Simpson s rule on [,b] To obtin n pproximtion of the integrl on the intervl [,b], we use stndrd technique. Suppose tht x nd y re relted by x = (b ) y (1.) so tht when y vries in the intervl [ 1,1], then x will vry in the intervl [,b]. We re going to use the reltion (1.) s substitution in n integrl, so we note tht d x = (b )d y/. We therefore hve where f (x)d x = b 1 ( ) b f (y + 1) + d y = b f (y)d y, (1.1) 1 ( ) f b (y) = f (y + 1) +. To determine n pproximtion to the integrl of f on the intervl [ 1,1], we use Simpson s rule (1.19). The result is 1 f (y)d y 1 ) f ( 1) + 4f () + f 1( (1) = f () + 4f (1/ ) + f (b) 3( ), 3 since the reltion in (1.) mps 1 to, the midpoint to 1/ = ( +b)/, nd the right endpoint b to 1. If we insert this in (1.1), we obtin Simpson s rule for the generl intervl [,b], see figure 1.5. Observtion Let f be n integrble function on the intervl [,b]. If f is interpolted by qudrtic polynomil p t the points, 1/ = ( + b)/ nd b, then the integrl of f cn be pproximted by the integrl of p, f (x)d x p (x)d x = b ( f () + 4f (1/ ) + f (b) ). (1.) 6 We my just s well derive this formul by doing the interpoltion directly on the intervl [,b], but then the lgebr becomes quite messy. 97

20 x 1 () x 1 4 (b) Figure 1.5. Simpson s rule with one subintervl () nd three subintervls (b). x 1 4 Figure 1.6. Simpson s rule with three subintervls Composite Simpson s rule In prctice, we will usully divide the intervl [, b] into smller subintervls nd use Simpson s rule on ech subintervl, see figure 1.5b. Note though tht Simpson s rule is not quite like the other numericl integrtion techniques we hve studied when it comes to splitting the intervl into smller pieces: The intervl over which f is to be integrted is split into subintervls, nd Simpson s rule is pplied on neighbouring pirs of intervls, see figure 1.6. In other words, ech prbol is defined over two subintervls which mens tht the totl number of subintervls must be even, nd the number of given vlues of f must be odd. If the prtition is {x i } n i= with x i = + i h, Simpson s rule on the intervl 98

21 [x i, x i ] is xi x i f (x)d x h 3 ( f (xi ) + 4f (x i 1 ) + f (x i ) ). The pproximtion of the totl integrl is therefore f (x)d x h 3 ( (f (xi ) + 4f (x i 1 ) + f (x i ) ). In this sum we observe tht the right endpoint of one subintervl becomes the left endpoint of the neighbouring subintervl to the right. Therefore, if this is implemented directly, the function vlues t the points with n even subscript will be evluted twice, except for the extreme endpoints nd b which only occur once in the sum. We cn therefore rewrite the sum in wy tht voids these redundnt evlutions. Observtion Suppose f is function defined on the intervl [, b], nd let {x i } n be uniform prtition of [,b] with step length h. The composite i= Simpson s rule pproximtes the integrl of f by f (x)d x I Simp (h) = h ( n 1 f () + f (b) + f (x i ) ) f (x i 1 ). With the midpoint rule, we computed sequence of pproximtions to the integrl by successively hlving the width of the subintervls. The sme is often done with Simpson s rule, but then cre should be tken to void unnecessry function evlutions since ll the function vlues computed t one step will lso be used t the next step. Exmple Let us test Simpson s rule on the sme exmple s the midpoint rule nd the trpezoidl rule, cos x d x = sin As in exmple 1.6, we split the intervl into k subintervls, for k = 1,,..., 1. 99

22 The result is where the error is defined by h I Simp (h) Error f (x)d x I Simp (h). When we compre this tble with exmples 1.6 nd 1.1, we note tht the error is now much smller. We lso note tht ech time the step length is hlved, the error is reduced by fctor of 16. In other words, by introducing one more function evlution in ech subintervl, we hve obtined method with much better ccurcy. This will be quite evident when we nlyse the error below The error An expression for the error in Simpson s rule cn be derived by using the sme technique s for the previous methods: We replce f (x), f () nd f (b) by cubic Tylor polynomils with reminders bout the point 1/, nd then collect nd simplify terms. However, these computtions become quite long nd tedious, nd s for the trpezoidl rule, the constnt in the error term is not the best possible. We therefore just stte the best possible error estimte here without proof. Lemm 1.16 (Locl error). If f is continuous nd hs continuous derivtives up to order 4 on the intervl [,b], the error in Simpson s rule is bounded by (b )5 E(f ) mx f (i v) (x). 88 x [,b] We note tht the error in Simpson s rule depends on (b ) 5, while the error in the midpoint rule nd trpezoidl rule depend on (b ) 3. This mens tht 3

23 the error in Simpson s rule goes to zero much more quickly thn for the other two methods when the width of the intervl [,b] is reduced. The globl error The pproch we used to deduce the globl error for the midpoint rule, see theorem 1.8, cn lso be used to derive the globl error in Simpson s rule. The following theorem sums this up. Theorem 1.17 (Globl error). Suppose tht f nd its first 4 derivtives re continuous on the intervl [, b], nd tht the integrl of f on [, b] is pproximted by Simpson s rule with n subintervls of equl width h. Then the error is bounded by E(f ) h 4 (b ) 88 mx f (i v) (x). (1.3) x [,b] The estimte (1.3) explins the behviour we noticed in exmple 1.15: Becuse of the fctor h 4, the error is reduced by fctor 4 = 16 when h is hlved, nd for this reson, Simpson s rule is very populr method for numericl integrtion. Exercises 1 Clculte n pproximtion to the integrl π/ sin x 1 + x d x = with Simpson s rule. Split the intervl into 6 subintervls. ) How mny function evlutions do you need to clculte the integrl d x 1 + x with the trpezoidl rule to mke sure tht the error is smller thn 1 1. b) How mny function evlutions re necessry to chieve the sme ccurcy with the midpoint rule? c) How mny function evlutions re necessry to chieve the sme ccurcy with Simpson s rule? 3 In this exercise you re going to progrm n lgorithm like lgorithm 1.5 for Simpson s rule. If you cnnot progrm, use Simpson s rule mnully with 1 subintervls, check the error, nd skip the second prt of (b). 31

24 ) Write progrm tht implements Simpson s rule s in lgorithm 1.5 nd test it on the integrl e x d x = e 1. b) Determine vlue of h tht gurntees tht the bsolute error is smller thn 1 1. Run your progrm nd check wht the ctul error is for this vlue of h. (You my hve to djust lgorithm 1.5 slightly nd print the bsolute error.) 4 ) Verify tht Simpson s rule is exct when f (x) = x i for i =, 1,, 3. b) Use () to show tht Simpson s rule is exct for ny cubic polynomil. c) Could you rech the sme conclusion s in (b) by just considering the error estimte (1.3)? 5 We wnt to design numericl integrtion method f (x)d x w 1 f () + w f ( 1/ ) + w 3 f (b). Determine the unknown coefficients w 1, w, nd w 3 by demnding tht the integrtion method should be exct for the three polynomils f (x) = x i for i =, 1,. Do you recognise the method? 1.5 Summry In this chpter we hve derived three methods for numericl integrtion. All these methods nd their error nlyses my seem rther overwhelming, but they ll follow common thred: Procedure The following is generl procedure for deriving numericl methods for integrtion of function f over the intervl [,b]: 1. Interpolte the function f by polynomil p t suitble points.. Approximte the integrl of f by the integrl of p. This mkes it possible to express the pproximtion to the integrl in terms of function vlues of f. 3. Derive n estimte for the locl error by expnding the function vlues in Tylor series with reminders bout the midpoint 1/ = ( + b)/. 4. Derive n estimte for the globl error by using the technique leding up to theorem

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical 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 information

Lecture 14: Quadrature

Lecture 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 information

Review of Calculus, cont d

Review 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 information

Lecture 20: Numerical Integration III

Lecture 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 information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 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 information

Math& 152 Section Integration by Parts

Math& 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 information

Numerical Integration

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 information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

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 information

Math 1B, lecture 4: Error bounds for numerical methods

Math 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 information

The Regulated and Riemann Integrals

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 information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: 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 information

Z 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...

Z 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 information

Chapter 0. What is the Lebesgue integral about?

Chapter 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 information

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Riemann 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 information

Riemann Sums and Riemann Integrals

Riemann 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 information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 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 information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced 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 information

7.2 The Definite Integral

7.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 information

Section 6.1 Definite Integral

Section 6.1 Definite Integral Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined

More information

P 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)

P 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 information

Riemann Sums and Riemann Integrals

Riemann 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 information

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.

Suppose 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 information

Definite integral. Mathematics FRDIS MENDELU

Definite 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 information

APPROXIMATE INTEGRATION

APPROXIMATE 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 information

NUMERICAL INTEGRATION

NUMERICAL 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 information

Chapter 5. Numerical Integration

Chapter 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 information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties 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 information

20 MATHEMATICS POLYNOMIALS

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

More information

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Definite 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 information

1 The Riemann Integral

1 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 information

Sections 5.2: The Definite Integral

Sections 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 information

Math 113 Exam 2 Practice

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

More information

Improper Integrals, and Differential Equations

Improper 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 information

Chapters 4 & 5 Integrals & Applications

Chapters 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 information

Math 131. Numerical Integration Larson Section 4.6

Math 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 information

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230 Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

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

More information

Review of basic calculus

Review 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 information

n 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

n 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 information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

Recitation 3: More Applications of the Derivative

Recitation 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 information

1 Part II: Numerical Integration

1 Part II: Numerical Integration Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble

More information

Chapter 6 Notes, Larson/Hostetler 3e

Chapter 6 Notes, Larson/Hostetler 3e Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn

More information

Best Approximation. Chapter The General Case

Best 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 information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

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

Chapter 7 Notes, Stewart 8e. 7.1 Integration by Parts Trigonometric Integrals Evaluating sin m x cos n (x) dx...

Chapter 7 Notes, Stewart 8e. 7.1 Integration by Parts Trigonometric Integrals Evaluating sin m x cos n (x) dx... Contents 7.1 Integrtion by Prts................................... 2 7.2 Trigonometric Integrls.................................. 8 7.2.1 Evluting sin m x cos n (x)......................... 8 7.2.2 Evluting

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE 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 information

Lecture 12: Numerical Quadrature

Lecture 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 information

UNIFORM 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 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 information

Numerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden

Numerical 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 information

5.7 Improper Integrals

5.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 information

The Riemann Integral

The Riemann Integral Deprtment of Mthemtics King Sud University 2017-2018 Tble of contents 1 Anti-derivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Anti-derivtive Function Definition Let f : I R be function

More information

Numerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1

Numerical 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 information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

ACCESS 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 information

1 Error Analysis of Simple Rules for Numerical Integration

1 Error Analysis of Simple Rules for Numerical Integration cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion

More information

Main topics for the First Midterm

Main topics for the First Midterm Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the

More information

3.4 Numerical integration

3.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 information

Lecture 1: Introduction to integration theory and bounded variation

Lecture 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 information

Anti-derivatives/Indefinite Integrals of Basic Functions

Anti-derivatives/Indefinite Integrals of Basic Functions Anti-derivtives/Indefinite Integrls of Bsic Functions Power Rule: In prticulr, this mens tht x n+ x n n + + C, dx = ln x + C, if n if n = x 0 dx = dx = dx = x + C nd x (lthough you won t use the second

More information

MATH 144: Business Calculus Final Review

MATH 144: Business Calculus Final Review MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives

More information

Overview of Calculus I

Overview 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 information

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

DOING 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 information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

MORE 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 information

III. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9

III. Lecture on 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 information

Chapter 3 Solving Nonlinear Equations

Chapter 3 Solving Nonlinear Equations Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,

More information

6.5 Numerical Approximations of Definite Integrals

6.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 information

INTRODUCTION TO INTEGRATION

INTRODUCTION 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 information

5.1 How do we Measure Distance Traveled given Velocity? Student Notes

5.1 How do we Measure Distance Traveled given Velocity? Student Notes . How do we Mesure Distnce Trveled given Velocity? Student Notes EX ) The tle contins velocities of moving cr in ft/sec for time t in seconds: time (sec) 3 velocity (ft/sec) 3 A) Lel the x-xis & y-xis

More information

LECTURE 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 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 information

Discrete Least-squares Approximations

Discrete Least-squares Approximations Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

f(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

f(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 information

Math 8 Winter 2015 Applications of Integration

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

More information

Calculus I-II Review Sheet

Calculus 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 information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Abstract inner product spaces

Abstract inner product spaces WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper 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 information

1. 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. 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 information

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

and 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 information

1 The Lagrange interpolation formula

1 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 information

4.4 Areas, Integrals and Antiderivatives

4.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 information

The practical version

The practical version Roerto s Notes on Integrl Clculus Chpter 4: Definite integrls nd the FTC Section 7 The Fundmentl Theorem of Clculus: The prcticl version Wht you need to know lredy: The theoreticl version of the FTC. Wht

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, 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 information

Interpreting Integrals and the Fundamental Theorem

Interpreting 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 information

A 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. 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 information

How 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?

How 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 information

The 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 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 information

Integration Techniques

Integration Techniques Integrtion Techniques. Integrtion of Trigonometric Functions Exmple. Evlute cos x. Recll tht cos x = cos x. Hence, cos x Exmple. Evlute = ( + cos x) = (x + sin x) + C = x + 4 sin x + C. cos 3 x. Let u

More information

different methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).

different 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 information

LECTURE. INTEGRATION AND ANTIDERIVATIVE.

LECTURE. INTEGRATION AND ANTIDERIVATIVE. ANALYSIS FOR HIGH SCHOOL TEACHERS LECTURE. INTEGRATION AND ANTIDERIVATIVE. ROTHSCHILD CAESARIA COURSE, 2015/6 1. Integrtion Historiclly, it ws the problem of computing res nd volumes, tht triggered development

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC 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 information

Riemann Integrals and the Fundamental Theorem of Calculus

Riemann Integrals and the Fundamental Theorem of Calculus Riemnn Integrls nd the Fundmentl Theorem of Clculus Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University September 16, 2013 Outline Grphing Riemnn Sums

More information

Lecture 23: Interpolatory Quadrature

Lecture 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 information

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

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

More information

10. AREAS BETWEEN CURVES

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

More information

Summer MTH142 College Calculus 2. Section J. Lecture Notes. Yin Su University at Buffalo

Summer 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 information

Numerical quadrature based on interpolating functions: A MATLAB implementation

Numerical 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 information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This

More information