MTODOS NUMRICOS

1.5 Serie de Taylor

Gustavo Rocha
2005-2
1.5 Serie de Taylor
1.5 Serie de Taylor
La serie de Taylor es, sin duda, el fundamento
matemtico ms importante para comprender, manejar y
formular mtodos numricos que se basan en la
aproximacin de funciones por medio de polinomios.
Aunque a veces no sea muy evidente, la mayora de los
mtodos numricos se basan en la aproximacin de
funciones por medio de polinomios.
1.5 Serie de Taylor
La expansin de Taylor de una funcin, es una serie infinita de potencias
que representa, de manera exacta, el comportamiento de la funcin en la
vecindad de un punto dado.
Si se ignoran todos los trminos de la serie de Taylor, excepto unos
cuantos, se obtiene un polinomio que aproxima a la funcin verdadera.
El error del mtodo numrico depende de la precisin con la que el
polinomio aproxima a a la funcin verdadera.
Los errores por truncamiento se evalan a travs de la comparacin del
desarrollo polinomial de la solucin numrica, con la serie de Taylor, de la
solucin exacta.
1.5.1 Expansin en serie de Taylor
Sea una funcin f(X) que tiene derivadas continuas hasta de orden
n en el punto X
i
, para el cual se conoce el valor de la funcin a
0
y el
de sus derivadas: a
1
, a
2
, a
3
, a
4
, a
n
,
f(x)
x
x
i
X
i+1

a
0

f(X
i+1
)
1.5.1 Expansin en serie de Taylor
Se trata de encontrar un polinomio de la forma:

_____ (1.13)

que permita predecir el valor de la funcin en un punto cualquiera X, en
trminos de la propia funcin y de sus derivadas en el punto X
i
.

El polinomio P(X) se hace coincidir con la funcin f(X), y las
primeras n derivadas del polinomio se hacen coincidir con las n
primeras derivadas de la funcin en el punto X
i
.


_____ (1.14)

i i
i
i i
(n) (n)
i i
P(X ) = f(X )
P'(Xi) = f'(X )
P''(X ) = f''(X )
...
P (X ) = f (X )
3 2 n
0 1 2 3 n
P(X) = a + a X + a X + a X + ... + a X + ...
1.5.1 Expansin en serie de Taylor
El valor de la funcin en un punto cualquiera X se puede evaluar a
travs de un polinomio equivalente al de la expresin (1.13):

____ (1.15)

Desarrollando la expresin (1.15) y comparndola con la expresin
(1.13), se obtiene:



_____(1.16)

2 3 n
0 1 i 2 i 3 i n i
f(X) = P(X) = b + b (X - X ) + b (X - X ) + b (X - X ) + ... + b (X - X ) + ...
2 3 4
0 0 1 i 2 i 3 i 4 i
2 3
1 1 2 i 3 i 4 i
2
2 2 3 4 i
n n
a = b - b X + b X - b X + b X - ...
a = b - 2b X + 3b X - 4b X + ...
a = b - 3b Xi + 6b X - ...
...
a = b - ...
1.5.1 Expansin en serie de Taylor
Las n primeras derivadas del polinomio son:


_____ (1.17)



Evaluando el polinomio y sus derivadas en el punto X
i
:



_____ (1.18)

