Save this PDF as:

Size: px
Start display at page:

## Transcription

1 Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl method is commonly referred to s qudrture In Clculus one discusses mny techniques for determining n ntiderivtive of fx) The integrl is then computed by evluting the ntiderivtive t the endpoints of the intervl However, mny integrnds of interest in science nd engineering do not hve known ntiderivtive Moreover, in some pplictions the integrnd is only known t few points in the intervl We would like to be ble to determine ccurte pproximtions of integrls lso in these situtions We will pproximte integrls 1) by sums fx j )w j ) These sums re referred to s qudrture rules The x j re the nodes nd the w j the weights of the qudrture rule We re interested in determining weights so tht the qudrture rules gives ccurte pproximtions of the integrl 1 for lrge clsses of integrnds fx) The nodes often cnnot be chosen freely, becuse the integrnd my only be known t certin points For instnce, the function fx) might not be explicitly known, only mesured vlues fx j ) t the nodes x j my be vilble We remrk tht you lredy encountered the pproximtion of integrls by sums in Clculus There, however, one ws less concerned with how smll the error E N f) = fx)dx fx j )w j 3) is for smll vlues of N; it ws sufficient tht the error E N f) converged to zero s N incresed to infinity In this lecture, we would like the error E N f) to be smll lredy for smll to modest number of terms N in the qudrture rule ) Method of undetermined coefficients Let the nodes x j, 1 j N, be given Throughout this lecture the nodes will be ordered so tht x 1 < x < < x N b 4) We would like to determine the weights w j so tht qudrture rule is exct for polynomils of s high degree s possible, ie, we would like the weights be such tht px)dx = for ll polynomils px) of s high degree s possible px j )w j 5) 1

2 The powers {1, x, x, x 3, } form bsis for polynomils Our objective therefore cn be expressed s follows: We would like the equlity 5) to hold for s mny powers px) = x j, j =, 1,,, s possible This requirement gives rise to the liner system of equtions for the weights, px) = 1 : w 1 + w + + w N = px) = x : x 1 w 1 + x w + + x N w N = px) = x : x 1w 1 + x w + + x N w N = px) = x N 1 : x N 1 1 w 1 + x N 1 w + + w N x N 1 N This system conveniently cn be expressed in the form x 1 x x N x N 1 1 x N 1 x N 1 N w 1 w w N = dx = b, xdx = 1 b ), x 1 dx = 3 b3 3 ), = xn 1 dx = 1 N bn N ) b 1 b ) 1 N bn N ) 6) 7) We recognize the mtrix s the trnspose of Vndermonde mtrix Vndermonde mtrices were encountered in the fll semester in connection with polynomil interpoltion, when we showed tht Vndermonde mtrices determined by distinct nodes re nonsingulr Determinnts re invrint under trnsposition This follows from the fct tht determinnts cn be expnded either by rows or by columns We conclude tht the mtrix in 7) is nonsingulr when the nodes x j re distinct Hence the bove liner system of equtions hs unique solution Thus, given n distinct nodes, we cn determine n unique weights w j, such tht the qudrture rule ) integrtes ll polynomils of degree strictly less thn n exctly Exmple 1 Consider the midpoint qudrture rule fx)dx fx 1 )w 1, x 1 = 1 + b) We only hve one weight to determine The liner system of equtions 6) reduces to px) = 1 : w 1 = dx = b This determines the weight w 1 Exmple Apply the midpoint rule to integrte the function fx) = exp 1 x/) ) 8) on the intervl [, 1] The midpoint rule pproximtes the integrnd by the constnt function f1/) = exp 7/8) on the intervl [, b] nd integrtes the ltter function exctly The blue grph of Figure 1 shows fx) nd the vlue of the integrl is the re below this grph The dshed red line displyes the constnt function tht is integrted exctly by the midpoint rule The vlue determined by the midpoint rule is the re below the dshed red line Sometimes N-node qudrture rules determined by solving the liner system of equtions 7) lso integrte powers x k for k N exctly The following exmple illustrtes this

