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

Similar documents
Chapter Trapezoidal Rule of Integration

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

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

TEACHERS ASSESS STUDENT S MATHEMATICAL CREATIVITY COMPETENCE IN HIGH SCHOOL

Integral Equations and their Relationship to Differential Equations with Initial Conditions

Numerical Differentiation and Integration

MTH 146 Class 7 Notes

Chapter 12-b Integral Calculus - Extra

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

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

Simpson s 1/3 rd Rule of Integration

STOCHASTIC CALCULUS I STOCHASTIC DIFFERENTIAL EQUATION

Calculation of Effective Resonance Integrals

Chapter Gauss-Seidel Method

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

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

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

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

Chapter Gauss Quadrature Rule of Integration

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

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

Trapezoidal Rule of Integration

Chapter 2. Review of Hydrodynamics and Vector Analysis

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

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.

Trapezoidal Rule of Integration

CS321. Introduction to Numerical Methods

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

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

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

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

CS321. Numerical Analysis

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

Observations on the transcendental Equation

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

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

Solution set Stat 471/Spring 06. Homework 2

Numerical Analysis Topic 4: Least Squares Curve Fitting

The Poisson Process Properties of the Poisson Process

Solution 2. Since. 1. Sums and Products

Modified Taylor's Method and Nonlinear Mixed Integral Equation

Conquering kings their titles take ANTHEM FOR CONGREGATION AND CHOIR

NEIGHBOURHOODS OF A CERTAIN SUBCLASS OF STARLIKE FUNCTIONS. P. Thirupathi Reddy. E. mail:

Chapter 2 Intro to Math Techniques for Quantum Mechanics

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

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

Derivation of the Metal-Semiconductor Junction Current

International Journal of Pure and Applied Sciences and Technology

X-Ray Notes, Part III

Expectation and Moments

Review for the Midterm Exam.

Mechanics of Materials CIVL 3322 / MECH 3322

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Numerical Integration

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

Matrix. Definition 1... a1 ... (i) where a. are real numbers. for i 1, 2,, m and j = 1, 2,, n (iii) A is called a square matrix if m n.

The Linear Regression Of Weighted Segments

A Class of Lobatto Methods of Order 2s

(1) Cov(, ) E[( E( ))( E( ))]

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

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Xidian University Liu Congfeng Page 1 of 49

The Z-Transform in DSP Lecture Andreas Spanias

Quantum Mechanics II Lecture 11 Time-dependent perturbation theory. Time-dependent perturbation theory (degenerate or non-degenerate starting state)

Calculating Exact Transitive Closure for a Normalized Affine Integer Tuple Relation

Isotropic Non-Heisenberg Magnet for Spin S=1

4.1 Schrödinger Equation in Spherical Coordinates

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

Invert and multiply. Fractions express a ratio of two quantities. For example, the fraction

INTERNATIONAL JOURNAL OF ADVANCED RESEARCH IN ENGINEERING AND TECHNOLOGY (IJARET) AN APPLIED TWO-DIMENSIONAL B-SPLINE MODEL FOR INTERPOLATION OF DATA

Dynamics of Structures

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

Density estimation III.

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers

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 ]

Introduction to Numerical Integration Part II

SOME REMARKS ON HORIZONTAL, SLANT, PARABOLIC AND POLYNOMIAL ASYMPTOTE

Shuai Dong. Isaac Newton. Gottfried Leibniz

ECE Microwave Engineering. Fall Prof. David R. Jackson Dept. of ECE. Notes 10. Waveguides Part 7: Transverse Equivalent Network (TEN)

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

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

Objective of curve fitting is to represent a set of discrete data by a function (curve). Consider a set of discrete data as given in table.

Boyce/DiPrima 9 th ed, Ch 7.6: Complex Eigenvalues

Through the fractional Riemann Liouville integral x

under the curve in the first quadrant.

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

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

Week 8 Lecture 3: Problems 49, 50 Fourier analysis Courseware pp (don t look at French very confusing look in the Courseware instead)

Introduction to Matrix Algebra

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

Chapter Unary Matrix Operations

On the Existence and uniqueness for solution of system Fractional Differential Equations

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

Interpolation. 1. What is interpolation?

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

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

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

CURVE FITTING LEAST SQUARES METHOD

Transcription:

Cper 7. Smpso s / Rule o Iegro Aer redg s per, you sould e le o. derve e ormul or Smpso s / rule o egro,. use Smpso s / rule o solve egrls,. develop e ormul or mulple-segme Smpso s / rule o egro,. use mulple-segme Smpso s / rule o egro o solve egrls, d. derve e rue error ormul or mulple-segme Smpso s / rule. 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 7.A-7.G. 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 Smpso s / rule o egrl ppromo, w mproves upo e ury o e rpezodl rule. Here, we wll dsuss e Smpso s / rule o ppromg egrls o e orm were I s lled e egrd, lower lm o egro upper lm o egro Smpso s / Rule Te rpezodl rule ws sed o ppromg e egrd y rs order polyoml, d e egrg e polyoml over ervl o egro. Smpso s / rule s 7..

