Métodos Numéricos modo virtual

(pandemia COVID-19)
Material Complementario

Interpolación - Error - versión 1.0

Este es material complementario de las diapos de la clase de interpolación usadas durante el dictado virtual
(pandemia COVID-19). En este documento deducimos la expresión del error del polinomio interpolante.

Sean f (x) una función definida en un intervalo [a, b] y pares ordenados (xi , f (xi )), xi ∈ [a, b] para i = 0, . . . , n.
Sabemos que existe un polinomio P (x) de grado ≤ n tal que P (xi ) = f (xi ) para todo i = 0, . . . , n.
Dado x̄ ∈ [a, b], x̄ 6= xi para todo i = 0, . . . , n estamos interesados en saber que error cometemos si aproxi-
mamos el valor de f (x̄) por P (x̄). En la próxima propiedad daremos respuesta a esto.

Proposición: Sea f (x) ∈ C n+1 [a, b], (xi , f (xi )), xi ∈ [a, b] para i = 0, . . . , n. Consideremos P (x) el
polinomio interpolante de grado ≤ n y x̄ ∈ [a, b]. Existe ξ(x̄) tal que

f n+1 (ξ(x̄))
f (x̄) = P (x̄) + (x̄ − x0 )(x̄ − x1 ) . . . (x̄ − xn )
(n + 1)!

Demostración:

ˆ Caso a: x̄ = xk para algún k ∈ {0, . . . , n}.

Sabemos que P (xk ) = f (xk ) porque P (x) es el polinomio interpolante en los puntos xi para i = 0, . . . , n.
Por otro lado (xk − x0 )(xk − x1 ) . . . (xk − xn ) se anula. Entonces ξ(x̄) puede elegirse en forma arbitraria
y la identidad es verdadera.
ˆ Caso b: x̄ 6= xk para todo k ∈ {0, . . . , n}.
n

Q (t−xi )
Definimos una función g(t) = f (t) − P (t) − f (x̄) − P (x̄) (x̄−xi ) para t ∈ [a, b].
i=0

Veamos que propiedades podemos deducir que cumple la función g(t). Sabemos que:
1. f (t) ∈ C n+1 [a, b] por hipótesis.
2. P (t) ∈ C n+1 [a, b] porque es un polinomio.
n
Q (t−xi ) n+1
3. (x̄−xi ) ∈ C [a, b] porque es un polinomio.
i=0

1
Métodos Numéricos modo virtual (pandemia COVID-19)
entonces podemos concluir que g(t) ∈ C n+1 [a, b].
¿Qué más podemos deducir? Evaluemos a g(t) en los puntos de interpolación:
n

Q (xk −xi )
g(xk ) = f (xk ) − P (xk ) − f (x̄) − P (x̄) (x̄−xi )
i=0
La última productoria se anula ya que k ∈ {0, . . . , n}. Además f (xk ) = P (xk ). Por lo tanto

g(xk ) = 0 para todo k ∈ {0, . . . , n}

Ahora evaluemos a g(t) en x̄:

n

Q (x̄−xi )
g(x̄) = f (x̄) − P (x̄) − f (x̄) − P (x̄) (x̄−xi ) .
i=0


La última productoria vale 1, entonces g(x̄) = f (x̄)−P (x̄)− f (x̄)−P (x̄) , lo que implica que g(x̄) = 0.
Resumiendo lo que sabemos de g(t) es que:
– g(t) ∈ C n+1 [a, b]
– g(t) se anula en x0 , . . . , xn y x̄.

Recordamos un resultado clásico del análisis (teorema de Rolle) que nos dice que si tenemos una
función h continua en [c, d] y diferenciable en (c, d) tal que h(c) = h(d), entonces existe ξ ∈ (a, b) tal
que h0 (ξ) = 0.
La función g(t) tiene al menos n + 2 puntos donde se anula. Si ordenamos x0 , x1 , . . . , xn , x̄ de menor a
mayor, podemos aplicar el teorema de Rolle a la función g(t) en cada intervalo definido por dos puntos
sucesivos (la función g(t) coincide en valor en los extremos de cada intervalo ya que vale cero en ambos
puntos). Entonces, podemos afirmar que g 0 (t) se anula en al menos un punto en cada intervalo. Por lo
tanto podemos afirmar que g 0 (t) se anula en al menos n + 1 puntos.