2 n-1
1 2 i 3 i n i
n-2
2 3 i n i
n-3
3 n i
(n)
P'(X) = b + 2b (X - X ) + 3b (X - X ) + ... + nb (X - X ) + ...
P''(X) = 2b + 3 2b (X - X ) + ... + n(n-1)b (X - X ) + ...
P'''(X) = 3 2b + ... + n(n-1)(n-2)b (X-X ) + ...
...
P (X) = n(n-

n n
1)(n-2) ... 3 2 1b + ... = n!b + ...
i 0 0
i 1 1
i 2 2
(n)
i n
P(X ) = b 0!b
P'(X ) = b 1!b
P''(X ) = 2b = 2!b
...
P (X ) = n!b
=
=
1.5.1 Expansin en serie de Taylor
Considerando simultneamente las expresiones (1.14) y (1.18):



_______ (1.19)




Sustituyendo los valores de los coeficientes dados en (1.19) en la expresin (1.15):

_____ (1.20)



que en forma sinttica se expresa:

_____ (1.20')

0 i
1 i
2 i
n i
b = f(X )
b = f'(X )/1!
b = f''(X )/2!
...
b = f(n)(X )/n!
2
i i i i i
3 (n) n
i i i i
f(X) = f(X ) + f'(X )(X - X ) + f''(X )(X - X ) /2!
+ f'''(X )(X - X ) /3! + ... + f (X )(X - X ) /n! + ...
j
i i
j=0
f(X) = f(j)(X )(X - X ) /j!

1.5.1 Expansin en serie de Taylor

Las expresiones (1.20) y (1.20') son equivalentes y representan la
expansin en serie de Taylor que permite evaluar el valor de la funcin en
cualquier punto X, en trminos de la propia funcin y de sus derivadas en el
punto X
i
. Se pueden presentar dos casos:

A) Cuando el valor de X se encuentra a la derecha de X
i
, se usa la
nomenclatura X
i+1
, con lo que se indica que es mayor que X
i
.




_____ (1.21)


donde h se denomina tamao del paso, tratndose en este caso de un
paso hacia adelante.
i+1 i i+1 i
(j) (j)
j j
i+1 i i+1 i i
j=0 j=0
X = X > X ; X - X = h > 0

f(X ) = f (X )(X - X ) /j! = f (X )h /j!


1.5.1 Expansin en serie de Taylor
B) Cuando el valor de X se encuentra a la izquierda de X
i
, se usa
la nomenclatura X
i-1
, con lo que se indica que es menor que X
i
.


_____ (1.22)


