Gauss Quadrature Rule of Integration

Similar documents
Gauss Quadrature Rule of Integration

Chapter Gauss Quadrature Rule of Integration

Simpson s 1/3 rd Rule of Integration

Trapezoidal Rule of Integration

Trapezoidal Rule of Integration

Numerical Methods. Lecture 5. Numerical integration. dr hab. inż. Katarzyna Zakrzewska, prof. AGH. Numerical Methods lecture 5 1

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

Chapter 3 Single Random Variables and Probability Distributions (Part 2)

Calculus AB. For a function f(x), the derivative would be f '(

Time in Seconds Speed in ft/sec (a) Sketch a possible graph for this function.

Integration. antidifferentiation

] dx (3) = [15x] 2 0

The Trapezoidal Rule

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Review Exercises for Chapter 4

University of Sioux Falls. MAT204/205 Calculus I/II

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

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

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

CHAPTER 6b. NUMERICAL INTERPOLATION

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

MA10207B: ANALYSIS SECOND SEMESTER OUTLINE NOTES

50 AMC Lectures Problem Book 2 (36) Substitution Method

The practical version

For a, b, c, d positive if a b and. ac bd. Reciprocal relations for a and b positive. If a > b then a ab > b. then

MAT 403 NOTES 4. f + f =

Formula for Trapezoid estimate using Left and Right estimates: Trap( n) If the graph of f is decreasing on [a, b], then f ( x ) dx

September 13 Homework Solutions

Chapter Direct Method of Interpolation More Examples Electrical Engineering

MA Exam 2 Study Guide, Fall u n du (or the integral of linear combinations

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

Polynomials. Polynomials. Curriculum Ready ACMNA:

Core 2 Logarithms and exponentials. Section 1: Introduction to logarithms

AP CALCULUS Test #6: Unit #6 Basic Integration and Applications

Logarithms LOGARITHMS.

Improper Integrals. Introduction. Type 1: Improper Integrals on Infinite Intervals. When we defined the definite integral.

Section 4: Integration ECO4112F 2011

Autar Kaw Benjamin Rigsby. Transforming Numerical Methods Education for STEM Undergraduates

Quadrature Rules for Evaluation of Hyper Singular Integrals

4.6 Numerical Integration

Chapter 6 Techniques of Integration

8.3 THE HYPERBOLA OBJECTIVES

The Ellipse. is larger than the other.

Lecture 14: Quadrature

Review of Gaussian Quadrature method

On the Scale factor of the Universe and Redshift.

Generalization of 2-Corner Frequency Source Models Used in SMSIM

Construction of Gauss Quadrature Rules

Chapter 9 Definite Integrals

More Properties of the Riemann Integral

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

Chapter 6 Notes, Larson/Hostetler 3e

APPM 1360 Exam 2 Spring 2016

The Double Integral. The Riemann sum of a function f (x; y) over this partition of [a; b] [c; d] is. f (r j ; t k ) x j y k

H (2a, a) (u 2a) 2 (E) Show that u v 4a. Explain why this implies that u v 4a, with equality if and only u a if u v 2a.

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Calculus AB Section I Part A A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION

Solutions to Assignment 1

Solids of Revolution

AP Calculus AB Unit 4 Assessment

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

MATH 409 Advanced Calculus I Lecture 22: Improper Riemann integrals.

The Trapezoidal Rule

LINEAR ALGEBRA APPLIED

Fundamental Theorem of Calculus

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

12.4 Similarity in Right Triangles

Mat 210 Updated on April 28, 2013

First Semester Review Calculus BC

are coplanar. ˆ ˆ ˆ and iˆ

Section 4.4. Green s Theorem

Lecture Summaries for Multivariable Integral Calculus M52B

Numerical Differentiation and Integration

Chapter Direct Method of Interpolation More Examples Civil Engineering

Section 7.1 Area of a Region Between Two Curves

1 Part II: Numerical Integration

Solving Radical Equations

Math 131. Numerical Integration Larson Section 4.6

QUADRATIC EQUATION. Contents

6.5 Improper integrals

PYTHAGORAS THEOREM WHAT S IN CHAPTER 1? IN THIS CHAPTER YOU WILL:

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

Applications of Definite Integral

B.Sc. in Mathematics (Ordinary)

T b a(f) [f ] +. P b a(f) = Conclude that if f is in AC then it is the difference of two monotone absolutely continuous functions.

Mathematics Number: Logarithms

Lecture 1 - Introduction and Basic Facts about PDEs

Type 2: Improper Integrals with Infinite Discontinuities

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?

1 This question is about mean bond enthalpies and their use in the calculation of enthalpy changes.

METHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS

Math 259 Winter Solutions to Homework #9

Chapter 8.2: The Integral

7.8 IMPROPER INTEGRALS

LESSON 11: TRIANGLE FORMULAE

MATH Final Review

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

Transcription:

Guss Qudrture Rule o Integrtion Mjor: All Engineering Mjors Authors: Autr Kw, Chrlie Brker http://numerilmethods.eng.us.edu Trnsorming Numeril Methods Edution or STEM Undergrdutes /0/00 http://numerilmethods.eng.us.edu

Guss Qudrture Rule o Integrtion http://numerilmethods.eng.us.edu

Wht is Integrtion? Integrtion The proess o mesuring the re under urve. y )d ) I )d Where: ) is the integrnd lower limit o integrtion upper limit o integrtion http://numerilmethods.eng.us.edu

