First Semester in Numerical Analysis with Julia, 2020a
First Semester in Numerical Analysis with Julia, 2020a
First Semester in Numerical Analysis with Julia, 2020a
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 3. INTERPOLATION 112<br />
there exists ξ ∈ (a, b) such that<br />
f[x 0 , ..., x n ]= f (n) (ξ)<br />
.<br />
n!<br />
To prove this theorem, we need the generalized Rolle’s theorem.<br />
Theorem 59 (Rolle’s theorem). Suppose f is a differentiable function on (a, b). If f(a) =<br />
f(b), then there exists c ∈ (a, b) such that f ′ (c) =0.<br />
Theorem 60 (Generalized Rolle’s theorem). Suppose f has n derivatives on (a, b). If f(x) =<br />
0 at (n +1) dist<strong>in</strong>ct numbers x 0 ,x 1 , ..., x n ∈ [a, b], then there exists c ∈ (a, b) such that<br />
f (n) (c) =0.<br />
Proof of Theorem 58 . Consider the function g(x) =p n (x)−f(x). Observe that g(x i )=0for<br />
i =0, 1, ..., n. From generalized Rolle’s theorem, there exists ξ ∈ (a, b) such that g (n) (ξ) =0,<br />
which implies<br />
p n (n) (ξ) − f (n) (ξ) =0.<br />
S<strong>in</strong>ce p n (x) =f[x 0 ]+f[x 0 ,x 1 ](x − x 0 )+... + f[x 0 , ..., x n ](x − x 0 ) ···(x − x n−1 ), p (n)<br />
n (x) equals<br />
n! times the lead<strong>in</strong>g coefficient f[x 0 , ..., x n ]. Therefore<br />
f (n) (ξ) =n!f[x 0 , ..., x n ].<br />
3.3 Hermite <strong>in</strong>terpolation<br />
In polynomial <strong>in</strong>terpolation, our start<strong>in</strong>g po<strong>in</strong>t has been the x and y-coord<strong>in</strong>ates of some<br />
data we want to <strong>in</strong>terpolate. Suppose, <strong>in</strong> addition, we know the derivative of the underly<strong>in</strong>g<br />
function at these x-coord<strong>in</strong>ates. Our new data set has the follow<strong>in</strong>g form.<br />
Data:<br />
x 0 ,x 1 , ..., x n<br />
y 0 ,y 1 , ..., y n ; y i = f(x i )<br />
y ′ 0,y ′ 1, ..., y ′ n; y ′ i = f ′ (x i )