Numerical Differentiation and Integration

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

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

Chapter 12-b Integral Calculus - Extra

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

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

Chapter Gauss-Seidel Method

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

CS321. Numerical Analysis

MTH 146 Class 7 Notes

CS321. Introduction to Numerical Methods

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

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

Numerical Analysis Topic 4: Least Squares Curve Fitting

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

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 ]

Chapter Trapezoidal Rule of Integration

Numerical differentiation (Finite difference)

Gauss Quadrature Rule of Integration

FORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS

Gauss Quadrature Rule of Integration

Chapter Gauss Quadrature Rule of Integration

under the curve in the first quadrant.

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

Numerical Analysis Third Class Chemical Engineering Department University of Technology Assoc. Prof. Dr. Zaidoon Mohsin shakor

GENERALIZED OPERATIONAL RELATIONS AND PROPERTIES OF FRACTIONAL HANKEL TRANSFORM

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

Numerical Integration

Chapter 5. Curve fitting

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

Implicit Runge-Kutta method for Van der pol problem

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

Computer Programming

19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in

20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2

Sequences and summations

ZETA REGULARIZATION METHOD APPLIED TO THE CALCULATION OF DIVERGENT DIVERGENT INTEGRALS

ICS141: Discrete Mathematics for Computer Science I

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

Introduction to Numerical Integration Part II

Section 5.1/5.2: Areas and Distances the Definite Integral

Interpolation. 1. What is interpolation?

Integration by Parts for D K

COMPLEX NUMBERS AND DE MOIVRE S THEOREM

The Z-Transform in DSP Lecture Andreas Spanias

Strategies for the AP Calculus Exam

Mathematics HL and further mathematics HL formula booklet

The definite Riemann integral

Diagonally Implicit Runge-Kutta Nystrom General Method Order Five for Solving Second Order IVPs

CHAPTER 5 Vectors and Vector Space

14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y

Analele Universităţii din Oradea, Fascicula: Protecţia Mediului, Vol. XIII, 2008

Mathematics HL and further mathematics HL formula booklet

Trapezoidal Rule of Integration

CHAPTER 6 CURVE FITTINGS

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

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.

Module 2: Introduction to Numerical Analysis

Asymptotic Dominance Problems. is not constant but for n 0, f ( n) 11. 0, so that for n N f

CURVE FITTING LEAST SQUARES METHOD

Numerical Solution of Higher Order Linear Fredholm Integro Differential Equations.

CHAPTER 7 SOLVING FUZZY LINEAR FRACTIONAL PROGRAMMING PROBLEM BY USING TAYLOR'S SERIES METHOD

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

Problem Set 4 Solutions

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS

Trapezoidal Rule of Integration

Solutions Manual for Polymer Science and Technology Third Edition

ME 501A Seminar in Engineering Analysis Page 1

On Several Inequalities Deduced Using a Power Series Approach

Mathematics HL and further mathematics HL formula booklet

SOLVING INITIAL VALUE PROBLEM USING RUNGE-KUTTA 6 th ORDER METHOD

Chapter Linear Regression

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Roberto s Notes on Integral Calculus Chapter 4: Definite integrals and the FTC Section 2. Riemann sums

ENGI 4430 Numerical Integration Page 5-01

The formulae in this booklet have been arranged according to the unit in which they are first

CS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department

5 - Determinants. r r. r r. r r. r s r = + det det det

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III

ON JENSEN S AND HERMITE-HADAMARD S INEQUALITY

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Revised by Prof. Jang, CAU

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:

Global Journal of Engineering and Technology Review

UNIT I. Definition and existence of Riemann-Stieltjes

Riemann Integral Oct 31, such that

Complex Variables. Chapter 19 Series and Residues. March 26, 2013 Lecturer: Shih-Yuan Chen

xl yl m n m n r m r m r r! The inner sum in the last term simplifies because it is a binomial expansion of ( x + y) r : e +.

Section 11.5 Notes Page Partial Fraction Decomposition. . You will get: +. Therefore we come to the following: x x

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ. Численные методы. Учебно-методическое пособие

Accuplacer Elementary Algebra Study Guide

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.

SOLVING SYSTEMS OF EQUATIONS, DIRECT METHODS

6.6 Moments and Centers of Mass

More Regression Lecture Notes CE 311K - McKinney Introduction to Computer Methods Department of Civil Engineering The University of Texas at Austin

Stats & Summary

Chapter 1. Infinite Sequences and Series. 1.1 Sequences. A sequence is a set of numbers written in a definite order

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

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

Lecture 3-4 Solutions of System of Linear Equations

Application: Work. 8.1 What is Work?

SOME REMARKS ON HORIZONTAL, SLANT, PARABOLIC AND POLYNOMIAL ASYMPTOTE

Transcription:

Numerl Deretto d Itegrto

Overvew Numerl Deretto Newto-Cotes Itegrto Formuls Trpezodl rule Smpso s Rules Guss Qudrture Cheyshev s ormul

Numerl Deretto Forwrd te dvded deree Bkwrd te dvded deree Ceter te dvded deree All sed o the Tylor Seres ' h h ''!...

Forwrd Fte Deree '' ' '...! '' ' O h h h O h h h h

Forwrd Dvded Deree ' O h Oh () ( +,y + ) (, y )

Bkwrd Deree Appromto o the Frst Dervtve Epd the Tylor seres kwrds '' h h '...! '' h h '...! ' h h The error s stll O(h)

Cetered Deree Appromto o the Frst Dervtve Sutrt kwrd deree ppromto rom orwrd Tylor seres epso ' 6 ''' ' '...! '' ' O h h h h h h