_____ (1.22')

donde h es el tamao del paso, tratndose en este caso de un paso
hacia atrs.
Para cada combinacin de puntos X
i
, X
i+1
en una funcin f(x), la
serie de Taylor es nica, es decir, no hay otra serie de potencias en
h = X
i+1
X
i
, para representar a f(X)
i-1 i i i-1
(j) (j)
j j
i-1 i i i-1 i i i-1
j par j impar
(j) (j)
j j
i-1 i i
j par j impar
X = X < X ; X - X = h > 0

f(X ) = f (X )(X - X ) /j! - f (X )(X - X ) /j!
f(X ) = f (X )h /j! - f (X )h /j!


1.5.1 Expansin en serie de Taylor
Ejemplo. En el punto X
i
= 1, la funcin f(X) y sus derivadas toman
los siguientes valores:
f(1) = 1; f'(1) = 6; f''(1) = 2; f'''(1) = 6.

A partir de estos datos y utilizando la expansin en serie de Taylor
dada en (1.21), encontrar el polinomio que permita predecir valor de
la funcin para cualquier valor de X, y, en particular, el valor de la
funcin para X
i+1
= 3.
f(X) = 1 + 6(X - 1) + 2(X - 1)
2
/2! + 6(X - 1)
3
/3!
= 1 + 6X - 6 + X
2
- 2X + 1 + X
3
- 3X
2
+ 3X - 1
= - 5 + 7X - 2X
2
+ X
3
h = X
i+1
- X
i
= 3 - 1 = 2
f(X
i+1
) = f(3) = 1 + 6(2) + 2(2)
2
/2! + 6(2)
3
/3!
= 1 + 12 + 4 + 8 = 25
1.5.1 Expansin en serie de Taylor
Vamos a repetir el ejercicio, pero ahora considerando la expansin
en serie de Taylor dada en (1.22) y obteniendo el valor de la funcin
para X
i-1
= 0.
f(X) = 1 - 6(1 - X) + 2(1 - X)
2
/2! - 6(1 - X)
3
/3!
= 1 - 6 + 6X + X
2
- 2X + 1 - 1 + 3X - 3X
2
+ X
3

= - 5 + 7X - 2X
2
+ X
3
h = X
i
- X
i-1
= 1 - 0 = 1
f(X
i-1
) = f(0) = 1 - 6(1) + 2(1)
2
/2! - 6(1)
3
/3!
= 1 - 6 + 1 - 1 = - 5
1.5.1 Expansin en serie de Taylor
En el ejercicio anterior, el polinomio obtenido se ajusta
perfectamente a la funcin, porque sta es algebraica, polinomial de
tercer grado; en este caso, las derivadas de orden superior al
tercero se anulan, por lo que los primeros cuatro trminos de la
expansin en serie de Taylor son suficientes para determinar, sin
error alguno, el comportamiento de la funcin, para cualquier valor
de X.
Pero no siempre es as; cuando se trata de funciones trascendentes
o mixtas, la expansin en serie de Taylor slo puede proporcionar
una aproximacin a la funcin de inters, porque, en ese caso, cada
uno de los trminos de la serie infinita tiene un valor absoluto
diferente de cero, con el que participa, as sea de manera mnima,
en el valor de la funcin. En virtud de que no es posible considerar
un nmero infinito de trminos, no hay ms remedio que truncar la
serie y considerar nicamente los n primeros.
1.5.1 Expansin en serie de Taylor
Una funcin f(x) es analtica en x
i
, si se puede representar por
medio de una serie de potencias en trminos de h = x
i+i
x
i
, dentro
de un radio de convergencia 0 < ,x
i+i
- x
i
,, y si todas sus derivadas
son continuas en la vecindad de x
i
. Los polinomios son funciones
analticas en todas partes.
Si la funcin f(x) es diferenciable en todas partes de la vecindad de
un punto x
0
, excepto en el mismo, el punto se denomina singular y
entonces la funcin no es analtica en x
0
. Algunas funciones
trascendentes tienen puntos singulares; por ejemplo, tan(x) es
analtica excepto en (n + )t.

1.5.1 Expansin en serie de Taylor
Ejemplo. Aproximar la funcin f(X) = cos X en 30, conociendo los
valores de la funcin y el de sus derivadas para 0 y considerando los
primeros siete trminos de la expansin en serie de Taylor. No
olvidemos trabajar en radianes:

X
i
= 0 = 0 ; X
i+1
= 30 = t/6 ; h = X
i+1
- X
i
= t /6 - 0 = t /6

f(X) = f(X
i
) + f'(X
i
)h + f''(X
i
)h
2
/2! + f'''(X
i
)h
3
/3! + f
iv
(X
i
)h
4
/4! + f
v
(X
i
)h
5
/5! + f
vi
(X
i
)h
6
/6!

f(X) = cos X f(0) = cos 0 = 1
f'(X) = - sen X f'(0) = - sen 0 = 0
f''(X) = - cos X f''(0) = - cos 0 = - 1
f'''(X) = sen X f'''(0) = sen 0 = 0
f
iv
(X) = cos X f
iv
(0) = cos 0 = 1
f
v
(X) = - sen X f
v
(0) = - sen 0 = 0
f
vi
(X) = - cos X f
vi
(0) = - cos 0 = - 1
1.5.1 Expansin en serie de Taylor
Ejemplos: Los desarrollos en serie de Taylor de e
-x
y de sen x, en la
vecindad de x = 1, son respectivamente:
El desarrollo en serie de Taylor de una funcin alrededor de x = 0 recibe el
nombre de serie de Maclaurin; por ejemplo: e
x
, cos x, y ln(x+1)
2 3 4
-x -1 -1 -1 -1 -1
2 3 4
h h h
e = e - he + e - e + e - ...
2! 3! 4!
h h h
sen(x) = sen(1) + h cos(1) sen(1) cos(1) sen(1) ...
2! 3! 4!
+ +
2 3 4
x
2 4 6 8
2 3 4
x x x
e = 1 + x + + + + ...
2! 3! 4!
x x x x
cos(x) = 1 ...
2! 4! 6! 8!
x x x
ln(x 1) x + + + + ...
2 3 4
+ +
+ =
1.5.1 Expansin en serie de Taylor
f(t/6) = 1 - 1(t/6)
2
/2! + 1(t/6)
4
/4! - 1(t/6)
6
/6!
= 1 - 0.1370778 + 0.0031317 - 0.0000286 = 0.8660252
Considerando como "verdadero" el valor que ofrece una calculadora
cientfica de 8 dgitos, que es: cos 30 = 0.8660254, se aprecia que
el truncamiento a siete trminos de la serie, conduce a un pequeo
error de 0.0000002
1.5.2 El residuo de la serie de Taylor
En la seccin 1.4 se esboz lo que era un error por truncamiento,
pero no qued lo suficientemente claro, porque para comprender
este concepto, faltaba conocer a detalle el comportamiento de la
expansin en serie de Taylor.
Ahora podemos entender con claridad qu es un truncamiento y
cmo repercute ste en un error, al aproximar el valor de una
funcin para un determinado valor de la variable, considerando
solamente los primeros n trminos de la serie infinita.
Los trminos de la serie que se desprecian constituyen un residuo
cuyo valor puede tener signo positivo, en detrimento del valor de la
funcin, o negativo, en profusin del valor de la funcin; en trminos
absolutos, este residuo puede ser significativo o insignificante
(como sucedi en el ejemplo anterior), lo cual depende de dos
factores:
1.5.2 El residuo de la serie de Taylor
1) El valor de n, es decir, el nmero de trminos de la serie,
considerados al aproximar el valor de la funcin; mientras mayor
sea el valor de n, menor ser el residuo y mejor ser la
aproximacin al valor de la funcin.
2) El valor de h, es decir, el tamao del paso o distancia entre el
valor de la variable para el cual se evala la funcin y el valor de la
variable para el que se conoce el valor de la funcin y el de sus
derivadas; mientras menor sea el valor de h, mayor ser la cercana
entre X
i
y X
i+1
y, por ende, mejor ser la aproximacin al valor de la
funcin.
1.5.2 El residuo de la serie de Taylor
En adelante, y en tanto no se indique lo contrario, usaremos nicamente la
expansin en serie de Taylor que considera el paso hacia adelante para
aproximar f(X
i+1
) a partir de f(X
i
) y sus derivadas, conforme a la expresin
(1.21), la que en forma explcita se escribe:
f(Xi+1) = f(Xi) + f'(Xi)h + f''(Xi)h
2
/2! + f'''(Xi)h
3
/3! + ... + f(n)(Xi)h
n
/n! + ...__ (1.21')
y en forma alternativa:
f(Xi+1) = f(Xi) + f'(Xi)h + f''(Xi)h
2
/2! + f'''(Xi)h
3
/3! + ... + f(n)(Xi)h
n
/n! + Rn (1.23)
Esta ltima expresin se conoce como expansin en serie de Taylor con
residuo, y es idntica a la expresin (1.21'), excepto porque los puntos
suspensivos se han sustituido por el trmino R
n
, que sintetiza los trminos
de la serie que se han despreciado y se conoce con el nombre de residuo
de la aproximacin al n-simo orden.
La serie se puede truncar en cualquier punto, de manera que el subndice n
indica que slo se han incluido en la aproximacin los primeros (n+1)
trminos de la serie.
1.5.2 El residuo de la serie de Taylor
Por ejemplo, podemos truncar la serie a un solo trmino (n = 0):
f(X
i+1
) ~ f(X
i
)
lo que implica suponer que la funcin que se va a aproximar es una
constante: P(X) = a
0
; si tal suposicin es cierta, la aproximacin resulta
perfecta y no hay error alguno, pero si no es as, existe un residuo R
0
tal
que se cumple:
f(X
i+1
) = f(X
i
) + R
0
R
0
= f'(X
i
)h + f''(X
i
)h
2
/2! + f'''(X
i
)h
3
/3! +...+ f
(n)
(X
i
)h
n
/n! +... _____ (1.24)
R
0
es el residuo de orden cero y representa una serie infinita idntica a la
de la expresin (1.21'), excepto por la exclusin del primer trmino.
Para simplificar, podramos truncar el residuo a solo un trmino: R
0
f'(X
i
)h,
despreciando todos los dems, pero esto obviamente no es exacto.
Conviene entonces encontrar una manera ms adecuada de valorar R
0
.
1.5.2 El residuo de la serie de Taylor
Auxilindonos en la siguiente figura, podemos ver fcilmente que la
recta que une los puntos [X
i
, f(X
i
)], [X
i+1
,f(X
i+1
)], tiene pendiente R
0
/h.
f(x)
x
x
i
X
i+1

f(x
i
)
f(X
i+1
)

f()
h
R
0
= f(X
i+1
) - f(x
i
)
P(X) = a
o

1.5.2 El residuo de la serie de Taylor
Invocando el teorema del valor medio, podemos asegurar que existe un
punto , entre X
i
y X
i+1
, para el cual el valor de la primera derivada f'(), es
decir, la pendiente de la tangente de la funcin en ese punto, es paralela a
la recta mencionada previamente: R
0
/h = f'(); y entonces:
R
0
= f'()h _____ (1.25)
De manera similar, si truncamos la serie a dos trminos (n=2):
f(Xi+1) f(X
i
) + f'(X
i
)h
estaremos suponiendo que la funcin que se va a aproximar es una recta:
P(X) = a
0
+ a
1
X; si la suposicin es correcta, la aproximacin es perfecta y
sin error, pero si no es as, existe un residuo R
1
tal que:
f(X
i+1
) = f(X
i
) + f'(X
i
)h + R
1
R
1
= f''(X
i
)h
2
/2! + f'''(X
i
)h
3
/3! + ... + f
(n)
(X
i
)h
n
/n! + ...
R
1
es un residuo de primer orden que, al igual que se hizo con R
0
, pero
ahora considerando el teorema extendido del valor medio, tambin se
puede evaluar de manera exacta mediante:
R1 = f''()h
2
/2!
Y as, sucesivamente:
1.5.2 El residuo de la serie de Taylor
f(x)
x
x
i
X
i+1

f(X
i+1
)
h
P(X) = a
o

P(X) = a
o
+ a
1
x
P(X) = a
o
+ a
1
x + a
2
x
2

f(x)
a
o

1.5.2 El residuo de la serie de Taylor
Un truncamiento a tres trminos (n = 2), supone que la funcin a aproximar
es una parbola P(X) = a
0
+ a
1
X + a
2
X
2
y un posible error dado por el
residuo de segundo orden R
2
.
Un truncamiento a cuatro trminos (n = 3), supone que la funcin a
aproximar es una parbola cbica P(X) = a
0
+ a
1
X + a
2
X
2
+ a
3
X
3
y un
posible error dado por el residuo de segundo orden R
3
.
En general, un truncamiento a (n+1) trminos de la serie, supone un
polinomio P(X) = a
0
+ a
1
X + a
2
X
2
+ a
3
X
3
+ ... + a
n
X
n
y un posible error dado
por el residuo de n-simo orden, que se expresa:
R
n
= f
(n+1)
()h
n+1
/(n+1)! _____ (1.26)
R
n
es el error por truncamiento al aproximar el valor de una funcin f(X
i+1
),
considerando solamente los (n+1) primeros trminos de la expansin en
serie de Taylor correspondiente a la funcin.
1.5.2 El residuo de la serie de Taylor
Ejemplo. Obtener una aproximacin al valor del nmero e, con mantisa de
ocho dgitos y considerando los primeros ocho trminos de la expansin en
serie de Taylor para la funcin f(X) = e
x
.
Sabemos que: e
0
= 1,
entonces: X
i
= 0 ; X
i+1
= 1 ; h = 1 - 0 = 1
f(0) = e
0
= 1 f(1) = e
f(1) = f(0) + f'(0)(1) + f''(0)(1)2/2! + f'''(0)(1)3/3! + fiv(0)(1)4/4! + ...
f'(X) = e
x
f'(0) = 1
f''(X) = e
x
f''(0) = 1
f'''(X) = e
x
f'''(0) = 1
...
f
'(n)
(X) = e
x
f
'(n)
(0) = 1
f(1) ~ 1 + 1 + 1/2! + 1/3! + 1/4! + 1/5! + 1/6! + 1/7!
e ~ 1 + 1 + 0.5 + 0.16666667 + 0.04166667 + 0.00833333 +
+ 0.00138889 + 0.00019841 = 2.71825397
El valor que arroja una calculadora de 9 dgitos es: e = 2.71828183
1.5.2 El residuo de la serie de Taylor
El error por truncamiento es:
R
7
= f
viii
()(1)8/8! = f
viii
()/40320 = 0.00002786
f
viii
() = e

