Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

Similar documents
OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS

Arithmetic Mean Derivative Based Midpoint Rule

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

QUADRATURE is an old-fashioned word that refers to

New Expansion and Infinite Series

Lab 11 Approximate Integration

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

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

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

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir

Interpolation. Gaussian Quadrature. September 25, 2011

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

A Generalized Inequality of Ostrowski Type for Twice Differentiable Bounded Mappings and Applications

Lecture 14: Quadrature

Undergraduate Research

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Review of Gaussian Quadrature method

1 The Lagrange interpolation formula

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

3.4 Numerical integration

Section 6.1 Definite Integral

Suppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.

Math& 152 Section Integration by Parts

METHODS OF APPROXIMATING THE RIEMANN INTEGRALS AND APPLICATIONS

NUMERICAL INTEGRATION

Lecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014

AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION

Numerical Analysis: Trapezoidal and Simpson s Rule

Research Article Composite Gauss-Legendre Formulas for Solving Fuzzy Integration

B.Sc. in Mathematics (Ordinary)

Review of Calculus, cont d

III. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9

Application Chebyshev Polynomials for Determining the Eigenvalues of Sturm-Liouville Problem

Hermite-Hadamard type inequalities for harmonically convex functions

The Shortest Confidence Interval for the Mean of a Normal Distribution

Lecture 12: Numerical Quadrature

Numerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden

Section 6.1 INTRO to LAPLACE TRANSFORMS

Construction of Gauss Quadrature Rules

Adomian Decomposition Method with Green s. Function for Solving Twelfth-Order Boundary. Value Problems

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

Definite integral. Mathematics FRDIS MENDELU

Composite Mendeleev s Quadratures for Solving a Linear Fredholm Integral Equation of The Second Kind

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Numerical Integration

Midpoint Approximation

Math 131. Numerical Integration Larson Section 4.6

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

Trapezoidal Rule, n = 1, x 0 = a, x 1 = b, h = b a. f (x)dx = h 2 (f (x 0) + f (x 1 )) h3

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Tangent Line and Tangent Plane Approximations of Definite Integral

Orthogonal Polynomials

An optimal 3-point quadrature formula of closed type and error bounds

THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS

Advanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.

Numerical Integration

New Integral Inequalities of the Type of Hermite-Hadamard Through Quasi Convexity

Numerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b

p(x) = 3x 3 + x n 3 k=0 so the right hand side of the equality we have to show is obtained for r = b 0, s = b 1 and 2n 3 b k x k, q 2n 3 (x) =

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University

Lecture 14 Numerical integration: advanced topics

Lecture 23: Interpolatory Quadrature

New Integral Inequalities for n-time Differentiable Functions with Applications for pdfs

Math 1B, lecture 4: Error bounds for numerical methods

APPROXIMATE INTEGRATION

Sections 5.2: The Definite Integral

Numerical integration

Journal of Inequalities in Pure and Applied Mathematics

Point-based methods for estimating the length of a parametric curve

A Bernstein polynomial approach for solution of nonlinear integral equations

Lecture 20: Numerical Integration III

On Second Derivative-Free Zero Finding Methods

Research Article The Modified Trapezoidal Rule for Computing Hypersingular Integral on Interval

COT4501 Spring Homework VII

1 Error Analysis of Simple Rules for Numerical Integration

Unit Six AP Calculus Unit 6 Review Definite Integrals. Name Period Date NON-CALCULATOR SECTION

Some estimates on the Hermite-Hadamard inequality through quasi-convex functions

8 J. Roumeliotis, P. Cerone, S.S. Drgomir In this pper, we extend the ove result nd develop n Ostrowski-type inequlity for weighted integrls. Applicti

Keywords : Generalized Ostrowski s inequality, generalized midpoint inequality, Taylor s formula.

Numerical quadrature based on interpolating functions: A MATLAB implementation

The Riemann Integral

Research Article Fejér and Hermite-Hadamard Type Inequalities for Harmonically Convex Functions

A Note on Feng Qi Type Integral Inequalities

f (a) + f (b) f (λx + (1 λ)y) max {f (x),f (y)}, x, y [a, b]. (1.1)

The Trapezoidal Rule

Numerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1

5: The Definite Integral

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

Czechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction

INEQUALITIES FOR TWO SPECIFIC CLASSES OF FUNCTIONS USING CHEBYSHEV FUNCTIONAL. Mohammad Masjed-Jamei

CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES

Fredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel

7. Numerical evaluation of definite integrals

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

4.6 Numerical Integration

