CÁLCULO

TEMA 4

SERIES DE TAYLOR
Y MacLAURIN

Apuntes preparados por:

Prof. Dr. Ignacio García-Juliá


 Series de Taylor y MacLaurin

1.- Conocimientos previos

Para poder seguir adecuadamente este tema, se requiere que el alumno repase:

• Límite de una función.

• Derivación y derivación n-ésima.

• Progresiones aritméticas y geométricas

• Aritmética de polinomios.

• Funciones exponenciales y logarítmicas.

2.- Introducción

En muchas ocasiones solemos encontrarnos con el problema de calcular el valor de una

función en un punto que puede ser, o bien difícil de calcular, o bien computacionalmente
costoso, o bien no disponemos de las herramientas analíticas o numéricas para hacerlo.

En estos casos, sería de gran ayuda el disponer de una herramienta matemática que fuera
capaz de “sustituir” la función que queremos calcular por un polinomio de fácil estructura
y, por tanto, de cálculo más fácil.

Normalmente, en los problemas del mundo real, no necesitamos conocer todo el

comportamiento de la función en todo su campo de existencia. Nos bastará conocer su
comportamiento en un determinado intervalo que es el que necesitamos para su cálculo o
estudio.

Un ejemplo: para la entrada en la atmósfera terrestre, una nave sigue una trayectoria
compleja, conocida, que tiene que ver con la solución a un conjunto de ecuaciones
diferenciales. Justo en el punto de entrada, punto crítico, cualquier desviación puede hacer
que la nave “rebote” y no entre en la atmósfera, según el ángulo de llegada, o que penetre
muy rápidamente dando lugar a un calentamiento excesivo. Ante cualquier desviación en
la llegada, los ordenadores tienen que ser capaces de corregir instantáneamente la
trayectoria de entrada, y para ello es fundamental que puedan recalcular esta lo más
rápida y precisamente posible. Esto es más fácil y rápido con un polinomio que con la
función original. Pero hay más, también necesitaremos conocer exactamente el error
cometido en dicha sustitución para determinar si es tolerable o no. Esto también será
posible calcularlo con la rapidez y la precisión necesarias para que el ángulo de entrada
sea el correcto.

Tema 4 -2- Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

3.- Aproximaciones polinómicas a funciones elementales

Para encontrar una función polinómica P que aproxime a otra función f, se empieza por
elegir un número c en el dominio de f en el que P y f tengan el mismo valor: P(c) = f(c);
esto es, ambas funciones pasan por el punto (c, f(c)).

Se dice entonces que la aproximación polinómica se expande alrededor de c o está

centrada en c. Geométricamente, el requisito de que P(c) = f(c) significa que la gráfica de P
debe pasar por el punto (c, f(c)). Hay muchos polinomios cuya gráfica pasa por ese punto;
la tarea es encontrar un polinomio cuya gráfica se parezca a la función f en la cercanía de
ese punto. Una manera de hacer esto es imponer el requisito adicional de que la pendiente
de la función polinómica sea la misma que la pendiente de la gráfica de f en el punto (c,
f(c)). Dicho de otro modo:

P '(c) = f '(c)

Ejemplo 1: Aproximación a f(x) = ex mediante un polinomio

Dada la función f(x) = ex, encontrar una función polinómica de primer grado en la forma

P1 (x) = a0 + a1 (x)

cuyo valor y pendiente en x = 0 coincidan con el valor y pendiente de f.

Solución

Como f(x) = ex y f ’(x) = ex , el valor y la pendiente de f en x = 0 están dados por:

!# f (0) = e 0 = 1
"
#$ f '(0) = e 0 = 1

Como

P1 (x) = a0 + a1 (x)

se puede usar la condición P1(0) = f(0) para concluir que a0 = 1. Es más, como P’(x) = a1, se
puede usar la condición P’(0) = f ‘(0), para concluir que a1 = 1 . Por tanto,

P1 (x) = 1+ x

En la siguiente figura se ven ambas gráficas:

Tema 4 -3- Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

y = ex
P1 (x) = 1+ x

