Module 11.4: nag quad util Numerical Integration Utilities. Contents
|
|
- Kevin Morrison
- 5 years ago
- Views:
Transcription
1 Qudrture Module Contents Module 11.4: ng qud util Numericl Integrtion Utilities ng qud util provides utility procedures for computtion involving integrtion of functions, e.g., the computtion of the weights nd bscisse for Gussin integrtion rules. Contents Introduction Procedures ng qud gs wt bsc Clcultion of weights nd bscisse for Gussin qudrture rules, generl choice of rule Exmples Exmple 1: Abscisse nd weights by Guss Legendre formul Additionl Exmples References [NP3245/3/pdf] Module 11.4: ng qud util
2 Module Contents Qudrture Module 11.4: ng qud util [NP3245/3/pdf]
3 Qudrture Module Introduction Introduction This module provides procedures which perform tsks in numericl integrtion. Currently, one procedure is vilble which computes the weights (norml or djusted) nd bscisse for six different types of Gussin integrtion rules. [NP3245/3/pdf] Module 11.4: ng qud util
4 Module Introduction Qudrture Module 11.4: ng qud util [NP3245/3/pdf]
5 Qudrture ng qud gs wt bsc Procedure: ng qud gs wt bsc 1 Description ng qud gs wt bsc computes the weights (norml or djusted) nd bscisse ssocited with six different Guss rules. These weights nd bscisse my be used for the pproximtion of definite integrls (see Stroud nd Secrest [2], or pges of Dvis nd Rbinowitz [1]). 2 Usge USE ng qud util CALL ng qud gs wt bsc(rule,, b, wt, x [, optionl rguments]) 3 A rguments Note. All rry rguments re ssumed-shpe rrys. The extent in ech dimension must be exctly tht required by the problem. Nottion such s x(n) is used in the rgument descriptions to specify tht the rry x must hve exctly n elements. This procedure derives the vlue of the following problem prmeter from the shpe of the supplied rrys. n 1 the number of Gussin weights nd bscisse 3.1 Mndtory Arguments rule integer, intent(in) Input: indictes the qudrture rule to be used s specified in Section 6.1: if rule = 0 for the Guss Legendre rule; if rule = 1 for the Guss Jcobi rule with norml weights; if rule = 1 for the Guss Jcobi rule with djusted weights; if rule = 2 for the Guss Exponentil rule with norml weights; if rule = 2 for the Guss Exponentil rule with djusted weights; if rule = 3 for the Guss Lguerre rule with norml weights; if rule = 3 for the Guss Lguerre rule with djusted weights; if rule = 4 for the Guss Hermite rule with norml weights; if rule = 4 for the Guss Hermite rule with djusted weights; if rule = 5 for the Guss Rtionl rule with norml weights; if rule = 5 for the Guss Rtionl rule with djusted weights. Constrints: 5 rule 5. rel(kind=wp), intent(in) b rel(kind=wp), intent(in) Input: the prmeters nd b which occur in the qudrture formule s specified in Section 6.1. Constrints: If rule =0,±1 or±2, then < b; if rule = ±3, then b 0; if rule = ±4, then b > 0; if rule = ±5, then + b 0. [NP3245/3/pdf] Module 11.4: ng qud util
6 ng qud gs wt bsc Qudrture wt(n) rel(kind=wp), intent(out) x(n) rel(kind=wp), intent(out) Output: the n weights nd bscisse of the Guss rule. Note: if rule = 2 or 4 nd the optionl rgument c 0.0, then n odd number of weights nd bscisse, n, my rise n error condition (see error%code = 204 in Section 4). 3.2 Optionl Arguments Note. Optionl rguments must be supplied by keyword, not by position. The order in which they re described below my differ from the order in which they occur in the rgument list. c rel(kind=wp), intent(in), optionl d rel(kind=wp), intent(in), optionl Input: the prmeters c nd d which occur in some of the qudrture formule (see Section 6.1). For some rules c nd d must not be too lrge (see Section 4). Note: If rule =0,thenc nd d should not be present; if rule = ±2, ±3 or±4, then d should not be present. Constrints: if rule = ±1, then c > 1 ndd > 1; if rule = ±2, ±3 or±4, then c > 1; if rule = ±5, then c > 1 ndd > c +1. error type(ngerror), intent(inout), optionl The NAG fl 90 error-hndling rgument. See the Essentil Introduction, or the module document ng error hndling (1.2). You re recommended to omit this rgument if you re unsure how to use it. If this rgument is supplied, it must be initilized by cll to ng set error before this procedure is clled. 4 Error Codes Ftl errors (error%level = 3): error%code Description 301 An input rgument hs n invlid vlue. 303 Arry rguments hve inconsistent shpes. 305 Invlid bsence of n optionl rgument. Filures (error%level = 2): error%code Description 201 Convergence not chieved. ng sym tridig eig ll hs filed to obtin convergence. 202 One or more of the weights re lrger thn HUGE(1.0 wp), the lrgest floting-point number on this mchine. The overflowing weights re returned s HUGE(1.0 wp). Possible solutions re to use smller vlue of n; or, if using djusted weights, to chnge to norml weights. To exmine the vlues of the weights fter this filure, you must supply the optionl rgument error nd initilise it using ng set error nd set hlt level to 3 (see the Essentil Introduction or the module ng error hndling for hndling error exits) Module 11.4: ng qud util [NP3245/3/pdf]
7 Qudrture ng qud gs wt bsc 203 One or more of the weights re too smll to be distinguished from zero on this mchine. The underflowing weights re returned s zero, which my be usble pproximtion. Possible solutions re to use smller vlue of n; or, if using norml weights, to chnge to djusted weights. To exmine the vlues of the weights fter this filure, you must supply the optionl rgument error nd initilise it using ng set error nd set hlt level to 3 (see the Essentil Introduction or the module ng error hndling for hndling error exits). 204 Guss Exponentil or Guss Hermite djusted weights with n odd nd c 0.0. Theoreticlly, in these cses: for c > 0.0, the centrl djusted weight is infinite, nd the exct function f(x) is zero t the centrl bsciss; for c < 0.0, the centrl djusted weight is zero, nd the exct function f(x) is infinite t the centrl bsciss. In either cse, indeterminte. the contribution of the centrl bsciss to the summtion is In prctice, the centrl weight my not hve overflowed or underflowed, if there is sufficient rounding error in the vlue of the centrl bsciss. If the soft fil option is used, the weights nd bsciss returned my be usble; the user must be prticulrly creful not to round the centrl bsciss to its true vlue without simultneously rounding the centrl weight to zero or s pproprite, or the summtion will suffer. It would be preferble to use norml weights, if possible. It is importnt to remember tht, when switching from norml weights to djusted weights or vice vers, redefinition of f(x) isinvolved. Wrnings (error%level = 1): error%code Description 101 Optionl rgument present but not used If rule =0,cnd d should not be present; if rule = ±2 or±3 or±4, d should not be present. 5 Exmples of Usge Acomplete exmple of the use of this procedure ppers in Exmple 1 of this module document. 6 Further Comments 6.1 Mthemticl Bckground This procedure computes the weights w i nd bscisse x i for use in the summtion S = n w i f(x i ) i=1 which pproximtes definite integrl (see Stroud nd Secrest [2], or pges of Dvis nd Rbinowitz [1]). The following types re provided: [NP3245/3/pdf] Module 11.4: ng qud util
8 ng qud gs wt bsc Qudrture 1. Guss Legendre: b Constrint: b>. 2. Guss Jcobi: norml weights: b djusted weights: b f(x) dx, exct for f(x) =P 2n 1 (x). (b x) c (x ) d f(x) dx, exct for f(x) =P 2n 1 (x), f(x) dx, exct for f(x) =(b x) c (x ) d P 2n 1 (x). Constrint: c> 1, d> 1, b>. 3. Guss Exponentil: norml weights: b x + b c 2 f(x) dx, exct for f(x) =P 2n 1 (x), djusted weights: b f(x) dx, exct for f(x) = x + b c 2 P 2n 1 (x). Constrint: c> 1, b>. 4. Guss Lguerre: norml weights: djusted weights: x c e bx f(x) dx (b >0), x c e bx f(x) dx (b <0), exct for f(x) =P 2n 1 (x), f(x) dx (b >0), Constrint: c> 1, b Guss Hermite: norml weights: + djusted weights: + Constrint: c> 1, b>0. f(x) dx (b <0), exct for f(x) = x c e bx P 2n 1 (x). x c e b(x )2 f(x) dx, exct for f(x) =P 2n 1 (x), f(x) dx, exct for f(x) = x c e b(x )2 P 2n 1 (x) Module 11.4: ng qud util [NP3245/3/pdf]
9 Qudrture ng qud gs wt bsc 6. Guss Rtionl: norml weights: S x c f(x) dx ( + b>0), x + b d x c ( ) 1 x + b d f(x) dx ( + b<0), exct for f(x) =P 2n 1, x + b djusted weights: S f(x) dx ( + b>0), f(x) dx ( + b<0), exct for f(x) = ( ) x c 1 x + b d P 2n 1. x + b Constrint: c> 1, d>c+1, + b 0. In the bove formule, P 2n 1 (x) stnds for ny polynomil of degree 2n 1orlessinx. 6.2 Algorithmic Detil The method used is to clculte the bscisse, using ng sym tridig eig ll to find the eigenvlues of the pproprite tridigonl mtrix (see Golub nd Welsch [3]). The weights re then determined by the formul: 1 n 1 w i = Pj (x i ) 2 j=0 where Pj (x) isthejth orthogonl polynomil with respect to the weight function over the pproprite intervl. 6.3 Accurcy The ccurcy depends minly on n, with incresing loss of ccurcy for lrger vlues of n. Typiclly, one or two deciml digits my be lost from mchine ccurcy with n 20, nd three or four deciml digits my be lost for n Timing The mjor portion of the time is tken up by ng sym tridig eig ll during the clcultion of the bscisse, where the time is roughly proportionl to n 3. [NP3245/3/pdf] Module 11.4: ng qud util
10 ng qud gs wt bsc Qudrture Module 11.4: ng qud util [NP3245/3/pdf]
11 Qudrture Exmple 1 Exmple 1: Abscisse nd weights by Guss Legendre formul Genertes bscisse nd weights using the Guss Legendre 7-point formul. 1 Progrm Text Note. The listing of the exmple progrm presented below is double precision. Single precision users re referred to Section 5.2 of the Essentil Introduction for further informtion. PROGRAM ng_qud_util_ex01! Exmple Progrm Text for ng_qud_util! NAG fl90, Relese 3. NAG Copyright 1997.!.. Use Sttements.. USE ng_exmples_io, ONLY : ng_std_out USE ng_qud_util, ONLY : ng_qud_gs_wt_bsc!.. Implicit None Sttement.. IMPLICIT NONE!.. Intrinsic Functions.. INTRINSIC KIND!.. Prmeters.. INTEGER, PARAMETER :: n = 7 INTEGER, PARAMETER :: wp= KIND(1.0D0)!.. Locl Sclrs.. INTEGER :: j, rule REAL (wp) ::, b!.. Locl Arrys.. REAL (wp) :: wt(n), x(n)!.. Executble Sttements.. WRITE (ng_std_out,*) Exmple Progrm Results for ng_qud_util_ex01 = -1.0_wp b = 1.0_wp rule = 0 CALL ng_qud_gs_wt_bsc(rule,,b,wt,x) WRITE (ng_std_out,*) WRITE (ng_std_out, (1X,A,I3,A) ) Legendre formul,, n, points WRITE (ng_std_out,*) WRITE (ng_std_out, (6X,A,8X,A) ) Abscisse, Weights WRITE (ng_std_out,*) WRITE (ng_std_out, (2F15.5) ) (x(j),wt(j),j=1,n) END PROGRAM ng_qud_util_ex01 2 Progrm Dt None. 3 Progrm Results Exmple Progrm Results for ng_qud_util_ex01 Legendre formul, Abscisse 7 points Weights [NP3245/3/pdf] Module 11.4: ng qud util
12 Exmple 1 Qudrture Module 11.4: ng qud util [NP3245/3/pdf]
13 Qudrture Additionl Exmples Additionl Exmples Not ll exmple progrms supplied with NAG fl 90 pper in full in this module document. The following dditionl exmples, ssocited with this module, re vilble. ng qud util ex02 Genertes bscisse with both weights nd djusted weights using the Guss Jcobi 7-point formul. ng qud util ex03 Genertes bscisse with both weights nd djusted weights using the Guss Exponentil 7-point formul. ng qud util ex04 Genertes bscisse with both weights nd djusted weights using the Guss Lguerre 7-point formul. ng qud util ex05 Genertes bscisse with both weights nd djusted weights using the Guss Hermite 7-point formul. ng qud util ex06 Genertes bscisse with both weights nd djusted weights using the Guss Rtionl 7-point formul. [NP3245/3/pdf] Module 11.4: ng qud util
14 References Qudrture References [1] Dvis P J nd Rbinowitz P (1975) Methods of Numericl Integrtion Acdemic Press [2] Stroud AH nd Secrest D (1966) Gussin Qudrture Formuls Prentice-Hll [3] Golub G H nd Welsch J H (1969) Clcultion of Guss qudrture rules Mth. Comput Module 11.4: ng qud util [NP3245/3/pdf]
D01BBF NAG Fortran Library Routine Document
D01 Qudrture D01BBF NAG Fortrn Librry Routine Document Note. Before using this routine, plese red the Users Note for your implementtion to check the interprettion of bold itlicised terms nd other implementtion-dependent
More informationNAG Fortran Library Routine Document D01BBF.1
D01 Qudrture NAG Fortrn Librry Routine Document Note: before using this routine, plese red the Users Note for your implementtion to check the interprettion of bold itlicised terms nd other implementtion-dependent
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationNumerical 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
More informationConstruction of Gauss Quadrature Rules
Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationNUMERICAL 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
More informationModule 10.2: nag nlin eqn Roots of a Single Nonlinear Equation. Contents
Nonlinear Equations Module Contents Module 10.2: nag nlin eqn Roots of a Single Nonlinear Equation nag nlin eqn provides a procedure for computing a single root of a continuous function in a given interval.
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
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
More informationLecture 14: Quadrature
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
More informationInterpolation. 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
More informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More information1. 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 ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
More informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More information63. 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 =
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationSuppose we want to find the area under the parabola and above the x axis, between the lines x = 2 and x = -2.
Mth 43 Section 6. Section 6.: 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
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationSection 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
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More informationNumerical Integration
Numericl Integrtion Wouter J. Den Hn London School of Economics c 2011 by Wouter J. Den Hn June 3, 2011 Qudrture techniques I = f (x)dx n n w i f (x i ) = w i f i i=1 i=1 Nodes: x i Weights: w i Qudrture
More informationCOSC 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
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationB.Sc. in Mathematics (Ordinary)
R48/0 DUBLIN INSTITUTE OF TECHNOLOGY KEVIN STREET, DUBLIN 8 B.Sc. in Mthemtics (Ordinry) SUPPLEMENTAL EXAMINATIONS 01 Numericl Methods Dr. D. Mckey Dr. C. Hills Dr. E.A. Cox Full mrks for complete nswers
More informationJDEP 384H: Numerical Methods in Business
BT 3.4: Solving Nonliner Systems Chpter 4: Numericl Integrtion: Deterministic nd Monte Crlo Methods Instructor: Thoms Shores Deprtment of Mthemtics Lecture 20, Februry 29, 2007 110 Kufmnn Center Instructor:
More informationTrapezoidal Rule, n = 1, x 0 = a, x 1 = b, h = b a. f (x)dx = h 2 (f (x 0) + f (x 1 )) h3
Trpezoidl Rule, n = 1, x 0 =, x 1 = b, h = b f (x)dx = h 2 (f (x 0) + f (x 1 )) h3 12 f (ξ). Simpson s Rule: n = 3, x 0 =, x 1 = +b 2, x 2 = b, h = b 2. Qudrture Rule, double node t x 1 P 3 (x)dx = f (x
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationModule 3.9: nag kelvin fun Kelvin Functions. Contents
Special Functions Module Contents Module 3.9: nag kelvin fun Kelvin Functions nag kelvin fun contains procedures for approximating Kelvin functions. Contents Introduction... 3.9.3 Procedures nag kelvin
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More informationINTRODUCTION TO INTEGRATION
INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationNOTES ON HILBERT SPACE
NOTES ON HILBERT SPACE 1 DEFINITION: by Prof C-I Tn Deprtment of Physics Brown University A Hilbert spce is n inner product spce which, s metric spce, is complete We will not present n exhustive mthemticl
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationDiscrete 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
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
More informationx = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is
Jnury 28, 2002 13. The Integrl The concept of integrtion, nd the motivtion for developing this concept, were described in the previous chpter. Now we must define the integrl, crefully nd completely. According
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationIII. 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
More informationModule 10.1: nag polynom eqn Roots of Polynomials. Contents
Nonlinear Equations Module Contents Module 10.1: nag polynom eqn Roots of Polynomials nag polynom eqn provides a procedure for computing the roots of a polynomial with real or complex coefficients. Contents
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationContinuous 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
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationLecture 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
More informationMath 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
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationSection 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
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationNumerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b
Numericl Integrtion Newton Cotes Formuls Given function f : R R nd two rel numbers, b R, < b, we clculte (pproximtely) the integrl I(f,, b) = f (x) dx K. Frischmuth (IfM UR) Numerics for CSE 08/09 8 /
More informationRead section 3.3, 3.4 Announcements:
Dte: 3/1/13 Objective: SWBAT pply properties of exponentil functions nd will pply properties of rithms. Bell Ringer: 1. f x = 3x 6, find the inverse, f 1 x., Using your grphing clcultor, Grph 1. f x,f
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationMATH 144: Business Calculus Final Review
MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives
More information4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve
Dte: 3/14/13 Objective: SWBAT pply properties of exponentil functions nd will pply properties of rithms. Bell Ringer: Use your clcultor to solve 4 7x =250; 5 3x =500; HW Requests: Properties of Log Equtions
More informationThe Fundamental Theorem of Calculus
The Fundmentl Theorem of Clculus MATH 151 Clculus for Mngement J. Robert Buchnn Deprtment of Mthemtics Fll 2018 Objectives Define nd evlute definite integrls using the concept of re. Evlute definite integrls
More informationp(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) =
Norwegin University of Science nd Technology Deprtment of Mthemticl Sciences Pge 1 of 5 Contct during the exm: Elen Celledoni, tlf. 73593541, cell phone 48238584 PLESE NOTE: this solution is for the students
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationTests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More informationJACOBI MATRICES FOR MEASURES MODIFIED BY A RATIONAL FACTOR. This paper decribes how, given the Jacobi matrix J for the measure d(t), it is possible to
JACOBI MATRICES FOR MEASURES MODIFIED BY A RATIONAL FACTOR SYLVAN ELHAY AND JAROSLAV KAUTSKY y Abstrct. This pper decribes how, given the Jcobi mtrix J for the mesure d(t), it is possible to produce the
More informationEfficient 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
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More information1 The Riemann Integral
The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re
More informationUniversity of Texas MD Anderson Cancer Center Department of Biostatistics. Inequality Calculator, Version 3.0 November 25, 2013 User s Guide
University of Texs MD Anderson Cncer Center Deprtment of Biosttistics Inequlity Clcultor, Version 3.0 November 5, 013 User s Guide 0. Overview The purpose of the softwre is to clculte the probbility tht
More informationWe divide the interval [a, b] into subintervals of equal length x = b a n
Arc Length Given curve C defined by function f(x), we wnt to find the length of this curve between nd b. We do this by using process similr to wht we did in defining the Riemnn Sum of definite integrl:
More informationMA 124 January 18, Derivatives are. Integrals are.
MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,
More information37 Kragujevac J. Math. 23 (2001) A NOTE ON DENSITY OF THE ZEROS OF ff-orthogonal POLYNOMIALS Gradimir V. Milovanović a and Miodrag M. Spalević
37 Krgujevc J. Mth. 23 (2001) 37 43. A NOTE ON DENSITY OF THE ZEROS OF ff-orthogonal POLYNOMIALS Grdimir V. Milovnović nd Miodrg M. Splević b Fculty of Electronic Engineering, Deprtment of Mthemtics, University
More informationThe 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
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationCAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.
Exmintion 1 Posted 23 October 2002. Due no lter thn 5pm on Mondy, 28 October 2002. Instructions: 1. Time limit: 3 uninterrupted hours. 2. There re four questions, plus bonus. Do not look t them until you
More informationOrthogonal 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
More information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
More informationUNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction
Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More information38 Riemann sums and existence of the definite integral.
38 Riemnn sums nd existence of the definite integrl. In the clcultion of the re of the region X bounded by the grph of g(x) = x 2, the x-xis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationModule 5.2: nag sym lin sys Symmetric Systems of Linear Equations. Contents
Linear Equations Module Contents Module 5.2: nag sym lin sys Symmetric Systems of Linear Equations nag sym lin sys provides a procedure for solving real or complex, symmetric or Hermitian systems of linear
More informationSection 7.1 Integration by Substitution
Section 7. Integrtion by Substitution Evlute ech of the following integrls. Keep in mind tht using substitution my not work on some problems. For one of the definite integrls, it is not possible to find
More informationDisclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.
Disclimer: This is ment to help you strt studying. It is not necessrily complete list of everything you need to know. The MTH 33 finl exm minly consists of stndrd response questions where students must
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationChapter 1. Basic Concepts
Socrtes Dilecticl Process: The Þrst step is the seprtion of subject into its elements. After this, by deþning nd discovering more bout its prts, one better comprehends the entire subject Socrtes (469-399)
More information