Efficient Computation of a Class of Singular Oscillatory Integrals by Steepest Descent Method

A Companion of Ostrowski Type Integral Inequality Using a 5-Step Kernel with Some Applications

Transcription:

IOSR Journl of Mthemtics (IOSR-JM) e-issn: - p-issn: 9-X. Volume Issue Ver. IV (My. - Jun. 0) PP - www.iosrjournls.org Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture T. Rmchndrn D.Udykumr nd R.Priml Deprtment of Mthemtics M.V.Muthih Government Arts College for Women Dindigul - 00 Tmil Ndu Indi. Deprtment of Mthemtics Government Arts College (Autonomous) Krur - 9 00 Tmil Ndu Indi. Astrct: A New method of evlution of Numericl integrtion y using Hrmonic Men derivtive - sed closed Newton cotes qudrture rule (HMDCNC) is presented in which the Hrmonic men vlue is used for Computing the function derivtive. It hs shown tht the proposed rule gives increse of single order of precision over the existing closed Newton cotes rule. The error terms re lso otined y using the concept of precision nd re compred with the existing methods. Finlly the ccurcy of the proposed rule is nlyzed using Numericl Exmples nd the results re compred with the existing methods. Keyword: Closed Newton-Cotes formul Error terms Hrmonic Men Derivtive Numericl Exmples Numericl Integrtion. I. Introduction Numericl integrtion is the study of how the numericl vlue of n integrl cn e found. It hs severl pplictions in the field of sttistics. In Sttistics it is used to evlute distriution functions nd other quntities. Mny recent Sttisticl Methods re dependent especilly on multiple integrtion possily in very high dimensions[].in the field of Mthemtics One of the most common method for the evlution of numericl integrtion is qudrture rule given y n w i f x i () i=0 Where there re (n+) distinct points =x 0 < x < < x n = x i =x 0 +ih i=0..n h = nd (n+) weights n w 0 w... w n.these w i cn e otined y using the method sed on the precision of qudrture formul. Select the vlues for w i i=0...n. so tht the error is zero tht is E n f = n w i f x i = 0 for f x = x j j = 0 n () i=0 Definition. An integrtion method of the form ( ) is sid to e of order P if it produces exct results (En[f] = 0) for ll polynomils of degree less thn or equl to P []. The list of Closed Newton cotes formuls (CNC) tht depend on the integer vlue of n re given elow When n= : Trpezoidl rule = f + f() When n =: Simpson's / rd rule f x dx = f + f + When n =: Simpson's / th rule = f + f + () f"(ξ) where ξ () + f + f + 0 f ξ where ξ () + f () 0 f ξ where ξ () DOI: 0.9/-00 www.iosrjournls.org Pge

When n= : Boole's rule = Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture f + f + + f + + f + + f f ξ where ξ () It is known tht the degree of precision is n+ for even vlue of n nd n for odd vlue of n. There re so mny methods re ville to increse the order of ccurcy in the closed Newton cotes formul. Dehghn et l. incresed the order of ccurcy in the closed Newton cotes formul[] y including the loction of oundries of the intervl s two dditionl prmeter nd rescling the originl integrl to fit the optiml oundry loctions. They hve pplied this technique to open semiopen Guss legendre nd Guss Cheyshev i n t e g r t i o n r ules[]. Burg h s p r o p o s e d d e r i v t i v e s e d c l o s e d o p e n n d M i d p o i n t q u d r t u r e r u l e s [ 9 0 ]. Weijing Zho nd Hongxing Li took different pproch y introducing Midpoint technique t the computtion of derivtive[]. Recently we proposed Midpoint derivtive sed open Newton cotes Qudrture rule[] nd Geometric men derivtive sed closed Newton cotes Qudrture rule[]. In this pper the use of Hrmonic men derivtive t the end points is investigted in closed Newton cotes qudrture formul. This new scheme gives single order of precision thn the existing formul. The error terms re otined nd compred with the existing method. The ccurcy of the proposed rule is oserved y illustrting the Numericl exmples. II. Hrmonic Men derivtive -sed Closed Newton Cotes Qudrture rule To evlute the definite integrl new method of evlution of Hrmonic men derivtive - sed Closed Newton Cotes formuls re derived. In this method the Hrmonic men derivtive is zero if either =0 or =0.Tht is this method is not pplicle if either =0 or =0. Theorem.. Closed T r p e z o i d l R u l e ( n = ) using H r mo n ic The precision of this method is. f + f() () f" + () M e n d erivtive is Since the rule () hs the degree of precision. Now u s e the rule () for f (x) = x. When f x = x x dx = ; n = + () =. It shows tht the solution is Exct. Therefore the precision of Closed T r p e z o i d l R u l e with H r mo n i c M e n d erivtive is. Theorem.. Closed S i m p s o n ' s / rd The precision of this method is. f + f + R u l e with H r mo n i c M e n d erivtive (n=) is () + f() 0 f + () Since the rule () hs the degree of precision.now u s e the rule () for f (x) = x. n = When f x = x x dx = ; + + + () 0 =. It shows tht the solution is Exct. Therefore the precision of Closed S i m p s o n ' s / rd R u l e with H r mo n i c M e n d erivtive is. DOI: 0.9/-00 www.iosrjournls.org Pge