Figura 4-1. Aproximación de

una función a una recta.

Vemos que la aproximación es buena, pero que si nos alejamos del punto (0,1), la precisión
de la aproximación disminuye. Para mejorar esta, podemos imponer otro requisito
todavía: que las segundas derivadas de P y f sean iguales en el punto x = 0. Tenemos
entonces que el polinomio que cumple las tres condiciones es:

1
P2 (x) = 1+ x + x 2
2

y = ex

P1 (x) = 1+ x
1
P2 (x) = 1+ x + x 2
2

Figura 4-2. Aproximación de

una función a una recta y a
un polinomio de segundo
grado

Y vemos que este nuevo polinomio se ajusta mejor a f que el polinomio de primer grado. Si
continuamos con este patrón de actuación obtenemos el polinomio de tercer grado:
Tema 4 -4- Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

y = ex

P1 (x) = 1+ x
1
P2 (x) = 1+ x + x 2
2

1 1
P3 (x) = 1+ x + x 2 + x 3
2 6 Figura 4-3. Aproximación de
una función a una recta, a un
polinomio de segundo grado y
a un polinomio de tercer
grado

Si continuamos este proceso, la aproximación polinómica de grado n será:

1 1 1
Pn (x) = 1+ x + x 2 + x 3 +... + x n ≈ e x
2 3! n!

4.- Polinomios de Taylor y McLaurin

Si f tiene n derivadas en c, entonces el polinomio:

f ''(c) f '''(c) f (n) (c)

Pn (x) = f (c) + f '(c)(x − c) + (x − c)2 + (x − c)3 +... + (x − c)n
2! 3! n!

se llama Polinomio de Taylor de grado n para la función f en el punto c.

Si c = 0, entonces

f ''(0) 2 f '''(0) 3 f (n) (0) n

Pn (x) = f (0) + f '(0)x + x + x +... + x
2! 3! n!

se llama Polinomio de McLaurin de grado n para f.

En el ejemplo 1 hemos calculado, sin haberlo nombrado todavía, el polinomio de

MacLaurin para la función f(x) = ex.

Ejemplo 2: Encontrar el polinomio de Taylor para ln x.

Encontrar los polinomios de Taylor P0, P1, P2, P3 y P4 para f(x) = ln x, centrado en c = 1.
Tema 4 -5- Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

Solución

Desarrollando se obtiene para c = 1 los siguientes valores de la función y sus derivadas:

#
% f (x) = ln x ⇒ f (1) = ln1 = 0
%
% f '(x) = 1 ⇒ f '(1) = 1 = 1
% x 1
%% 1
$ f ''(x) = − 2 ⇒ f ''(1) = −1
% x
% 2!
% f '''(x) = x 3 ⇒ f '''(1) = 2
%
% f (4) (x) = − 3! ⇒ f (4) (1) = −6
%& x4

Por consiguiente, los polinomios de Taylor quedan como sigue:

P0 (x) = f (1) = 0
P1 (x) = f (1) + f '(1)(x −1) = (x −1)
f ''(1) 1
P2 (x) = f (1) + f '(1)(x −1) + (x −1)2 = (x −1) − (x −1)2
2! 2
f ''(1) f '''(1) 1 1
P3 (x) = f (1) + f '(1)(x −1) + (x −1)2 + (x −1)3 = (x −1) − (x −1)2 + (x −1)3
2! 3! 2 3
(4)
f ''(1) f '''(1) f (1) 1 1 1
P4 (x) = f (1) + f '(1)(x −1) + (x −1)2 + (x −1)3 + (x −1)4 = (x −1) − (x −1)2 + (x −1)3 − (x −1)4
2! 3! 4! 2 3 4

Ejemplo 3: Encontrar los polinomios de MacLaurin para y = cos x

Encontrar los polinomios de MacLaurin P0, P2, P4 y P6 para f(x) = cos x. Use P6(x) para aproximar
el valor de cos(0.1).

Solución

Desarrollando respecto de c = 0 se obtiene:

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

Tema 4 -6- Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