Two-Point Gussin Qudrture Rule 4 http://numerilmethods.eng.us.edu

http://numerilmethods.eng.us.edu 5 Bsis o the Gussin Qudrture Rule Previously, the Trpezoidl Rule ws developed y the method o undetermined oeiients. The result o tht development is summrized elow. ) ) ) ) ) d

Bsis o the Gussin Qudrture Rule The two-point Guss Qudrture Rule is n etension o the Trpezoidl Rule pproimtion where the rguments o the untion re not predetermined s nd ut s unknowns nd. In the two-point Guss Qudrture Rule, the integrl is pproimted s I )d ) ) 6 http://numerilmethods.eng.us.edu

http://numerilmethods.eng.us.edu 7 Bsis o the Gussin Qudrture Rule The our unknowns,, nd re ound y ssuming tht the ormul gives et results or integrting generl third order polynomil,. ) 0 Hene ) d )d 0 4 4 0 ) 4 4 4 0

http://numerilmethods.eng.us.edu 8 Bsis o the Gussin Qudrture Rule It ollows tht ) ) 0 0 )d Equting Equtions the two previous two epressions yield ) 4 4 4 0 ) ) 0 0 ) ) ) ) 0

http://numerilmethods.eng.us.edu 9 Bsis o the Gussin Qudrture Rule Sine the onstnts 0,,, re ritrry 4 4 4

http://numerilmethods.eng.us.edu 0 Bsis o Guss Qudrture The previous our simultneous nonliner Equtions hve only one eptle solution,

http://numerilmethods.eng.us.edu Bsis o Guss Qudrture Hene Two-Point Gussin Qudrture Rule ) ) ) d

Higher Point Gussin Qudrture Formuls http://numerilmethods.eng.us.edu

http://numerilmethods.eng.us.edu Higher Point Gussin Qudrture Formuls ) ) ) ) d is lled the three-point Guss Qudrture Rule. The oeiients,, nd, nd the untionl rguments,, nd re lulted y ssuming the ormul gives et epressions or ) d 5 5 4 4 0 Generl n-point rules would pproimte the integrl )....... ) ) )d n n integrting ith order polynomil

Arguments nd Weighing Ftors or n-point Guss Qudrture Formuls In hndooks, oeiients nd rguments given or n-point Guss Qudrture Rule re given or integrls g )d n i s shown in Tle. i g i ) Tle : Weighting tors nd untion rguments used in Guss Qudrture Formuls. Points Weighting Ftors.000000000.000000000 0.555555556 0.888888889 0.555555556 4 0.47854845 0.654555 0.654555 4 0.47854845 Funtion Arguments -0.5775069 0.5775069-0.774596669 0.000000000 0.774596669-0.866-0.998044 0.998044 4 0.866 4 http://numerilmethods.eng.us.edu