Si este mismo razonamiento lo aplicamos ahora a la función g 0 (t) en los intervalos definidos por los
n + 1 puntos donde se anula , llegaremos a la conclusión que g 00 (t) se anula en al menos n puntos.

Repitiendo el proceso, llegaremos a que g n+1 (t) se anula en al menos 1 punto. Este punto depende de
los valores x0 , x1 , . . . , xn , x̄. Llamemos ξ(x̄) a dicho punto.

Volvamos ahora a la expresión de g(t):

n

Y (t − xi )
g(t) = f (t) − P (t) − f (x̄) − P (x̄)
i=0
(x̄ − xi )

Desde aqui, derivando término a término, podemos obtener la expresión de la derivada de orden n + 1:

n

Y (t − xi ) n+1
g n+1 (t) = f n+1 (t) − P n+1 (t) − f (x̄) − P (x̄) ( )
i=0
(x̄ − xi )

Sabemos que P (t) es un polinomio de grado ≤ n, por lo tanto la deriva de orden n + 1 es cero. Además
n
Q (t−xi )
(x̄−xi ) es un polinomio de grado n + 1, por lo tanto la deriva de orden n + 1 es igual al coeficiente
i=0
n
1
Q
que acompaña a la potencia de orden n + 1 (que vale (x̄−xi ) ), multiplicada por (n + 1)!
i=0
De estas observaciones, deducimos que:

n
Y 1
g n+1 (t) = f n+1 (t) − f (x̄) − P (x̄) (n + 1)!(


)
i=0
(x̄ − xi )
Métodos Numéricos modo virtual (pandemia COVID-19)

Si ahora evaluamos la expresión anterior en ξ(x̄), tendremos que

n
Y 1
g n+1 (ξ(x̄)) = 0 = f n+1 (ξ(x̄)) − f (x̄) − P (x̄) (n + 1)!(


)
i=0
(x̄ − xi )

n
Y 1
f n+1 (ξ(x̄)) = f (x̄) − P (x̄) (n + 1)!(


)
i=0
(x̄ − xi )

n
f n+1 (ξ(x̄)) Y

(x̄ − xi ) = f (x̄) − P (x̄)
(n + 1)! i=0

n
f n+1 (ξ(x̄)) Y
P (x̄) + (x̄ − xi ) = f (x̄)
(n + 1)! i=0


 Métodos Numéricos modo virtual
(pandemia COVID-19)
Material Complementario

Interpolación - Diferencias Divididas - versión 1.0

Este es material complementario de las diapos de la clase de interpolación usadas durante el dictado virtual
(pandemia COVID-19). En este documento deducimos la expresión del polinomio interpolante mediante el
uso de diferencias divididas.

Sean f (x) una función definida en un intervalo [a, b] y pares ordenados (xi , f (xi )), xi ∈ [a, b] para i = 0, . . . , n.
Sabemos que existe un polinomio P (x) de grado ≤ n tal que P (xi ) = f (xi ) para todo i = 0, . . . , n.
n n
P Q (x−xi )
La expresion para este (único!) polinomio es P (x) = yk Lnk (x) donde Lnk = (xk −xi ) .
k=0 i=0,i6=k
En el caso que se agregara un punto más al conjunto de los puntos de interpolación, deberı́amos rehacer la
expresión de cada término. ¿Cómo podremos evitar este trabajo?
Definimos las Diferencias Divididas como

ˆ Orden 0 : f [xi ] = f (xi )

f [xi+1 ]−f [xi ]
ˆ Orden 1 : f [xi , xi+1 ] = xi+1 −xi

f [xi+1 ,...,xi+k ]−f [xi ,...xi+k−1 ]

