Simpson s 1/3 rd Rule of Integration

Similar documents
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

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

Gauss Quadrature Rule of Integration

Gauss Quadrature Rule of Integration

Numerical Integration

Approximate Integration

Limit of a function:

Numerical Solutions of Fredholm Integral Equations Using Bernstein Polynomials

1 Tangent Line Problem

AP Calculus Notes: Unit 6 Definite Integrals. Syllabus Objective: 3.4 The student will approximate a definite integral using rectangles.

Approximations of Definite Integrals

1. (25 points) Use the limit definition of the definite integral and the sum formulas to compute. [1 x + x2

Multiplicative Versions of Infinitesimal Calculus

Interpolation. 1. What is interpolation?

0 otherwise. sin( nx)sin( kx) 0 otherwise. cos( nx) sin( kx) dx 0 for all integers n, k.

BC Calculus Review Sheet

Crushed Notes on MATH132: Calculus

Surds, Indices, and Logarithms Radical

The Definite Integral

INTEGRATION TECHNIQUES (TRIG, LOG, EXP FUNCTIONS)

Week 13 Notes: 1) Riemann Sum. Aim: Compute Area Under a Graph. Suppose we want to find out the area of a graph, like the one on the right:

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

Topic 4 Fourier Series. Today

Fast Fourier Transform 1) Legendre s Interpolation 2) Vandermonde Matrix 3) Roots of Unity 4) Polynomial Evaluation

lecture 16: Introduction to Least Squares Approximation

Chapter Gauss Quadrature Rule of Integration

Important Facts You Need To Know/Review:

Options: Calculus. O C.1 PG #2, 3b, 4, 5ace O C.2 PG.24 #1 O D PG.28 #2, 3, 4, 5, 7 O E PG #1, 3, 4, 5 O F PG.

Eigenfunction Expansion. For a given function on the internal a x b the eigenfunction expansion of f(x):

 n. A Very Interesting Example + + = d. + x3. + 5x4. math 131 power series, part ii 7. One of the first power series we examined was. 2!

Graphing Review Part 3: Polynomials

Chapter Simpson s 1/3 Rule of Integration. ( x)

1.1 The FTC and Riemann Sums. An Application of Definite Integrals: Net Distance Travelled

Fourier Series and Applications

Unit 1 Chapter-3 Partial Fractions, Algebraic Relationships, Surds, Indices, Logarithms

UNIVERSITY OF BRISTOL. Examination for the Degrees of B.Sc. and M.Sci. (Level C/4) ANALYSIS 1B, SOLUTIONS MATH (Paper Code MATH-10006)

f ( x) ( ) dx =

Numerical Integration by using Straight Line Interpolation Formula

4. When is the particle speeding up? Why? 5. When is the particle slowing down? Why?

10.5 Power Series. In this section, we are going to start talking about power series. A power series is a series of the form

Math1242 Project I (TI 84) Name:

Numerical Integration - (4.3)

1.3 Continuous Functions and Riemann Sums

( a n ) converges or diverges.

(200 terms) equals Let f (x) = 1 + x + x 2 + +x 100 = x101 1

Mathacle. PSet Stats, Concepts In Statistics Level Number Name: Date:

Homework 2 solutions

5.3. The Definite Integral. Limits of Riemann Sums

12.2 The Definite Integrals (5.2)

Closed Newton-Cotes Integration

y udv uv y v du 7.1 INTEGRATION BY PARTS

the midpoint of the ith subinterval, and n is even for

FOURIER SERIES PART I: DEFINITIONS AND EXAMPLES. To a 2π-periodic function f(x) we will associate a trigonometric series. a n cos(nx) + b n sin(nx),

General properties of definite integrals

Pre-Calculus - Chapter 3 Sections Notes

Orthogonal functions - Function Approximation

PhysicsAndMathsTutor.com

EVALUATING DEFINITE INTEGRALS

A GENERAL METHOD FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS: THE FROBENIUS (OR SERIES) METHOD

GRAPHING LINEAR EQUATIONS. Linear Equations. x l ( 3,1 ) _x-axis. Origin ( 0, 0 ) Slope = change in y change in x. Equation for l 1.