3 Figure 1: Integrnd of Exmple nd pproximtion used by the midpoint rule Exmple 3 Consider the midpoint qudrture rule of Exmple 1 Appliction of this rule to fx) = x yields x 1 w 1 = 1 + b)b ) = 1 b ) The right-hnd side equls xdx, ie, the midpoint rule integrtes not only constnts but lso liner polynomils exctly Exercise 1 Determine the weights of the -node qudrture rule fx)dx fx 1 )w 1 + fx )w, x 1 =, x = b This rule is known s the trpezoidl rule By construction, it integrtes liner polynomils exctly Does it integrte polynomils of higher degree exctly? Justify your nswer Exercise Determine the weights of the 3-node qudrture rule 1 fx)dx fx 1 )w 1 + fx )w + fx 3 )w 3, x 1 =, x = 1/, x 3 = 1 This rule is known s Simpson s rule It integrtes qudrtic polynomils exctly by construction Does it integrte polynomils of higher degree exctly? Justify your nswer Wht is the nlogous qudrture rule for the intervl [1, ]? Hint: Do not recompute the qudrture rule, just mke chnge of vribles Integrting Lgrnge polynomils When discussing polynomil interpoltion lst fll, we considered different polynomil bses, such s monomils, Lgrnge polynomils, nd Newton polynomils We cn lso in the present context use bses different 3

4 from the monomil one For instnce, the nodes x 1, x,,x N determine the Lgrnge polynomils l k x) = N j k x x j x k x j, k = 1,,, N, 9) which form bsis for ll polynomils of degree t most N 1 We therefore my require tht n N-node qudrture rule integrtes the Lgrnge polynomils 9) exctly Substituting px) = l k x) into 5) yields l k x)dx = l k x j )w j = l k x k )w k = w k, k = 1,,,N, where the simplifictions of the right-hnd side follow from the fct tht { 1, k = j, l k x j ) =, k j Thus, the weights cn be obtined by integrting the Lgrnge polynomils However, the evlution of these integrls is tedious unless N is smll Composite qudrture rules nd Tylor expnsion Assume tht the midpoint rule of Exmple does not yield sufficiently ccurte pproximtion of the integrl We my then subdivide the intervl [, 1] into subintervls nd pply the midpoint rule on ech subintervl This defines the composite midpoint rule Figure : Integrnd of Exmple 4 nd pproximtion used by the composite midpoint rule obtined by dividing the intervl [, 1] into two subintervls of equl length 4

5 Exmple 4 Divide the intervl [, 1] into the subintervls [, 1/] nd [1/, 1] nd pply the midpoint rule on ech subintervl to the integrnd fx) defined by 8) This yields the composite midpoint rule 1 fx)dx fx 1 )w 1 + fx )w with x 1 = 1/4, x = 3/4, nd w 1 = w = 1/ The integrl is pproximted by the re below the dshed curve of Figure Compring the grphs of Figures 1 nd suggests tht subdivision of the intervl should increse the ccurcy of the computed pproximtion While the computtion of the weights by solving the liner system of equtions 7) is esily done in MATLAB or Octve, this pproch does not shed ny light on how the error behves when we increse the number of nodes Further insight on the behvior of qudrture rules cn be gined by expnding the integrnd into Tylor series Consider the pproximtion of the integrl fx)dx by the midpoint rule nd use the Tylor expnsion of fx) t x = h/ Here x h/ should be thought of s firly smll, nd we ssume tht fx) is continuously differentible s mny times s required Then ) ) h h fx) = f + f x h ) + f ) h! x h ) + f ) h x h 3! Integrting the left-hnd nd right-hnd sides from to h yields ) h ) h fx)dx = f dx + f x h + f ) h x h ) 3 dx + f ) h 3! 4! ) dx + f h! ) 3 + f ) h x h 4 + 4! ) ) x h ) 4 dx + x h ) dx The integrl over odd powers of x h/ vnishes due to symmetry The right-hnd side therefore simplifies to ) h fx)dx = f h + f ) h x h ) dx + f ) h x h 4 dx +! 4! ) ) h = f h + f ) h h 3 + f ) h h The midpoint rule pplied to the integrl on the left-hnd side gives the first term in the right-hnd side The remining terms express the qudrture error Hence, this error is given by ) h fx)dx f h = f ) h h 3 + f ) h h 5 + 1) 1 8 Since ll derivtives of order nd higher of liner function vnish, the right-hnd side vnishes for such functions It follows tht the midpoint rule is exct for liner functions We know this lredy, nd here it is consequence of the expnsion of the integrl in powers of h Exercise 3 Use Tylor expnsions round x = nd x = h to determine how the pproximtion of the integrl fx)dx determined by the trpezoidl rule depends on h Thus, express the qudrture error fx)dx h f) + fh)) 5