ˆ Orden k : f [xi , xi+1 , . . . , xi+k ] = xi+k −xi

Veamos que el polinomio interpolante se puede expresar en función de estas diferencias

Proposición: Dada f (x) un función definida en [a, b] y pares ordenados (xi , f (xi )), xi ∈ [a, b] para
i = 0, . . . , n, el polinomio interpolante se puede expresar como

P (x) = f [x0 ] + f [x0 , x1 ] (x − x0 ) + f [x0 , x1 , x2 ] (x − x0 )(x − x1 ) + . . . + f [x0 , . . . , xn ] (x − x0 ) . . . (x − xn−1 )

Demostración:
Haremos la demostración por inducción en n.

ˆ Caso base: n=1 Los puntos de interpolación son x0 y x1 y el polinomio interpolante tiene grado ≤ 1.
Por la expresión del polinomio en función de los Lnk , tenemos que

(x − x1 ) (x − x0 )
P (x) = f (x0 ) + f (x1 )
(x0 − x1 ) (x1 − x0 )

1
Métodos Numéricos modo virtual (pandemia COVID-19)
Sumando y restando x0 en el primer término

(x − x0 + x0 − x1 ) (x − x0 )
P (x) = f (x0 ) + f (x1 )
(x0 − x1 ) (x1 − x0 )

(x0 − x1 ) + (x − x0 ) (x − x0 )
P (x) = f (x0 ) + f (x1 )
(x0 − x1 ) (x1 − x0 )

(x0 − x1 ) (x − x0 ) (x − x0 )
P (x) = f (x0 ) + f (x0 ) + f (x1 )
(x0 − x1 ) (x0 − x1 ) (x1 − x0 )

Simplificando en el primer término y sacando factor común (x − x0 ) entre los dos últimos, obtenemos:

(f (x1 ) − f (x0 ))
P (x) = f (x0 ) + (x − x0 )
(x1 − x0 )

Usando las definiciones de las diferencias divididas obtenemos la expresión de P (x) en función de ellas:

P (x) = f [x0 ] + f [x0 , x1 ](x1 − x0 )

ˆ Paso inductivo
Sea Pn (x) el polinomio interpolante en los puntos x0 , . . . , xn , es decir Pn (xi ) = f (xi ). Por hipótesis
inductiva,

Pn (x) = f [x0 ]+f [x0 , x1 ] (x−x0 )+f [x0 , x1 , x2 ] (x−x0 )(x−x1 )+. . .+f [x0 , . . . , xn ] (x−x0 ) . . . (x−xn−1 )

Sea Qn (x) el polinomio interpolante en los puntos x1 , . . . , xn+1 , es decir Qn (xi ) = f (xi ). Por hipótesis
inductiva,

Qn (x) = f [x1 ]+f [x1 , x2 ] (x−x1 )+f [x1 , x2 , x3 ] (x−x1 )(x−x1 )+. . .+f [x1 , . . . , xn+1 ] (x−x1 ) . . . (x−xn )

Sea Pn+1 (x) el polinomio interpolante en los puntos x0 , . . . , xn+1 . Queremos ver que

Pn+1 (x) = f [x0 ] + . . . + f [x0 , . . . , xn ] (x − x0 ) . . . (x − xn−1 ) + f [x0 , . . . , xn+1 ] (x − x0 ) . . . (x − xn )

Nos construimos el polinomio P (x) = Pn (x) + a(x − x0 ) . . . (x − xn ). Veamos que propiedades tiene
P (x).
Claramente P (x) es un polinomio de grado ≤ n + 1 y además es fácil ver que P (xi ) = Pn (xi ) = f (xi )
para i = 0, . . . , n.
Por otro lado, eligiendo convenientemente a podemos conseguir que P (xn+1 ) = f (xn+1 ). ¿Cómo
hacemos esto? Si queremos que P (xn+1 ) = f (xn+1 ), entonces debe cumplirse que Pn (xn+1 ) + a(xn+1 −
f (xn+1 )−Pn (xn+1 )
x0 ) . . . (xn+1 − xn ) = f (xn+1 ). Basta tomar a = (xn+1 −x0 )...(xn+1 −xn ) que siempre está definido.

