Numerical Integration


 Oswald Hunter
 1 years ago
 Views:
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 roundoff error, integrtion behves differently from differentition: Numericl integrtion is very insensitive to roundoff errors, so we will ignore roundoff 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 roundoff 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 soclled 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 soclled 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 righthnd 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 roundoff error, which ws mjor concern in our study of numericl differentition. The good news is tht roundoff error is not usully problem in numericl integrtion. The only sitution where roundoff 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
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 relvlues nd smooth The pproximtion of n integrl by numericl
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 200910 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 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 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& 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 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 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 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 informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More 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 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 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 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 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 informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixedpoint itertion to converge when solving the eqution
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 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 informationSection 6.1 Definite Integral
Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined
More 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 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 informationSuppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = 2.
Mth 43 Section 6. Section 6.: Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot
More 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 informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose ntiderivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
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 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 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 information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More 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 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 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 informationMath 113 Exam 2 Practice
Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number
More informationImproper Integrals, 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 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 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 informationPolynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230
Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given
More 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 informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More informationSection 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 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 information1 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 informationChapter 6 Notes, Larson/Hostetler 3e
Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn
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 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 informationChapter 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 informationTHE EXISTENCEUNIQUENESS THEOREM FOR FIRSTORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCEUNIQUENESS THEOREM FOR FIRSTORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrdlindeloftheorem/ This document is proof of the existenceuniqueness theorem
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 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 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 ThreePoint
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 informationThe Riemann Integral
Deprtment of Mthemtics King Sud University 20172018 Tble of contents 1 Antiderivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Antiderivtive Function Definition Let f : I R be function
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 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 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 informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 23, Sections 4.14.8, nd Sections 5.15.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of 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 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 informationAntiderivatives/Indefinite Integrals of Basic Functions
Antiderivtives/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 informationMATH 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 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 informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OEDIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More informationW. 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 informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecturenotes/numericalinter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecturenotes/numericalinter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More informationChapter 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 noninteger. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More information6.5 Numerical Approximations of Definite Integrals
Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 6.5 Numericl Approximtions of Definite Integrls Sometimes the integrl of function cnnot be expressed with elementry functions, i.e., polynomil,
More 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 information5.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 xxis & yxis
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 informationDiscrete Leastsquares Approximations
Discrete Lestsqures Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More 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 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 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 informationCalculus III Review Sheet
Clculus III 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 informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06073 HIKARI Ltd, www.mhikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationAbstract 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 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. GaussJacobi 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. GussJcobi 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 informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
More 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 information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationThe 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 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 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 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 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 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 informationIntegration 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 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 informationLECTURE. 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 informationpadic Egyptian Fractions
padic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Setup 3 4 pgreedy Algorithm 5 5 pegyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationRiemann 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 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 informationWe know that if f is a continuous nonnegative function on the interval [a, b], then b
1 Ares Between Curves c 22 Donld Kreider nd Dwight Lhr We know tht if f is continuous nonnegtive function on the intervl [, b], then f(x) dx is the re under the grph of f nd bove the intervl. We re going
More information10. AREAS BETWEEN CURVES
. AREAS BETWEEN CURVES.. Ares etween curves So res ove the xxis re positive nd res elow re negtive, right? Wrong! We lied! Well, when you first lern out integrtion it s convenient fiction tht s true in
More 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 informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More informationMath 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