6 s function of h in similrly mnner s 1) Consider the pproximtion of the integrl 1) by the N-point composite midpoint rule with the nodes nd weights x j = + j 1 ) h, h = b N 1, w j = h, j = 1,,,N 11) Thus, the qudrture rule is given by M h f) = h fx j ) Anlogously to 1), we obtin for ech subintervl [x j h/, x j + h/] of length h the expression xj+h/ x j h/ nd summing over j = 1,,,N yields fx) M h f) = fx)dx fx j )h = f x j ) h 3 +, 1) 1 f x j ) h 3 + = 1 1 N f x j ) Nh ) The verge of the second derivtive vlues 1 N N f x j ) is not smller thn the smllest of the vlues f x j ) Similrly, the verge is not lrger thn lrgest of the vlues f x j ) This is expressed by the inequlities min 1 j N f x j ) 1 N f x j ) mx 1 j N f x j ) We ssumed f x) to be continuous Therefore there is vlue ξ in [, b], such tht f ξ) = 1 N f x j ) Substituting this expression into the right-hnd side of 13) yields Multiplying the expression for h in 11) by the denomintor gives fx) M h f) = f ξ) Nh ) which when substituted into 14) yields Nh = b + h, fx) M h f) = f ξ) b 1 h + 15) The right-hnd side indictes tht we cn expect the error to decrese by fctor 4 when h is hlved nd the number of nodes x j is doubled) 6

7 Generlly, one does not know in dvnce how mny nodes to use in order to chieve desired ccurcy We therefore my be interested in hlving h until consecutive pproximtions M h f), M h/ f), M h/4 f),, of the integrl do not vry much when h is further reduced Exercise 4 In pplictions with complicted functions, the evlution of the function vlues my dominte the rithmetic work required to evlute qudrture rules The computtion of M h f) requires the evlution of N function vlues How mny of these function vlues cn be used gin when evluting M h/ f)? Exercise 5 Consider the N-point composite trpezoidl rule for the pproximtion of 1), 1 T h f) = h fx 1) + fx ) + fx 3 ) + + fx N 1 ) + 1 ) fx N) with ) Use the representtion x j = + j 1)h, T h f) = N 1 h = b, j = 1,,, N N 1 h fx j) + fx j+1 )) nd the result of Exercise 3 to determine n expnsion of the qudrture error similr to 13) b) Wht is the nlog of formul 15)? How much is the error reduced when h is hlved? c) Assume tht T h f) hs been evluted, nd we would like to compute T h/ f) How mny dditionl function evlutions re required? Singulr integrnds, infinite intervls, nd dptive qudrture rules The composite midpoint nd trpezoidl rules cn be used lso for functions tht re not twice continuously differentible However, the error then my converge to zero slower when h is reduced nd mny nodes might be required to give smll qudrture error We discuss n pproch to remedy this sitution Consider the evlution of the integrl 1 fx) lnx)dx 16) where fx) is ssumed to be smooth function We note tht the trpezoidl rule cnnot be used since the integrnd is infinite t x = Nevertheless, the integrl exists for mny smooth functions fx) We therefore my consider using composite midpoint rule; see Exercise 6 Alterntively, the method of undetermined coefficients cn be pplied to the function fx) only Thus, we would like to determine qudrture rule of the form fx j )w j with given nodes, nd we seek to determine weights w j so tht the reltion px) lnx)dx = px j )w j 17) 7