MA123, Chapter 9: Computing some integrals (pp )

Inner Product Spaces (Chapter 5)

B. Examples 1. Finite Sums finite sums are an example of Riemann Sums in which each subinterval has the same length and the same x i

HIGHER SCHOOL CERTIFICATE EXAMINATION MATHEMATICS 3 UNIT (ADDITIONAL) AND 3/4 UNIT (COMMON) Time allowed Two hours (Plus 5 minutes reading time)

Lecture 4 Recursive Algorithm Analysis. Merge Sort Solving Recurrences The Master Theorem

REVIEW OF CHAPTER 5 MATH 114 (SECTION C1): ELEMENTARY CALCULUS

Math 3B Midterm Review

Feedback & Assessment of Your Success. 1 Calculus AP U5 Integration (AP) Name: Antiderivatives & Indefinite Integration (AP) Journal #1 3days

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

Solutions to Problem Set 7

Theorem 5.3 (Continued) The Fundamental Theorem of Calculus, Part 2: ab,, then. f x dx F x F b F a. a a. f x dx= f x x

The Reimann Integral is a formal limit definition of a definite integral

MATH 118 HW 7 KELLY DOUGAN, ANDREW KOMAR, MARIA SIMBIRSKY, BRANDEN LASKE

THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK SUMMER EXAMINATION 2005 FIRST ENGINEERING

Linford 1. Kyle Linford. Math 211. Honors Project. Theorems to Analyze: Theorem 2.4 The Limit of a Function Involving a Radical (A4)

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

[Q. Booklet Number]

Fourier Series. Topic 4 Fourier Series. sin. sin. Fourier Series. Fourier Series. Fourier Series. sin. b n. a n. sin

f(t)dt 2δ f(x) f(t)dt 0 and b f(t)dt = 0 gives F (b) = 0. Since F is increasing, this means that

Name of the Student:

Infinite Series Sequences: terms nth term Listing Terms of a Sequence 2 n recursively defined n+1 Pattern Recognition for Sequences Ex:

Review of the Riemann Integral

For students entering Honors Precalculus Summer Packet

M3P14 EXAMPLE SHEET 1 SOLUTIONS

Homework 3 solutions

Assessment Center Elementary Algebra Study Guide for the ACCUPLACER (CPT)

Limits and an Introduction to Calculus

Mathematical Notation Math Calculus & Analytic Geometry I

where i is the index of summation, a i is the ith term of the sum, and the upper and lower bounds of summation are n and 1.

Chapter 25 Sturm-Liouville problem (II)

F x = 2x λy 2 z 3 = 0 (1) F y = 2y λ2xyz 3 = 0 (2) F z = 2z λ3xy 2 z 2 = 0 (3) F λ = (xy 2 z 3 2) = 0. (4) 2z 3xy 2 z 2. 2x y 2 z 3 = 2y 2xyz 3 = ) 2

CITY UNIVERSITY LONDON

Mathematical Notation Math Calculus & Analytic Geometry I

MAS221 Analysis, Semester 2 Exercises

Chapter 5. The Riemann Integral. 5.1 The Riemann integral Partitions and lower and upper integrals. Note: 1.5 lectures

n 2 + 3n + 1 4n = n2 + 3n + 1 n n 2 = n + 1

SOLUTION OF SYSTEM OF LINEAR EQUATIONS. Lecture 4: (a) Jacobi's method. method (general). (b) Gauss Seidel method.

Particle in a Box. and the state function is. In this case, the Hermitian operator. The b.c. restrict us to 0 x a. x A sin for 0 x a, and 0 otherwise

Transcription:

Simpso s 1/3 rd Rule o Itegrtio Mjor: All Egieerig Mjors Authors: Autr Kw, Chrlie Brker Trsormig Numericl Methods Eductio or STEM Udergrdutes 1/10/010 1

Simpso s 1/3 rd Rule o Itegrtio

Wht is Itegrtio? Itegrtio The process o mesurig the re uder curve. y ( )d () I ( )d Where: () is the itegrd lower limit o itegrtio upper limit o itegrtio 3

