Chapter Trapezoidal Rule of Integration

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

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

Mathematically, integration is just finding the area under a curve from one point to another. It is b

TEACHERS ASSESS STUDENT S MATHEMATICAL CREATIVITY COMPETENCE IN HIGH SCHOOL

4. Runge-Kutta Formula For Differential Equations

4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula

Integral Equations and their Relationship to Differential Equations with Initial Conditions

Numerical Differentiation and Integration

Chapter 12-b Integral Calculus - Extra

Laplace Transform. Definition of Laplace Transform: f(t) that satisfies The Laplace transform of f(t) is defined as.

Trapezoidal Rule of Integration

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

MTH 146 Class 7 Notes

STOCHASTIC CALCULUS I STOCHASTIC DIFFERENTIAL EQUATION

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

Calculation of Effective Resonance Integrals

ASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE EQUATIONS ON DISCRETE REAL TIME SCALES

Trapezoidal Rule of Integration

3/20/2013. Splines There are cases where polynomial interpolation is bad overshoot oscillations. Examplef x. Interpolation at -4,-3,-2,-1,0,1,2,3,4

International Journal of Pure and Applied Sciences and Technology

Observations on the transcendental Equation

Simpson s 1/3 rd Rule of Integration

Chapter Gauss-Seidel Method

Modeling and Predicting Sequences: HMM and (may be) CRF. Amr Ahmed Feb 25

Chapter Gauss Quadrature Rule of Integration

Numerical Methods using the Successive Approximations for the Solution of a Fredholm Integral Equation

X-Ray Notes, Part III

Calculating Exact Transitive Closure for a Normalized Affine Integer Tuple Relation

The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

Mechanics of Materials CIVL 3322 / MECH 3322

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X

4.1 Schrödinger Equation in Spherical Coordinates

Advanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University

Eurasian International Center of Theoretical Physics, Eurasian National University, Astana , Kazakhstan

Numerical Analysis Topic 4: Least Squares Curve Fitting

Modified Taylor's Method and Nonlinear Mixed Integral Equation

The z-transform. LTI System description. Prof. Siripong Potisuk

Xidian University Liu Congfeng Page 1 of 49

Representation of Solutions of Linear Homogeneous Caputo Fractional Differential Equations with Continuous Variable Coefficients

Approximate Integration. Left and Right Endpoint Rules. Midpoint Rule = 2. Riemann sum (approximation to the integral) Left endpoint approximation

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

Solution set Stat 471/Spring 06. Homework 2

Solution 2. Since. 1. Sums and Products

SIMULATING THE SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL DIFFERENTIAL EQUATIONS BY BLOCK-PULSE WAVELETS

P-Convexity Property in Musielak-Orlicz Function Space of Bohner Type

SOME REMARKS ON HORIZONTAL, SLANT, PARABOLIC AND POLYNOMIAL ASYMPTOTE

Gauss Quadrature Rule of Integration

The Poisson Process Properties of the Poisson Process

Review for the Midterm Exam.

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

Numerical Integration

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]

Optimality of Strategies for Collapsing Expanded Random Variables In a Simple Random Sample Ed Stanek

1 4 6 is symmetric 3 SPECIAL MATRICES 3.1 SYMMETRIC MATRICES. Defn: A matrix A is symmetric if and only if A = A, i.e., a ij =a ji i, j. Example 3.1.

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

Accuplacer Elementary Algebra Study Guide

Appendix A Examples for Labs 1, 2, 3 1. FACTORING POLYNOMIALS

Some Unbiased Classes of Estimators of Finite Population Mean

under the curve in the first quadrant.

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

Differential Equation of Eigenvalues for Sturm Liouville Boundary Value Problem with Neumann Boundary Conditions

( ) 2 3 ( ) I. Order of operations II. Scientific Notation. Simplify. Write answers in scientific notation. III.

ME 501A Seminar in Engineering Analysis Page 1

UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM

CHAPTER 4 RADICAL EXPRESSIONS

Exercise # 2.1 3, 7, , 3, , -9, 1, Solution: natural numbers are 3, , -9, 1, 2.5, 3, , , -9, 1, 2 2.5, 3, , -9, 1, , -9, 1, 2.

CS321. Numerical Analysis