8 holds for ll polynomils of s high degree s possible This yields the liner system of equtions, which is nlogous to 6), px) = 1 : w 1 + w + + w N = px) = x : x 1 w 1 + x w + + x N w N = px) = x : x 1w 1 + x w + + x N w N = px) = x N 1 : x N 1 1 w 1 + x N 1 w + + w N x N 1 N lnx)dx, b xlnx)dx, b x lnx)dx, = xn 1 lnx)dx The right-hnd side expressions cn be evluted by integrtion by prts Note tht the lloction of nodes x j is quite rbitrry; in prticulr, we my put node t the origin Exercise 6 Compute n pproximtion of 16) with fx) = expx ) by the -node composite midpoint rule Use the MATLAB function qud to determine the exct vlue Note tht the MATLAB function for the integrnd hs to written to llow vector rguments, eg, f=expx^)*logx); How lrge is the error? The MATLAB function qud pplies composite Simpson rule The MATLAB function qud is n exmple of n dptive qudrture rule Adptive rules re composite rules tht estimte the qudrture error by compring the result obtined with different mesh sizes h New nodes re llocted in subintervls [x j, x j+1 ] in which the qudrture error is deemed to be lrger thn prescribed tolernce Exercise 7 Determine -point qudrture rule with the sme nodes s in Exercise 6 by solving liner system of equtions of the form 18) Apply it to the function fx) of Exercise 6 Is this rule more ccurte thn the one in Exmple 6? Exercise 8 Evlute the integrl exp x 3 )dx with t lest 5 correct deciml digits using the MATLAB function qud This function requires the intervl of integrtion to be finite Therefore the integrl hs to be split, exp x 3 )dx = c exp x 3 )dx + c exp x 3 )dx, where the constnt c > is chosen lrge enough so tht the second integrl cn be neglected nd the first integrl is evluted with the function qud For instnce, we my choose c lrge enough so tht c exp x 3 )dx Determine such vlue of c nd justify your choice The error bound for the integrl evluted with qud then should not be lrger thn In this lecture, we hve fixed the nodes nd then determined the weights so tht we integrte polynomils of s high degree s possible A clever choice of nodes s well s weights in n N-point qudrture rule mkes it possible to integrte polynomils of degree up to N 1 exctly The midpoint rule is n exmple These rules re known s Gussin qudrture rules The nodes nd weights generlly re not known in closed form, but they cn be computed quite efficiently by solving n eigenvlue problem symmetric N N tridigonl mtrix This is often done with the QR-lgorithm or modifiction thereof 18) 8

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

NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

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

III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the

Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the

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

3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

### Numerical Integration

Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl

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

THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS CARLOS SUERO, MAURICIO ALMANZAR CONTENTS 1 Introduction 1 2 Proof of Gussin Qudrture 6 3 Iterted 2-Dimensionl Gussin Qudrture 20 4

### Numerical integration

2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

### Section 6.1 Definite Integral

Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined

### Math 8 Winter 2015 Applications of Integration

Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

### Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions

### Discrete Least-squares Approximations

Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve

### Orthogonal Polynomials and Least-Squares Approximations to Functions

Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny

### We know that if f is a continuous nonnegative function on the interval [a, b], then b

1 Ares Between Curves c 22 Donld Kreider nd Dwight Lhr We know tht if f is continuous nonnegtive function on the intervl [, b], then f(x) dx is the re under the grph of f nd bove the intervl. We re going

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

COSC 336 Numericl Anlysis I Numericl Integrtion nd Dierentition III - Guss Qudrture nd Adptive Qudrture Edgr Griel Fll 5 COSC 336 Numericl Anlysis I Edgr Griel Summry o the lst lecture I For pproximting

### Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given

### 1 Error Analysis of Simple Rules for Numerical Integration

cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion

### a a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.

Section 9 The Lplce Expnsion In the lst section, we defined the determinnt of (3 3) mtrix A 12 to be 22 12 21 22 2231 22 12 21. In this section, we introduce generl formul for computing determinnts. Rewriting

### The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

Mth 104: finl informtion The finl exm will tke plce on Fridy My 11th from 8m 11m in Evns room 60. The exm will cover ll prts of the course with equl weighting. It will cover Chpters 1 5, 7 15, 17 21, 23

### Math 554 Integration

Mth 554 Integrtion Hndout #9 4/12/96 Defn. A collection of n + 1 distinct points of the intervl [, b] P := {x 0 = < x 1 < < x i 1 < x i < < b =: x n } is clled prtition of the intervl. In this cse, we

### 20 MATHEMATICS POLYNOMIALS

0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of

### Math 113 Exam 2 Practice

Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number

### Improper Integrals. Introduction. Type 1: Improper Integrals on Infinite Intervals. When we defined the definite integral.

Improper Integrls Introduction When we defined the definite integrl f d we ssumed tht f ws continuous on [, ] where [, ] ws finite, closed intervl There re t lest two wys this definition cn fil to e stisfied:

### Week 10: Riemann integral and its properties

Clculus nd Liner Algebr for Biomedicl Engineering Week 10: Riemnn integrl nd its properties H. Führ, Lehrstuhl A für Mthemtik, RWTH Achen, WS 07 Motivtion: Computing flow from flow rtes 1 We observe the