7.. Cper 7. eeso o Trpezodl rule were e egrd s ppromed y seod order polyoml. Fgure Iegro o uo Meod : Hee I were s seod order polyoml gve y. Coose,,,, d, s e ree pos o e uo o evlue, d. Solvg e ove ree equos or ukows,, d gve

Smpso s / Rule o Iegro 7.. Te I Susug vlues o, d gve Se or Smpso / rule, e ervl [ ], s roke o segmes, e segme wd Hee e Smpso s / rule s gve y Se e ove orm s / s ormul, s lled Smpso s / rule. Meod : Smpso s / rule lso e derved y ppromg y seod order polyoml usg Newo s dvded deree polyoml s were

7.. Cper 7. Iegrg Newo s dvded deree polyoml gves us Susug vlues o,, d o s equo yelds e sme resul s eore Meod : Oe ould eve use e Lgrge polyoml o derve Smpso s ormul. Noe y meod o ree-po qudr erpolo e used o ompls s sk. I s se, e erpolg uo eomes Iegrg s uo ges

Smpso s / Rule o Iegro 7.. Beleve or o, smplyg d org s lrge epresso yelds you e sme resul s eore. Meod : Smpso s / rule lso e derved y e meod o oees. Assume Le e rg-d sde e e epresso or e egrls,, d. Ts mples e rg d sde wll e e epressos or egrls o y ler omo o e ree egrls or geerl seod order polyoml. Now

7.. Cper 7. Solvg e ove ree equos or, d gve Ts gves Te egrl rom e rs meod e vewed s e re uder e seod order polyoml, wle e equo rom Meod e vewed s e sum o e res o ree regles. Emple I emp o udersd e mesm o e depolrzo proess uel ell, elero-ke model or med oyge-meol urre o plum ws developed e lorory FAMU. A very smpled model o e reo developed suggess uol relo egrl orm. To d e me requred or % o e oyge o e osumed, e me, s T s gve y T.. 7...7 Use Smpso s / rule o d e me requred or % o e oyge o e osumed.

Smpso s / Rule o Iegro 7..7 Fd e rue error, E, or pr. Fd e solue relve rue error,, or pr. Soluo T...9 7.7.. 7.7....8.. 7.7.......7.9..9..9 7.89 T....9.. [.8.89. ] 9 s [ ] Te e vlue o e ove egrl s, 7..7. T...9 s so e rue error s E True Vlue Approme Vlue.9 9.

7..8 Cper 7. Te solue relve rue error,, would e e True Vlue True Error.9..7 % Mulple-segme Smpso s / Rule Jus lke mulple-segme rpezodl rule, oe sudvde e ervl [ ], o segmes d pply Smpso s / rule repeedly over every wo segmes. Noe eeds o e eve. Dvde ervl [ ], o equl segmes, so e segme wd s gve y. Now were... Apply Smpso s /rd Rule over e ervl,... Se...,,, e...

Smpso s / Rule o Iegro 7..9 [ {... } {... } ] odd eve odd eve Emple I emp o udersd e mesm o e depolrzo proess uel ell, elero-ke model or med oyge-meol urre o plum ws developed e lorory FAMU. A very smpled model o e reo developed suggess uol relo egrl orm. To d e me requred or % o e oyge o e T s s gve y osumed, e me, T..7... 7 Use our-smpso s / Rule o d e me requred or % o e oyge o e osumed. Fd e rue error, E, or pr. Fd e solue relve rue error,, or pr. Soluo T odd eve.....

7.. Cper 7..7.. 7 So. 7.7....8.....7 7.7.7..7.799..7.7..9 7.7.9..9.89..9.9..7 7.7.7..7..7.9. 7.7...... T odd eve... odd. eve [.. ].

Smpso s / Rule o Iegro 7......7.7.9..8.799.9.89. 9 s Te e vlue o e ove egrl s..7. T.. 7.9 s so e rue error s E True Vlue Approme Vlue.9 9.988 Te solue relve rue error,, would e e True Error True Vlue.988.9. % Tle Vlues o Smpso s / Rule or Emple w mulple segmes. Approme Vlue E 9. 9.988 9. 8 9. 9. %.7..9 7.8.98

7.. Cper 7. Error Mulple-segme Smpso s / rule Te rue error sgle pplo o Smpso s /rd Rule s gve y E ζ, < ζ < 88 I mulple-segme Smpso s / rule, e error s e sum o e errors e pplo o Smpso s / rule. Te error e segmes Smpso s /rd Rule s gve y E ζ, < ζ < 88 ζ 9 E ζ, < ζ < 88 ζ 9 : E ζ, < ζ < 88 ζ 9 : E, < < 88 ζ ζ ζ 9 E ζ, < ζ < 88 Hee, e ol error e mulple-segme Smpso s / rule s ζ 9 E E Te e rue error epresso sds or e our dervve o e uo.

Smpso s / Rule o Iegro 7.. 9 9 ζ ζ ζ 9 ζ Te erm s pprome verge vlue o E 9 were ζ, < <. Hee INTEGRATION Top Smpso s / rule Summry Teook oes o Smpso s / rule Mjor Ceml Egeerg Auors Aur Kw, Mel Keels De Novemer 8, We Se p://umerlmeods.eg.us.edu