A través de repetida derivación vemos que el patrón 1, 0, -1, 0 se repite y se obtienen los polinomios
de MacLaurin siguientes:

P0 (x) = 1
1 2
P2 (x) = 1− x
2!
1 1
P4 (x) = 1− x 2 + x 4
2! 4!
1 2 1 4 1 6
P6 (x) = 1− x + x − x
2! 4! 6!

Usando P6(x), se obtiene una aproximación de cos(0.1) ≈ 0.995004165. La gráfica de ambas

funciones, en el entorno de x = 0 es la siguiente:

f (x) = cos x

1 2 1 4 1 6 Figura  4-­‐4.  Aproximación  de  la  

P6 (x) = 1− x + x − x función  y  =  cos  x  con  un  
2! 4! 6!
polinomio  P6(x)

Ejemplo 4: Encontrar un polinomio de Taylor para y = sen x

Encontrar un polinomio de Taylor para f(x) = sen x, desarrollado en el entorno de c = ∏/6

Solución

Desarrollando respecto de c = ∏/6, se obtiene los siguiente:

Tema 4 -7- Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

"π % "π % 1
f (x) = sen x ⇒ f $ ' = sen $ ' =
#6& #6& 2

"π % "π % 3
f '(x) = cos x ⇒ f ' $ ' = cos $ ' =
#6& #6& 2
"π % "π % 1
f ''(x) = −sen x ⇒ f '' $ ' = −sen $ ' = −
#6& #6& 2
"π % "π % 3
f '''(x) = −cos x ⇒ f ''' $ ' = −cos $ ' =-
#6& #6& 2

Así, el tercer polinomio de Taylor para la función y = sen x desarrollado en el entorno del
punto c = ∏/6, es

!π $ !π $
f '' # & 2 f ''' # & 3 2 3
!π $ ! π $! π$ "6% ! π $ " 6 %! π$ 1 3! π$ 1 ! π$ 3 ! π$
P3 (x) = f # & + f ' # &# x − & + # x − & + # x − & = + # x − & − # x − & − # x − &
"6% " 6 %" 6% 2! " 6% 3! " 6% 2 2 " 6 % 2(2!) " 6 % 2(3!) " 6%

Y su representación gráfica con respecto a y = sen x es la siguiente:

f (x) = sen x

1 3" π % 1 " π%
2
3 " π%
3
Figura 4-5. Aproximación de la
P3 (x) = + $x − '− $x − ' − $x − '
2 2 # 6 & 2(2!) # 6 & 2(3!) # 6& función y = sen x con un
polinomio P3(x)

5.- Resto de un polinomio de Taylor

Una técnica de aproximación es de poco valor sin alguna idea de su precisión. Para medir
la precisión de una aproximación al valor de una función f(x) mediante un polinomio de
Taylor Pn(x), se puede usar el concepto de resto Rn(x), definido como sigue:
Tema 4 -8- Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

f (x) = Pn (x) + Rn (x)

Así, al valor absoluto de Rn(x) se le llama error de la aproximación. Es decir

Rn (x) = f (x) − Pn (x)

El siguiente teorema da un procedimiento general para estimar el resto de un polinomio

de Taylor. Este teorema es conocido como Teorema de Taylor, y el resto dado en el
teorema se denomina Resto de Lagrange:

Teorema de Taylor: si una función f es derivable hasta el orden n +1 en un intervalo I que

contiene a c, entonces, para toda x en I, existe z entre x y c tal que

f ''(c) 2 f (n) (c)

f (x) = f (c) + f '(c)(x − c) + (x − c) +... + (x − c)n + Rn (x)
2! n!

donde

f (n+1) (z)
Rn (x) = (x − c)n+1
(n +1)!

Veremos esto con unos ejemplos:

Ejemplo 5: Determinar la precisión de una aproximación

El polinomio de MacLaurin de tercer grado para la función y = sen x está dado por:

x3
P3 (x) = x −
3!

Usar el teorema de Taylor para aproximar sen (0.1) mediante P3(0.1) y determinar la precisión de la
aproximación.

Solución