### If u = g(x) is a differentiable function whose range is an interval I and f is continuous on I, then f(g(x))g (x) dx = f(u) du

Integrtion by Substitution: The Fundmentl Theorem of Clculus demonstrted the importnce of being ble to find nti-derivtives. We now introduce some methods for finding ntiderivtives: If u = g(x) is differentible

### Lecture 1. Functional series. Pointwise and uniform convergence.

1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

### The Trapezoidal Rule

_.qd // : PM Pge 9 SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion

### Convex Sets and Functions

B Convex Sets nd Functions Definition B1 Let L, +, ) be rel liner spce nd let C be subset of L The set C is convex if, for ll x,y C nd ll [, 1], we hve 1 )x+y C In other words, every point on the line

### 7. Numerical evaluation of definite integrals

7. Numericl evlution of definite integrls Tento učení text yl podpořen z Operčního progrmu Prh - Adptilit Hn Hldíková Numericl pproximtion of definite integrl is clled numericl qudrture, the formuls re

### Numerical Methods I. Olof Widlund Transcribed by Ian Tobasco

Numericl Methods I Olof Widlund Trnscribed by In Tobsco Abstrct. This is prt one of two semester course on numericl methods. The course ws offered in Fll 011 t the Cournt Institute for Mthemticl Sciences,

### The Definite Integral

CHAPTER 3 The Definite Integrl Key Words nd Concepts: Definite Integrl Questions to Consider: How do we use slicing to turn problem sttement into definite integrl? How re definite nd indefinite integrls

### 2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

AP Clculus BC Review Chpter 8 Prt nd Chpter 9 Things to Know nd Be Ale to Do Know everything from the first prt of Chpter 8 Given n integrnd figure out how to ntidifferentite it using ny of the following

### Introduction to Numerical Analysis

Introduction to Numericl Anlysis Doron Levy Deprtment of Mthemtics nd Center for Scientific Computtion nd Mthemticl Modeling (CSCAMM) University of Mrylnd June 14, 2012 D. Levy CONTENTS Contents 1 Introduction

### Distance And Velocity

Unit #8 - The Integrl Some problems nd solutions selected or dpted from Hughes-Hllett Clculus. Distnce And Velocity. The grph below shows the velocity, v, of n object (in meters/sec). Estimte the totl

### 31.2. Numerical Integration. Introduction. Prerequisites. Learning Outcomes

Numericl Integrtion 3. Introduction In this Section we will present some methods tht cn be used to pproximte integrls. Attention will be pid to how we ensure tht such pproximtions cn be gurnteed to be

### MATH 573 FINAL EXAM. May 30, 2007

MATH 573 FINAL EXAM My 30, 007 NAME: Solutions 1. This exm is due Wednesdy, June 6 efore the 1:30 pm. After 1:30 pm I will NOT ccept the exm.. This exm hs 1 pges including this cover. There re 10 prolems.

### Numerical Methods 2007

Physics Mster Course: Numericl Methods 2007 Hns Mssen Rdboud Universiteit Nijmegen Onderwijsinstituut Wiskunde, Ntuur- en Sterrenkunde Toernooiveld 1 6525 ED Nijmegen September 2007 1 Introduction In these

### MATH 174A: PROBLEM SET 5. Suggested Solution

MATH 174A: PROBLEM SET 5 Suggested Solution Problem 1. Suppose tht I [, b] is n intervl. Let f 1 b f() d for f C(I; R) (i.e. f is continuous rel-vlued function on I), nd let L 1 (I) denote the completion

### (0.0)(0.1)+(0.3)(0.1)+(0.6)(0.1)+ +(2.7)(0.1) = 1.35

7 Integrtion º½ ÌÛÓ Ü ÑÔÐ Up to now we hve been concerned with extrcting informtion bout how function chnges from the function itself. Given knowledge bout n object s position, for exmple, we wnt to know

### Anonymous Math 361: Homework 5. x i = 1 (1 u i )

Anonymous Mth 36: Homewor 5 Rudin. Let I be the set of ll u (u,..., u ) R with u i for ll i; let Q be the set of ll x (x,..., x ) R with x i, x i. (I is the unit cube; Q is the stndrd simplex in R ). Define

### The Wave Equation I. MA 436 Kurt Bryan

