Lecture 19: Continuous Least Squares Approximation


 Sherilyn Nichols
 8 months ago
 Views:
Transcription
1 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 some m n This exmple motivted our study of discrete lest squres problems ( subject with mny other diverse pplictions), but the choice of the m points is somewht rbitrry Suppose we simply wish for the pproximting polynomil to represent f throughout ll of [, b] Wht vlue should m tke? How should one pick the points {x k }? Suppose we uniformly distribute these pproximtion points over [, b]: set h m := (b )/m nd let x k = + kh m The lest squres error formul, when scled by h m, tkes the form of Riemnn sum tht, in the m limit, pproximtes n integrl: lim m h m m (f(x k ) p(x k )) 2 = (f(x) p(x)) 2 dx Tht is, in the limit of infinitely mny uniformly spced pproximtion points, we re ctully minimizing n integrl, rther thn sum In this lecture, we will see how to pose such problems s mtrix problem of dimension (n + 1)by(n + 1), insted of discrete lest squres problem with mtrix of dimension by(n + 1) 331 Inner products for function spces To fcilitte the development of continuous lest squres pproximtion theory, we introduce forml structure for C[, b] First, recognize tht C[, b] is liner spce: ny liner combintion of continuous functions on [, b] must itself be continuous on [, b] Definition The inner product of the functions f, g C[, b] is given by f, g = f(x)g(x) dx This inner product stisfies the following bsic xioms: αf + g, h = α f, h + g, h for ll f, g, h C[, b] nd ll α R; f, g = g, f for ll f, g C[, b]; f, f for ll f C[, b] Just s the vector 2norm nturlly follows from the vector inner product ( x 2 = x x), so we hve ( f L 2 := f, f 1/2 b 1/2 = f(x) dx) 2 Here the superscript 2 in L 2 refers to the fct tht the integrnd involves the squre of the function f; the L stnds for Lebesgue, coming from the fct tht this inner product cn be generlized from R b If we wnted to consider complexvlued functions f nd g, the inner product would be generlized to f, g = f(x)g(x) dx, giving f, g = g, f 18 Jnury M Embree, Rice University
2 C[, b] to the set of ll functions tht re squreintegrble in the sense of Lebesgue integrtion By restricting our ttention to continuous functions, we dodge the mesuretheoretic complexities (The Lebesgue theory gives more robust definition of the integrl thn the conventionl Riemnn pproch; for detils, consult MATH 425) 332 Lest squres minimiztion vi clculus Given some f C[, b], the bsic L 2 pproximtion problem seeks the polynomil p P n tht minimizes the error f p in the L 2 norm In symbols: min p P n f p L 2 We shll denote the polynomil tht ttins this minimum by p We cn solve this minimiztion problem using bsic clculus Consider this exmple for n = 1, where we optimize the error over polynomils of the form p(x) = c + c 1 x Note tht f p L 2 will be minimized by the sme polynomil s f p 2 L 2 Thus for ny given p P 1, the error function is given by E(c, c 1 ) := f(x) (c + c 1 x) 2 L 2 = = = (f(x) c c 1 x) 2 dx ( ) f(x) 2 2f(x)(c + c 1 x) + (c 2 + 2c c 1 x + c 2 1x 2 ) dx f(x) 2 dx 2c f(x) dx 2c 1 xf(x) dx + c 2 (b ) + c c 1 (b 2 2 ) c2 1(b 3 3 ) To find the optiml polynomil, p, we need to optimize E over c nd c 1, ie, we must find the vlues of c nd c 1 for which = = c c 1 First, compute c c 1 = 2 = 2 Setting these prtil derivtives equl to zero yields f(x) dx + 2c (b ) + c 1 (b 2 2 ) xf(x) dx + c (b 2 2 ) c 1(b 3 3 ) 2c (b ) + c 1 (b 2 2 ) = 2 c (b 2 2 ) c 1(b 3 3 ) = 2 f(x) dx xf(x) dx These equtions, liner in the unknowns c nd c 1, cn be written in the mtrix form [ 2(b ) b 2 2 ] [ ] [ ] c 2 b (b3 3 = f(x) dx ) c 1 2 xf(x) dx 18 Jnury M Embree, Rice University
3 When b this system lwys hs unique solution The resulting c nd c 1 re the coefficients for the monomilbsis expnsion of the lest squres pproximtion p P 1 to f on [, b] Exmple: f(x) = e x We pply this result to the function f(x) = e x for x [, 1] Since 1 e x dx = e 1, 1 xe x dx = [e x (x 1)] 1 x= = 1, we must solve the system [ ] [ ] [ ] 2 1 c 2e 2 2 = 1 3 c 1 2 The desired solution is c = 4e 1, c 1 = 18 6e Below we show plot of this pproximtion (left), nd the error f(x) p (x) 3 25 f(x) p(x) f(x) p * (x) x x We cn see from these pictures tht the pproximtion looks decent to the eye, but the error is not terribly smll (In fct, f p L 2 = 6277 ) We cn decrese tht error by incresing the degree of the pproximting polynomil Just s we used 2by2 liner system to find the best liner pproximtion, generl (n + 1)by(n + 1) liner system cn be constructed to yield the L 2 optiml degreen pproximtion 333 Generl polynomil bses Note tht we performed the bove minimiztion in the monomil bsis: p(x) = c + c 1 x is liner combintion of 1 nd x Our experience with interpoltion suggests tht different choices for the bsis my yield pproximtion lgorithms with superior numericl properties Thus, we develop the form of the pproximting polynomil in n rbitrry bsis Suppose {φ k } n is bsis for P n Then ny p P n cn be written s p(x) = c k φ k (x) 18 Jnury M Embree, Rice University
4 The error expression tkes the form E(c,, c n ) := f(x) p(x) 2 L 2 = ( f(x) 2 c k φ k (x)) dx As before, compute / c j for j =,, n: = f, f 2 c j = 2 f, φ j + Setting / c j = gives the n + 1 equtions f, φ j = c k f, φ k + 2c k φ k, φ j c k φ k, φ j l= c k c l φ k, φ l This is simply system of liner lgebric equtions, which cn be written in the mtrix form φ, φ φ, φ 1 φ, φ n c f, φ φ 1, φ φ 1, φ 1 c 1 = f, φ 1, φ n, φ φ n, φ 1 φ n, φ n c n f, φ n which we shll denote s Hc = b Suppose we pply this method on the intervl [, b] = [, 1] with the monomil bsis, φ k (x) = x k In tht cse, 1 φ k, φ j = x k, x j = x j+k 1 dx = j + k + 1, nd the coefficient mtrix hs n elementry structure In fct, this is form of the notorious Hilbert mtrix It is exceptionlly difficult to obtin ccurte solutions with this mtrix in floting point rithmetic, reflecting the fct tht the monomils re poor bsis for P n on [, 1] Let H denote the n + 1dimensionl Hilbert mtrix, nd suppose b is constructed so tht the exct solution to the system Hc = b is c = (1, 1,, 1) T Let ĉ denote computed solution to the system in MATLAB Idelly the forwrd error c ĉ 2 will be nerly zero (if the rounding errors incurred while constructing b nd solving the system re smll) Unfortuntely, this is not the cse entirely consistent with our nlysis of the sensitivity of liner systems, studied in Section 142 n κ(h) c ĉ See MD Choi, Tricks or trets with the Hilbert mtrix, Americn Mth Monthly 9 (1983) Jnury M Embree, Rice University
5 Clerly these errors re not cceptble! The lst few 2norm condition numbers re in fct smller thn they ought to be, consequence of the fct tht MATLAB is not computing the singulr vlue decomposition of the Hilbert mtrix exctly (MATLAB computes the condition number s the rtio of the mximum nd minimum singulr vlues) The stndrd lgorithm for computing singulr vlues obtins nswers with smll bsolute ccurcy, but not smll reltive ccurcy Thus we expect tht singulr vlues smller thn bout 1 16 H 2 my not even be computed to the correct order of mgnitude In the next lecture, we will see how betterconditioned bses for P n yield mtrices H for which we cn solve Hx = b much more ccurtely 334 Connection to discrete lest squres Why did the continuous lest squres pproximtion problem studied bove directly led to squre (n + 1) (n + 1) liner system, while the discrete lest squres problem introduced in Lecture 16 led to n (m + 1) (n + 1) lest squres problem? In the discrete cse, we seek to minimize c Af 2, where (using the monomil bsis) 1 x x 2 x n c f(x ) 1 x 1 x 2 1 x n 1 c 1 f(x 1 ) A = 1 x 2 x 2 2 x n 2, c = c 2, f = f(x 2 ) 1 x m x 2 m x n m c n f(x m ) We hve seen tht this discrete problem cn be solved vi the norml equtions A Ac = A f Now compute n f(x k) n x kf(x k ) A f = n x2 k f(x k) C n+1 n xn k f(x k) Notice tht if m + 1 pproximtion points re uniformly spced over [, b], x k h m = (b )/m, we hve f(x) dx f, 1 xf(x) dx f, x lim h ma f = m x2 f(x) dx = f, x 2, xn f(x) dx f, x n = + kh m for 18 Jnury M Embree, Rice University
6 which is precisely the right hnd side vector b C n+1 obtined for the continuous lest squres problem Similrly, the (j + 1, k + 1) entry of the mtrix A A C (n+1) (n+1) for the discrete problem cn be formed s m m (A A) j+1,k+1 = x j l xk l = x j+k l, l= nd thus for uniform grids, we hve in the limit tht Thus in ggregte we hve lim h m(a A) j+1,k+1 = m lim h ma A = H, m l= x j+k dx = x j, x k where H is the mtrix tht rose in the continuous lest squres problem We rrive t the following beutiful conclusion: The norml equtions A Ac = A f formed for polynomil pproximtion by discrete lest squres converges to exctly the sme (n + 1) (n + 1) system Hc = b s we independently derived for polynomil pproximtion by continuous lest squres In the ltter cse, clculus led us directly to the norml eqution form of the solution 18 Jnury M Embree, Rice University
Lecture 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 relvlues nd smooth The pproximtion of n integrl by numericl
More informationDiscrete Leastsquares Approximations
Discrete Lestsqures 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 informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationOrthogonal Polynomials and LeastSquares Approximations to Functions
Chpter Orthogonl Polynomils nd LestSqures Approximtions to Functions **4/5/3 ET. Discrete LestSqures Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationLecture 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
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 informationSections 5.2: The Definite Integral
Sections 5.2: The Definite Integrl In this section we shll formlize the ides from the lst section to functions in generl. We strt with forml definition.. The Definite Integrl Definition.. Suppose f(x)
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 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 informationThe 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
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 xxis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationMath 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
More informationUNIFORM 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
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 relvlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
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 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 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 informationMapping the delta function and other Radon measures
Mpping the delt function nd other Rdon mesures Notes for Mth583A, Fll 2008 November 25, 2008 Rdon mesures Consider continuous function f on the rel line with sclr vlues. It is sid to hve bounded support
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 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at UrbanaChampaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t UrbnChmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
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 informationMath 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:
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 informationThinPlate Splines. Contents
ThinPlte Splines Dvid Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Cretive Commons Attribution 4.0 Interntionl License. To view copy of this
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecturenotes/numericalinter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecturenotes/numericalinter03.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 informationTHE 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 2Dimensionl Gussin Qudrture 20 4
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06073 HIKARI Ltd, www.mhikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
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 informationWe 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
More informationLecture 3. Limits of Functions and Continuity
Lecture 3 Limits of Functions nd Continuity Audrey Terrs April 26, 21 1 Limits of Functions Notes I m skipping the lst section of Chpter 6 of Lng; the section bout open nd closed sets We cn probbly live
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 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
More information31.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
More informationCzechoslovak 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 HENSTOCKKURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When relvlued
More information1 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
More informationWeek 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
More informationdifferent methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).
Mth 1A with Professor Stnkov Worksheet, Discussion #41; Wednesdy, 12/6/217 GSI nme: Roy Zho Problems 1. Write the integrl 3 dx s limit of Riemnn sums. Write it using 2 intervls using the 1 x different
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 informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationQuantum Physics II (8.05) Fall 2013 Assignment 2
Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.
More informationSummer MTH142 College Calculus 2. Section J. Lecture Notes. Yin Su University at Buffalo
Summer 6 MTH4 College Clculus Section J Lecture Notes Yin Su University t Bufflo yinsu@bufflo.edu Contents Bsic techniques of integrtion 3. Antiderivtive nd indefinite integrls..............................................
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the righthnd side limit equls to the lefthnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationEulerMaclaurin Summation Formula 1
Jnury 9, EulerMclurin Summtion Formul Suppose tht f nd its derivtive re continuous functions on the closed intervl [, b]. Let ψ(x) {x}, where {x} x [x] is the frctionl prt of x. Lemm : If < b nd, b Z,
More informationMath 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,
More informationIf 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 ntiderivtives. We now introduce some methods for finding ntiderivtives: If u = g(x) is differentible
More informationChapter 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
More informationCalculus and linear algebra for biomedical engineering Week 11: The Riemann integral and its properties
Clculus nd liner lgebr for biomedicl engineering Week 11: The Riemnn integrl nd its properties Hrtmut Führ fuehr@mth.rwthchen.de Lehrstuhl A für Mthemtik, RWTH Achen Jnury 9, 2009 Overview 1 Motivtion:
More informationSeries: Elementary, then Taylor & MacLaurin
Lecture 3 Series: Elementry, then Tylor & McLurin This lecture on series strts in strnge plce, by revising rithmetic nd geometric series, nd then proceeding to consider other elementry series. The relevnce
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
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 informationPolynomial 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
More informationConvex 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
More informationBig idea in Calculus: approximation
Big ide in Clculus: pproximtion Derivtive: f (x) = df dx f f(x +h) f(x) =, x h rte of chnge is pproximtely the rtio of chnges in the function vlue nd in the vrible in very short time Liner pproximtion:
More informationNumerical 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
More informationChapter 4. Lebesgue Integration
4.2. Lebesgue Integrtion 1 Chpter 4. Lebesgue Integrtion Section 4.2. Lebesgue Integrtion Note. Simple functions ply the sme role to Lebesgue integrls s step functions ply to Riemnn integrtion. Definition.
More informationThe 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
More informationMathematical Analysis: Supplementary notes I
Mthemticl Anlysis: Supplementry notes I 0 FIELDS The rel numbers, R, form field This mens tht we hve set, here R, nd two binry opertions ddition, + : R R R, nd multipliction, : R R R, for which the xioms
More informationChapter 6. Riemann Integral
Introduction to Riemnn integrl Chpter 6. Riemnn Integrl WonKwng Prk Deprtment of Mthemtics, The College of Nturl Sciences Kookmin University Second semester, 2015 1 / 41 Introduction to Riemnn integrl
More informationFINALTERM EXAMINATION 2011 Calculus &. Analytical GeometryI
FINALTERM EXAMINATION 011 Clculus &. Anlyticl GeometryI 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...
More informationPractice 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
More informationMATH 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 relvlued function on I), nd let L 1 (I) denote the completion
More informationMidpoint Approximation
Midpoint Approximtion Sometimes, we need to pproximte n integrl of the form R b f (x)dx nd we cnnot find n ntiderivtive in order to evlute the integrl. Also we my need to evlute R b f (x)dx where we do
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 information(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
More informationMA123, 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. 27233, Gootmn) Chpter Gols: Assignments: Understnd the sttement of the Fundmentl Theorem of Clculus.
More information2 Sturm Liouville Theory
2 Sturm Liouville Theory So fr, we ve exmined the Fourier decomposition of functions defined on some intervl (often scled to be from π to π). We viewed this expnsion s n infinite dimensionl nlogue of expnding
More informationNecessary and Sufficient Conditions for Differentiating Under the Integral Sign
Necessry nd Sufficient Conditions for Differentiting Under the Integrl Sign Erik Tlvil 1. INTRODUCTION. When we hve n integrl tht depends on prmeter, sy F(x f (x, y dy, it is often importnt to know when
More informationIntegrals along Curves.
Integrls long Curves. 1. Pth integrls. Let : [, b] R n be continuous function nd let be the imge ([, b]) of. We refer to both nd s curve. If we need to distinguish between the two we cll the function the
More informationCLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTONCOTES 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
More information20 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
More informationMATH 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.
More information13.3 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS
33 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS As simple ppliction of the results we hve obtined on lgebric extensions, nd in prticulr on the multiplictivity of extension degrees, we cn nswer (in
More informationMA Handout 2: Notation and Background Concepts from Analysis
MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,
More information5.2 Volumes: Disks and Washers
4 pplictions of definite integrls 5. Volumes: Disks nd Wshers In the previous section, we computed volumes of solids for which we could determine the re of crosssection or slice. In this section, we restrict
More informationLine and Surface Integrals: An Intuitive Understanding
Line nd Surfce Integrls: An Intuitive Understnding Joseph Breen Introduction Multivrible clculus is ll bout bstrcting the ides of differentition nd integrtion from the fmilir single vrible cse to tht of
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 informationMATH 409 Advanced Calculus I Lecture 22: Improper Riemann integrals.
MATH 409 Advned Clulus I Leture 22: Improper Riemnn integrls. Improper Riemnn integrl If funtion f : [,b] R is integrble on [,b], then the funtion F(x) = x f(t)dt is well defined nd ontinuous on [,b].
More informationNWI: Mathematics. Various books in the library with the title Linear Algebra I, or Analysis I. (And also Linear Algebra II, or Analysis II.
NWI: Mthemtics Literture These lecture notes! Vrious books in the librry with the title Liner Algebr I, or Anlysis I (And lso Liner Algebr II, or Anlysis II) The lecture notes of some of the people who
More informationa 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
More informationIntroduction 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
More informationSection 14.3 Arc Length and Curvature
Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in
More informationA Convergence Theorem for the Improper Riemann Integral of Banach Spacevalued Functions
Interntionl Journl of Mthemticl Anlysis Vol. 8, 2014, no. 50, 24512460 HIKARI Ltd, www.mhikri.com http://dx.doi.org/10.12988/ijm.2014.49294 A Convergence Theorem for the Improper Riemnn Integrl of Bnch
More information10. AREAS BETWEEN CURVES
. AREAS BETWEEN CURVES.. Ares etween curves So res ove the xxis 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
More informationRiemann Integrals and the Fundamental Theorem of Calculus
Riemnn Integrls nd the Fundmentl Theorem of Clculus Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University September 16, 2013 Outline Grphing Riemnn Sums
More informationTopic 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?
More informationVariational Techniques for SturmLiouville Eigenvalue Problems
Vritionl Techniques for SturmLiouville Eigenvlue Problems Vlerie Cormni Deprtment of Mthemtics nd Sttistics University of Nebrsk, Lincoln Lincoln, NE 68588 Emil: vcormni@mth.unl.edu Rolf Ryhm Deprtment
More informationu(t)dt + i a f(t)dt f(t) dt b f(t) dt (2) With this preliminary step in place, we are ready to define integration on a general curve in C.
Lecture 4 Complex Integrtion MATHGA 2451.001 Complex Vriles 1 Construction 1.1 Integrting complex function over curve in C A nturl wy to construct the integrl of complex function over curve in the complex
More informationFor a continuous function f : [a; b]! R we wish to define the Riemann integral
Supplementry Notes for MM509 Topology II 2. The Riemnn Integrl Andrew Swnn For continuous function f : [; b]! R we wish to define the Riemnn integrl R b f (x) dx nd estblish some of its properties. This
More informationNOTES AND PROBLEMS: INTEGRATION THEORY
NOTES AND PROBLEMS: INTEGRATION THEORY SAMEER CHAVAN Abstrct. These re the lecture notes prepred for prticipnts of AFSI to be conducted t Kumun University, Almor from 1st to 27th December, 2014. Contents
More informationAnonymous 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
More information4.1. Probability Density Functions
STT 1 4.14. 4.1. Proility Density Functions Ojectives. Continuous rndom vrile  vers  discrete rndom vrile. Proility density function. Uniform distriution nd its properties. Expected vlue nd vrince of
More information1 Ordinary Differential Equations
1 Ordinry Differentil Equtions 1. Mthemticl Bckground 1..1 Smoothness Definition 1.1 A function f defined on [, b] is continuous t x [, b] if lim x x f(x) = f(x ). Remrk Note tht this implies existence
More informationMath 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
More informationSection 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40
Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since
More informationThe 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
More informationSummary 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 numberline nd its elements re often referred
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 informationA. 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
More informationMATH 409 Advanced Calculus I Lecture 18: Darboux sums. The Riemann integral.
MATH 409 Advnced Clculus I Lecture 18: Drboux sums. The Riemnn integrl. Prtitions of n intervl Definition. A prtition of closed bounded intervl [, b] is finite subset P [,b] tht includes the endpoints
More information