En conclusión, P (x) es un polinomio de grado ≤ n+1 que interpola en los puntos x0 , . . . , xn , xn+1 . Como
ya sabemos que el polinomio interpolante en un conjunto de puntos es único, entonces P (x) = Pn+1 (x).
Si demostramos que a = f [x0 , . . . , xn+1 ] entonces tendremos la propiedad requerida.
(x−xn+1 )
Consideremos un nuevo polinonio Q(x) = Qn (x) + (x n+1 −x0 )
(Qn (x) − Pn (x)). Por la expresión de Q(x)
deducimos que es un polinomio de grado ≤ n + 1

Vamos a evaluar a Q(x) en los puntos xi para todo i = 0, . . . , n + 1.

Métodos Numéricos modo virtual (pandemia COVID-19)

(xi − xn+1 )
Q(xi ) = Qn (xi ) + (Qn (xi ) − Pn (xi ))
(xn+1 − x0 )

Si i = 1, . . . , n, sabemos que xi es un punto de interpolación tanto para Qn (x) como para Pn (x).
Entonces Qn (xi ) − Pn (xi ) = 0, de donde se deduce que Q(xi ) = Qn (xi ) = f (xi ).
Si evaluamos en xi = xn+1 , el segundo término se anula y resulta que Q(xn+1 ) = Qn (xn+1 ) = f (xn+1 )
ya que xn+1 es un punto de interpolación para Qn (x).
0 −xn+1 )
Finalmente, si evaluamos en x0 , Q(x0 ) = Qn (x0 ) + (x
(xn+1 −x0 ) (Qn (x0 ) − Pn (x0 )) = Qn (x0 ) − (Qn (x0 ) −
Pn (x0 )) = Pn (x0 ). Como x0 es punto de interpolación para Pn (x), sabemos que Pn (x0 ) = f (x0 ) por lo
tanto resulta Q(x0 ) = f (x0 )

En resumen: Q(xi ) = f (xi ) para todo i = 0, . . . , n + 1. Pero entonces Q(x) = Pn+1 (x) ya que sabemos
que el polinomio interpolante es único.
Si dos polinomios son iguales, entonces los coeficientes que acompañan a cada potencia deben coincidir.
Recordemos la expresion de Pn+1 (x) y Q(x):

Pn+1 (x) = Pn (x) + a(x − x0 ) . . . (x − xn )

(x − xn+1 )
Q(x) = Qn (x) + (Qn (x) − Pn (x))
(xn+1 − x0 )

El coeficiente que acompaña a la potencia n + 1 de Pn+1 (x) es a que es el que queremos demostrar que
vale f [x0 , . . . , xn+1 ].
El coeficiente que acompaña a la potencia n + 1 de Q(x) es el coeficiente de la potencia n de Qn (x),
menos el coeficiente de la potencia n de Pn (x), dividido por (xn+1 − x0 ).
El coeficiente de la potencia n de Qn (x), por hipótesis inductiva es f [x1 , . . . , xn+1 ].
El coeficiente de la potencia n de Pn (x), por hipótesis inductiva es f [x0 , . . . , xn ].
f [x1 ,...,xn+1 ]−f [x0 ,...,xn ]
Entonces a = (xn+1 −x0 ) que es la definición de f [x0 , . . . , xn+1 ].


 Métodos Numéricos modo virtual
(pandemia COVID-19)
Material Complementario

Interpolación - Proceso recursivo (Neville) - versión 1.0

Este es material complementario de las diapos de la clase de interpolación usadas durante el dictado virtual
(pandemia COVID-19). En este documento deducimos la expresión del polinomio interpolante mediante el
uso de la recursión e Neville.

Sean f (x) una función definida en un intervalo [a, b] y pares ordenados (xi , f (xi )), xi ∈ [a, b] para i = 0, . . . , n.
Sabemos que existe un polinomio P (x) de grado ≤ n tal que P (xi ) = f (xi ) para todo i = 0, . . . , n.
n n
P Q (x−xi )
La expresion para este (único!) polinomio es P (x) = yk Lnk (x) donde Lnk = (xk −xi ) .
k=0 i=0,i6=k
Vamos a mostrar otra manera de expresar a este polinomio. La idea es ver si podemos expresar a un
polinomio en función de otros dos polinomios que interpolan en un punto menos y como usar esto en un
proceso recursivo.

Definición:
Notamos Pm1 ,m2 ,...mk (x) al polinomio interpolante en los puntos xm1 , xm2 , . . . , xmk ∈ [a, b].

Proposición: Sean x0 , . . . , xk ∈ [a, b], i, j ∈ {0, . . . , k}. El polinomio interpolante P0,...,k (x) puede expre-
sarse como:

(x − xj )P0,...,j−1,j+1,...,k (x) − (x − xi )P0,...,i−1,i+1,...,k (x)

P0,...,k (x) =
(xi − xj )

Demostración:
Observar que P0,...,j−1,j+1,...,k (x) interpola en los mismo puntos que P0,...,k (x) salvo xj y P0,...,i−1,i+1,...,k (x)
interpola en los mismo puntos que P0,...,k (x) salvo xi .
(x−xj )P0,...,j−1,j+1,...,k (x)−(x−xi )P0,...,i−1,i+1,...,k (x)
Debemos ver P (x) = (xi −xj ) , polinomio de grado ≤ k, interpola en
los puntos x0 , . . . , xk y dado que el polinomio interpolante es único, esto asegura la validez de la expresión.
Sea xr con r ∈ {1, . . . , k}, r 6= i, j. Por lo tanto xr es un punto interpolante tanto para P0,...,j−1,j+1,...,k (x)
como para P0,...,i−1,i+1,...,k (x).

(xr − xj )P0,...,j−1,j+1,...,k (xr ) − (xr − xi )P0,...,i−1,i+1,...,k (xr )

P (xr ) =
(xi − xj )

1
Métodos Numéricos modo virtual (pandemia COVID-19)

(xr − xj )f (xr ) − (xr − xi )f (xr )

P (xr ) = = f (xr )
(xi − xj )
Consideremos ahora xi . Sabemos que xi es punto interpolante para P0,...,j−1,j+1,...,k (x). Entonces:

(xi − xj )P0,...,j−1,j+1,...,k (xi ) − (xi − xi )P0,...,i−1,i+1,...,k (xi )

P (xi ) =
(xi − xj )
(xi − xj )P0,...,j−1,j+1,...,k (xi )
P (xi ) = = f (xi )
(xi − xj )
Por último, consideremos xj . Sabemos xj es punto interpolante para P0,...,i−1,i+1,...,k (x). Entonces:

(xj − xj )P0,...,j−1,j+1,...,k (xj ) − (xj − xi )P0,...,i−1,i+1,...,k (xj )

P (xj ) =
(xi − xj )
(xj − xi )P0,...,i−1,i+1,...,k (xj )
P (xj ) = − = f (xi )
(xi − xj )
Entonces P (x) = P0,...,k (x)

Esto nos permite un proceso recursivo. Definimos Qij (x) = Pi−j,i−j+1,...,i (x) para i ≥ j. Tenemos que

(x − x0 )Qnn−1 (x) − (x − xn )Qn−1n−1 (x)

P0,...,n (x) = Qnn (x) = .
(xn − x0 )

Por la relación establecida, tenemos que:

(x − xi−j )Qij−1 (x) − (x − xi )Qi−1j−1 (x)

Qij =
(xi − xi−j )
lo que nos permite una recursión

Por ejemplo:

x0 f (x0 ) = Q00 (x)

Q11 (x)
x1 f (x1 ) = Q10 (x) Q22 (x)
Q21 (x) Q33 (x)
x2 f (x2 ) = Q20 (x) Q32 (x) Q44 (x)
Q31 (x) Q43 (x)
x3 f (x3 ) = Q30 (x) Q42 (x)
Q41 (x)
x4 f (x4 ) = Q40 (x)
 Métodos Numéricos modo virtual
(pandemia COVID-19)
Material Complementario

Interpolación segmentaria usando trazadores cúbicos - versión 1.0

Este es material complementario de las diapos de la clase de interpolación usadas durante el dictado virtual
(pandemia COVID-19). En este documento deducimos la existencia y unicidad de un trazador cúbico.

Sean f (x) una función definida en un intervalo [a, b] y pares ordenados (xi , f (xi )), xi ∈ [a, b] para i = 0, . . . , n.
Una trazador cúbico es un función S(x) tal que verifica la siguientes propiedades:

1. S(x) = Si (x) para x ∈ [xi , xi+1 ] con Si (x) = ai +bi (x−xi )+ci (x−xi )2 +di (x−xi )3 para i = 0, . . . , n−1

2. S(xi ) = f (xi ) para i = 0, . . . , n

3. Si (xi+1 ) = Si+1 (xi+1 ) para i = 0, . . . , n − 2
4. Si0 (xi+1 ) = Si+1
0
(xi+1 ) para i = 0, . . . , n − 2
5. Si00 (xi+1 ) = Si+1
00
(xi+1 ) para i = 0, . . . , n − 2

6. S 00 (x0 ) = S 00 (xn ) = 0 ó S 0 (x0 ) = f 0 (x0 ), S 0 (xn ) = f 0 (xn )

El objetivo del desarrollo que haremos a continuación es mostrar porque podemos asegurar que existe una
función que cumple con todas estas condiciones.
Notemos en primer lugar que tenemos 4 coeficientes a determinar para cada Si (x), lo que nos da un total de
4n coeficientes. La segunda propiedad nos impone n + 1 condiciones. La tercera, cuarta y quinta propiedad
imponen n − 1 condiciones cada una. Tenemos entonces un total de n + 1 + n − 1 + n − 1 + n − 1 = 4n − 2
condiciones. La última propiedad aporta 2 condiciones. Por lo tanto, tenemos tantas condiciones como
coeficientes a determinar. Debemos ver que existen coeficientes que satisfacen todas estas condiciones.
Analicemos cada una de estas condiciones. Comenzamos con S(xi ) = f (xi ) para i = 0, . . . , n. Como
S(x) = Si (x) para x ∈ [xi , xi+1 ], entonces tendremos que:

S(xi ) = ai + bi (xi − xi ) + ci (xi − xi )2 + di (xi − xi )3 = f (xi ) ∀i = 0, . . . , n − 1

S(xn ) = an−1 + bn−1 (xn − xn−1 ) + cn−1 (xn − xn−1 )2 + dn−1 (xn − xn−1 )3 = f (xn )
De aqui derivamos que
ai = f (xi ) ∀i = 0, . . . , n − 1
an−1 + bn−1 (xn − xn−1 ) + cn−1 (xn − xn−1 )2 + dn−1 (xn − xn−1 )3 = f (xn )

1
Métodos Numéricos modo virtual (pandemia COVID-19)

bn−1 (xn − xn−1 ) + cn−1 (xn − xn−1 )2 + dn−1 (xn − xn−1 )3 = f (xn ) − f (xn−1 )
La próxima condición es Si (xi+1 ) = Si+1 (xi+1 ) para i = 0, . . . , n − 2. Considerando la expresión de cada
Si (x), tenemos la siguiente relación:

ai +bi (xi+1 −xi )+ci (xi+1 −xi )2 +di (xi+1 −xi )3 = ai+1 +bi+1 (xi+1 −xi+1 )+ci+1 (xi+1 −xi+1 )2 +di+1 (xi+1 −xi+1 )3

ai + bi (xi+1 − xi ) + ci (xi+1 − xi )2 + di (xi+1 − xi )3 = ai+1 para i = 0, . . . , n − 2

f (xi ) + bi (xi+1 − xi ) + ci (xi+1 − xi )2 + di (xi+1 − xi )3 = f (xi+1 ) para i = 0, . . . , n − 2

0
La cuarta condición es Si0 (xi+1 ) = Si+1 (xi+1 ) para i = 0, . . . , n − 2

bi + 2ci (xi+1 − xi ) + 3di (xi+1 − xi )2 = bi+1 + 2ci+1 (xi+1 − xi+1 ) + 3di+1 (xi+1 − xi+1 )2

bi + 2ci (xi+1 − xi ) + 3di (xi+1 − xi )2 = bi+1

La quinta condición es Si00 (xi+1 ) = Si+1
00
(xi+1 ) para i = 0, . . . , n − 2

2ci + 6di (xi+1 − xi ) = 2ci+1 + 6di+1 (xi+1 − xi+1 )

2ci + 6di (xi+1 − xi ) = 2ci+1

Finalmente, analicemos una de las dos últimas alternativas: S 00 (x0 ) = S 00 (xn ) = 0 (la otra alternativa es
similar)
S 00 (x0 ) = S000 (x0 ) = 2c0 = 0
S 00 (xn ) = Sn−1
00
(xn ) = 2cn−1 + 6dn−1 (xn − xn−1 ) = 0
Veamos entonces todas las condiciones que nos quedaron:

1. ai = f (xi ) para i = 0, . . . , n − 1
2. bn−1 (xn − xn−1 ) + cn−1 (xn − xn−1 )2 + dn−1 (xn − xn−1 )3 = f (xn ) − f (xn−1 )

3. f (xi ) + bi (xi+1 − xi ) + ci (xi+1 − xi )2 + di (xi+1 − xi )3 = f (xi+1 ) para i = 0, . . . , n − 2

4. bi + 2ci (xi+1 − xi ) + 3di (xi+1 − xi )2 = bi+1 para i = 0, . . . , n − 2
5. 2ci + 6di (xi+1 − xi ) = 2ci+1 para i = 0, . . . , n − 2

6. c0 = 0
7. 2cn−1 + 6dn−1 (xn − xn−1 ) = 0

La idea de lo que vamos a hacer a continuación es tratar de poner a todas las variables en función de los
coeficientes ai que ya tenemos determinados y de los ci . Notamos hi = (xi − xi−1 ) para i = 1, . . . , n.
cn−1
De (7) podemos despejar dn−1 → dn−1 = − 3hn
.
(f (xn )−f (xn−1 )−cn−1 h2n −dn−1 h3n )
De (2) podemos despejar bn−1 → bn−1 = hn . Reemplazando la expresión que
cn−1 3
(f (xn )−f (xn−1 )−cn−1 h2n + 3h h ) (f (xn )−f (xn−1 ))
ya tenemos de dn−1 , obtenemos bn−1 = hn
n n
, bn−1 = hn
− 23 cn−1 hn
(2ci+1 −2ci )
De (5) podemos despejar di → di = 6hi+1
para i = 0, . . . , n − 2
Métodos Numéricos modo virtual (pandemia COVID-19)

(f (xi+1 )−f (xi )−ci h2i+1 −di h3i+1 )

De (3) podemos despejar bi → bi = hi+1 para i = 0, . . . , n − 2. Reemplazando la
expresión de di , obtenemos
(2ci+1 −2ci ) 3
(f (xi+1 ) − f (xi ) − ci h2i+1 − 6hi+1 hi+1 )
bi =
hi+1

(f (xi+1 ) − f (xi )) (2ci+1 − 2ci )

bi = − ci hi+1 − hi+1
hi+1 6
(f (xi+1 ) − f (xi )) 2 ci+1
bi = − ci hi+1 − hi+1
hi+1 3 3
Finalmente, vamos a usar (4). Por un lado lo hacemos para i = 0, . . . , n − 3

bi + 2ci (xi+1 − xi ) + 3di (xi+1 − xi )2 = bi+1

Reemplazamos la expresión que tenemos de bi , bi+1 y di

(f (xi+1 ) − f (xi )) 2 ci+1 (2ci+1 − 2ci ) (f (xi+2 ) − f (xi+1 )) 2 ci+2

− ci hi+1 − hi+1 +2ci hi+1 +3h2i+1 = − ci+1 hi+2 − hi+2
hi+1 3 3 6hi+1 hi+2 3 3

2 1 2 1 (f (xi+2 ) − f (xi+1 )) (f (xi+1 ) − f (xi ))

ci (− hi+1 +2hi+1 −hi+1 )+ci+1 (− hi+1 +hi+1 + hi+2 )+ci+2 ( hi+2 ) = −
3 3 3 3 hi+2 hi+1

1 2 1 (f (xi+2 ) − f (xi+1 )) (f (xi+1 ) − f (xi ))

ci ( hi+1 ) + ci+1 ( (hi+1 + hi+2 )) + ci+2 ( hi+2 ) = −
3 3 3 hi+2 hi+1
Nos queda el caso i = n − 2:

bn−2 + 2cn−2 (xn−1 − xn−2 ) + 3dn−2 (xn−1 − xn−2 )2 = bn−1

Reemplazamos la expresión que tenemos de bn−2 ,bn−1 y dn−2

(f (xn−1 ) − f (xn−2 )) 2 cn−1 (2cn−1 − 2cn−2 ) 2 (f (xn ) − f (xn−1 )) 2

− cn−2 hn−1 − hn−1 +2cn−2 hn−1 +3 hn−1 = − cn−1 hn
hn−1 3 3 6hn−1 hn 3

2 1 2 (f (xn ) − f (xn−1 )) (f (xn−1 ) − f (xn−2 ))

cn−2 (− hn−1 + 2hn−1 − hn−1 ) + cn−1 (− hn−1 + hn−1 + hn ) = −
3 3 3 hn hn−1

1 2 (f (xn ) − f (xn−1 )) (f (xn−1 ) − f (xn−2 ))

cn−2 ( hn−1 ) + cn−1 (hn + hn−1 ) = −
3 3 hn hn−1
En definitiva tenemos las siguientes n ecuaciones que involucran a los n coeficientes c0 , . . . , cn−1

c0 = 0
(f (xi+2 )−f (xi+1 )) (f (xi+1 )−f (xi ))
ci ( 13 hi+1 ) + ci+1 23 (hi+1 + hi+2 ) + ci+2 ( 31 hi+2 ) = hi+2
− hi+1
para i=0,. . . ,n-3
(f (xn )−f (xn−1 )) (f (xn−1 )−f (xn−2 ))
cn−2 ( 13 hn−1 ) + cn−1 23 (hn + hn−1 ) = hn
− hn−1

La matriz asociada al sistema es

Métodos Numéricos modo virtual (pandemia COVID-19)

 
c0 c1 c2 c3 ... ci ci+1 ci+2 ... cn−2 cn−1
 1 0 0 0 ... 0 0 0 ... 0 0 
1 2 1

 h1 (h1 + h2 ) 3 h2 0 ... 0 0 0 ... 0 0 
3 3
1 2 1

 0
3 h2 (h2 + h3 ) 3 h3 ... 0 0 0 ... 0 0 
3
 .. .. .. .. .. .. .. .. ..
 

 . . . . ... . . . ... . . 
 1 2 1

 0 0 0 0 ... 3 hi+1 3 (hi+1 + hi+2 ) 3 hi+2 ... 0 0 
 
 . .. .. .. .. .. .. .. ..
 ..

. . . ... . . . ... . . 
1 2
0 0 0 0 ... 0 0 0 ... 3 hn−1 3 (hn + hn−1 )

que resulta ser estrictamente diagonal dominante, por lo cual existe solución del sistema y la solución es
única. De esta manera obtenemos en forma única los coeficientes c0 , . . . , cn−1 . Dado que el resto de los
coeficientes se encuentran expresados en función de c0 , . . . , cn−1 , podemos afirmar que el trazador cúbico
existe y es único.