Simpso s 1/3 rd Rule 4

Bsis o Simpso s 1/3 rd Rule Trpezoidl rule ws sed o pproimtig the itegrd y irst order polyomil, d the itegrtig the polyomil i the itervl o itegrtio. Simpso s 1/3rd rule is etesio o Trpezoidl rule where the itegrd is pproimted y secod order polyomil. Hece I ( )d ( )d Where ( ) is secod order polyomil. ( ) 0 1 5

Bsis o Simpso s 1/3 rd Rule Choose (, ( )),,, d (, ( )) s the three poits o the uctio to evlute 0, 1 d. ( ) ( ) 0 1 0 1 ( ) ( ) 0 1 6

Bsis o Simpso s 1/3 rd Rule 7 0 1 Solvig the previous equtios or 0, 1 d give ( ) ( ) ( ) 4 4 ( ) 3 ( ) ( ) 3 ( ) ( ) ( ) 4 ( )

8 Bsis o Simpso s 1/3 rd Rule The )d ( I ( ) d 1 0 3 3 1 0 3 3 3 1 0 ) (

Bsis o Simpso s 1/3 rd Rule Sustitutig vlues o 0, 1, give ( )d 6 ( ) 4 ( ) Sice or Simpso s 1/3rd Rule, the itervl [, ] is roke ito segmets, the segmet width h 9

Bsis o Simpso s 1/3 rd Rule Hece ( )d h 3 ( ) 4 ( ) Becuse the ove orm hs 1/3 i its ormul, it is clled Simpso s 1/3rd Rule. 10

Emple 1 The distce covered y rocket rom t8 to t30 is give y 30 140000 000l 9. 8t dt 8 140000 100t ) Use Simpso s 1/3rd Rule to id the pproimte vlue o ) Fid the true error, Et c) Fid the solute reltive true error, t 11

Solutio ) 30 8 ( t) dt ( ) 4 ( ) 6 30 8 6 [ ( 8 ) 4 ( 19 ) ( 30 )] 1 6 11065.7 [ 177. 667 4( 484. 7455 ) 901. 6740] m

Solutio (cot) ) The ect vlue o the ove itegrl is 30 140000 000l 9. 8t dt 8 140000 100t 11061.34 m True Error E t 11061. 34 11065. 7 4.38 m 13

Solutio (cot) )c) Asolute reltive true error, 11061. 34 11065. 7 t 100% 11061. 34 0.0396% 14

Multiple Segmet Simpso s 1/3rd Rule 15

Multiple Segmet Simpso s 1/3 rd Rule Just like i multiple segmet Trpezoidl Rule, oe c sudivide the itervl [, ] ito segmets d pply Simpso s 1/3rd Rule repetedly over every two segmets. Note tht eeds to e eve. Divide itervl [, ] ito equl segmets, hece the segmet width h ( )d 0 ( )d where 0 16

Multiple Segmet Simpso s 1/3 rd Rule () ( )d 4 ( )d 0... 4 ( )d ( )d... ( )d... 0 - Apply Simpso s 1/3rd Rule over ech itervl, ( )d ( 0 ( 0 ) 4 ( 1 ) ( ) 6 )... 4 4 ( ) ( 3 ) ( ) ( 4 ) 6... 17

Multiple Segmet Simpso s 1/3 rd Rule ( 4 ) ( 3 ) (... ( 4 ) 6 4 )... ( ) 4 ( 1 ) ) 6 ( ( ) Sice h i, 4,..., i i 18

Multiple Segmet Simpso s 1/3 rd Rule The ( )d ( 0 ) 4 ( 1 ) h 6 ( )... ( ) 4 ( 3 ) ( ) h 6 4... ( 4 ) 4 ( 3 ) h 6 ( ) ( ) 4 ( ) h 6 1 ( )... 19

Multiple Segmet Simpso s 1/3 rd Rule 0 h ( )d ( 0 ) 4{ ( 1 ) ( 3 )... ( 1 )} 3 [...] { ( ) ( )... ( )} ( )}]... 4 h 1 ( ) 4 ( ) ( ) ( ) 0 3 1 i i i i i i odd eve 1 ( 0 ) 4 ( ) ( ) ( ) 3 1 i i i i i i odd eve