1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

### Math 360: A primitive integral and elementary functions

Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:

### 10. AREAS BETWEEN CURVES

. AREAS BETWEEN CURVES.. Ares etween curves So res ove the x-xis re positive nd res elow re negtive, right? Wrong! We lied! Well, when you first lern out integrtion it s convenient fiction tht s true in

### Lecture Notes: Orthogonal Polynomials, Gaussian Quadrature, and Integral Equations

18330 Lecture Notes: Orthogonl Polynomils, Gussin Qudrture, nd Integrl Equtions Homer Reid My 1, 2014 In the previous set of notes we rrived t the definition of Chebyshev polynomils T n (x) vi the following

### Notes on Calculus II Integral Calculus. Miguel A. Lerma

Notes on Clculus II Integrl Clculus Miguel A. Lerm November 22, 22 Contents Introduction 5 Chpter. Integrls 6.. Ares nd Distnces. The Definite Integrl 6.2. The Evlution Theorem.3. The Fundmentl Theorem

### CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES

Filomt 27:4 (2013) 649 658 DOI 10.2298/FIL1304649M Published by Fculty of Sciences nd Mthemtics University of Niš Serbi Avilble t: http://www.pmf.ni.c.rs/filomt CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED

### Math 4200: Homework Problems

Mth 4200: Homework Problems Gregor Kovčič 1. Prove the following properties of the binomil coefficients ( n ) ( n ) (i) 1 + + + + 1 2 ( n ) (ii) 1 ( n ) ( n ) + 2 + 3 + + n 2 3 ( ) n ( n + = 2 n 1 n) n,

### Problem Set 3

14.102 Problem Set 3 Due Tuesdy, October 18, in clss 1. Lecture Notes Exercise 208: Find R b log(t)dt,where0

### MA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.

MA123, Chpter 1: Formuls for integrls: integrls, ntiderivtives, nd the Fundmentl Theorem of Clculus (pp. 27-233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.

### MATH 222 Second Semester Calculus. Fall 2015

MATH Second Semester Clculus Fll 5 Typeset:August, 5 Mth nd Semester Clculus Lecture notes version. (Fll 5) This is self contined set of lecture notes for Mth. The notes were written by Sigurd Angenent,

### STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

### 5: The Definite Integral

5: The Definite Integrl 5.: Estimting with Finite Sums Consider moving oject its velocity (meters per second) t ny time (seconds) is given y v t = t+. Cn we use this informtion to determine the distnce

### Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

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

### 2

1 Notes for Numericl Anlysis Mth 5466 by S. Adjerid Virgini Polytechnic Institute nd Stte University (A Rough Drft) 2 Contents 1 Differentition nd Integrtion 5 1.1 Introduction............................

### Test , 8.2, 8.4 (density only), 8.5 (work only), 9.1, 9.2 and 9.3 related test 1 material and material from prior classes

Test 2 8., 8.2, 8.4 (density only), 8.5 (work only), 9., 9.2 nd 9.3 relted test mteril nd mteril from prior clsses Locl to Globl Perspectives Anlyze smll pieces to understnd the big picture. Exmples: numericl

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

Generliztions of the Ostrowski s inequlity K. S. Anstsiou Aristides I. Kechriniotis B. A. Kotsos Technologicl Eductionl Institute T.E.I.) of Lmi 3rd Km. O.N.R. Lmi-Athens Lmi 3500 Greece Abstrct Using

### Summary of Elementary Calculus

Summry of Elementry Clculus Notes by Wlter Noll (1971) 1 The rel numbers The set of rel numbers is denoted by R. The set R is often visulized geometriclly s number-line nd its elements re often referred

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

IJRET: Interntionl Journl of Reserch in Engineering nd Technology eissn: 9-6 pissn: -78 OPEN NEWTON - COTES QUADRATURE WITH MIDPOINT DERIVATIVE FOR INTEGRATION OF ALGEBRAIC FUNCTIONS T. Rmchndrn R.Priml

### MAC-solutions of the nonexistent solutions of mathematical physics

Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

### 7 - Continuous random variables

7-1 Continuous rndom vribles S. Lll, Stnford 2011.01.25.01 7 - Continuous rndom vribles Continuous rndom vribles The cumultive distribution function The uniform rndom vrible Gussin rndom vribles The Gussin