= 1.1233152 ; = 0.11628431
Observamos que efectivamente se localiza entre X
i
y X
i+1
: 0 < < 1,
aunque bastante ms cerca de X
i
que de X
i+1
Si hubisemos truncado a solo tres trminos: e ~ 2.5,
R
2
= f'''()(1)
3
/3! = f'''()/6 = 0.21828183
f'''() = e

= 1.30969098 ; = 0.26979122
Vemos tambin que el valor de es distinto para residuos de diferente
orden, pero siempre cumple con localizarse entre X
i
y X
i+1
.
1.5.2 El residuo de la serie de Taylor
Los valores dados en la siguiente tabla, referidos a este ejemplo,
pueden verificarse fcilmente:
n e Rn f(n+1)()
0 1 1.71828183 1.71828183 0.54132486
1 2 0.71828183 1.43656366 0.36225391
2 2.5 0.21828183 1.30969098 0.26979172
3 2.66666667 0.05161516 1.23876384 0.21411398
4 2.70833334 0.00994849 1.19381880 0.17715724
5 2.71666667 0.00161516 1.16291520 0.15092996
6 2.71805556 0.00022627 1.14040080 0.13137978
7 2.71825397 0.00002786 1.12331520 0.11628431

1.5.2 El residuo de la serie de Taylor
En este ejemplo, tuvimos manera de consultar el valor "exacto" de e, echando mano
de una calculadora, igual que lo pudimos haber consultado en un libro; el nmero e
es conocido por toda la comunidad cientfica, por eso su valor es accesible a
cualquiera.
Pero no ocurre lo mismo cuando estamos estimando el valor de una funcin
complicada, ligada a un experimento en el que apenas tenemos idea de su
comportamiento y del orden de magnitud que puede tomar la funcin, para cada
determinado valor de la variable. En tal caso, no hay manera de calcular con
exactitud los residuos y solo habr que conformarse con una estimacin burda de
ellos.
Para el efecto, y siempre que sea factible derivar analticamente la funcin de
inters, se sugiere considerar como valor estimado de el punto medio entre X
i
y X
i+1
,
es decir:
* = (X
i
+ X
i+1
)/2 _____ (1.27)
con la seguridad de que los residuos estimados a partir de este valor y, por ende, los
errores asociados a ellos, siempre sern superiores a los verdaderos.
R
n
= f
(n+1)
(*)h
n+1
/(n+1)! _____ (1.26')