Arguments nd Weighing Ftors or n-point Guss Qudrture Formuls Tle ont.) : Weighting tors nd untion rguments used in Guss Qudrture Formuls. Points Weighting Ftors 5 0.696885 0.47868670 0.568888889 4 0.47868670 5 0.696885 6 0.7449 0.607657 0.467995 4 0.467995 5 0.607657 6 0.7449 Funtion Arguments -0.90679846-0.584690 0.000000000 4 0.584690 5 0.90679846-0.946954-0.660986-0.869860 4 0.869860 5 0.660986 6 0.946954 5 http://numerilmethods.eng.us.edu

Arguments nd Weighing Ftors or n-point Guss Qudrture So i the tle is given or g Formuls )d integrls, how does one solve )d? The nswer lies in tht ny integrl with limits o [, ] n e onverted into n integrl with limits [,] mt Let I, then t I then t, Suh tht: 6 m http://numerilmethods.eng.us.edu

Arguments nd Weighing Ftors or n-point Guss Qudrture Formuls Then Hene t d dt Sustituting our vlues o, nd d into the integrl gives us ) d t dt 7 http://numerilmethods.eng.us.edu

Emple For n integrl Rule. )d, derive the one-point Gussin Qudrture Solution The one-point Gussin Qudrture Rule is )d ) 8 http://numerilmethods.eng.us.edu

Solution The two unknowns, nd re ound y ssuming tht the ormul gives et results or integrting generl irst order polynomil, ) 0. ) d ) 0 d 0 0 ) 9 http://numerilmethods.eng.us.edu

Solution It ollows tht ) ) d Equting Equtions, the two previous two epressions yield 0 ) 0 ) 0 ) ) 0 0 http://numerilmethods.eng.us.edu

Bsis o the Gussin Qudrture Rule Sine the onstnts 0, nd re ritrry giving http://numerilmethods.eng.us.edu

Solution Hene One-Point Gussin Qudrture Rule ) d ) ) http://numerilmethods.eng.us.edu

Emple ) Use two-point Guss Qudrture Rule to pproimte the distne overed y roket rom t8 to t0 s given y 0 40000 000ln 9. 8t dt 8 40000 00t ) Find the true error, E t or prt ). ) Also, ind the solute reltive true error, or prt ). http://numerilmethods.eng.us.edu

Solution First, hnge the limits o integrtion rom [8,0] to [-,] y previous reltions s ollows 0 t )dt 0 8 0 8 0 8 8 d 9) d 4 http://numerilmethods.eng.us.edu

Solution ont) Net, get weighting tors nd untion rgument vlues rom Tle or the two point rule,. 000000000 0. 5775069. 000000000 0. 5775069 5 http://numerilmethods.eng.us.edu

Solution ont.) Now we n use the Guss Qudrture ormul 9) d 9) 9) 0. 57750 ) 9) 0. 57750 ) 9). 6495 ) 5. 5085 ) 96. 87 ) 708. 48) 058.44 m 6 http://numerilmethods.eng.us.edu

Solution ont) sine 40000. 6495 ) 000ln 9. 8. 6495 ) 40000 00. 6495 ) 96.87 40000 5. 5085 ) 000ln 9. 8 5. 5085 ) 40000 00 5. 5085 ) 708.48 7 http://numerilmethods.eng.us.edu

Solution ont) ) The true error, t ) E t E, is True Vlue Approimte Vlue 06.4 058.44.9000 m The solute reltive true error, t, is Et vlue 06.4m) 06. 4 058. 44 t 00% 06. 4 8 0.06% http://numerilmethods.eng.us.edu

Additionl Resoures For ll resoures on this topi suh s digitl udiovisul letures, primers, tetook hpters, multiple-hoie tests, worksheets in MATLAB, MATHEMATICA, MthCd nd MAPLE, logs, relted physil prolems, plese visit http://numerilmethods.eng.us.edu/topis/guss_qu drture.html

THE END http://numerilmethods.eng.us.edu