Chapter 2. Review of Hydrodynamics and Vector Analysis

Suggested Solution for Pure Mathematics 2011 By Y.K. Ng (last update: 8/4/2011) Paper I. (b) (c)

A Class of Lobatto Methods of Order 2s

1. Consider an economy of identical individuals with preferences given by the utility function

Pre-Calculus - Chapter 3 Sections Notes

AML710 CAD LECTURE 12 CUBIC SPLINE CURVES. Cubic Splines Matrix formulation Normalised cubic splines Alternate end conditions Parabolic blending

A L A BA M A L A W R E V IE W

8. Queueing systems lect08.ppt S Introduction to Teletraffic Theory - Fall

Section 1.3 Triangles

Density estimation III.

Gauss Quadrature Rule of Integration

Project 3: Using Identities to Rewrite Expressions

ON THE EXTENSION OF WEAK ARMENDARIZ RINGS RELATIVE TO A MONOID

Problem Set 4 Solutions

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Data Compression Techniques (Spring 2012) Model Solutions for Exercise 4

Through the fractional Riemann Liouville integral x

Introduction to Matrix Algebra

The Linear Regression Of Weighted Segments

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Numerical differentiation (Finite difference)

Lecture 02: Bounding tail distributions of a random variable

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Differential Entropy 吳家麟教授

Introduction to Neural Networks Computing. CMSC491N/691N, Spring 2001

Mathematics. Area under Curve.

02/15/04 INTERESTING FINITE AND INFINITE PRODUCTS FROM SIMPLE ALGEBRAIC IDENTITIES

CONTRIBUTIONS TO THE STUDY OF THE PASSING THROUGH THE RESONANCE OF THE LINEAR SYSTEMS HAVING A FINITE NUMBER OF DEGREES OF FREEDOM

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

G8-11 Congruence Rules

Chapter Unary Matrix Operations

LOCUS 1. Definite Integration CONCEPT NOTES. 01. Basic Properties. 02. More Properties. 03. Integration as Limit of a Sum

Transcription:

Cper 7 Trpezodl Rule o Iegro Aer redg s per, you sould e le o: derve e rpezodl rule o egro, use e rpezodl rule o egro o solve prolems, derve e mulple-segme rpezodl rule o egro, 4 use e mulple-segme rpezodl rule o egro o solve prolems, d 5 derve e ormul or e rue error e mulple-segme rpezodl rule o egro W s egro? Iegro s e proess o mesurg e re uder uo ploed o grp Wy would we w o egre uo? Amog e mos ommo emples re dg e veloy o ody rom elero uo, d dspleme o ody rom veloy uo Trougou my egeerg elds, ere re w somemes seems lke ouless pplos or egrl lulus You red ou some o ese pplos Cpers 7A-7G Somemes, e evluo o epressos volvg ese egrls eome dug, o deerme For s reso, wde vrey o umerl meods s ee developed o smply e egrl Here, we wll dsuss e rpezodl rule o ppromg egrls o e orm were I d s lled e egrd, lower lm o egro upper lm o egro W s e rpezodl rule? Te rpezodl rule s sed o e Newo-Coes ormul oe ppromes e egrd y order polyoml, e e egrl o e uo s ppromed y 7

7 Cper 7 e egrl o order polyoml Iegrg polyomls s smple d s sed o e lulus ormul Fgure Iegro o uo d, So we w o pprome e egrl I d o d e vlue o e ove egrl, oe ssumes were 4 were s order polyoml Te rpezodl rule ssumes, s, ppromg e egrl y ler polyoml srg le, d d Dervo o e Trpezodl Rule Meod : Derved rom Clulus d d d 5

Trpezodl Rule 7 Bu w s d? Now oe ooses,, d, s e wo pos o pprome y srg le rom o, 6 7 Solvg e ove wo equos or d, 8 Hee rom Equo 5, d 8 9 Meod : Also Derved rom Clulus lso e ppromed y usg Newo s dvded deree polyoml s Hee d d d

74 Cper 7 Ts gves e sme resul s Equo euse ey re jus dere orms o wrg e sme polyoml Meod : Derved rom Geomery Te rpezodl rule lso e derved rom geomery Look Fgure Te re uder e urve s e re o rpezod Te egrl rpezod o Are d Sum o leg o prllel sdesperpedulr dse ewee prllel sdes Fgure Geomer represeo o rpezodl rule Meod 4: Derved rom Meod o Coees Te rpezodl rule lso e derved y e meod o oees Te ormul d were