Theorem.. Closed S i m p s o n ' s / rd f + f + The precision of this method is. Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture + f + R u l e with H r mo n i c M e n d erivtive ( n=) is () + f 0 f + (9) Since the rule () hs the degree of precision.now u s e the rule (9) for f (x) = x. When f x = x x dx = ; n = + + + + + () = 0. It shows tht the solution is Exct. Therefore the precision of Closed S i m p s o n ' s / rd R u l e with H r mo n i c M e n d erivtive is. Theorem.. Closed B o o l e ' s R u l e with H r mo n i c M e n d erivtive (n=) is The precision of this method is. f + f + + f + + f + + f () f + Since the rule () hs the degree of precision.now u s e the rule (0) for f (x) = x. n = + =. When f x = x x dx = ; + + + DOI: 0.9/-00 www.iosrjournls.org Pge + + 0( + ) It shows tht the solution is Exct. Therefore the precision of Closed B o o l e ' s R u l e with M e n d erivtive is III. The Error terms of Hrmonic Men derivtive -sed Closed Newton Cotes Qudrture rule The error of pproximtion for the method sed on the precision of qudrture formul is otined y using the difference etween the qudrture formul for the monomil xp + p+! nd the exct result p+! x p+ (0) dx where p is the precision of the qudrture formul. Theorem.. H r mo n i c M e n d erivtive-sed Closed T r p e z o i d l R u l e ( n=)with the error term is f + f f" + where ξ ( ).This is f i f t h order ccurte with the error term E f = ( + ) f (ξ). Let f x = x!! x ( + ) f (ξ) () dx = ; () f + f() f" + = ( + ) ( + )( + ) () () Therefore + ( + )( + ) () = +.

Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture Therefore the Error term is E f = ( + ) f (ξ). Theorem.. H r mo n i c M e n d erivtive-sed Closed S i m p s o n ' s / rd the error term is f + f + where ξ ( ).This is s e v e n t h order ccurte with the error term E f = 0( + ) f (ξ). f + f + Let f x = x! + f Therefore 0 Therefore the Error term is 0( + ) R u l e (n=)with () + f() 0 f + 0 + f ξ ()! x 0 f + = 0( + ) dx = 0 ; + + + ( + ) () + + + ( + ) () = 0( + ). E f = 0 + f ξ. Theorem.. H r mo n i c M e n d erivtive-sed Closed S i m p s o n ' s /th R u l e ( n=) with the error term is f + f + + f + + f Where ξ ( ).This is s e v e n t h order ccurte with the error term E f = ( + ) f (ξ) Therefore f + f + = Let f x = x! + f + + f! x dx = 0 f + 0 ; 0 f + ( + ) f (ξ) () + (( + + + + + )( + ) ) 0 + (( + + + + + )( + ) ) = ( + ). Therefore the Error term is E f = ( + ) f (ξ). DOI: 0.9/-00 www.iosrjournls.org 9 Pge

Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture Theorem.. H r mo n i c M e n d erivtive-sed Closed B o o l e ' s r u l e ( n=) with the error term is f + f + + f + + f + + f f + 9 00 + f ξ () Where ξ ( ).This is n i n t h order ccurte with the error term 9 E f = 00 + f ξ. + f + f = Therefore 00 Let f x = x! + f +! x + f + dx = 00 ; + f f +!. ((9 + 9 + 0 + 9 + 9 + 0 + 9 + 9 )( + )!. ((9 + 9 + 0 + 9 + 9 + 0 + 9 +9 )( + ) () ) () ) 9 = 00 +. Therefore the Error term is 9 E f = 00 + f ξ. The summry of Precision the orders nd the error terms for H r m o n i c men derivtive sed Closed Newton- Cotes Qudrture re shown in Tle. Tle : Comprison of Error terms Rules Precision Order Error terms Trpezoidl rule (n=) ( + ) f (ξ) Simpson's / rd rule (n=) 0 + f ξ Simpson's / th rule (n=) ( + ) f (ξ). Boole's rule (n=) 9 9 00 + f ξ IV. Numericl Results In this section The vlues of e x dx dx nd re estimted using the Hrmonic Men derivtive - sed closed Newton cotes formul nd the results re compred with the existing closed Newton-Cotes qudrture formul. The comprisons re shown in Tle nd. we know tht Error = Exct vlue - Approximte vlue +x Exmple.:solve e x dx nd compre the solutions with the CNC nd HMDCNC rules. DOI: 0.9/-00 www.iosrjournls.org 0 Pge

Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture Solution: Exct vlue of e x dx=.0 Tle : Comprison of CNC nd HMDCNC rules vlue CNC HMDCNC of n Approximte vlue Error Approximte vlue Error n=.09 0.99.99 0.00 n=. 0.00.0 0.0009 n=.0 0.000000.0 0.000 n=.00 0.00000.0 0.000000 Exmple : solve dx +x nd compre the solutions with the CNC nd HMDCNC rules. +x Solution: Exct vlue of dx =0.00 Tle : Comprison of CNC nd HMDCNC rules Vlue of CNC HMDCNC n Approximte vlue Error Approximte vlue Error n= 0. 0.009 0.0 0.0099 n= 0.0 0.0000 0.009 0.0000009 n= 0.009 0.00000 0.00 0.00000 n= 0.0 0.0000000 0.00 0.000000 V. Conclusion In this pper new scheme of Hrmonic men derivtive - sed Closed Newton- Cotes qudrture formuls were presented long with their error terms. This Hrmonic men derivtive vlue is included in the existing formul nd the error terms were derived y using the difference etween the qudrture formul for the monomils nd the exct results. This proposed scheme gives single order of precision thn the existing formul. Finlly the ccurcy of the proposed scheme is illustrted with Numericl exmples. References [] M.Evn nd T.Swrtz "Methods for pproximting integrls in sttistics with specil Emphsis on Byesin integrtion prolem"sttisticl Science vol.0 pp. - 99. [] M.K.JinS.R.K.Iyengr nd R.K.Jin Numericl methods for Scientific nd Computtion New Age Interntionl (P) limited Fifth Edition 00. [] M.Dehghn M.Msjed-Jmei nd M.R.Eslhchi On numericl improvement of closed Newton-Cotes qudrture rules Applied Mthemtics nd Computtions vol. pp. -0 00. [] M.Dehghn M.Msjed-Jmei nd M.R.Eslhchi On numericl improvement of open Newton-Cotes qudrture rules Applied Mthemtics nd Computtions vol. pp. - 00. [] M.Dehghn M.Msjed-Jmei nd M.R.Eslhchi The semi-open Newton-Cotes qudrture rule nd its numericl improvement Applied Mthemtics nd Computtions vol. pp. 9-0 00. [] E.Bolin M.Msjed-Jmei nd M.R.Eslhchi On numericl improvement of Guss - legendre Qudrture rules Applied Mthemtics nd Computtions vol.0 pp. 9-9 00. [] M.R.Eslhchi M.Dehghn nd M.Msjed-Jmei On numericl improvement of the first kind Guss-Cheyshev Qudrture rules Applied Mthemtics nd Computtionsvol. pp. - 00. [ ] Clrence O.E.Burg Derivtive-sed closed Newton-cotes numericl qudrture Applied Mthemtics nd Computtions vol. pp. 0-0 0. [9] Fiz ZfrSir Sleem nd Clrence O.E.Burg New Derivive sed open Newton- cotes qudrture rules Astrct nd Applied AnlysisVolume 0 Article ID 09 pges0. [ 0 ] Clrence O.E.Burg nd Ezechiel Degny Derivtive-sed midpoint qudrture rule Applied Mthemtics nd Computtions vol. pp. - 0. [] Weijing Zho nd Hongxing Midpoint Derivtive-Bsed Closed Newton-Cotes Qudrture Astrct nd Applied Anlysis vol.0 Article ID 0 pges0. [] T.Rmchndrn nd R.Priml"Open Newton cotes qudrture with midpoint derivtive for integrtion of Algeric furnctions"interntionl Journl of Reserch in Engineering nd TechnologyVol. pp. 0-0. [] T.RmchndrnD.Udykumr nd R.Priml"Geometric men derivtive -sed closed Newton cotes qudrture " Interntionl Journl of Pure & Engineering. MthemticsVol. pp. 0-0. DOI: 0.9/-00 www.iosrjournls.org Pge