() ' h O h (,y ) ( +,y + ) ( -,y - )

() () true dervtve orwrd te dvded deree ppro. () () kwrd te dvded deree ppro. etered te dvded deree ppro.

Newto-Cotes Itegrto Commo umerl tegrto sheme Bsed o the strtegy o replg omplted uto or tulted dt wth some ppromtg uto tht s esy to tegrte I d d...

Newto-Cotes Itegrto Commo umerl tegrto sheme Bsed o the strtegy o replg omplted uto or tulted dt wth some ppromtg uto tht s esy to tegrte I d d... () s th order polyoml

() () The ppromto o tegrl y the re uder - rst order polyoml - seod order polyoml 5 5 5 5

() () The ppromto o tegrl y the re uder - rst order polyoml - seod order polyoml 5 5 5 5 We lso ppromted the tegrl y usg seres o polyomls ppled pee wse.

() 5 5 6 7 8 9 A ppromto o tegrl y the re uder strght le segmets.

() 5 5 6 7 8 9 A ppromto o tegrl y the re uder strght le segmets.

Trpezodl Rule Frst o the Newto-Cotes losed tegrto ormuls Correspods to the se where the polyoml s rst order I d d A strght le e represeted s:

Trpezodl Rule I d d d Itegrte ths equto. Results the trpezodl rule. I Multple Applto o the Trpezodl Rule I ( ) d ( ) d ( ) d I h h h I

Multple Applto o the Trpezodl Rule I } } wdth verge heght The verge heght represets weghted verge o the uto vlues Note tht the teror pots re gve twe the weght o the two ed pots

Smpso s / Rule Correspods to the se where the uto s seod order polyoml I d d Desgte d s d, d estmte () s seod order Lgrge polyoml I d d... d

Smpso s / Rule Ater tegrto d lger mpulto, we get the ollowg equtos I h 6 } } wdth verge heght Multple Applto o Smpso s / Rule I d d d I j j ) ( ) ( ) (,,5..,,6..

Smpso s /8 Rule Correspods to the se where the uto s thrd order polyoml 8 h I d d I

Guss Qudrture Rule Prevously, the Trpezodl Rule ws developed y the method o udetermed oeets. The result o tht developmet s summrzed elow. ( ) d ( ) ( ) ( ) ( ) The two-pot Guss Qudrture Rule s eteso o the Trpezodl Rule ppromto where the rgumets o the uto re ot predetermed s d ut s ukows d. I the two-pot Guss Qudrture Rule, the tegrl s ppromted s I ( )d ( ) ( )

Bss o the Guss Qudrture Rule The our ukows,, d re oud y ssumg tht the ormul gves et results or tegrtg geerl thrd order polyoml,. ) ( Hee d )d (

It ollows tht )d ( Equtg Equtos the two prevous two epressos yeld

Se the ostts,,, re rtrry The our smulteous oler Equtos hve oly oe eptle soluto, ) ( d Hee Two-Pot Guss Qudrture Rule

6 Hgher Pot Guss Qudrture Formuls ) ( ) ( ) ( ) ( d s lled the three-pot Guss Qudrture Rule. The oeets,, d, d the utol rgumets,, d re lulted y ssumg the ormul gves et epressos or d 5 5 Geerl -pot rules would ppromte the tegrl ) (....... ) ( ) ( )d ( tegrtg th order polyoml

Argumets d Weghg Ftors or -pot Guss Qudrture Formuls I hdooks, oeets d rgumets gve or -pot Guss Qudrture Rule re gve or tegrls g( )d s show Tle. g( ) Tle : Weghtg tors d uto rgumets used Guss Qudrture Formuls. Pots Weghtg Ftors =. =. =.555555556 =.888888889 =.555555556 =.78585 =.65555 =.65555 =.78585 Futo Argumets = -.577569 =.577569 = -.77596669 =. =.77596669 = -.866 = -.998 =.998 =.866 7

Argumets d Weghg Ftors or -pot Guss Qudrture Formuls Tle (ot.) : Weghtg tors d uto rgumets used Guss Qudrture Formuls. Pots Weghtg Ftors 5 =.696885 =.786867 =.568888889 =.786867 5 =.696885 6 =.79 =.67657 =.67995 =.67995 5 =.67657 6 =.79 Futo Argumets = -.967986 = -.5869 =. =.5869 5 =.967986 = -.9695 = -.66986 = -.86986 =.86986 5 =.66986 6 =.9695 8

Argumets d Weghg Ftors or -pot Guss Qudrture Formuls So the tle s gve or g( )d tegrls, how does oe solve ( )d? The swer les tht y tegrl wth lmts o, e overted to tegrl wth lmts, Let mt I, the t I the t, Suh tht: 9 m

Argumets d Weghg Ftors or -pot Guss Qudrture Formuls The Hee t dt d Susttutg our vlues o, d d to the tegrl gves us ) ( dt t d

Cheyshev's Formul Cheyshev s Qudrture -pot rules would ppromte the tegrl wth ll equl weghts. ( ) d w[ ( ) ( )....... ( )] Tle : Weghtg tors d uto rgumets used Guss Qudrture Formuls. Pots Futo Argumets = -.577569 =.577569 = -.77678 =. =.77678 = -.79657 = -.87597 =. 87597 =. 79657 Pots Futo Argumets 5 = -.89787 = -.7596 =. =.7596 5 =.89787 6 = -.866688 = -.58658 = -.66655 =.66655 5 =.58658 6 =.866688

THANK YOU Sumtted y: Rkesh Kumr, (Deptt. O Mthemts), P.G.G.C.G, Se-, Chdgrh.