Trpezodl Rule 75 Fgure Are y meod o oees Te erpreo s s evlued pos d, d e uo evluo s gve weg o Geomerlly, Equo s looked s e re o rpezod, wle Equo s vewed s e sum o e re o wo regles, s sow Fgure How oe derve e rpezodl rule y e meod o oees? Assume d 4 Le e rg d sde e e epresso or egrls o d d d, s, e ormul wll e lso e e or ler omos o d, s, or d 5 d 6 Solvg e ove wo equos gves 7 Hee

76 Cper 7 d 8 Meod 5: Aoer ppro o e Meod o Coees Te rpezodl rule lso e derved y e meod o oees y oer ppro d Assume d 9 Le e rg d sde e e or egrls o e orm d So d Bu we w d o gve e sme resul s Equo or d Hee rom Equos d, Se d re rrry or geerl srg le Ag, solvg e ove wo equos gves 4

Trpezodl Rule 77 Tereore d 5 Emple Hum vso s e remrkle ly o er D spes rom D mges Te rgug queso s: we reple some o ese les o ompuer? Yes, e doe d o do s, egro o veor elds s requred Te ollowg egrl eeds o egred Were, I d, < < 6 9688, 7 < < 796 8487 9 6778, Use sgle segme Trpezodl rule o d e vlue o e egrl Fd e rue error, E, or pr Fd e solue relve rue error or pr 7 Soluo I, were, < < 6 9688 796, 7 < < 7 7 I 85 8487 9 6778, 7 6 9688 796 8487 9 6778 Te e vlue o e ove egrl s oud usg Mple or lulg e rue error d relve rue error

78 Cper 7 I d 679 so e rue error s E True Vlue Approme Vlue 679 85 664 Te solue relve rue error,, would e e True Error % True Vlue 679 679 4 % 859 % Mulple-Segme Trpezodl Rule I Emple, e rue error usg sgle segme rpezodl rule ws lrge We dvde e ervl [,] o [,5] d [ 5,] ervls d pply e rpezodl rule over e segme Hee d 5 d d d 5 5 5 5 5 6 5 9688 5 796 5 8487 7845 6 9688 796 8487 7 5 9 6778 9 6778 7845 7845 7 d 5 5 64975 Te rue error, E s E 679 64975 745