Aplicando el teorema de Taylor se tiene que,

x3 x 3 f (4) (z) 4
sen x = x − + R3 (x) = x − + x
3! 3! (4)!

Tema 4 -9- Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

donde 0 < z < 0.1, es decir, el error es menor que el primer término suprimido. Por
consiguiente:

(0.1)3
sen (0.1) ≈ 0.1− ≈ 0.1− 0.000167 = 0.099833
3!

Como

f (4) (z) = sen z

se sigue que el error

R3 (0.1)

puede acotarse como sigue:1

sen z 0.0001
0 < R3 (0.1) = (0.1)4 < ≈ 0.000004
4! 4!

Esto implica que

0.099833 < sen (0.1) = 0.099833 + R3 (x) < 0.099833+ 0.000004

0.09983 < sen (0.1) < 0.099837

6.- Cálculo de los polinomios de Taylor con MATLAB

MATLAB nos proporciona las funciones necesarias para trabajar con polinomios de Taylor.
Recordemos que la función diff nos proporciona las derivadas de una función. En el
siguiente ejemplo calculamos las cinco primeras derivadas de la función dada:

>> f = 1/x;
>> [diff(f), diff(f,2), diff(f,3), diff(f,4), diff(f,5)]
ans = [ -1/x^2, 2/x^3, -6/x^4, 24/x^5, -120/x^6]

Recordemos que podemos arreglar la solución para que quede más legible con pretty:
>> pretty(ans)

1 Tener en cuenta que sen z es menor que la unidad, por tanto, si se elimina del producto la cantidad será
siempre mayor.
Tema 4 - 10 - Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

+- -+
| 1 2 6 24 120 |
| - --, --, - --, --, - --- |
| 2 3 4 5 6 |
| x x x x x |
+- -+

MATLAB nos proporciona la función taylor para calcular directamente los polinomios
de Taylor de una función dada:

A continuación se realiza el desarrollo del polinomio de Taylor de grado 5 para la función

f(x) = 1/(2-x) en el entorno del punto x =1:

>> f = 1/(2-x);
>> pretty(taylor(f, x, 1))
2 3 4 5
x + (x - 1) + (x - 1) + (x - 1) + (x - 1)

Si queremos el polinomio de MacLaurin, basta con no especificar el punto de entorno, y

MATLAB tomará por defecto el punto x = 0:

>> pretty(taylor(f, x))

5 4 3 2
x x x x x
-- + -- + -- + -- + - + 1/2
64 32 16 8 4

Por defecto, MATLAB calcula el polinomio de Taylor (o MacLaurin) de orden 5. Si

queremos un orden superior, debemos especificarlo:

>> pretty(taylor(f, x, 1, 'Order', 8))

2 3 4 5 6 7
x + (x - 1) + (x - 1) + (x - 1) + (x - 1) + (x - 1) + (x - 1)

Otra función de MATLAB que utilizaremos con frecuencia es subs cuya misión es, como
hemos visto en temas anteriores, sustituir una variable por un valor cualquier expresión.
Por ejemplo:

>> f = sin(x);
>> r = taylor(f, x, pi/6, ‘Order’, 4)

Tema 4 - 11 - Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

r =
(3^(1/2)*(pi/6 - x)^3)/12 - (pi/6 - x)^2/4 - (3^(1/2)*(pi/6 - x))/2 +
+ 1/2

>> p = subs(r, x, 0.1)

p =
(3^(1/2)*(pi/6 - 1/10)^3)/12 - (pi/6 - 1/10)^2/4 - (3^(1/2)*…
(pi/6 - 1/10))/2 + 1/2
>> eval(p)
ans = 0.099264649295008
Aquí debemos detenernos un momento y hacer una serie de comentarios. Primero, cuando
llamamos a la función taylor, establecemos el parámetro ‘Order’ que nos indica el
número de términos que queremos obtener. También, observemos la forma de invocar a la
función subs: en la expresión r, polinomio de Taylor que acabamos de obtener,
cambiamos la variable x por el valor 0.1. Como estamos trabajando con variables
simbólicas, para obtener el valor real de la variable p deberemos “evaluar” la expresión
obtenida. Evaluar es equivalente a escribir la expresión en la línea de comando y pulsar
[Intro]. El resultado numérico obtenido tiene 15 decimales. MATLAB, por defecto,
proporciona el resultado con 4 decimales (format short). Para tener más precisión,
como en este caso, deberemos escribir previamente format long.