### 1 i n x i x i 1. Note that kqk kp k. In addition, if P and Q are partition of [a, b], P Q is finer than both P and Q.

Chpter 6 Integrtion In this chpter we define the integrl. Intuitively, it should be the re under curve. Not surprisingly, fter mny exmples, counter exmples, exceptions, generliztions, the concept of the

### USA Mathematical Talent Search Round 1 Solutions Year 21 Academic Year

1/1/21. Fill in the circles in the picture t right with the digits 1-8, one digit in ech circle with no digit repeted, so tht no two circles tht re connected by line segment contin consecutive digits.

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

Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued

### Continuous Random Variables

STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

### Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

Improper Integrls The First Fundmentl Theorem of Clculus, s we ve discussed in clss, goes s follows: If f is continuous on the intervl [, ] nd F is function for which F t = ft, then ftdt = F F. An integrl

### A basic logarithmic inequality, and the logarithmic mean

Notes on Number Theory nd Discrete Mthemtics ISSN 30 532 Vol. 2, 205, No., 3 35 A bsic logrithmic inequlity, nd the logrithmic men József Sándor Deprtment of Mthemtics, Bbeş-Bolyi University Str. Koglnicenu

### MATH 409 Advanced Calculus I Lecture 19: Riemann sums. Properties of integrals.

MATH 409 Advnced Clculus I Lecture 19: Riemnn sums. Properties of integrls. Drboux sums Let P = {x 0,x 1,...,x n } be prtition of n intervl [,b], where x 0 = < x 1 < < x n = b. Let f : [,b] R be bounded

### S. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:

FACTA UNIVERSITATIS NIŠ) Ser Mth Inform 9 00) 6 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Dedicted to Prof G Mstroinni for his 65th birthdy

### Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Chpter 3 MTRIX In this chpter: Definition nd terms Specil Mtrices Mtrix Opertion: Trnspose, Equlity, Sum, Difference, Sclr Multipliction, Mtrix Multipliction, Determinnt, Inverse ppliction of Mtrix in

### Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero

### Continuous probability distributions

Chpter 1 Continuous probbility distributions 1.1 Introduction We cll x continuous rndom vrible in x b if x cn tke on ny vlue in this intervl. An exmple of rndom vrible is the height of dult humn mle, selected

### Practice final exam solutions

University of Pennsylvni Deprtment of Mthemtics Mth 26 Honors Clculus II Spring Semester 29 Prof. Grssi, T.A. Asher Auel Prctice finl exm solutions 1. Let F : 2 2 be defined by F (x, y (x + y, x y. If

### Part III. Calculus of variations

Prt III. Clculus of vritions Lecture notes for MA342H P. Krgeorgis pete@mths.tcd.ie 1/43 Introduction There re severl pplictions tht involve expressions of the form Jy) = Lx,yx),y x))dx. For instnce, Jy)

### 7.6 The Use of Definite Integrals in Physics and Engineering

Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 7.6 The Use of Definite Integrls in Physics nd Engineering It hs been shown how clculus cn be pplied to find solutions to geometric problems

### Fundamental Methods of Numerical Extrapolation With Applications

Fundmentl Methods of Extrpoltion 1 Fundmentl Methods of Numericl Extrpoltion With Applictions Eric Hung-Lin Liu ehliu@mit.edu Keywords: numericl nlysis, extrpoltion, richrdson, romberg, numericl differentition,

### SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics

SCHOOL OF ENGINEERING & BUIL ENVIRONMEN Mthemtics An Introduction to Mtrices Definition of Mtri Size of Mtri Rows nd Columns of Mtri Mtri Addition Sclr Multipliction of Mtri Mtri Multipliction 7 rnspose

### FINALTERM EXAMINATION 2011 Calculus &. Analytical Geometry-I

FINALTERM EXAMINATION 011 Clculus &. Anlyticl Geometry-I Question No: 1 { Mrks: 1 ) - Plese choose one If f is twice differentible function t sttionry point x 0 x 0 nd f ''(x 0 ) > 0 then f hs reltive...

### First variation. (one-variable problem) January 14, 2013

First vrition (one-vrible problem) Jnury 14, 2013 Contents 1 Sttionrity of n integrl functionl 2 1.1 Euler eqution (Optimlity conditions)............... 2 1.2 First integrls: Three specil cses.................

### Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions

### Week 7 Riemann Stieltjes Integration: Lectures 19-21

Week 7 Riemnn Stieltjes Integrtion: Lectures 19-21 Lecture 19 Throughout this section α will denote monotoniclly incresing function on n intervl [, b]. Let f be bounded function on [, b]. Let P = { = 0

### Interpolation. Gaussian Quadrature. September 25, 2011

Gussin Qudrture September 25, 2011 Approximtion of integrls Approximtion of integrls by qudrture Mny definite integrls cnnot be computed in closed form, nd must be pproximted numericlly. Bsic building

### ECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance

Generl structure ECO 37 Economics of Uncertinty Fll Term 007 Notes for lectures 4. Stochstic Dominnce Here we suppose tht the consequences re welth mounts denoted by W, which cn tke on ny vlue between

### Matrix Eigenvalues and Eigenvectors September 13, 2017

Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

### Problem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:

(x + y ) = y + (x + y ) = x + Problem Set 9 Discussion: Nov., Nov. 8, Nov. (on probbility nd binomil coefficients) The nme fter the problem is the designted writer of the solution of tht problem. (No one

### The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.

ADVANCED CALCULUS PRACTICE PROBLEMS JAMES KEESLING The problems tht follow illustrte the methods covered in clss. They re typicl of the types of problems tht will be on the tests. 1. Riemnn Integrtion

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

Applied Mthemticl Sciences, Vol. 8, 214, no. 31, 1535-1542 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/1.12988/ms.214.43166 Efficient Computtion of Clss of Singulr Oscilltory Integrls by Steepest Descent

### D i k B(n;D)= X LX ;d= l= A (il) d n d d+ (A ) (kl)( + )B(n;D); where c B(:::;n c;:::) = B(:::;n c ;:::), in prticulr c n = n c. Using the reltions [n

Explicit solutions of the multi{loop integrl recurrence reltions nd its ppliction? P. A. BAKOV ; nstitute of ucler Physics, Moscow Stte University, Moscow 9899, Russi The pproch to the constructing explicit

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS +. The solution of the eqution will e (), () 0,, 5, 5. The roots of the given eqution ( p q) ( q r) ( r p) 0 + + re p q r p (), r p p q, q r p q (), (d), q r p q.

### Calculus MATH 172-Fall 2017 Lecture Notes

Clculus MATH 172-Fll 2017 Lecture Notes These notes re concise summry of wht hs been covered so fr during the lectures. All the definitions must be memorized nd understood. Sttements of importnt theorems

### UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE

UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE 1. Pointwise Convergence of Sequence Let E be set nd Y be metric spce. Consider functions f n : E Y for n = 1, 2,.... We sy tht the sequence

### Matrix & Vector Basic Linear Algebra & Calculus

Mtrix & Vector Bsic Liner lgebr & lculus Wht is mtrix? rectngulr rry of numbers (we will concentrte on rel numbers). nxm mtrix hs n rows n m columns M x4 M M M M M M M M M M M M 4 4 4 First row Secon row

### Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?

Synoptic Meteorology I: Finite Differences 16-18 September 2014 Prtil Derivtives (or, Why Do We Cre About Finite Differences?) With the exception of the idel gs lw, the equtions tht govern the evolution

### Not for reproduction

AREA OF A SURFACE OF REVOLUTION cut h FIGURE FIGURE πr r r l h FIGURE A surfce of revolution is formed when curve is rotted bout line. Such surfce is the lterl boundry of solid of revolution of the type

### UniversitaireWiskundeCompetitie. Problem 2005/4-A We have k=1. Show that for every q Q satisfying 0 < q < 1, there exists a finite subset K N so that

Problemen/UWC NAW 5/7 nr juni 006 47 Problemen/UWC UniversitireWiskundeCompetitie Edition 005/4 For Session 005/4 we received submissions from Peter Vndendriessche, Vldislv Frnk, Arne Smeets, Jn vn de

### 1 Online Learning and Regret Minimization

2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

### Topic 1 Notes Jeremy Orloff

Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble

### A. Limits - L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.

A. Limits - L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c

### C1M14. Integrals as Area Accumulators

CM Integrls s Are Accumultors Most tetbooks do good job of developing the integrl nd this is not the plce to provide tht development. We will show how Mple presents Riemnn Sums nd the ccompnying digrms

### Topic 6b Finite Difference Approximations

/8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?