Emple Use 4-segmet Simpso s 1/3rd Rule to pproimte the distce covered y rocket rom t 8 to t30 s give y 30 140000 000l 9.8t dt 8 140000 100t 1 ) Use our segmet Simpso s 1/3rd Rule to id the pproimte vlue o. ) Fid the true error, E t or prt (). c) Fid the solute reltive true error, or prt ().

Solutio ) Usig segmet Simpso s 1/3rd Rule, 30 8 h 5. 5 4 So ( t 0 ) (8) ( t 1 ) (8 5.5) (13.5) ( t ) (13.5 5.5) (19) ( t 3 ) (19 5.5) (4.5) ( t 4 ) (30)

Solutio (cot.) 1 ( t0 ) 4 ( t ) ( ) ( ) 3 1 i i i i i i t t odd eve 30 8 3 (8) 4 ( t ) ( ) 3(4) 1 i i i i i i t odd eve 1 (30) [ (8) 4 ( t ) 4 ( t ) ( t ) (30)] 1 3 3

Solutio (cot.) cot. 11 6 [ (8) 4 (13.5) 4 (4.5) (19) (30)] 11 6 [ 177.667 4(30.469) 4(676.0501) (484.7455) 901.6740] 11061.64 m 4

Solutio (cot.) ) I this cse, the true error is E t 11061.34 11061.64 0.30 m c) The solute reltive true error t 11061.34 11061.64 11061.34 100% 0.007% 5

Solutio (cot.) Tle 1: Vlues o Simpso s 1/3rd Rule or Emple with multiple segmets Approimte Vlue E t Є t 4 6 8 10 11065.7 11061.64 11061.40 11061.35 11061.34 4.38 0.30 0.06 0.01 0.00 0.0396% 0.007% 0.0005% 0.0001% 0.0000% 6

Error i the Multiple Segmet Simpso s 1/3 rd Rule The true error i sigle pplictio o Simpso s 1/3rd Rule is give s 5 ( ) (4) E t ( ζ), < ζ < 880 I Multiple Segmet Simpso s 1/3rd Rule, the error is the sum o the errors i ech pplictio o Simpso s 1/3rd Rule. The error i segmet Simpso s 1/3rd Rule is give y 7 E E 1 5 h ) 5 ( 0 ) ( 4 ) ( 4 ( ζ1 ) ( ζ1 ), 880 90 5 5 ( 4 ) ( 4 ) h ( 4 ) ( ζ ) ( ζ ), 880 90 < 0 < ζ1 < < ζ 4

E i ( Error i the Multiple Segmet Simpso s 1/3 rd Rule i 880 ) 5 ( i 1 ) ( 4 ) ( ζ i 5 ) h 90 ( 4 ) ( ζ 5 ( 4 ) 5 ( 4 ) E ζ h ( ) 1 880 1 4, ζ 90 1... i ), ( i 1 ) < ζi < i 4 < ζ < 1 8 ( E 880 ) 5 4 ζ h ζ 5 (4), < ζ < 90

Error i the Multiple Segmet Simpso s 1/3 rd Rule Hece, the totl error i Multiple Segmet Simpso s 1/3rd Rule is E t E i i 1 5 h 90 i 1 (4) ( ζ i ) ( ) 5 90 5 i 1 (4) ( ζ i ) (4) 5 ( ζi ) ( ) i 1 4 90 9

Error i the Multiple Segmet Simpso s 1/3 rd Rule The term i 1 (4) ( ζ i ) is pproimte verge vlue o (4) ( ), < < Hece E t ( ) 4 90 5 (4) 30 where (4) i 1 (4) ( ζ i )

Additiol Resources For ll resources o this topic such s digitl udiovisul lectures, primers, tetook chpters, multiple-choice tests, worksheets i MATLAB, MATHEMATICA, MthCd d MAPLE, logs, relted physicl prolems, plese visit /topics/simpsos_ 13rd_rule.html

THE END