Trpezodl Rule 79 Te rue error ow s redued rom 64 6 o 745 Eedg s proedure o dvdg ], [ o equl segmes d pplyg e rpezodl rule over e segme, e sum o e resuls oed or e segme s e pprome vlue o e egrl Dvde o equl segmes s sow Fgure 4 Te e wd o e segme s 6 Te egrl I e roke o egrls s d I d d d d 7 Fgure 4 Mulple 4 segme rpezodl rule Applyg rpezodl rule Equo 7 o e segme gves [ ] d [ ] [ ] [ ]

7 Cper 7 8 Emple Hum vso s e remrkle ly o er D spes rom D mges Te rgug queso s: we reple some o ese les o ompuer? Yes, e doe d o do s, egro o veor elds s requred Te ollowg egrl eeds o egred d I were 7, 7 9 6778, 8487 796 9688 6 < < < <, Use -segme Trpezodl rule o d e vlue o e egrl Fd e rue error, E, or pr Fd e solue relve rue error,,or pr Soluo I 5

Trpezodl Rule 7 Se I, < < 6 9688, 7 < < 796 5 4 [ 5 ] 4 784 7 4 6497 [ ] 8487 9 6778, 7 6 9688 796 8487 9 6778 7 9688 784 6 5 5 796 5 8487 5 9 6778 Te e vlue o e ove egrl s oud usg Mple or lulg e rue error d relve rue error I d 679 so e rue error s E True Vlue Approme Vlue 679 6497 749 Te solue relve rue error,, would e e True Error % True Vlue 679 6497 679 44494 % %

7 Cper 7 Tle Vlues oed usg mulple-segme Trpezodl rule or, < < 6 9688 796 8487 9 6778, 7, 7 < < Vlue E % % 85 664 4 --- 6498 749 44494 4 6 5465 8 497 4 66 47 468 85 5 5847 65 896 879 6 77769 6977 795 487 7 458 865 44 8866 8 55754 588 8885 7 Emple Use e mulple-segme rpezodl rule o d e re uder e urve e rom o Soluo Usg wo segmes, we ge 5 e 5 5 9 5 e 6 e I 5 [ 5 ] 4 [ 9 6] 5 57 4

Trpezodl Rule 7 So w s e rue vlue o s egrl? d 4659 e Mkg e solue relve rue error 4659 555 4659 7956% Wy s e rue vlue so r wy rom e pprome vlues? Jus ke look Fgure 5 As you see, e re uder e rpezods ye, ey relly look lke rgles ow overs smll poro o e re uder e urve As we dd more segmes, e ppromed vlue qukly pproes e rue vlue Fgure 5 -segme rpezodl rule ppromo Tle Vlues oed usg mulple-segme rpezodl rule or d e Approme E Vlue 68 459 9974% 555 965 7955% 4 76 75978 8% 8 74 9546 797% 6 47 4887 98% 457 495% 64 468 5 4%

74 Cper 7 Emple 4 Use mulple-segme rpezodl rule o d I d Soluo We o use e rpezodl rule or s egrl, s e vlue o e egrd s e However, s kow dsouy urve wll o ge e re uder We ssume y vlue or e uo Te lgorm o dee e uo so we use e mulple-segme rpezodl rule s gve elow Fuo I Te I Te ^ 5 Ed Fuo Bslly, we re jus ssgg e uo vlue o zero Everywere else, e uo s ouous Ts mes e rue vlue o our egrl wll e jus rue Le s see w ppes usg e mulple-segme rpezodl rule Usg wo segmes, we ge 77 I [ ] 4 [ 77] 4 56 So w s e rue vlue o s egrl? d 884 Tus mkg e solue relve rue error

Trpezodl Rule 75 884 56 884 545% Tle Vlues oed usg mulple-segme rpezodl rule or Approme E Vlue 54 474 54% 4 79 6 664% 8 97 7 585% 6 56 86% 46 65 9% 64 57 58 98% 8 646 8 6454% 56 699 9 4564% 5 77 9 7% 4 764 64 8% 48 78 45 6% 496 796 4% d Error Mulple-segme Trpezodl Rule Te rue error or sgle segme Trpezodl rule s gve y E "ζ, < ζ < Were ζ s some po [, ] W s e error e e mulple-segme rpezodl rule? I wll e smply e sum o e errors rom e segme, were e error e segme s o e sgle segme rpezodl rule Te error e segme s [ ] E " ζ, < ζ < " ζ [ ] E " ζ, < ζ < " ζ [ ] E " ζ, < ζ <

76 Cper 7 " ζ [{ } { } ] E " ζ, < ζ < " ζ [ { } ] E " ζ, < ζ < " ζ Hee e ol error e mulple-segme rpezodl rule s E E " ζ " ζ " ζ " ζ Te erm dervve ", < < Hee s pprome verge vlue o e seod " ζ E I Tle 4, e pprome vlue o e egrl 4 l 98 d 4 8 s gve s uo o e umer o segmes You vsulze s e umer o segmes re douled, e rue error ges ppromely qurered

Trpezodl Rule 77 Tle 4 Vlues oed usg mulple-segme rpezodl rule or 4 l 98 d 4 8 Approme E % Vlue % 66-5 85 54 4-5 47 594 8 74-75 56 6 65-4 66 4 For emple, or e -segme rpezodl rule, e rue error s -5, d qurer o error s -55 T s lose o e rue error o -48 or e 4-segme rpezodl rule C you swer e queso wy s e rue error o ely -55? How does s ormo elp us umerl egro? You wll d ou s orms e ss o Romerg egro sed o e rpezodl rule, were we use e rgume rue error ges ppromely qurered we e umer o segmes s douled Romerg egro sed o e rpezodl rule s ompuolly more ee usg e rpezodl rule y sel developg uom egro seme INTEGRATION Top Trpezodl Rule Summry Teook oes o e rpezodl rule o egro Mjor Compuer Egeerg Auors Aur Kw, Mel Keels De Novemer, We Se p://umerlmeodsegusedu