Taylor PDF
Taylor PDF
Taylor PDF
Practica 8
Polinomio de Taylor
1. Polinomio de Taylor
El analisis completo de una funcion puede resultar muy difcil. Una forma de abordar este
problema es aproximar la funcion por una mas sencilla. En este caso vamos a aproximar
las funciones por polinomios. Dicha aproximacion se hace cerca de un valor concreto y solo
servira para valores cercanos. A medida que nos alejemos, la aproximacion sera menos con-
fiable y es posible que el polinomio se aleje mucho de la funcion bajo estudio.
En la version mas sencilla podemos aproximar la funcion a estudiar por una constante. Es
claramente mucho mas facil operar con un numero (la constante) que con una funcion que
puede ser complicada. Esta aproximacion se usa mucho en las aplicaciones. Por ejemplo, la
aceleracion de la gravedad cerca de la superficie de la tierra se aproxima por la constante
m m
g = 9,8 2 (o a veces directamente g 10 2 ), por mas que en realidad es una funcion que
s s
g
depende de la altura. Otro ejemplo es la densidad del agua que se aproxima por = 1 3
cm
por mas que dependa de la presion y la temperatura.
Para otras aplicaciones se necesitan aproximaciones mas precisas, como una funcion lineal.
Por ejemplo, dados una funcion f y un valor de x podemos tomar la recta tangente al grafico
de f , que ya hemos calculado. En fsica, un caso tpico es aproximar la longitud de un ob-
jeto por un valor constante mas una correccion lineal pequena que depende linealmente de
la temperatura. Para mejorar la precision se agregan coeficientes cuadraticos o cubicos, a
medida que es necesario.
Vamos a analizar otro ejemplo. Supongamos que queremos calcular aproximadamente el
valor de 65. Como primera aproximacion, podemos tomar que 65 64 = 8. Muchas
veces con ese nivel de precision alcanza. En realidad, el valor es un poco mas grande que
8, informalmente decimos que es 8 y pico. Para poder aproximar este numero con mas
precision podemos usar la recta tangente a x en el punto correspondiente a x = 64. Es es-
perable que cerca de ese punto la recta tangente sea una mejor aproximacion que la funcion
constante. Esto no solo nos da un aproximacion de 65, sino que tenemos una aproximacion
de x cuando x vale aproximadamente 64.
Luego vamos a ver como extender esta idea a polinomios de grado mas alto. Tambien vamos
a estimar el error que se comete en cada aproximacion, de manera que podamos saber si la
1
P1 ( x ) = 8 + ( x 64) .
16
1
x + 4, pero si vamos
Distribuyendo, esta funcion tambien se podra escribir como P1 ( x ) =
16
a utilizar valores de x cercanos a 64 la primera expresion hace que las cuentas sean mas
faciles.
Por ejemplo, para calcular aproximadamente f (65) = 65 podemos calcular el valor de
1
P1 (65) = 8 + (65 64) = 8,0625.
16
Este valor es muy cercano a 65 8,0622...
1
P1 ( x ) = 8 + ( x 64)
16
estamos pidiendo que f (64) = P(64) para que ambas coincidan en el punto. Ademas pedi-
mos que f 0 (64) = P0 (64) para que las respectivas tangentes tengan la misma pendiente.
Para extenderlo a funciones cuadraticas en las que aproximamos el grafico de la funcion
por una parabola, podemos pedir que ambas tengan la misma curvatura en ese punto. La
curvatura es difcil de calcular, pero esta condicion es equivalente a pedir que tengan la
misma derivada segunda, o sea que f 00 (64) = P00 (64).
Si tomamos
P2 ( x ) = a + b( x 64) + c( x 64)2
f 00 (64)
y derivamos dos veces nos queda que P200 ( x ) = 2c, por lo que c = . Ademas, como
2
tenamos antes, a = f (64) y b = f 0 (64).
1 1 1
En nuestro ejemplo f 00 ( x ) = , por lo que f 00 (64) = yc= , por lo tanto
4 x3 2048 4096
1 1
P2 ( x ) = 8 + ( x 64) ( x 64)2 .
16 4096
Entonces, al estimar f (65) = 65 por
1 1 2
P2 (65) = 8 + 1 1 = 8,062256...
16 4096
tenemos una mejor aproximacion.
(Notemos que el valor de 65 es aproximadamente 8,062258....)
Lo bueno de escribir a
P2 ( x ) = a + b( x 64) + c( x 64)2
en vez de
P2 ( x ) = u + vx + wx2
10
P0 ( x )
8
P1 ( x )
5
P2 ( x )
f (x)
0
0 20 40 60 64 80 100
que cumpla las igualdades anteriores y que ademas f 000 (64) = P000 (64). Las formulas para a,
b y c no cambian. Al igualar las derivadas terceras tenemos que f 000 (64) = 3.2.1d = 6d. (Al
derivar baja primero un 3, despues un 2 y despues un 1.)
1 1 1
P3 ( x ) = 8 + ( x 64) ( x 64)2 + ( x 64)3
16 4096 524288
y reemplazando queda que
1 1 1
(65 64)
P3 (65) = 8 + (65 64)2 + (65 64)3
16 4096 524288
1 1 2 1
= 8+ 1 1 + 13 = 8,06225777...
16 4096 524288
(Notemos que en realidad 65 8,06225775....)
f ( x0 ) f 0 ( x0 ) f 00 ( x0 ) f 000 ( x0 ) f ( n ) ( x0 )
Pn ( x ) = + ( x x0 )+ ( x x0 )2 + ( x x0 )3 + + ( x x0 ) n
0! 1! 2! 3! n!
En la expresion, n! es el factorial de n, o sea el producto de todos los numeros naturales de 1
a n, o sea n! = 1.2.3. . . . .n y 0! = 1.
Tambien podemos escribir los primeros terminos directamente, usando que 0! = 1, 1! = 1,
2! = 2 y 3! = 1.2.3 = 6 y queda
f 00 ( x0 ) f 000 ( x0 ) f ( n ) ( x0 )
Pn ( x ) = f ( x0 )+ f 0 ( x0 )( x x0 )+ ( x x0 )2 + ( x x0 )3 + + ( x x0 ) n
2 6 n!
que es una expresion mas simple y mas parecida a la que vimos antes.
Para ver como usarlo, resolvamos otro ejemplo.
Solucion
f 0 ( x ) = 2e2x ,
f 00 ( x ) = 2.2e2x ,
f 000 ( x ) = 2.2.2e2x ,
...
f (5) ( x ) = 25 e2x .
4 8 16 32 5
P5 ( x ) = 1 + 2x + x2 + x3 + x4 + x
2 6 24 120
4 2 4
P5 ( x ) = 1 + 2x + 2x2 + x3 + x4 + x5 .
3 3 15
1
Al aproximar f = e0,2 obtenemos
10
2 3 4 5
1 1 1 4 1 2 1 4 1
P5 = 1+2 +2 + + +
10 10 10 3 10 3 10 15 10
229013
= 1,22140266...
187500
Sabiendo que el valor exacto es e0,2 1,22140275... vemos que los valores son similares.
En cambio al aproximar f (1) = e2 obtenemos
4 2 4 109
P5 (1) = 1 + 2.1 + 2.12 + .13 + .14 + .15 = 7,266...
3 3 15 15
Sabiendo que el valor exacto es e2 7,389... vemos que la aproximacion tiene una mayor
diferencia. 2
Para poder estimar estas diferencias y ver si la aproximacion nos resulta util, vamos a desa-
rrollar una expresion para el error que se comete al usar el polinomio de Taylor para apro-
ximar la funcion.
R3 ( x ) R3 ( x ) R3 (0) R30 ( )
= =
x4 x4 04 4. 3
en donde es un punto intermedio, que esta en el intervalo (0; x ). (En general no se puede
determinar facilmente. Solo sabemos que existe.)
Con esta misma idea, tomando la funcion 4x3 tenemos que
en donde ahora esta en el intervalo (0; ). Para simplificar esta expresion, podemos usar
que el intervalo (0; ) esta incluido en el intervalo (0; x ) y decir directamente que esta en
el intervalo (0; x ).
Repetimos una vez mas el razonamiento con la funcion 4.3.x2 y entonces
en donde esta en el intervalo (0; ), o para simplificar tomamos que esta en el intervalo
(0; x ).
Repetimos una ultima vez mas el procedimiento usando ahora la funcion 4.3.2.x y entonces
(4)
R3000 ( ) R3000 ( ) R3000 (0) R (c)
= = 3
4.3.2. 4.3.2. ( 0) 4.3.2.1
en donde c esta en el intervalo (0; ), o para simplificar, tomamos que c esta en el intervalo
(0; x ).
Pero como P3 es un polinomio de grado 3 a lo sumo, entonces al derivarlo cuatro veces se
(4)
anula. As que P3 ( x ) = 0 para todo x y por ello
(4) (4) (4) (4)
R3 (c) = f 3 (c) P3 (c) = f 3 (c).
Solucion
Calculemos las derivadas
f ( x ) = ln(1 + x ),
1
f 0 (x) = ,
1+x
1
f 00 ( x ) = ,
(1 + x )2
2
f 000 ( x ) = ,
(1 + x )3
2.3
f (4) ( x ) = .
(1 + x )4
As que f (0) = 0, f 0 (0) = 1, f 00 (0) = 1, f 000 (0) = 2. Entonces,
1 1 137
P3 (0,2) = 0,2 0,22 + 0,23 = 0,18267...
2 3 750
1
la diferencia es menor que = 0,0004 . 2
2500
Generalizando el ejemplo anterior:
Solucion
Para el primer caso, tenemos que
6
f (6) ( c ) 6 26 e2c 6 26 e2c e2c
1
R5 (0,1) = 0,1 = 0,1 = =
6! 6! 6! 10 11250000
En esta acotacion, e0,2 es mucho menor que 3, pero la acotamos de esta manera para que el
resultado sea mas facil de manejar.
Entonces,
2c
e 3 1
| R5 (0,1)| = < = .
11250000 11250000 3750000
1
As que | R5 (0,1)| < = 0,000000267... por lo que demostramos que los valores de
3750000
f (0,1) y P5 (0,1) son muy cercanos, como habamos visto numericamente antes.
Para el segundo caso, tenemos que
f (6) ( c ) 26 e2c 6 4 2c
| R5 (1)| = 16 = 1 = e
6! 6! 45
y entonces
4 2c 4 4
| R5 (1)| = e < 32 = 0,8.
45 45 5
As que tenemos una acotacion del error grande y probablemente sea una mala idea usar
P5 (1) como una aproximacion de f (1). Esto es coherente con los valores que obtuvimos
antes numericamente para comparar. 2
Solucion
f ( x ) = sen( x ),
f 0 ( x ) = cos( x ),
f 00 ( x ) = sen( x ),
f 000 ( x ) = cos( x ),
f (4) ( x ) = sen( x ).
4 (4) 4
R3 ( ) = f (c) = sen(c)
4! 24
en donde c esta en el intervalo (0; ). Como sen( x ) no es creciente ni decreciente en este
intervalo, hay que tener mas cuidado.
Si cometemos el error de evaluar solo en ambos extremos tenemos que sen(0) = 0 y que
sen( ) = 0. Es importante recordar que no hay que evaluar en los extremos, sino acotar la
funcion. Por suerte, sen( x ) esta acotada entre 1 y 1. Por ello,
4 4 44 32
| R3 ( )| = |sen(c)| 1 < =
24 24 24 3
en donde usamos la cota de < 4 para obtener una expresion mas sencilla. 2
1
P3 ( ) = 3 2,0261
6
la diferencia claramente no es 0. Es mas, es una mala aproximacion. La acotacion del resto en
32
este caso da | R3 ( )| < 10,667..., lo cual es coherente con la diferencia que obtuvimos.
3
En general, estas aproximaciones son utiles cuando se puede dar una cota chica del error.
Solo para comparar, analicemos la aproximacion en otro valor facil de calcular, pero que este
mas cercano a x0 = 0.
1
Si utilizamos el mismo polinomio para estimar f = sen = queda
6 6 2
1 3 3
P3 = = 0,49967
6 6 6 6 6 1296
2
| R3 ( )| < 0,0082...
243
Solucion
Calculemos las derivadas
f ( x ) = xe x1 ,
f 0 ( x ) = e x1 + xe x1 ,
f 00 ( x ) = 2e x1 + xe x1 ,
f 000 ( x ) = 3e x1 + xe x1 .
f 00 (1)
P2 ( x ) = f (1) + f 0 (1)( x 1) + ( x 1)2
2
3
= 1 + 2 ( x 1) + ( x 1)2 .
2
Por un lado, como c < 1 tenemos que ec1 < e11 = e0 = 1. Ademas 3 + c < 3 + 1 = 4.
Entonces,
3ec1 + cec1 = (3 + c) ec1 < (3 + 1)e11 = 4.
Por ello,
1 3 1 3
10 10 1
| R3 (0,9)| = 3ec1 + cec1 < 4= .
6 6 1500
Con esto vemos que si aproximamos f (0,9) = 0,9e0,91 por
1 2
3 2 1 3 163
P2 (0,9) = 1 + 2 (0,9 1) + (0,9 1) = 1 + 2 + = ,
2 10 2 10 200
1
la acotacion del resto da | R2 (0,9)| < . 2
1500
3. Problemas varios
A veces es util encontrar una cota del error para todos los x en un intervalo cercano a x0 .
Solucion
La funcion f y x0 son iguales a los que utilizamos anteriromente. Por ello,
3
P2 ( x ) = 1 + 2 ( x 1) + ( x 1)2
2
y
| x 1 | 3 c 1 c 1
| x 1|3
| R3 ( x )| = 3e + ce = (3 + c ) e c 1
6 6
con x en el intervalo (0,9; 1,1) y c en el intervalo ( x; 1) o (1; x ) , dependiendo de que x < 1
o x > 1. En todos los casos, para cualquier valor de x, tenemos que c esta en el intervalo
(0,9; 1,1).
Para acotar, tenemos que por un lado si x esta en el intervalo (0,9; 1,1) entonces
1
| x 1| <
10
1
y entonces | x 1|3 < . Por otro lado, tenemos que
103
41
3 + c < 3 + 1,1 < 4,1 =
10
y que como e x es creciente
En esta ultima parte, e0,1 es mucho mas chico que 3, pero lo acotamos por 3 para simplificar
la expresion. Entonces,
1
| x 1|3 103 41 41
| R3 ( x )| = (3 + c ) e c 1 < 3 = .
6 6 10 20000
2
En algunos casos, puede ser util simplificar aun mas la expresion del error, por ejemplo
tomando
41 50 1
| R3 ( x )| =
< = .
20000 20000 400
Con cada acotacion de la cota perdemos un poco de informacion sobre la precision del
calculo, pero a cambio ganamos en simplicidad. Si la acotacion que utilizamos es muy burda
podemos perder de vista la calidad de la aproximacion.
Vimos en algunos ejemplos que no podamos acotar el error por un numero chico. En la
mayora de los casos para obtener aproximaciones mas precisas alcanza con aumentar el
orden del polinomio de Taylor utilizado.
Solucion
Primero, recordemos que f 0 ( x ) = 2e2x , f 00 ( x ) = 2.2e2x , ..., f (5) ( x ) = 25 e2x . En este caso, es
posible encontrar una formula general para la derivada n-esima de f , que es
f (n) ( x ) = 2n e2x .
De la misma manera, f (0) = 1, f 0 (0) = 2, f 00 (0) = 2.2, ..., f (5) (0) = 25 y obtenemos que
f ( n ) (0 ) = 2n .
Entonces,
4 8 2n
Pn ( x ) = 1 + 2x + x2 + x3 + + x n
2 6 n!
210+1 32 8 1
R10 (1) < = 0,00046... <
(10 + 1)! 17325 1000
Cintia Buxton, Lisi DAlfonso, Flora Gutierrez, Gabriela Jeronimo, Gustavo Massaccesi, Juan Carlos Pedraza y
Juan Sabia (2015), Polinomio de Taylor, Teoricas de Analisis Matematico (28).