Ahora ya estamos en condiciones de repetir el ejemplo 5 usando MATLAB.

Ejemplo 5(bis): Determinar la precisión de una aproximación (MATLAB)

El polinomio de MacLaurin de tercer grado para la función y = sen x está dado por:

x3
P3 (x) = x −
3!

Usando MATLAB, utilizar el teorema de Taylor para aproximar sen (0.1) mediante P3(0.1) y
determinar la precisión de la aproximación.

Solución

% Cálculo de sen(0.1) con polinomio de Taylor de grado 3 y su

% aproximación con el Resto de Lagrange

% Se declaran las variables simbólicas, el punto de evaluación y la

% función

Tema 4 - 12 - Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

syms x f r
punto = 0.1
f = sin(x);

% Cálculo del polinomio de Taylor de cuatro términos

r = taylor(f, x, 'Order', 4)

% Hacemos la sustitución en el polinomio de x por 0.1 y obtenemos su

% valor numérico

p = subs(r, x, punto)
u = eval(p)

% Calculamos el grado del Resto de Lagrange

dim = sym2poly(r);
grado = size(dim);

% Calculamos el valor del resto de Lagrange

R = (punto)^grado(2)/factorial(grado(2))

% El máximo error cometido al quedarnos con el polinomio de grado 3 será

% el valor del polinomio más el resto de Lagrange

max_error = u + R

% Presentamos los resultados obtenidos

fprintf('El sen(0.1) está comprendido entre %1.6f y %1.6f \n',...

u, max_error)

% ---------------------------- FIN ------------------------------------

También podemos utilizar el polinomio de Taylor para aproximar el valor de una función
con una precisión determinada. Esto es verdaderamente útil cuando se debe evaluar la
decisión de sustituir la función por el polinomio atendiendo a su carga de computación:

Ejemplo 6: Aproximar un valor con una precisión determinada

Tema 4 - 13 - Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

Determine el grado del polinomio de Taylor Pn(x) desarrollado respecto de c = 1 que debe usarse
para aproximar ln(1.2) de manera que el error sea menor de 0.001.

Solución

Podemos empezar por hallar la forma que tiene la derivada enésima de la función
propuesta. Para ello, veremos la estructura de las cinco primeras derivadas:

>> f = log(x);
>> diff(f,x)
ans = 1/x
>> diff(f,x,2)
ans = -1/x^2
>> diff(f,x,3)
ans = 2/x^3
>> diff(f,x,4)
ans = -6/x^4
>> diff(f,x,5)
ans = 24/x^5

Vemos que el denominador tiene la variable x elevada a la misma potencia que la

derivada, por tanto podemos aventurar que el término general será xn. Por otro lado, el
signo de la derivada va alternando, positivo para los valores impares de n y negativo para
los valores pares de n. Esto puede expresarse como (-1)n-1. En cuanto al numerador,
observamos que es el desarrollo de un factorial de grado inferior al orden de la derivada,
es decir, para una derivada de orden n, el numerador sería (n-1)!. Por tanto:

(n −1)!
f (n) (x) = (−1)n−1
xn

Para nuestro propósito, y atendiendo al Resto de Lagrange, es conveniente trabajar con el

término n+1:

n!
f (n+1) (x) = (−1)n
x n+1

Usando el teorema de Taylor sabemos que:

Tema 4 - 14 - Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

f (n+1) (z) n+1 n! " 1 % n+1 (0.2)n+1

Rn (1.2) = (1.2 −1) = n+1 $ '(0.2) = n+1
(n +1)! z # (n +1)!& z (n +1)

