Lecture 12: Numerical Quadrature
|
|
- Charlene Black
- 5 years ago
- Views:
Transcription
1 Lecture 12: Numericl Qudrture J.K. WI3097TU Delft Institute of Applied Mthemtics Delft University of Technology 5 December 2012 () Numericl Qudrture 5 December / 46
2 Outline 1 Review Newton-Rphson Systems of Nonliner Equtions Riemnn integrtion 2 Errors for Riemnn integrtion 3 Newton-Cotes Qudrture () Numericl Qudrture 5 December / 46
3 Newton-Rphson method We need to determine the solution to the problem: Find x such tht f (x) = 0 Assumptions f (x) C 2 [, b] (two continuous derivtives) x is n pproximtion to the root, α. f ( x) 0 x α = ɛ is smll. () Numericl Qudrture 5 December / 46
4 Newton-Rphson method By expnding f (x) in Tylor series round x, nd plugging in x = α, we hve: 0 = f (α) = f ( x)+(α x)f ( x)+ 1 2 (α x)2 f (ξ(α)), ξ(α) (α, x). This gives us the Newton-Rphson method: α = x f ( x) f ( x) () Numericl Qudrture 5 December / 46
5 Newton-Rphson method Algorithm 1 Given α 0 2 Define α n = α n 1 f (α n 1) f (α n 1 ) 3 If α n α n 1 ɛ then STOP! α = α n. 4 Else go to Step 2. () Numericl Qudrture 5 December / 46
6 Vritions to Newton s method Secnt Method (Chord Newton) Newton-Rphson, α n = α n 1 f (α n 1) f (α n 1 ) Replce f (α n 1 ) by n pproximtion, So tht α n = α n 1 f (α n 1 ) f (α n 1) f (α n 2 ) α n 1 α n 2 f (α n 1 ) f (α n 1 ) f (α n 2 ) (α n 1 α n 2 ) This is the Secnt Method () Numericl Qudrture 5 December / 46
7 Vritions to Newton s method Method of flse position Combines the ides of Chord Newton nd Bisection methods. Given two pproximtions α 0, α 1 such tht f (α 0 )f (α 1 ) < 0 : While n N 0 If qq 1 < 0 then else q 0 = f (α 0 ) q 1 = f (α 1 ), n = 2 α = α 1 q 1 (α 1 α 0 ) q 1 q 0 q = f (α), n = n + 1 α 0 = α q 0 = q α 1 = α q 1 = q () Numericl Qudrture 5 December / 46
8 Systems of Nonliner Equtions We cn write the system of equtions f 1 (x 1, x 2, x 3,..., x n ) = 0, f 2 (x 1, x 2, x 3,..., x n ) = 0, f 3 (x 1, x 2, x 3,..., x n ) = 0, f n (x 1, x 2, x 3,..., x n ) = 0, s F(x) = 0. () Numericl Qudrture 5 December / 46
9 Systems of Nonliner Equtions The Newton-Rphson method is then x (p) = x (p 1) J 1 (x (p 1) )F(x (p 1) ) where J(x) = Jcobin = f 1 f 1 x 1 x 2 f 2 f 2 x 1 x 2 f n 1 x n 1 f n x n 1 f n 1 x n f n x n () Numericl Qudrture 5 December / 46
10 Riemnn integrtion Recll from Clculus the formul for Riemnn integrtion: 1 Define prtition, P of intervl [, b] : = x 0 < x 1 < x 2 < < x n = b, where the length of the subintervls is given by h k = x k x k 1, k 1. 2 Define scttering, T, to be set of intermedite points: x k 1 t k x k, k = 1,..., n. 3 The Riemnn sum is then given by R n (f ) = n f (t k )h k. k=1 () Numericl Qudrture 5 December / 46
11 Riemnn integrtion Rectngle rule Recll tht b f (x) dx = lim n R n(f ) = lim n n f (t k )h k. k=1 where h = h k = b n. To obtin the Rectngle rule, we tke t k to be n endpoint of the subintervl. Tht is, t k = x k 1 (left endpoint) OR t k = x k (right endpoint). Then I R = h n f (x k 1 ) OR I R = h k=1 n f (x k ). k=1 () Numericl Qudrture 5 December / 46
12 Riemnn integrtion Midpoint rule For the Midpoint rule, we tke t k to be the midpoint of the subintervls: t k = 1 2 (x k 1 + x k ). The Midpoint rule is then n I M = h f (t k ). k=1 () Numericl Qudrture 5 December / 46
13 Riemnn integrtion Errors The errors for these methods re Rectngle rule: b first order method. Midpoint rule: b which is second order. f (x) dx I R 1 2 M 1(b )h, f (x) dx I M 1 24 M 2(b )h 2, But! Don t forget there is lso rounding error. () Numericl Qudrture 5 December / 46
14 Riemnn integrtion Errors Wht is the totl ffect of the errors? Let ˆf (x) = f (x) + ɛ(x) be the perturbed exct solution. Then b b (f (x) ˆf (x)) dx f (x) ˆf (x) dx ɛ mx (b ). () Numericl Qudrture 5 December / 46
15 Riemnn integrtion Errors Consider the rectngle rule, I R (f ) = h n f (x k ). k=1 In this cse, the error is b (f (x) ˆf (x)) dx b f (x) dx I R (f ) + I R (f ) I R (ˆf ) }{{}}{{} Rounding error Approximtion error 1 n 2 M 1(b )h + h f (x k ) ˆf (x k ) }{{} k=1 ɛ(x k ) () Numericl Qudrture 5 December / 46
16 Riemnn integrtion Errors The totl error is given by b (f (x) ˆf (x)) dx 1 2 M 1(b )h + h n ɛ(x k ) k=1 If we ssume ɛ(x) ɛ mx, then b (f (x) ˆf (x)) dx 1 n 2 M 1(b )h + h ɛ mx k=1 So tht b 1 2 M 1(b )h + }{{} hn ɛ mx =(b ) ( ) 1 (f (x) ˆf (x)) dx 2 M 1h + ɛ mx (b ) () Numericl Qudrture 5 December / 46
17 Qudrture Errors Remember: The pproximtion errors should dominte. 1 2 M 1h > ɛ mx h > 2 M 1 ɛ mx () Numericl Qudrture 5 December / 46
18 Qudrture Errors Assume Exct Approximte Exct f (x) ˆf (x) f (x) ɛ. Recll, the reltive error is given by b = f (x) dx b b f (x) dx ˆf (x) dx b f (x) dx b f (x) dx ɛ = Reltive error. () Numericl Qudrture 5 December / 46
19 Qudrture Errors Definition (Condition number) For numericl qudrture, the condition number is given by K I = b f (x) dx b f (x) dx. If K I >> 1, then the problem is ill-posed. () Numericl Qudrture 5 December / 46
20 Qudrture Exmple: Condition Number Exmple The profits or losses per dy of cr mnufcturer depend on the seson. Assume the following profit formuls for the first nd third qurters: w spr (t) = sin(π(t 1/2)), t [0, 1], w fll (t) = sin(π(t + 1/2)). Totl first qurter profit is given by 1 0 w spr (t) dt. Is this well-posed problem? Wht re the first nd third qurter profits if we use the rectngle rule with n = 50? () Numericl Qudrture 5 December / 46
21 Qudrture Exmple: Condition Number To find out if the problem is well-posed, we need to find the condition number: K I = w spr(t) dt 0 w spr(t) dt where 1 0 w spr (t) dt = 1 ( sin(π(t 1/2))) dt 0 = 0.01t π cos(π(t 1/2)) = 0.01 () Numericl Qudrture 5 December / 46
22 Qudrture Exmple: Condition Number w spr (t) dt = = sin(π(t 1/2)) dt 0 1/ ( sin(π(t 1/2))) dt ( sin(π(t 1/2))) dt 1/2 ( = ) ( ) π π = π () Numericl Qudrture 5 December / 46
23 Qudrture Exmple: Condition Number The condition number is then K I = w spr(t) dt 0 w spr(t) dt = Lrge condition number so ill-posed. If we use the rectngle rule with n = 50, w spr = 0.01 w fll = Probbly need bigger n. () Numericl Qudrture 5 December / 46
24 Newton-Cotes Qudrture Let us look t nother method for evluting b f (x) dx. Ide: Replce f (x) by n pproximtion nd then integrte. b f (x) dx b p(x) dx. For exmple, using liner interpolnt, we cn replce f (x) by p 1 (x) : (b x)f () + (x )f (b) p 1 (x) =. b () Numericl Qudrture 5 December / 46
25 Newton-Cotes Qudrture Then we hve the Trpezoidl rule over one domin: T 1 (f ) = b where x 1 = nd x 2 = b. p 1 (x) dx = 1 (b )(f () + f (b)). 2 () Numericl Qudrture 5 December / 46
26 Newton-Cotes Qudrture However, we cn obtin better pproximtion by breking the intervl [, b] into smller subintervls. () Numericl Qudrture 5 December / 46
27 Newton-Cotes Qudrture Define n to be the number of subintervls nd h = b n. Define the nodes to be x j = + jh, j = 0, 1,..., n. Then I(f ) = b f (x) dx = x1 x 0 f (x) dx + + n = k=1 xn 1 xk x2 x n 2 f (x) dx + x k 1 f (x) dx. x 1 f (x) dx + xn x n 1 x3 x 2 f (x) dx f (x) dx () Numericl Qudrture 5 December / 46
28 Newton-Cotes Qudrture We then use liner interpoltion for ech integrl: xk x k 1 f (x) dx x k x k 1 (f (x k 1 ) + f (x k )) = h 2 2 (f (x k 1) + f (x k )). So tht I(f ) = b f (x) dx = = n k=1 n k=1 xk x k 1 f (x) dx h 2 (f (x k 1) + f (x k )) T (f )= h n 1 2 (f (x 0) + f (x n )) + h k=1 f (x k ) which is the Trpezoidl Rule. () Numericl Qudrture 5 December / 46
29 Newton-Cotes Qudrture The error is given by b f (x) dx T (f ) (b ) h 2 f (ξ). 12 Rectngle rule is O(h). Trpezoidl rule is O(h 2 ). But require the sme mount of work, but the Trpezoidl rule is higher order. () Numericl Qudrture 5 December / 46
30 Newton-Cotes Qudrture Exmple Evlute 4 0 dx 1 + x 2 using the trpezoidl rule with two subintervls. Compre it to the exct solution. () Numericl Qudrture 5 December / 46
31 Newton-Cotes Qudrture The exct solution is given by 4 0 dx 1 + x 2 = 2 0 dx 1 + x = tn 1 (4) dx 1 + x 2 () Numericl Qudrture 5 December / 46
32 Newton-Cotes Qudrture Now, let us evlute the integrl using the Trpezoidl rule with n = 2. Then Then = 0, b = 4, h = = 2. x 0 = = 0, x 1 = x 0 + h = 2, x 2 = b = 4. () Numericl Qudrture 5 December / 46
33 Newton-Cotes Qudrture Then T (f ) = h n 1 2 (f (x 0) + f (x n )) + h k=1 = h 2 (f (x 0) + f (x 2 )) + hf (x 1 ) f (x k ) = 2 2 ( ) = = () Numericl Qudrture 5 December / 46
34 Newton-Cotes Qudrture The error is bounded by b f (x) dx T (f ) (b ) h 2 f (ξ) 12 ξ [0,4] mx = 8 3 2(3ξ 2 1) (1 + x 2 ) 3 }{{} = 2 t ξ=0 The ctul error is b f (x) dx T (f ) n=1 = tn 1 (4) = () Numericl Qudrture 5 December / 46
35 Qudrture Summry There re two pproches to pproximting b f (x) dx Approximting the entire intervl using the ide of re under curve: Riemnn integrtion: Rectngle nd Midpoint rule Approximting f (x) using liner interpoltion nd then integrting Trpezoidl rule. () Numericl Qudrture 5 December / 46
36 Qudrture Exmple Exmple We wnt to determine the integrl 1 1 ((10x) ) dx. 1 The reltive rounding error in the function vlues is less thn ɛ. Determine the reltive error in the integrl due to rounding errors. 2 Use the composite midpoint rule s the numericl qudrture method nd ɛ < Give resonble vlue for the step size, h. () Numericl Qudrture 5 December / 46
37 Qudrture Exmple 1. Reltive error due to rounding errors. In this exmple, = 1, b = 1, nd f (x) = (10x) Assume f (x) ˆf (x) f (x) ɛ. Then, the reltive error is b f (x) dx b ˆf (x) dx b b f (x) dx b f (x) ˆf (x) dx b f (x) dx f (x) dx ɛ b f (x) dx = Kɛ () Numericl Qudrture 5 December / 46
38 Qudrture Exmple Reltive error = Kɛ where K = b f (x) dx b f (x) dx. This is just the condition number. Let us determine the vlue of K. () Numericl Qudrture 5 December / 46
39 Qudrture Exmple First, the numertor: b 1 f (x) dx = = (10x) dx ((10x) ) dx 0.01 = ( x x =500 ((10x) ) dx ) ( ) x x 0.0 () Numericl Qudrture 5 December / 46
40 Qudrture Exmple Now, the denomintor: b f (x) dx = 1 ((10x) ) dx 1 ( 1 1 = x x) = () Numericl Qudrture 5 December / 46
41 Qudrture Exmple Therefore, the reltive error is given by b f (x) dx b ˆf (x) dx b b f (x) dx f (x) ˆf (x) dx b f (x) dx = Kɛ = ɛ = ɛ We cn see from the condition number tht this problem is ill-posed. () Numericl Qudrture 5 December / 46
42 Qudrture Exmple 2. Determining the step size. Let the rounding error be given by ɛ < We need for the pproximtion error to be less thn the rounding error where f (x) ˆf (x) f (x) ɛ. Therefore the bsolute error is given through f (x) dx ˆf (x) dx 1 1 f (x) dx ˆf (x) dx 1 ɛ 1 1 f (x) dx () Numericl Qudrture 5 December / 46
43 Qudrture Exmple So tht the roundoff error is f (x) dx ˆf (x) dx 1 1 ɛ (10x) dx 1 500ɛ = if ɛ = () Numericl Qudrture 5 December / 46
44 Qudrture Exmple The pproximtion error is given by where M = b 24 h2 M mx f (x) = mx 6000x. x [ 1,1] x [ 1,1] Thus, the pproximtion error is bounded by 1 f (x) dx I M (f ) 2 24 h = 500h 2. 1 () Numericl Qudrture 5 December / 46
45 Qudrture Exmple We need to choose the step size, h, so tht the pproximtion error is less thn the rounding error: 500h 2 < 500ɛ h < ɛ = if ɛ = () Numericl Qudrture 5 December / 46
46 Mteril ddressed 1 Review Newton-Rphson Systems of Nonliner Equtions Riemnn integrtion 2 Errors for Riemnn integrtion 3 Newton-Cotes Qudrture Mteril in book: Chpter 4 Useful exercises: Ch4: 1-5 () Numericl Qudrture 5 December / 46
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 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 informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More 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 informationThe Fundamental Theorem of Calculus
The Fundmentl Theorem of Clculus MATH 151 Clculus for Mngement J. Robert Buchnn Deprtment of Mthemtics Fll 2018 Objectives Define nd evlute definite integrls using the concept of re. Evlute definite integrls
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More 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 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 informationMidpoint Approximation
Midpoint Approximtion Sometimes, we need to pproximte n integrl of the form R b f (x)dx nd we cnnot find n ntiderivtive in order to evlute the integrl. Also we my need to evlute R b f (x)dx where we do
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 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 informationThe 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 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 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 informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More 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 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 informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More 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 informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More 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 informationLecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),
More 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 informationB.Sc. in Mathematics (Ordinary)
R48/0 DUBLIN INSTITUTE OF TECHNOLOGY KEVIN STREET, DUBLIN 8 B.Sc. in Mthemtics (Ordinry) SUPPLEMENTAL EXAMINATIONS 01 Numericl Methods Dr. D. Mckey Dr. C. Hills Dr. E.A. Cox Full mrks for complete nswers
More informationBig idea in Calculus: approximation
Big ide in Clculus: pproximtion Derivtive: f (x) = df dx f f(x +h) f(x) =, x h rte of chnge is pproximtely the rtio of chnges in the function vlue nd in the vrible in very short time Liner pproximtion:
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More 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 informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
More 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 informationIII. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More informationWe divide the interval [a, b] into subintervals of equal length x = b a n
Arc Length Given curve C defined by function f(x), we wnt to find the length of this curve between nd b. We do this by using process similr to wht we did in defining the Riemnn Sum of definite integrl:
More informationNumerical Integration
Numericl Integrtion Wouter J. Den Hn London School of Economics c 2011 by Wouter J. Den Hn June 3, 2011 Qudrture techniques I = f (x)dx n n w i f (x i ) = w i f i i=1 i=1 Nodes: x i Weights: w i Qudrture
More informationDefinite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +
Definite Integrls --5 The re under curve cn e pproximted y dding up the res of rectngles. Exmple. Approximte the re under y = from x = to x = using equl suintervls nd + x evluting the function t the left-hnd
More informationSection 5.4 Fundamental Theorem of Calculus 2 Lectures. Dr. Abdulla Eid. College of Science. MATHS 101: Calculus 1
Section 5.4 Fundmentl Theorem of Clculus 2 Lectures College of Science MATHS : Clculus (University of Bhrin) Integrls / 24 Definite Integrl Recll: The integrl is used to find re under the curve over n
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More 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 x-xis & y-xis
More information1 The fundamental theorems of calculus.
The fundmentl theorems of clculus. The fundmentl theorems of clculus. Evluting definite integrls. The indefinite integrl- new nme for nti-derivtive. Differentiting integrls. Tody we provide the connection
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationConstruction of Gauss Quadrature Rules
Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture
More 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 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 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 information5: The Definite Integral
5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce
More 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 informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationMain 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 informationTangent Line and Tangent Plane Approximations of Definite Integral
Rose-Hulmn Undergrdute Mthemtics Journl Volume 16 Issue 2 Article 8 Tngent Line nd Tngent Plne Approximtions of Definite Integrl Meghn Peer Sginw Vlley Stte University Follow this nd dditionl works t:
More informationAP Calculus Multiple Choice: BC Edition Solutions
AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this
More 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 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 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 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 informationf(a+h) f(a) x a h 0. This is the rate at which
M408S Concept Inventory smple nswers These questions re open-ended, nd re intended to cover the min topics tht we lerned in M408S. These re not crnk-out-n-nswer problems! (There re plenty of those in the
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More information7.2 Riemann Integrable Functions
7.2 Riemnn Integrble Functions Theorem 1. If f : [, b] R is step function, then f R[, b]. Theorem 2. If f : [, b] R is continuous on [, b], then f R[, b]. Theorem 3. If f : [, b] R is bounded nd continuous
More 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 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 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, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationChapter 6. Riemann Integral
Introduction to Riemnn integrl Chpter 6. Riemnn Integrl Won-Kwng Prk Deprtment of Mthemtics, The College of Nturl Sciences Kookmin University Second semester, 2015 1 / 41 Introduction to Riemnn integrl
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 information1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
More 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 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 informationAdvanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.
Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationIntegrals - Motivation
Integrls - Motivtion When we looked t function s rte of chnge If f(x) is liner, the nswer is esy slope If f(x) is non-liner, we hd to work hrd limits derivtive A relted question is the re under f(x) (but
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationWeek 10: Riemann integral and its properties
Clculus nd Liner Algebr for Biomedicl Engineering Week 10: Riemnn integrl nd its properties H. Führ, Lehrstuhl A für Mthemtik, RWTH Achen, WS 07 Motivtion: Computing flow from flow rtes 1 We observe the
More 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 informationSection 6: Area, Volume, and Average Value
Chpter The Integrl Applied Clculus Section 6: Are, Volume, nd Averge Vlue Are We hve lredy used integrls to find the re etween the grph of function nd the horizontl xis. Integrls cn lso e used to find
More informationAB Calculus Review Sheet
AB Clculus Review Sheet Legend: A Preclculus, B Limits, C Differentil Clculus, D Applictions of Differentil Clculus, E Integrl Clculus, F Applictions of Integrl Clculus, G Prticle Motion nd Rtes This is
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 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 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 informationMTH 122 Fall 2008 Essex County College Division of Mathematics Handout Version 10 1 October 14, 2008
MTH 22 Fll 28 Essex County College Division of Mthemtics Hndout Version October 4, 28 Arc Length Everyone should be fmilir with the distnce formul tht ws introduced in elementry lgebr. It is bsic formul
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 information( ) as a fraction. Determine location of the highest
AB Clculus Exm Review Sheet - Solutions A. Preclculus Type prolems A1 A2 A3 A4 A5 A6 A7 This is wht you think of doing Find the zeros of f ( x). Set function equl to 0. Fctor or use qudrtic eqution if
More informationTrapezoidal Rule, n = 1, x 0 = a, x 1 = b, h = b a. f (x)dx = h 2 (f (x 0) + f (x 1 )) h3
Trpezoidl Rule, n = 1, x 0 =, x 1 = b, h = b f (x)dx = h 2 (f (x 0) + f (x 1 )) h3 12 f (ξ). Simpson s Rule: n = 3, x 0 =, x 1 = +b 2, x 2 = b, h = b 2. Qudrture Rule, double node t x 1 P 3 (x)dx = f (x
More information( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).
AB Clculus Exm Review Sheet A. Preclculus Type prolems A1 Find the zeros of f ( x). This is wht you think of doing A2 A3 Find the intersection of f ( x) nd g( x). Show tht f ( x) is even. A4 Show tht f
More informationMath 100 Review Sheet
Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s
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 informationLecture 21: Order statistics
Lecture : Order sttistics Suppose we hve N mesurements of sclr, x i =, N Tke ll mesurements nd sort them into scending order x x x 3 x N Define the mesured running integrl S N (x) = 0 for x < x = i/n for
More informationMath 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8
Mth 3 Fll 0 The scope of the finl exm will include: Finl Exm Review. Integrls Chpter 5 including sections 5. 5.7, 5.0. Applictions of Integrtion Chpter 6 including sections 6. 6.5 nd section 6.8 3. Infinite
More 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 informationArithmetic Mean Derivative Based Midpoint Rule
Applied Mthemticl Sciences, Vol. 1, 018, no. 13, 65-633 HIKARI Ltd www.m-hikri.com https://doi.org/10.1988/ms.018.858 Arithmetic Men Derivtive Bsed Midpoint Rule Rike Mrjulis 1, M. Imrn, Symsudhuh Numericl
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 informationCOT4501 Spring Homework VII
COT451 Spring 1 Homework VII The ssignment is due in clss on Thursdy, April 19, 1. There re five regulr problems nd one computer problem (using MATLAB). For written problems, you need to show your work
More 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 informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued
More informationAPPLICATIONS OF THE DEFINITE INTEGRAL
APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More information. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. sin 2 (θ) =
Review of some needed Trig. Identities for Integrtion. Your nswers should be n ngle in RADIANS. rccos( 1 ) = π rccos( - 1 ) = 2π 2 3 2 3 rcsin( 1 ) = π rcsin( - 1 ) = -π 2 6 2 6 Cn you do similr problems?
More 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 informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
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 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 information