formulas with complex poles
|
|
- Mavis Nelson
- 5 years ago
- Views:
Transcription
1 Computing rational Gauss-Chebyshev quadrature formulas with complex poles Karl Deckers Joris Van Deun Adhemar Bultheel Department of Computer Science K.U.Leuven Augustus 7, 2006
2 Rational Gauss-Chebyshev quadrature Algorithm to compute the nodes and weights for rational Gauss-Chebyshev quadrature formulas. Gauss quadrature: 1 1 f(x)w(x)dx n λ nk f(x nk ) k=1 Chebyshev weight functions: w(x) = (1 x) a (1 + x) b, a, b { ± 1 } 2
3 Notations X I C Riemann sphere C { } I interval [ 1, 1] complement of I with respect to a set X A n sequence of poles {α 1,..., α n } C I L n space of rational functions with poles in A n
4 Back to the quadrature formula Theorem There exist a set of nodes x nk and weights λ nk, k = 1,..., n so that the quadrature formula 1 1 f(x)w(x)dx n λ nk f(x nk ) k=1 is exact for f L n 1 L n 1. In the special case in which α n is real, this quadrature formula is exact for f L n L n 1.
5 Nodes and weights nodes x nk = cos θ nk I satisfy F n (θ nk ) = πk, k = 1,..., n F n (θ) is strictly increasing with increasing θ [0, π] the nodes have to be computed numerically, e.g. using Newton s method
6 Nodes and weights weights the weights are given by λ nk = G n (x nk ), k = 1,..., n the weights can be computed straightforwardly
7 Introduction Asymptotic inflection point distribution Two methods for determining a set of initial values for Newton s method: Asymptotic Zero Distribution (AZD) Asymptotic Inflection Point Distribution (AIPD)
8 Asymptotic inflection point distribution Theorem Assume the sequence of poles is bounded away from I and the asymptotic distribution of the poles is given by a measure ν on (a subset of) C I, then the asymptotic distribution of the nodes is given by an absolutely continuous measure µ and the density of the nodes on [ 1, x] is given by t(x) = x 1 µ (u)du.
9 Asymptotic inflection point distribution distribution of the poles is known lim n α n = α R I θ (0) n,k = f AZD(t n,k ) f AZD (t) is the inverse of t(θ) {t n,k } n k=1 is a set of n equally distributed points in [0, 1] distribution of the poles is unknown t(θ) can be approximated by a finite sum t n (θ) we can use the cubic interpolating spline s AZD (t) to approximate the inverse of t n (θ) θ (0) n,k = s AZD(t n,k )
10 Asymptotic inflection point distribution
11 Asymptotic inflection point distribution Problem : does not work well for poles close to the boundary introducing large local maxima of df n(θ) dθ Example: a = [-.5+i*1e-3*ones(1,2),.75+i*1e-2*ones(1,4), 1.01,-2] w = 2 %w(x) = \sqrt((1-x)/(1+x))
12 Asymptotic inflection point distribution Asymptotic inflection point distribution define θ (0) b j l j = F n = f AIP D (α j ) θ bj θ (0) b j π approximate from the left θ (0) n,k = θ(0) b j, with j = arg max j (l j k) approximate from the right θ (0) n,k = θ(0) b j, with j = arg min j (l j k)
13 Asymptotic inflection point distribution Asymptotic inflection point distribution
14 Syntax for gqcorf [x,l,err,fail] = gqcorf(a,w) where x vector with the resulting nodes L vector with the resulting weights err (optional) to check whether the computations succeeded fail (optional) vector with indices of nodes/weights for which the computations failed (if any) a vector with poles w (optional) choice of weight function
15 Introduction Example: a = [-.5+i*1e-3*ones(1,2),.75+i*1e-2*ones(1,4), 1.01,-2] w = 2 %w(x) = \sqrt((1-x)/(1+x)) method x λ π n k=1 λ k i total found bisection AZD AIPD 17 2
16 Introduction
17 Introduction The complexity of the algorithm is of order ϑ(m n). If m << n the complexity is of order ϑ(n). Example: m = 5 n t n t n t
18 Introduction Example: n = 8192 m t m t m t
Jim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationLearning Target: I can sketch the graphs of rational functions without a calculator. a. Determine the equation(s) of the asymptotes.
Learning Target: I can sketch the graphs of rational functions without a calculator Consider the graph of y= f(x), where f(x) = 3x 3 (x+2) 2 a. Determine the equation(s) of the asymptotes. b. Find the
More informationRational Gauss-Radau and rational Szegő-Lobatto quadrature on the interval and the unit circle respectively
Rational Gauss-Radau and rational Szegő-Lobatto quadrature on the interval and the unit circle respectively Francisco Perdomo-Pio Departament of Mathematical Analysis. La Laguna University. 38271 La Laguna.
More informationMath 112 (Calculus I) Midterm Exam 3 KEY
Math 11 (Calculus I) Midterm Exam KEY Multiple Choice. Fill in the answer to each problem on your computer scored answer sheet. Make sure your name, section and instructor are on that sheet. 1. Which of
More informationCalculus 221 worksheet
Calculus 221 worksheet Graphing A function has a global maximum at some a in its domain if f(x) f(a) for all other x in the domain of f. Global maxima are sometimes also called absolute maxima. A function
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationi x i y i
Department of Mathematics MTL107: Numerical Methods and Computations Exercise Set 8: Approximation-Linear Least Squares Polynomial approximation, Chebyshev Polynomial approximation. 1. Compute the linear
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationExam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20
Exam 2 Average: 85.6 Median: 87.0 Maximum: 100.0 Minimum: 55.0 Standard Deviation: 10.42 Fall 2011 1 Today s class Multiple Variable Linear Regression Polynomial Interpolation Lagrange Interpolation Newton
More information4.1 Analysis of functions I: Increase, decrease and concavity
4.1 Analysis of functions I: Increase, decrease and concavity Definition Let f be defined on an interval and let x 1 and x 2 denote points in that interval. a) f is said to be increasing on the interval
More informationx 2 x n r n J(x + t(x x ))(x x )dt. For warming-up we start with methods for solving a single equation of one variable.
Maria Cameron 1. Fixed point methods for solving nonlinear equations We address the problem of solving an equation of the form (1) r(x) = 0, where F (x) : R n R n is a vector-function. Eq. (1) can be written
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More information1. Introduction. 2. Outlines
1. Introduction Graphs are beneficial because they summarize and display information in a manner that is easy for most people to comprehend. Graphs are used in many academic disciplines, including math,
More informationFINAL EXAMINATION, MAT 2010 December 12, Cell phones are strictly prohibited!
FINAL EXAMINATION, MAT 00 December, 03 Write your solutions in a blue book. To receive full credit you must show all work. You are allowed to use an approved graphing calculator unless otherwise indicated.
More informationCalculation of Highly Oscillatory Integrals by Quadrature Method
Calculation of Highly Oscillatory Integrals by Quadrature Methods Krishna Thapa Texas A&M University, College Station, TX May 18, 211 Outline Purpose Why Filon Type Quadrature Methods Our Expectations
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationA numerical solution of the constrained weighted energy problem and its relation to rational Lanczos iterations
A numerical solution of the constrained weighted energy problem and its relation to rational Lanczos iterations Karl Deckers, Andrey Chesnokov, and Marc Van Barel Department of Computer Science, Katholieke
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationMATH section 4.4 Concavity and Curve Sketching Page 1. is increasing on I. is decreasing on I. = or. x c
MATH 0100 section 4.4 Concavity and Curve Sketching Page 1 Definition: The graph of a differentiable function y = (a) concave up on an open interval I if df f( x) (b) concave down on an open interval I
More informationMath 121 Test 3 - Review 1. Use differentials to approximate the following. Compare your answer to that of a calculator
Math Test - Review Use differentials to approximate the following. Compare your answer to that of a calculator.. 99.. 8. 6. Consider the graph of the equation f(x) = x x a. Find f (x) and f (x). b. Find
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
More informationDefinitions & Theorems
Definitions & Theorems Math 147, Fall 2009 December 19, 2010 Contents 1 Logic 2 1.1 Sets.................................................. 2 1.2 The Peano axioms..........................................
More informationGeometric Interpolation by Planar Cubic Polynomials
1 / 20 Geometric Interpolation by Planar Cubic Polynomials Jernej Kozak, Marjeta Krajnc Faculty of Mathematics and Physics University of Ljubljana Institute of Mathematics, Physics and Mechanics Avignon,
More informationFinal Examination 201-NYA-05 May 18, 2018
. ( points) Evaluate each of the following limits. 3x x + (a) lim x x 3 8 x + sin(5x) (b) lim x sin(x) (c) lim x π/3 + sec x ( (d) x x + 5x ) (e) lim x 5 x lim x 5 + x 6. (3 points) What value of c makes
More informationName: Instructor: 1. a b c d e. 15. a b c d e. 2. a b c d e a b c d e. 16. a b c d e a b c d e. 4. a b c d e... 5.
Name: Instructor: Math 155, Practice Final Exam, December The Honor Code is in effect for this examination. All work is to be your own. No calculators. The exam lasts for 2 hours. Be sure that your name
More informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
More informationEdexcel GCE Further Pure Mathematics (FP1) Required Knowledge Information Sheet. Daniel Hammocks
Edexcel GCE Further Pure Mathematics (FP1) Required Knowledge Information Sheet FP1 Formulae Given in Mathematical Formulae and Statistical Tables Booklet Summations o =1 2 = 1 + 12 + 1 6 o =1 3 = 1 64
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationSolutions 2. January 23, x x
Solutions 2 January 23, 2016 1 Exercise 3.1.1 (a), p. 149 Let us compute Lagrange polynomials first for our case x 1 = 0, x 2 = 2, x 3 = 3: In [14]: from sympy import * init_printing(use_latex=true) x=symbols(
More informationMIDTERM 2. Section: Signature:
MIDTERM 2 Math 3A 11/17/2010 Name: Section: Signature: Read all of the following information before starting the exam: Check your exam to make sure all pages are present. When you use a major theorem (like
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationMathematical Economics: Lecture 3
Mathematical Economics: Lecture 3 Yu Ren WISE, Xiamen University October 7, 2012 Outline 1 Example of Graphing Example 3.1 Consider the cubic function f (x) = x 3 3x. Try to graph this function Example
More informationf (x) = 2x x = 2x2 + 4x 6 x 0 = 2x 2 + 4x 6 = 2(x + 3)(x 1) x = 3 or x = 1.
F16 MATH 15 Test November, 016 NAME: SOLUTIONS CRN: Use only methods from class. You must show work to receive credit. When using a theorem given in class, cite the theorem. Reminder: Calculators are not
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationSections 4.1 & 4.2: Using the Derivative to Analyze Functions
Sections 4.1 & 4.2: Using the Derivative to Analyze Functions f (x) indicates if the function is: Increasing or Decreasing on certain intervals. Critical Point c is where f (c) = 0 (tangent line is horizontal),
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationMIDTERM 2 REVIEW: ADDITIONAL PROBLEMS. 1 2 x + 1. y = + 1 = x 1/ = 1. y = 1 2 x 3/2 = 1. into this equation would have then given. y 1.
MIDTERM 2 REVIEW: ADDITIONAL PROBLEMS ) If x + y =, find y. IMPLICIT DIFFERENTIATION Solution. Taking the derivative (with respect to x) of both sides of the given equation, we find that 2 x + 2 y y =
More informationAPPLICATIONS OF DIFFERENTIATION
4 APPLICATIONS OF DIFFERENTIATION APPLICATIONS OF DIFFERENTIATION Many applications of calculus depend on our ability to deduce facts about a function f from information concerning its derivatives. APPLICATIONS
More informationMath 104A - Homework 3
Math 104A - Homework 3 Due 7/7 2.4.6 Show that the following sequences converge linearly to p = 0. How large must n be before we have p n p 5 10 2? a p n = 1/n. Since p n+1 0 p n 0 = 1/(n + 1) 1/n = n
More informationVariable. Peter W. White Fall 2018 / Numerical Analysis. Department of Mathematics Tarleton State University
Newton s Iterative s Peter W. White white@tarleton.edu Department of Mathematics Tarleton State University Fall 2018 / Numerical Analysis Overview Newton s Iterative s Newton s Iterative s Newton s Iterative
More informationSpring 2015 Sample Final Exam
Math 1151 Spring 2015 Sample Final Exam Final Exam on 4/30/14 Name (Print): Time Limit on Final: 105 Minutes Go on carmen.osu.edu to see where your final exam will be. NOTE: This exam is much longer than
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationA.P. Calculus BC Test Three Section Two Free-Response No Calculators Time 45 minutes Number of Questions 3
A.P. Calculus BC Test Three Section Two Free-Response No Calculators Time 45 minutes Number of Questions 3 Each of the three questions is worth 9 points. The maximum possible points earned on this section
More informationIn 1854, Karl Weierstrauss gave an example of a continuous function which was nowhere di erentiable: cos(3 n x) 2 n. sin(3 n x), 2
Why non-pictured analysis? CHAPTER 1 Preliminaries f is continuous at x if lim f(x + h) = f(x) h!0 and f(x + h) f(x) f is di erentiable at x if lim h!0 h Then but Di erentiability =) continuity, continuity
More informationThe Condition of Linear Barycentric Rational Interpolation from Equispaced Samples
The Condition of Linear Barycentric Rational Interpolation from Equispaced Samples Georges Klein University of Fribourg (Switzerland) georges.klein@unifr.ch homeweb.unifr.ch/kleing/pub Vancouver, ICIAM,
More informationApplications of Derivatives
Applications of Derivatives Extrema on an Interval Objective: Understand the definition of extrema of a function on an interval. Understand the definition of relative extrema of a function on an open interval.
More informationConformal Mapping Lecture 20 Conformal Mapping
Let γ : [a, b] C be a smooth curve in a domain D. Let f (z) be a function defined at all points z on γ. Let C denotes the image of γ under the transformation w = f (z). The parametric equation of C is
More information= π + sin π = π + 0 = π, so the object is moving at a speed of π feet per second after π seconds. (c) How far does it go in π seconds?
Mathematics 115 Professor Alan H. Stein April 18, 005 SOLUTIONS 1. Define what is meant by an antiderivative or indefinite integral of a function f(x). Solution: An antiderivative or indefinite integral
More informationMath 180, Final Exam, Fall 2007 Problem 1 Solution
Problem Solution. Differentiate with respect to x. Write your answers showing the use of the appropriate techniques. Do not simplify. (a) x 27 x 2/3 (b) (x 2 2x + 2)e x (c) ln(x 2 + 4) (a) Use the Power
More informationc) xy 3 = cos(7x +5y), y 0 = y3 + 7 sin(7x +5y) 3xy sin(7x +5y) d) xe y = sin(xy), y 0 = ey + y cos(xy) x(e y cos(xy)) e) y = x ln(3x + 5), y 0
Some Math 35 review problems With answers 2/6/2005 The following problems are based heavily on problems written by Professor Stephen Greenfield for his Math 35 class in spring 2005. His willingness to
More informationOn Gauss-type quadrature formulas with prescribed nodes anywhere on the real line
On Gauss-type quadrature formulas with prescribed nodes anywhere on the real line Adhemar Bultheel, Ruymán Cruz-Barroso,, Marc Van Barel Department of Computer Science, K.U.Leuven, Celestijnenlaan 2 A,
More informationMATH 2053 Calculus I Review for the Final Exam
MATH 05 Calculus I Review for the Final Exam (x+ x) 9 x 9 1. Find the limit: lim x 0. x. Find the limit: lim x + x x (x ).. Find lim x (x 5) = L, find such that f(x) L < 0.01 whenever 0 < x
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationMATH 151 Engineering Mathematics I
MATH 151 Engineering Mathematics I Fall, 2016, WEEK 4 JoungDong Kim Week4 Section 2.6, 2.7, 3.1 Limits at infinity, Velocity, Differentiation Section 2.6 Limits at Infinity; Horizontal Asymptotes Definition.
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationSuppose that f is continuous on [a, b] and differentiable on (a, b). Then
Lectures 1/18 Derivatives and Graphs When we have a picture of the graph of a function f(x), we can make a picture of the derivative f (x) using the slopes of the tangents to the graph of f. In this section
More informationChapter 3 : Likelihood function and inference
Chapter 3 : Likelihood function and inference 4 Likelihood function and inference The likelihood Information and curvature Sufficiency and ancilarity Maximum likelihood estimation Non-regular models EM
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationInterpolation Theory
Numerical Analysis Massoud Malek Interpolation Theory The concept of interpolation is to select a function P (x) from a given class of functions in such a way that the graph of y P (x) passes through the
More informationKatholieke Universiteit Leuven Department of Computer Science
Separation of zeros of para-orthogonal rational functions A. Bultheel, P. González-Vera, E. Hendriksen, O. Njåstad Report TW 402, September 2004 Katholieke Universiteit Leuven Department of Computer Science
More informationMA 113 Calculus I Fall 2012 Exam 3 13 November Multiple Choice Answers. Question
MA 113 Calculus I Fall 2012 Exam 3 13 November 2012 Name: Section: Last 4 digits of student ID #: This exam has ten multiple choice questions (five points each) and five free response questions (ten points
More informationCS412: Introduction to Numerical Methods
CS412: Introduction to Numerical Methods MIDTERM #1 2:30PM - 3:45PM, Tuesday, 03/10/2015 Instructions: This exam is a closed book and closed notes exam, i.e., you are not allowed to consult any textbook,
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More information10550 PRACTICE FINAL EXAM SOLUTIONS. x 2 4. x 2 x 2 5x +6 = lim x +2. x 2 x 3 = 4 1 = 4.
55 PRACTICE FINAL EXAM SOLUTIONS. First notice that x 2 4 x 2x + 2 x 2 5x +6 x 2x. This function is undefined at x 2. Since, in the it as x 2, we only care about what happens near x 2 an for x less than
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationAP Calculus BC Chapter 4 AP Exam Problems. Answers
AP Calculus BC Chapter 4 AP Exam Problems Answers. A 988 AB # 48%. D 998 AB #4 5%. E 998 BC # % 5. C 99 AB # % 6. B 998 AB #80 48% 7. C 99 AB #7 65% 8. C 998 AB # 69% 9. B 99 BC # 75% 0. C 998 BC # 80%.
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationFundamentals of Numerical Mathematics. Charles University Prague Faculty of Mathematics and Physics Czech Republic
Fundamentals of Numerical Mathematics Vít Dolejší Charles University Prague Faculty of Mathematics and Physics Czech Republic dolejsi@karlin.mff.cuni.cz December 18, 2017 Preface These lecture notes more
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationEXAM 3 MAT 167 Calculus I Spring is a composite function of two functions y = e u and u = 4 x + x 2. By the. dy dx = dy du = e u x + 2x.
EXAM MAT 67 Calculus I Spring 20 Name: Section: I Each answer must include either supporting work or an explanation of your reasoning. These elements are considered to be the main part of each answer and
More informationQuick Review Sheet for A.P. Calculus Exam
Quick Review Sheet for A.P. Calculus Exam Name AP Calculus AB/BC Limits Date Period 1. Definition: 2. Steps in Evaluating Limits: - Substitute, Factor, and Simplify 3. Limits as x approaches infinity If
More informationNumerical Methods I: Numerical Integration/Quadrature
1/20 Numerical Methods I: Numerical Integration/Quadrature Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 umerical integration 2/20 We want to approximate the definite integral
More informationMath 241 Final Exam, Spring 2013
Math 241 Final Exam, Spring 2013 Name: Section number: Instructor: Read all of the following information before starting the exam. Question Points Score 1 5 2 5 3 12 4 10 5 17 6 15 7 6 8 12 9 12 10 14
More informationD(f/g)(P ) = D(f)(P )g(p ) f(p )D(g)(P ). g 2 (P )
We first record a very useful: 11. Higher derivatives Theorem 11.1. Let A R n be an open subset. Let f : A R m and g : A R m be two functions and suppose that P A. Let λ A be a scalar. If f and g are differentiable
More informationCHAPTER Order points and then develop a divided difference table:
1 CHAPTER 17 17.1 Order points and then develop a divided difference table: x y First Second Third Fourth Fifth Sixth 5 5.75 -.45 0.75 -.E-17 1.51788E-18 7.54895E-19 -.5719E-19 1.8 16.415 -.075 0.75-1.7E-17
More informationSOLUTIONS OF SELECTED PROBLEMS
SOLUTIONS OF SELECTED PROBLEMS Problem 36, p. 63 If µ(e n < and χ En f in L, then f is a.e. equal to a characteristic function of a measurable set. Solution: By Corollary.3, there esists a subsequence
More information1.18 Multiple Choice Questions on Limits
24 The AP CALCULUS PROBLEM BOOK 3x 4 2x + 33. lim x 7x 8x 5 =.8 Multiple Choice Questions on Limits A) B) C) 0 D) 3 7 E) 3 8 34. lim x 0 x = A) B) C) 0 D) E) does not exist 9x 2 35. lim x /3 3x = A) B)
More informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationPhysics 331 Introduction to Numerical Techniques in Physics
Physics 331 Introduction to Numerical Techniques in Physics Instructor: Joaquín Drut Lecture 12 Last time: Polynomial interpolation: basics; Lagrange interpolation. Today: Quick review. Formal properties.
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one
More informationComputational Methods. Solving Equations
Computational Methods Solving Equations Manfred Huber 2010 1 Solving Equations Solving scalar equations is an elemental task that arises in a wide range of applications Corresponds to finding parameters
More informationComputation of Rational Szegő-Lobatto Quadrature Formulas
Computation of Rational Szegő-Lobatto Quadrature Formulas Joint work with: A. Bultheel (Belgium) E. Hendriksen (The Netherlands) O. Njastad (Norway) P. GONZÁLEZ-VERA Departament of Mathematical Analysis.
More information18.01 EXERCISES. Unit 3. Integration. 3A. Differentials, indefinite integration. 3A-1 Compute the differentials df(x) of the following functions.
8. EXERCISES Unit 3. Integration 3A. Differentials, indefinite integration 3A- Compute the differentials df(x) of the following functions. a) d(x 7 + sin ) b) d x c) d(x 8x + 6) d) d(e 3x sin x) e) Express
More informationReview for the Final Exam
Math 171 Review for the Final Exam 1 Find the limits (4 points each) (a) lim 4x 2 3; x x (b) lim ( x 2 x x 1 )x ; (c) lim( 1 1 ); x 1 ln x x 1 sin (x 2) (d) lim x 2 x 2 4 Solutions (a) The limit lim 4x
More informationDIFFERENTIAL GEOMETRY 1 PROBLEM SET 1 SOLUTIONS
DIFFERENTIAL GEOMETRY PROBLEM SET SOLUTIONS Lee: -4,--5,-6,-7 Problem -4: If k is an integer between 0 and min m, n, show that the set of m n matrices whose rank is at least k is an open submanifold of
More informationSOLUTIONS TO EXAM 2, MATH 10550
SOLUTIONS TO EXAM 2, MATH 0550. Find the critical numbers of f(x) = 6 x2 x /3. We have f (x) = 3 x 3 x 2/3 = [ x 5/3 ] 3 x 2/3. So x = 0 is a critical point. For x 0, the equation f (x) = 0 can be written
More informationINTEGRATION WORKSHOP 2003 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 23 COMPLEX ANALYSIS EXERCISES DOUGLAS ULMER 1. Meromorphic functions on the Riemann sphere It s often useful to allow functions to take the value. This exercise outlines one way to
More informationConvergence rates of derivatives of a family of barycentric rational interpolants
Convergence rates of derivatives of a family of barycentric rational interpolants J.-P. Berrut, M. S. Floater and G. Klein University of Fribourg (Switzerland) CMA / IFI, University of Oslo jean-paul.berrut@unifr.ch
More informationSMOOTH APPROXIMATION OF DATA WITH APPLICATIONS TO INTERPOLATING AND SMOOTHING. Karel Segeth Institute of Mathematics, Academy of Sciences, Prague
SMOOTH APPROXIMATION OF DATA WITH APPLICATIONS TO INTERPOLATING AND SMOOTHING Karel Segeth Institute of Mathematics, Academy of Sciences, Prague CONTENTS The problem of interpolating and smoothing Smooth
More informationMath 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More informationChow groups of Zero-Cycles Relative to Hyperplane Arrangements. Karl-Dieter Crisman
Chow groups of Zero-Cycles Relative to Hyperplane Arrangements Karl-Dieter Crisman October 5, 2004 Notation for hyperplane arrangements: We will denote an arrangement by A. Hyperplanes {H i } m i=1 α i.
More informationMA102: Multivariable Calculus
MA102: Multivariable Calculus Rupam Barman and Shreemayee Bora Department of Mathematics IIT Guwahati Differentiability of f : U R n R m Definition: Let U R n be open. Then f : U R n R m is differentiable
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationDESIGN OF THE QUESTION PAPER
DESIGN OF THE QUESTION PAPER MATHEMATICS - CLASS XI Time : 3 Hours Max. Marks : 00 The weightage of marks over different dimensions of the question paper shall be as follows:. Weigtage of Type of Questions
More information