06.09.2021 Views

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

SHOW MORE
SHOW LESS

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 )

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!