Polynomial Interpolation with n + 1 nodes
|
|
- Jonathan Small
- 6 years ago
- Views:
Transcription
1 Polynomial Interpolation with n + 1 nodes Given n + 1 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )),, (x n, f (x n )), Interpolating polynomial of degree n P(x) = f (x 0 )L 0 (x) + f (x 1 )L 1 (x) + + f (x n )L n (x), with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ),, P(x n ) = f (x n ), and Lagrangian polynomials L i (x) = j i x x j x i x j. For all i, j, L i (x j ) = { 1, if i = j, 0, if i j,
2 Polynomial Interpolation with n + 1 nodes Analysis: Estimate errors with polynomial interpolation Computation: How to compute P(x) numerically.
3 Polynomial Interpolation: Analysis Theorem: Assume that the nodal points x 0,, x n are mutually distinct, then the interpolating polynomial P(x) of degree n exists and is unique.
4 Polynomial Interpolation: Analysis Lemma: Let f (x) C n+1 [a, b] be such that f (x 1 ) = 0, f (x 2 ) = 0,, f (x n ) = 0, where x 1, x 2,, x n [a, b] are mutually distinct. Then there exists a ξ [min(x 1, x 2,, x n ), max(x 1, x 2,, x n )] such that f (n 1) (ξ) = 0.
5 Polynomial Interpolation Error Theorem: Suppose x 0, x 1,, x n are distinct numbers in the interval [a, b] and f C n+1 [a, b]. Then, for each x [a, b], a number ξ(x) between, x 0, x 1,, x n (hence (a, b)) exists with f (x) = P(x) + f (n+1) (ξ(x)) (n + 1)! where P(x) is the interpolating polynomial. (x x 0 )(x x 1 ) (x x n ),
6 Polynomial Interpolation Error: Proof If x = x 0, x 1,, x n, then error = 0 and theorem is true. Now let x be not equal to any node. Define function g for t [a, b] g(t) def = (f (t) P(t)) (f (x) P(x)) (t x 0)(t x 1 ) (t x n ) (x x 0 )(x x 1 ) (x x n ) = (f (t) P(t)) (f (x) P(x)) n j=0 (t x j ) (x x j ) C n+1 [a, b]. Then g(t) vanishes at n + 2 distinct points: g(x) = 0, g(x k ) = 0,, for k = 0, 1,, n. There must be a ξ between x and nodal points such that g (n+1) (ξ) = 0.
7 Polynomial Interpolation Error: Proof Since g (n+1) (ξ) = (f (t) P(t)) (n+1) t=ξ (f (x) P(x))( = f (n+1) (ξ) (f (x) P(x)) = 0 Therefore f (x) = P(x) + f (n+1) (ξ(x)) (n + 1)! (n + 1)! n j=0 (x x j) n j=0 (t x j ) (x x j ) )(n+1) t=ξ (x x 0 )(x x 1 ) (x x n ),
8 Polynomial Interpolation: Corollary Corollary: Suppose x 0, x 1,, x n are distinct numbers in the interval [a, b] and f is polynomial of degree at most n, then P(x) = f (x).
9 Secant Method: Order of Convergence Assume that the secant method p k+1 = p k f (p k)(p k p k 1 ), k = 1, 2,, f (p k ) f (p k 1 ) converges to the root p (f (p) = 0.) It follows that p k+1 p = p k p (f (p k) f (p))(p k p k 1 ) f (p k ) f (p k 1 ) = (p k p)(f (p k ) f (p k 1 )) (f (p k ) f (p))(p k p k 1 ) f (p k ) f (p k 1 ) ( p k p k 1 = f (p) f (p k ) (p ) k p)(f (p k ) f (p k 1 )) f (p k ) f (p k 1 ) (p k p k 1 )
10 Secant Method: Order of Convergence But the last expression P(p) def = f (p k ) + (p k p)(f (p k ) f (p k 1 )) (p k p k 1 ) is linear interpolation on nodes p k and p k 1, so there exists ξ k between p, p k and p k 1 with and therefore p k+1 p = f (p) P(p) = f (2) (ξ k ) 2 (p p k )(p p k 1 ), p k p k 1 f (p k ) f (p k 1 ) f (2) (ξ k ) (p p k )(p p k 1 ). 2
11 Secant Method: Order of Convergence = α Let α = > 1 2 be the golden ratio. Then α(α 1) = 1, and p k+1 p p k p α ( ) pk p α 1 p k 1 p α = p k p k 1 f (p k ) f (p k 1 ) f (2) (ξ k ) 2 f (2) (p) 2f (p), provided f (p) 0, or p k+1 p lim k p k p α = f (2) (p) 2f (p) 1 α. (the limit does exist.)
12 Polynomial Interpolation, Error Bounds: Estimate maximum error in second order polynomial interpolation of function f (x) = 1 x over interval [2, 4] using nodes x 0 = 2, x 1 = 2.75, and x 2 = 4.
13 Polynomial Interpolation, Error Bounds: f (x) = P(x) + f (3) (ξ(x)) (x 2)(x 2.75)(x 4), 3! where P(x) is the interpolating polynomial on [2, 4]. f (3) (x) = 6x 4, f (3) (ξ(x)) 6 2 4, def g(x) = (x 2)(x 2.75)(x 4) = x 3 35x x 22. d g(x) d x = 3x x = 1 (3x 7)(2x 7). 2 g( 7 2 ) = 9 16 = max x [2,4] g(x). ( g( 7 25 ) = 3 f (x) P(x) 1 3! = < 9 16 )
14 Polynomial Interpolation: Experiments: Easy function: f (x) = e x on [ 1, 1]. f (n) (ξ) e for all ξ ( 1, 1). 7 nodal points (n = 6); random x points.
15 Polynomial Interpolation: Experiments: Easy function: f (x) = e x on [ 1, 1]. f (n) (ξ) e for all ξ ( 1, 1). 21 nodal points (n = 20); random x points.
16 Polynomial Interpolation: Experiments: Hard function: f (x) = x f (n) (ξ) = on [ 1, 1]. n! (1.1 + ξ) n+1 10n+1 n! for all ξ ( 1, 1). 7 nodal points (n = 6); random x points.
17 Polynomial Interpolation: Experiments: Hard function: f (x) = x f (n) (ξ) = on [ 1, 1]. n! (1.1 + ξ) n+1 10n+1 n! for all ξ ( 1, 1). 21 nodal points (n = 20); random x points.
18 Polynomial Interpolation: Experiments: Hardest function: f (x) = x 2 on [ 1, 1]. f (n) (ξ) can be very large for some ξ ( 1, 1). 7 nodal points (n = 6); random x points.
19 Polynomial Interpolation: Experiments: Hardest function: f (x) = x 2 on [ 1, 1]. f (n) (ξ) can be very large for some ξ ( 1, 1). 21 nodal points (n = 20); random x points.
20 Polynomial Interpolation: Experiments: Hardest function: f (x) = x 2 on [ 1, 1]. f (n) (ξ) can be very large for some ξ ( 1, 1). 81 nodal points (n = 80); random x points.
21 Polynomial Interpolation: Analysis and Computation Analysis: Estimate errors with polynomial interpolation Computation: How to compute P(x) numerically
22 Polynomial Interpolation: Neville s Method Let Q(x) interpolate f (x) at x 0, x 1,, x k, Let Q(x) interpolate f (x) at x 1,, x k, x k+1. Then P(x) def = (x x k+1)q(x) (x x 0 ) Q(x) x 0 x k+1 interpolates f (x) at x 0, x 1,, x k, x k+1 Build higher degree interpolating polynomial from lower degreed ones.
23 Neville s Method: Proof Let Q(x) interpolate f (x) at x 0, x 1,, x k, Let Q(x) interpolate f (x) at x 1,, x k, x k+1. for 1 j k, P(x) def = (x x k+1)q(x) (x x 0 ) Q(x) x 0 x k+1 P(x j ) = (x j x k+1 )Q(x j ) (x j x 0 ) Q(x j ) x 0 x k+1 = (x j x k+1 )f (x j ) (x j x 0 )f (x j ) x 0 x k+1 = f (x j ).
24 Neville s Method: Proof Let Q(x) interpolate f (x) at x 0, x 1,, x k, Let Q(x) interpolate f (x) at x 1,, x k, x k+1. for j = 0, k + 1, P(x) def = (x x k+1)q(x) (x x 0 ) Q(x) x 0 x k+1 P(x 0 ) = (x 0 x k+1 )Q(x 0 ) (x 0 x 0 ) Q(x 0 ) x 0 x k+1 = f (x 0 ), P(x k+1 ) = (x k+1 x k+1 )Q(x k+1 ) (x k+1 x 0 ) Q(x k+1 ) x 0 x k+1 = f (x k+1 ). Thus P(x) interpolates f (x) at x 0, x 1,, x k, x k+1.
25 Neville s Method, general case Let Q(x) interpolate f (x) at all of x 0, x 1,, x k, x k+1 but x j, Let Q(x) interpolate f (x) at all of x 0, x 1,, x k, x k+1 but x i, i j. Then P(x) def = (x x j)q(x) (x x i ) Q(x) x i x j interpolates f (x) at all of x 0, x 1,, x k, x k+1
26 Neville s Method, recursive Step 1: compute Q i,i+1 (x), linear polynomial interpolating f (x) at x i, x i+1, i = 0, 1,, n 1. Step 2: compute Q i,i+1,,j+1 (x), polynomial of degree j i + 1, interpolating f (x) at x i, i + 1,, x j+1, from Q i,i+1,,j (x) and Q i+1,i+2,,j+1 (x), for i = 0, 1,, n 1, and j = i + 1, n.
27 Neville s Method, pseudo-code
28 Divided Differences Given n + 1 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )),, (x n, f (x n )), Interpolating polynomial of degree n It follows P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ) + + a n (x x 0 )(x x 1 ) (x x n 1 ), with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ),, P(x n ) = f (x n ). a 0 = P(x 0 ) = f (x 0 ), P(x) f (x 0 ) x x 0 = a 1 + a 2 (x x 1 ) + + a n (x x 1 ) (x x n 1 ).
29 Divided Differences Define f [x i ] = f (x i ), f [x i, x i+1 ] = f (x i+1) f (x i ) x i+1 x i, for all i Let x = x 1 in P(x) f (x 0 ) x x 0 = a 1 + a 2 (x x 1 ) + + a n (x x 1 ) (x x n 1 ). It follows that a 1 = P(x 1) f (x 0 ) = f [x 0, x 1 ] x 1 x 0 P(x) P(x 0 ) x x 0 = f [x 0, x 1 ] + a 2 (x x 1 ) + +a n (x x 1 ) (x x n 1 ) P[x 0, x] P[x 0, x 1 ] x x 1 = a a n (x x 2 ) (x x n 1 ).
30 Divided Differences recursive, f [x i, x i+1,, x j+1 ] = f [x i+1,, x j+1 ] f [x i, x i+1,, x j ] x j+1 x i, then in we have P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ) + + a n (x x 0 )(x x 1 ) (x x n 1 ), a 0 = f [x 0 ] a 1 = f [x 0, x 1 ] a 2 = f [x 0, x 1, x 2 ].. a n 1 = f [x 0, x 1,, x n ]
31 Newton Divided Difference Example
32 Newton Divided Difference Example
33 Coefficients are numbers on top of all f [ ] columns
34 Newton Divided Difference Theorem: Suppose that f C n [a, b] and x 0, x 1,, x n are distinct nodes in [a, b]. Then a number ξ (a, b) exists such that f [x 0, x 1,, x n ] = f (n) (ξ). n! Proof: Let g(x) = f (x) P(x), where P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ) + + a n (x x 0 )(x x 1 ) (x x n 1 ). Since g(x j ) = f (x j ) P(x j ) = 0, j = 0, 1,, n, it follows that a number ξ (a, b) exists such that g (n) (ξ) = 0. But g (n) (ξ) = f (n) (ξ) P (n) (ξ) = f (n) (ξ) a n n! = f (n) (ξ) f [x 0, x 1,, x n ] n! = 0.
35 Polynomial Interpolation: Experiments: Easy function: f (x) = e x on [ 1, 1]. f (n) (ξ) e for all ξ ( 1, 1). 7 nodal points (n = 6); random x points.
36 Divided difference coefficients for e x on [ 1, 1] 10 0 Divided Difference Coefficients, semilog scale
37 Polynomial Interpolation: Experiments: Hardest function: f (x) = x 2 on [ 1, 1]. f (n) (ξ) can be very large for some ξ ( 1, 1). 21 nodal points (n = 20); random x points.
38 Divided difference coefficients for 1/(0.2 + x 2 ) on [ 1, 1] 10 4 Divided Difference Coefficients, semilog scale
39 Divided difference, book version
40 Divided difference, O(n) memory
41 Interpolation via NDD
42 Forward Differences Let x 1 = x 0 + h, x 2 = x 0 + 2h f [x 0, x 1 ] = f [x 0, x 0 + h] = f (x 0 + h) f (x 0 ) h f [x 0, x 1, x 2 ] = f [x 0, x 0 + h, x 0 + 2h] = f [x 1, x 2 ] f [x 0, x 1 ] 2h = f (x 2) + f (x 0 ) 2f (x 1 ) 2h 2 (first order FD), (second order FD)
43 Forward Differences Let x 1 = x 0 + h, x 2 = x 0 + 2h f [x 0, x 1 ] = f [x 0, x 0 + h] = f (x 0 + h) f (x 0 ) h f [x 0, x 1, x 2 ] = f [x 0, x 0 + h, x 0 + 2h] = f [x 1, x 2 ] f [x 0, x 1 ] 2h = f (x 2) + f (x 0 ) 2f (x 1 ) 2h 2 (first order FD), (second order FD)
44 Backward Differences Let x n 1 = x n h, x n 2 = x n 2h f [x n 1, x n ] = f [x n h, x n ] = f (x n) f (x n h) (1st order BD), h f [x n 2, x n 1, x n ] = f [x n 2h, x n h, x n ] = f [x n h, x n ] f [x n 2h, x n h] 2h = f (x n) + f (x n 2 ) 2f (x n 1 ) 2h 2 (2nd order BD)
45 Backward Differences Let x n 1 = x n h, x n 2 = x n 2h f [x n 1, x n ] = f [x n h, x n ] = f (x n) f (x n h) (1st order BD), h f [x n 2, x n 1, x n ] = f [x n 2h, x n h, x n ] = f [x n h, x n ] f [x n 2h, x n h] 2h = f (x n) + f (x n 2 ) 2f (x n 1 ) 2h 2 (2nd order BD)
46 Double nodes (I) Given 2 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )) Interpolating polynomial of degree 1 P(x) = a 0 + a 1 (x x 0 ) with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ), Where a 0 = f (x 0 ), a 1 = f (x 1) f (x 0 ) x 1 x 0.
47 Double nodes (I) Given 2 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )) Interpolating polynomial of degree 1 Where P(x) = a 0 + a 1 (x x 0 ) with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ), a 0 = f (x 0 ), a 1 = f (x 1) f (x 0 ) x 1 x 0. Now let x 1 x 0, we obtain (Taylor expansion) P(x) = a 0 + a 1 (x x 0 ), where a 0 = f (x 0 ), a 1 = f (x 0 ) def = f [x 0, x 0 ]. P(x 0 ) = f (x 0 ), P (x 0 ) = f (x 0 ).
48 Double nodes (II) Given 3 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )), (x 2, f (x 2 )), Interpolating polynomial of degree 2 P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ), with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ), P(x 2 ) = f (x 2 ), where we have a 0 = f [x 0 ], a 1 = f [x 0, x 1 ], a 2 = f [x 0, x 1, x 2 ] = f [x 1, x 2 ] f [x 0, x 1 ] x 2 x 0.
49 Double nodes (II) Given 3 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )), (x 2, f (x 2 )), Interpolating polynomial of degree 2 P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ), with P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ), P(x 2 ) = f (x 2 ), where we have a 0 = f [x 0 ], a 1 = f [x 0, x 1 ], a 2 = f [x 0, x 1, x 2 ] = f [x 1, x 2 ] f [x 0, x 1 ] x 2 x 0. Now let x 2 = x 1, we obtain (mixed interpolation) f [x 1, x 2 ] = f (x 1 ) = f [x 1, x 1 ], a 2 = f [x 1, x 1 ] f [x 0, x 1 ] x 1 x 0 def = f [x 0, x 1, x 1 ] P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 )(x x 1 ) P(x 0 ) = f (x 0 ), P(x 1 ) = f (x 1 ) P (x 1 ) = f (x 1 ).
50 Double nodes (III) Given 4 distinct points (x 0, f (x 0 )), (x 1, f (x 1 )), (x 2, f (x 2 )), (x 3, f (x 3 )), Interpolating polynomial of degree 3 P(x) = a 0 +a 1 (x x 0 )+a 2 (x x 0 )(x x 1 )+a 3 (x x 0 )(x x 1 )(x x 2 ). where a 0 = f [x 0 ], a 1 = f [x 0, x 1 ], It follows a 2 = f [x 0, x 1, x 2 ] = f [x 1, x 2 ] f [x 0, x 1 ] x 2 x 0 a 3 = f [x 0, x 1, x 2, x 3 ] = f [x 1, x 2, x 3 ] f [x 0, x 1, x 2 ] x 3 x 0.
51 Double nodes (IV) Let x 1 = x 0, and x 3 = x 2. It follows that a 1 = f [x 0, x 1 ] = f [x 0, x 0 ] a 2 = f [x 1, x 2 ] f [x 0, x 1 ] x 2 x 0 = f [x 0, x 2 ] f [x 0, x 0 ] x 2 x 0 a 3 = f [x 0, x 1, x 2, x 3 ] = f [x 0, x 2, x 2 ] f [x 0, x 0, x 2 ] x 2 x 0. Hermite Interpolation: P(x) = a 0 + a 1 (x x 0 ) + a 2 (x x 0 ) 2 + a 3 (x x 0 ) 2 (x x 2 ) P(x 0 ) = f (x 0 ), P (x 0 ) = f (x 0 ), P(x 2 ) = f (x 2 ), P (x 2 ) = f (x 2 ).
Order of convergence
Order of convergence Linear and Quadratic Order of convergence Computing square root with Newton s Method Given a > 0, p def = a is positive root of equation Newton s Method p k+1 = p k p2 k a 2p k = 1
More informationBisection Method. and compute f (p 1 ). repeat with p 2 = a 2+b 2
Bisection Method Given continuous function f (x) on the interval [a, b] with f (a) f (b) < 0, there must be a root in (a, b). To find a root: set [a 1, b 1 ] = [a, b]. set p 1 = a 1+b 1 2 and compute f
More informationJim 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 informationA first order divided difference
A first order divided difference For a given function f (x) and two distinct points x 0 and x 1, define f [x 0, x 1 ] = f (x 1) f (x 0 ) x 1 x 0 This is called the first order divided difference of f (x).
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
More informationChapter 3 Interpolation and Polynomial Approximation
Chapter 3 Interpolation and Polynomial Approximation Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Polynomial Interpolation
More informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationNeville s Method. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Neville s Method
Neville s Method MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have learned how to approximate a function using Lagrange polynomials and how to estimate
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
More informationHermite Interpolation
Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationSolving Differential Equations Using Power Series
LECTURE 8 Solving Differential Equations Using Power Series We are now going to employ power series to find solutions to differential equations of the form () y + p(x)y + q(x)y = 0 where the functions
More informationSolving Differential Equations Using Power Series
LECTURE 25 Solving Differential Equations Using Power Series We are now going to employ power series to find solutions to differential equations of the form (25.) y + p(x)y + q(x)y = 0 where the functions
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 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 informationUNIT-II INTERPOLATION & APPROXIMATION
UNIT-II INTERPOLATION & APPROXIMATION LAGRANGE POLYNAMIAL 1. Find the polynomial by using Lagrange s formula and hence find for : 0 1 2 5 : 2 3 12 147 : 0 1 2 5 : 0 3 12 147 Lagrange s interpolation formula,
More informationInterpolation & Polynomial Approximation. Hermite Interpolation I
Interpolation & Polynomial Approximation Hermite Interpolation I Numerical Analysis (th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationMA 3021: Numerical Analysis I Numerical Differentiation and Integration
MA 3021: Numerical Analysis I Numerical Differentiation and Integration Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationMath 471. Numerical methods Root-finding algorithms for nonlinear equations
Math 471. Numerical methods Root-finding algorithms for nonlinear equations overlap Section.1.5 of Bradie Our goal in this chapter is to find the root(s) for f(x) = 0..1 Bisection Method Intermediate value
More informationNotes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University. (A Rough Draft)
Notes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) 1 2 Contents 1 Error Analysis 5 2 Nonlinear Algebraic Equations 7 2.1 Convergence
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
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 informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
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 informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More information2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13
2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13 REVIEW Lecture 12: Classification of Partial Differential Equations (PDEs) and eamples with finite difference discretizations Parabolic PDEs Elliptic
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
More information3 Interpolation and Polynomial Approximation
CHAPTER 3 Interpolation and Polynomial Approximation Introduction A census of the population of the United States is taken every 10 years. The following table lists the population, in thousands of people,
More informationConvergence Rates on Root Finding
Convergence Rates on Root Finding Com S 477/577 Oct 5, 004 A sequence x i R converges to ξ if for each ǫ > 0, there exists an integer Nǫ) such that x l ξ > ǫ for all l Nǫ). The Cauchy convergence criterion
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
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 informationLagrange Interpolation and Neville s Algorithm. Ron Goldman Department of Computer Science Rice University
Lagrange Interpolation and Neville s Algorithm Ron Goldman Department of Computer Science Rice University Tension between Mathematics and Engineering 1. How do Mathematicians actually represent curves
More informationChapter 4. Unconstrained optimization
Chapter 4. Unconstrained optimization Version: 28-10-2012 Material: (for details see) Chapter 11 in [FKS] (pp.251-276) A reference e.g. L.11.2 refers to the corresponding Lemma in the book [FKS] PDF-file
More informationNumerical Methods I: Polynomial Interpolation
1/31 Numerical Methods I: Polynomial Interpolation Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 16, 2017 lassical polynomial interpolation Given f i := f(t i ), i =0,...,n, we would
More information, applyingl Hospital s Rule again x 0 2 cos(x) xsinx
Lecture 3 We give a couple examples of using L Hospital s Rule: Example 3.. [ (a) Compute x 0 sin(x) x. To put this into a form for L Hospital s Rule we first put it over a common denominator [ x 0 sin(x)
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationEquations with regular-singular points (Sect. 5.5).
Equations with regular-singular points (Sect. 5.5). Equations with regular-singular points. s: Equations with regular-singular points. Method to find solutions. : Method to find solutions. Recall: The
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
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 information1 Review of Interpolation using Cubic Splines
cs412: introduction to numerical analysis 10/10/06 Lecture 12: Instructor: Professor Amos Ron Cubic Hermite Spline Interpolation Scribes: Yunpeng Li, Mark Cowlishaw 1 Review of Interpolation using Cubic
More information8.3 Numerical Quadrature, Continued
8.3 Numerical Quadrature, Continued Ulrich Hoensch Friday, October 31, 008 Newton-Cotes Quadrature General Idea: to approximate the integral I (f ) of a function f : [a, b] R, use equally spaced nodes
More informationNovember 20, Interpolation, Extrapolation & Polynomial Approximation
Interpolation, Extrapolation & Polynomial Approximation November 20, 2016 Introduction In many cases we know the values of a function f (x) at a set of points x 1, x 2,..., x N, but we don t have the analytic
More informationMath 578: Assignment 2
Math 578: Assignment 2 13. Determine whether the natural cubic spline that interpolates the table is or is not the x 0 1 2 3 y 1 1 0 10 function 1 + x x 3 x [0, 1] f(x) = 1 2(x 1) 3(x 1) 2 + 4(x 1) 3 x
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 informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationEVALUATING A POLYNOMIAL
EVALUATING A POLYNOMIAL Consider having a polynomial p(x) = a + a 1 x + a 2 x 2 + + a n x n which you need to evaluate for many values of x. How do you evaluate it? This may seem a strange question, but
More informationHomework 8 Solutions to Selected Problems
Homework 8 Solutions to Selected Problems June 7, 01 1 Chapter 17, Problem Let f(x D[x] and suppose f(x is reducible in D[x]. That is, there exist polynomials g(x and h(x in D[x] such that g(x and h(x
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
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 informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationNumerical Methods I: Interpolation (cont ed)
1/20 Numerical Methods I: Interpolation (cont ed) Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 Interpolation Things you should know 2/20 I Lagrange vs. Hermite interpolation
More informationJuly 21 Math 2254 sec 001 Summer 2015
July 21 Math 2254 sec 001 Summer 2015 Section 8.8: Power Series Theorem: Let a n (x c) n have positive radius of convergence R, and let the function f be defined by this power series f (x) = a n (x c)
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationMA 8019: Numerical Analysis I Solution of Nonlinear Equations
MA 8019: Numerical Analysis I Solution of Nonlinear Equations Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationPower Series Solutions to the Legendre Equation
Power Series Solutions to the Legendre Equation Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationPart 2: Linesearch methods for unconstrained optimization. Nick Gould (RAL)
Part 2: Linesearch methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
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 informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationA rational interpolation scheme with super-polynomial rate of convergence
Center for Turbulence Research Annual Research Briefs 2008 31 A rational interpolation scheme with super-polynomial rate of convergence By Q. Wang, P. Moin AND G. Iaccarino 1. Motivation and objectives
More informationAPPLIED MATHEMATICS ADVANCED LEVEL
APPLIED MATHEMATICS ADVANCED LEVEL INTRODUCTION This syllabus serves to examine candidates knowledge and skills in introductory mathematical and statistical methods, and their applications. For applications
More informationA Note on Multiplicity Weight of Nodes of Two Point Taylor Expansion
Applied Mathematical Sciences, Vol, 207, no 6, 307-3032 HIKARI Ltd, wwwm-hikaricom https://doiorg/02988/ams2077302 A Note on Multiplicity Weight of Nodes of Two Point Taylor Expansion Koichiro Shimada
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 informationTHE SECANT METHOD. q(x) = a 0 + a 1 x. with
THE SECANT METHOD Newton s method was based on using the line tangent to the curve of y = f (x), with the point of tangency (x 0, f (x 0 )). When x 0 α, the graph of the tangent line is approximately the
More informationChapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn
Chapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide n c: Solve the problem
More informationInterpolating Accuracy without underlying f (x)
Example: Tabulated Data The following table x 1.0 1.3 1.6 1.9 2.2 f (x) 0.7651977 0.6200860 0.4554022 0.2818186 0.1103623 lists values of a function f at various points. The approximations to f (1.5) obtained
More informationCaculus 221. Possible questions for Exam II. March 19, 2002
Caculus 221 Possible questions for Exam II March 19, 2002 These notes cover the recent material in a style more like the lecture than the book. The proofs in the book are in section 1-11. At the end there
More informationMAT137 Calculus! Lecture 45
official website http://uoft.me/mat137 MAT137 Calculus! Lecture 45 Today: Taylor Polynomials Taylor Series Next: Taylor Series Power Series Definition (Power Series) A power series is a series of the form
More informationMathematical Olympiad Training Polynomials
Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationPower Series Solutions to the Legendre Equation
Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial
More informationMATH 409 Advanced Calculus I Lecture 16: Mean value theorem. Taylor s formula.
MATH 409 Advanced Calculus I Lecture 16: Mean value theorem. Taylor s formula. Points of local extremum Let f : E R be a function defined on a set E R. Definition. We say that f attains a local maximum
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationInfinite Continued Fractions
Infinite Continued Fractions 8-5-200 The value of an infinite continued fraction [a 0 ; a, a 2, ] is lim c k, where c k is the k-th convergent k If [a 0 ; a, a 2, ] is an infinite continued fraction with
More informationDownloaded from
Question 1: Exercise 2.1 The graphs of y = p(x) are given in following figure, for some polynomials p(x). Find the number of zeroes of p(x), in each case. (i) (ii) (iii) Page 1 of 24 (iv) (v) (v) Page
More informationFunctions and Equations
Canadian Mathematics Competition An activity of the Centre for Education in Mathematics and Computing, University of Waterloo, Waterloo, Ontario Euclid eworkshop # Functions and Equations c 006 CANADIAN
More informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
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 informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationPolynomial Interpolation Part II
Polynomial Interpolation Part II Prof. Dr. Florian Rupp German University of Technology in Oman (GUtech) Introduction to Numerical Methods for ENG & CS (Mathematics IV) Spring Term 2016 Exercise Session
More informationNewton's forward interpolation formula
Newton's Interpolation Formulae Interpolation is the process of approximating a given function, whose values are known at N + 1 tabular points, by a suitable polynomial, P N (x) of degree N which takes
More informationThe Mean Value Theorem and the Extended Mean Value Theorem
The Mean Value Theorem and the Extended Mean Value Theorem Willard Miller September 21, 2006 0.1 The MVT Recall the Extreme Value Theorem (EVT) from class: If the function f is defined and continuous on
More information