donde 1 < z < 1.2. En este intervalo, (0.2)n+1/zn+1(n+1) es menor que (0.2)n+1/(n+1), ya que en
el denominador tenemos zn+1 que siempre será mayor que 1, y este término, a su vez,
deberá ser menor que la precisión requerida. En otras palabras:

(0.2)n+1
< 0.001 ⇒ 1000 < (n +1)5n+1
n +1

Para resolver esta ecuación podemos proceder de dos formas: por prueba y error2, dando
valores a n hasta ver qué valor satisface la ecuación, o utilizando MATLAB:

>> solve((x+1)*5^(x+1)-1000)
ans = lambertw(0, 1000*log(5))/log(5) - 1

En este caso, MATLAB no nos da directamente la solución, pero nos dice que esta es la
solución a una ecuación de Lambert dividida por el logaritmo neperiano de 5 menos 1. No
tenemos que saber lo que es una ecuación Lambert3, bastará con que MATLAB la resuelva:

>> lambertw(0, 1000*log(5))

ans = 5.651686311474081

Si dividimos por ln 5 y restamos 1, obtenemos el resultado de la ecuación, es decir, el valor

de n : 2.5115, que nos llevaría a un polinomio de grado 3.

7.- Problemas propuestos

Ejercicio 1

Considerar la función f(x) = x2ex.

(a) Encontrar los polinomios de MacLaurin de orden 2, 3 y 4 para f.

(b) Representar gráficamente f y los polinomios calculados en el apartado anterior.

2 (2+1)5(2+1) = 375; (3+1)5(3+1) = 2500. Por tanto, el grado que satisface la ecuación es 3.
3 Una ecuación de Lambert es aquella que tiene la forma x = wew donde w es a su vez una función de x.
Tema 4 - 15 - Prof. Dr. Ignacio García Juliá
 Series de Taylor y MacLaurin

(c) Evaluar y comparar los valores de f(n)(0) y Pn(n)(0) para n = 2, 3 y 4.

(d) Usar los resultados del apartado anterior para hacer una conjetura sobre f(n)(0) y Pn(n)(0)

Ejercicio 2

En los apartados siguientes, encontrar los polinomios de Taylor para el grado que se indica
y centrado en c:

1
(a) f (x) = , n = 4, c = 1
x

2
(b) f (x) = , n = 4, c = 2
x2

(c) f (x) = x, n = 4, c = 1

(d) f (x) = 3 x , n = 4, c = 8

(e) f (x) = ln x, n = 4, c = 1

( f ) f (x) = x 2 cos x, n = 2, c = π

(g) f (x) = cosh x, n = 2, c = 2

(h) f (x) = senh x, n = 3, c = 1

(i) f (x) = tanh x, n = 4, c = 1

Ejercicio 3

En los apartados siguientes, encontrar los polinomios de MacLaurin para el grado que se
indica:

(a) f (x) = e− x , n = 3

(b) f (x) = e− x , n = 5

(c) f (x) = e 2 x , n = 4

Tema 4 - 16 - Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

(d) f (x) = e3x , n = 4

(e) f (x) = sen x, n = 5

( f ) f (x) = sen π x, n = 3

(g) f (x) = xe x , n = 4

(h) f (x) = x 2 e− x , n = 4

1
(i) f (x) = ,n = 4
1+ x

x
( j) f (x) = ,n = 4
1+ x

(k) f (x) = sec x, n = 2

(l) f (x) = tan x, n = 3

Ejercicio 4

En los apartados siguientes, determinar el grado del polinomio de MacLaurin requerido

para que el error en la aproximación de la función en el valor indicado de x sea menor de
una milésima.

(a) sen (0.3)

(b) e 0.6

(c) cos(0.1)

(d) e 0.3

Ejercicio 5 (avanzado)

Considérese la siguiente función:

(x − 3)2 + (y −1)2 = 4

Tema 4 - 17 - Prof. Dr. Ignacio García Juliá

 Series de Taylor y MacLaurin

(a) Encontrar el polinomio de Taylor centrado en el punto c = 3 de grado 2

(b) Analizar los resultados.

o0o

Tema 4 - 18 - Prof. Dr. Ignacio García Juliá