5.

Integración en R
5.1. Definición y propiedades
Sea f acotada en [a, b] . Dividimos [a, b] en n subintervalos de
la misma longitud ∆x por medio de los n+1 puntos:
M4
a = x0 < x1 < · · · < xn = b con xk+1 −xk = b−a
n ≡ ∆x .
m4
Para cada n llamamos suma inferior Ln y superior Un a:
a=x0 x 1 D x x2 x3 b=x 4
n n m = ínf{ f (x) : x ∈ [xk−1 , xk ]}
Ln = ∑ mk ∆x , Un = ∑ Mk ∆x , con k
k=1 k=1 Mk = sup{ f (x) : x ∈ [xk−1 , xk ]}

Si ambas sucesiones {Ln } y {Un } convergen hacia un

mismo límite, decimos que f es integrable en [a, b] ,
Def.
representamos ese límite común por ab f ó ab f (x)dx
R R

y le llamamos integral de f en [a, b] .

[Esta no es la definición de ‘integral de Riemann’ habitual (ver Spivak), pero es mucho más corta].
El significado geométrico es claro: si f ≥0 , la integral (≥0) representa el área A de la región
limitada por la gráfica de f , el eje x y las rectas x=a y x=b : A es para todo n mayor que
la suma Ln de las áreas de los rectángulos pequeños y menor que la suma Un de los grandes;
al crecer n , ambas sumas tienden hacia A . Si f ≤ 0 , Ln y Un son negativas. La integral (≤ 0)
en valor absoluto es el área de la región (situada bajo el eje x ) limitada por el eje x , la gráfica
de f y las rectas x = a y x = b . Si f es positiva y negativa, la integral
Rb
+
a f será la diferencia entre las áreas de las regiones que queden por +
encima y las áreas de las que queden por debajo del eje x : a – b

Con los teoremas que probaremos, para saber si f es integrable y calcular la integral no se
necesitará usar la definición casi nunca. Por ahora, sólo con lo visto, unos ejemplos:
n 2 1
Ln = ∑ (k−1) 1 1
 2
0 + · · · + (n−1)2

n2 n = n3
Ej. f (x) = x2 , x ∈ [0, 1] . k=1
n 2
Un = ∑ nk2 1n = n13 12 + · · · + n2
 
k=1
Usando el resultado que vimos en un problema de sucesiones: 0 1
n[n+1][2n+1] [n−1]n[2n−1] n[n+1][2n+1] 1 R1
12 + · · · + n2 = 6 ; Ln = 6n3
, Un = 6n3
; Ln , Un → 3 = 0 f.

−1 si x = a
(
Rb 1
Ej. g(x) = 0 si a < x < b . Ln = − b−a b−a
n , Un = n ⇒ a g = 0 .
1 si x = b a
Rb
g es discontinua, pero integrable. Seguiría siendo a g = 0 si cambiamos b
–1 __
b–a
n
el valor 0 por cualquier otro en un número finito de puntos de (a, b).
Veremos pronto que funciones acotadas con un número finito de discontinuidades son siempre integra-
bles, así que las funciones no integrables tienen que ser tan patológicas como la siguiente.

1, x∈Q 1
Ej. h(x) = , x ∈ [a, b]. En cada [xk−1 , xk ] hay puntos de Q y de
0 , x ∈ R−Q
n n
R−Q ⇒ Ln = ∑ 0 b−a b−a
n = 0 , Un = ∑ 1 n = b − a ∀n ⇒ h no integrable.
k=1 k=1 a b

91
Las siguientes propiedades son intuitivamente claras a la vista del significado geométrico de la
integral y se demuestran mecánicamente usando de las definiciones (las dos primeras se resumen
diciendo que ‘la integral es lineal’, como la derivada):
Teorema:
Sean f y g integrables en [a, b] . Entonces:
M
Rb Rb Rb Rb Rb
a cf = c a f , c∈R ; a [ f + g] = a f+ a g.
m Rb
Si m ≤ f ≤ M en [a, b] ⇒ m(b − a) ≤ a f ≤ M(b − a) . g
R  R +
 b
f  ≤ ab | f | . Si f ≤ g en [a, b] ⇒ ab f ≤ ab g .
R R
f

 a
Ra Ra Ra
Si f es impar −a f =0. Si f es par −a f =2 0 f.

La siguiente sigue siendo intuitiva, pero es pesada de demostrar con nues-

tra definición. [Para f continua será trivial usando los teoremas de 5.2,
pero la propiedad es cierta también para f integrable y discontinua]. a c b
Teorema:
Rb Rc Rb
Si a<c<b ; f integrable en [a, b] ⇒ f integrable en [a, c] y [c, b] , e a f= a f+ c f.
Ra Rb Ra
Definiendo a f =0 e a f =− b f , la igualdad es válida para a, b, c cualesquiera.

Los dos siguientes teoremas dicen que las f (acotadas) no integrables son extrañas:
Teorema: f continua en [a, b] ⇒ f integrable en [a, b] .
Idea de la demostración. Se prueba que {Ln } y {Un } siempre convergen (sus límites se llaman
integral inferior y superior de f ) y así f es integrable si Un − Ln → 0 . Como f es uniformemente
continua, la diferencia entre dos valores cualesquiera de f en dos x, y∈[a, b] es tan pequeña como
queramos para x e y suficientemente próximos; en particular, lo es la diferencia entre los valores
máximo y mínimo de f en [xk−1 , xk ] , si el intervalo es muy pequeño. Así, si n es suficientemente
grande tenemos que para cada k :
n
ε
Mk −mk < b−a ∀ε ⇒ Un −Ln = ∑ (Mk −mk )∆x < ε para n grande ⇒ Un −Ln → 0 .
k=1
f continua a trozos
Ya vimos que funciones discontinuas podían ser integrables. Una (integrable)
f se dice continua a trozos en [a, b] si es continua salvo en un
número finito de puntos y en ellos posee límites laterales.
1
o sen ( 1x ) en [0, 1] , definámoslas como las definamos en x = 0 .
 
No lo son las funciones x

Teorema: f continua a trozos en [a, b] ⇒ f integrable en [a, b] .

[dividimos en subintervalos de modo que f sólo tenga
= + discontinuidades en los extremos; en cada intervalo es
fácil ver que es integrable por ser suma de una función
continua y de otra integrable como la g del ejemplo].
Como es 0 el valor de la integral de una función como g se ve que cambiando el valor de una
f integrable en un número finito de puntos, la nueva función h continúa siendo integrable
y el valor de la integral es el mismo, pues h = f +g con una g de esas y la integral es lineal.

Hay funciones integrables con infinitas discontinuidades; por ejemplo, una f creciente y acotada
es integrable pues Mk−1 =mk , k=1, ..., n ⇒ Un −Ln = f (b)−n f (a) , aunque tenga infinitos saltos .



92
5.2. Teoremas fundamentales
Estos teoremas fundamentales del cálculo infinitesimal relacionan las derivadas y las integra-
les y nos permitirán hallar muchísimas integrales prescindiendo de la definición.
Rx
Sea f acotada e integrable en [a, b] ; para cada x ∈ [a, b] la a f existe (y es un número).
'área'=F(x)
Rx
Podemos, pues, definir una nueva función: F(x) = a f , x ∈ [a, b]
b
a x f
Primer teorema fundamental del cálculo infinitesimal:
f integrable en [a, b] ⇒ F continua en [a, b] .
Si además f es continua en un c∈(a, b) entonces F es derivable en c y F 0 (c)= f (c) .
Y por tanto, si f es continua en todo [a, b] entonces F 0 (x) = f (x) ∀x ∈ [a, b] .

Como f es acotada en [a, b] , existen supremo e ínfimo de f en cada intervalo ⊂ [a, b] .

Sea c ∈ (a, b) y sea h > 0 . Llamemos: c+h
c f
 
Mh = sup f (x) : x ∈ [c, c+h] , mh = ínf f (x) : x ∈ [c, c+h] .
R c+h R c+h Rc
Entonces: mh h ≤ c f= a f− a f = F(c+h)−F(c) ≤ Mh h
⇒ [F(c+h) − F(c)] → 0 , si h → 0+ . a c c+h b

Así pues, F es continua en c (cambiando detalles se vería para h → 0− ).

Sea ahora f continua en c . Si h > 0 se deduce que: mh ≤ F(c+h)−F(c)
h ≤ Mh
(y a la misma igualdad se llegaría, de forma análoga, para h < 0 ).
F(c+h)−F(c)
Como f continua en c , Mh , mh → f (c) si h → 0 , y por tanto h −→ f (c) .
h→0
x x [El teorema nos dice que, al contrario que al derivarla,
!0 f !0 F
f F la F obtenida integrando f es ‘más suave’ que ella. Si
f es discontinua en c , F es continua en c (aunque F
d/dx d/dx tenga un ‘pico’ en c ); y si f tiene picos, desaparecen
al integrarla. En general, f ∈Cn ⇒ F ∈Cn+1 ].

Ej. f (x) = esen (ch x) continua ∀x ⇒ F(x) = 0x f , G(x) = 7x f , . . . tienen por derivada f (x) ∀x .
R R
 0
F (x)= f (x) también para x<a (si f continua en x ), pues si c<x con f integrable en [c, a] :
F(x) = ax f = cx f − ca f .
R R R 

Segundo teorema fundamental del cálculo infinitesimal:

Rb b
f es continua en [a, b] y f = g0 para alguna función g ⇒ a f = g(b)−g(a) ≡ g a

Como F 0 = f = g0 ⇒ F(x) = g(x)+k para algún número k . Como 0 = F(a) = g(a)+k

⇒ k = −g(a) ⇒ F(x) = g(x)−g(a) . En particular: F(b) = ab f = g(b)−g(a) .
R

Def. Dada una función f , una g cuya derivada sea f se llama primitiva de f .
El segundo teorema dice que para calcular la integral de una f continua basta hallar una
primitiva de f (y no es necesario utilizar las sumas superiores e inferiores). Si g es primitiva
de f , es claro que cualquier otra primitiva de f es de la forma g + K .

93
 R
Def. El conjunto de todas las primitivas se designa por f (x) dx .
Rb
(Son funciones y no un número como a f ; a veces se llama integral definida
de f entre a y b a esta última, e integral indefinida al conjunto de primitivas).
En algunos casos, hallar la primitiva de una f es inmediato y, por tanto, lo es calcular algunas
integrales. Por ejemplo, es ahora ya trivial calcular la primera integral vista en 5.1:
R1 2 x3 1 1 x3 2 d x3 2

Ej. 0 x dx = 3 0 = 3 pues 3 es una primitiva de x ya que dx 3 = x ;
2
R 2 1 3
(todas las primitivas de x son x dx = 3 x +K ; si para el cálculo de esta integral
tomásemos otro valor de la K 6= 0 , llegaríamos, desde luego, al mismo resultado).

De hecho, cada derivada conocida nos proporciona una fórmula de integración:

Z Z Z
dx dx 1 2x dx
= log |x2 −1| , . . .
R
Ej. cos2 x
= tan x , sh x = ch x , [4−x]2
= 4−x , x2 −1

(más exacto sería escribir tan x+K , ch x+K , ... ; nosotros no lo haremos pero tengámoslo en cuenta).
A menudo al integrando le faltarán constantes que se calculan derivando de cabeza o tanteando
(como con la x2 de arriba: derivando x3 se tiene 3x2 , luego falta 31 en la primitiva):
Z Z Z
dx dx
dx 1
R −5/3 3 −2/3 1 x
Ej. √ = 3 arc sen (3x) , x dx = − 2 x , = 4[1+(x/2)2 ] = 2 arctan 2 , . . .
1−9x2 4+x2

Pero en muchísimas ocasiones calcular primitivas puede ser largo o muy complicado (a ello
nos dedicaremos en la próxima sección). Más aún, hay funciones de apariencia sencilla para
las que se demuestra que no tienen primitivas que puedan escribirse como suma, producto,
composición,... de funciones elementales, como:
R 2 R x R dx R √ R√
senx2 dx , ex dx , senx x dx , ex dx , log 1+x3 dx , 3 1+x2 dx , . . .
R R
x,

Si f es continua una primitiva suya es la F(x) de los teoremas fundamentales (pero esto
no sirve para calcular una integral concreta). Así F(x) = 0x sent 2 dt , F ∗ (x) = −1
Rx
sent 2 dt , ...
R

son todas primitivas de f (x) = sen x2 ; es decir, sen x2 dx = ax sent 2 dt + K .

R R

[Las variables x, t, ... son mudas, pero no se repite la letra del límite de integración en el
integrando porque podría dar lugar a errores: F(1) es 01 sent 2 dt , pero no es 01 sen 1 dx
R R

y a esto nos podría llevar una incorrecta notación como 0x sen x2 dx ].

R

También hay funciones integrables sin primitivas (claramente no pueden ser continuas):

1 si x = 0 1 Rx

Ej. f (x) = no tiene primitiva F(x) = a f = 0 ∀x no lo es .
0 si x 6= 0 0

R x dt
De los TFCI se deducen las propiedades del log x = 1 t que habíamos adelantado:

f (x) = 1x continua si x > 0 ⇒ F(x) = log x derivable (y continua) si x > 0 y F 0 (x) = 1x .

De la definición también salrían el resto de sus propiedades. Probemos una de ellas

(aunque utilizando un cambio de variable de los que veremos en la próxima sección):
R ab dt R a dt R ab dt
para a, b > 0 , log(ab) = 1 t = 1 t + a t = log a + log b
↑
haciendo t = as en la segunda integral

94
El segundo TFCI permite también probar con facilidad algunas de las propiedades generales de
las integrales vistas en 5.1, en el caso particular de que el integrando sea continuo; por ejemplo,
si F y G son primitivas de f y g se tiene:
Rb Rb Rb
a [ f + g] = [F +G](b) − [F +G](a) = F(b)−F(a) + G(b)−G(a) = a f+ a g,
Rb Rc Rb
a f = F(b) − F(a) = F(b)−F(c) + F(c)−F(a) = a f+ c f , ...
Pero recordemos que también son ciertas estas propiedades para las funciones que sólo son
continuas a trozos. De hecho, sabemos hallar ya fácilmente integrales de muchas f de ese tipo,
dividiendo el intervalo y aplicando los TFCI en cada subintervalo:
n cos x , 0 ≤ x ≤ π
1 cos x
Z π
Ej. Hallemos f , si f (x) = 2 .
π
0 −1 , 2 <x≤π π/2 π
Rπ R π/2 Rπ π/2 π 0
0 f = 0 cos x dx + π/2 [−1] dx = 1− .
= [sen x]0 + [−x]ππ/2 2
 Rπ Rπ  –1
pues π/2 f = π/2 [−1] dx , ya que coinciden salvo en x = π2 .
Rx 1 F
También sabemos hallar para todo x ∈ [0, π] la primitiva F(x) = 0 f.
0 π
Para x > π2 hay dos expresiones de la función: π/2
(R x π
1-π/2
0 cost dt , 0 ≤ x ≤ 2 sen x , 0 ≤ x ≤ π2
Z x 
f = R π/2 Rx =
0
0 cost dt + π/2 [−1] dt , π2 ≤ x ≤ π 1+ π2 −x , π2 ≤ x ≤ π
[función que, como nos aseguraba el primer TFCI, es continua también en x = π/2 ,
aunque no es derivable en ese punto: F 0 (π/2− ) = 0 6= −1 = F 0 (π/2+ ) ].

Como sabemos hallar derivadas de funciones definidas por integrales, sabemos hacer con ellas
todo lo que hemos visto en cálculo diferencial: rectas tangentes, crecimiento y decrecimiento,
extremos, límites indeterminados,...
Z x
t3
Ej. Hallemos la ecuación de la recta tangente a la gráfica de la función F(x)= dt en x = 1 :
−1 t 4 −4
3 R1 t3
F 0 (x) = x4x−4 , F 0 (1) = 31 ; F(1) = = 0 (integrando impar) ⇒ tangente: y = − x−1
−1 t 4 −4 dt 3 .
1
Aunque no se necesita, podríamos (primitiva inmediata) hallar la F(x) = 4 log |x4 −4| − log 3 .
 

Rx
Ej. Sea F(x) = 1 arctan(et ) dt . Estudiar dónde F es inyectiva. Precisar si F(0) > 0 ó F(0) < 0 .
TFC
Integrando continuo ∀x ⇒ F 0 (x) = arctan(ex ) > 0 ∀x , pues ex > 0 ∀x
⇒ F es estrictamente creciente en R ⇒ F es inyectiva en todo R .
R1
Será F(0) < 0 porque F es estrictamente creciente y F(1) = 1 =0 ,
R0 R1 R1
o porque F(0) = 1 =− 0 y es 0 > 0 (integrando positivo en todo R).

| cost 3 |
Z x
Ej. Determinemos, si existe, el límite de G(x) = 1x t 2 +1
dt cuando x → 0 y cuando x → ∞ .
0
Rx R0
El numerador F = 0 es continuo y derivable ∀x (integrando continuo) y es F(0) = 0 = 0.
Cuando x → 0 tenemos indeterminación 0/0 . Por L’Hôpital,
F 0 (x) | cos x3 |
lı́m G(x) = lı́m 1 = lı́m x2 +1 = 1 .
x→0 x→0 x→0

Si x → ∞ , tal vez no valga L’Hôpital (¿tenderá F a ∞?). De hecho, no hay indeterminación,

pues vamos a ver (aunque la primitiva sea no calculable) que F está acotada. En efecto:
| cos x3 | 1 R x dt
0≤ x2 +1
≤ x2 +1
∀x ⇒ 0 ≤ F(x) ≤ 0 t 2 +1 = arctan x ∀x ⇒
arctan x
0 ≤ G(x) ≤ x ⇒ 0 ≤ lı́m G ≤ lı́m arctan
x
x
= 0 ⇒ G −→ 0 .
x→∞ x→∞ x→0

95
En ocasiones se trabaja con funciones similares a la F(x) , definidas por integrales de funciones
f continuas, pero con límites de integración que también son funciones (derivables) de x .
Los TFCI también nos permiten derivarlas:
R b(x)
Si H(x) = a(x) f entonces H 0 (x) = f [b(x)] b0 (x) − f [a(x)] a0 (x) .


Para los x tales que f sea continua en [a(x), b(x)] , o en [b(x), a(x)] si a(x) > b(x) .
R b(x) R a(x)
H(x) = 0 f − 0 f = F[b(x)]−F[a(x)] , con F(x) = 0x f , y regla de la cadena .
 R 

R 3x 2
Ej. Utilicemos la fórmula anterior para hallar dos derivadas de la función: H(x) = x 2x e−t dt .
R 3x −t 2 2 2 2 2 2 2
H 0 (x) = 2x e dt + x e−9x · 3− e−4x · 2 → H 00 (x) = 2 3e−9x −2e−4x +2x2 8e−4x − 27e−9x
   

[expresiones válidas ∀x , tanto si es positivo como negativo].

R 1+x
Ej. Estudiemos el crecimiento y decrecimiento de L(x) = 1−x logt dt en el intervalo I = 0, 21 .
 

Como log x es continua en 21 , 32 valores donde integramos si x ∈ 0, 12 podemos derivar L :

   

L0 (x) = log(1+x)− log(1−x) · [−1] = log(1−x2 ) < 0 si x ∈ I ⇒ decrece.

1
Era esperable: las áreas negativas que aparecen son mayores que las positivas.
En este caso la primitiva sí sería calculable (por partes, como veremos), pero
log x
es un rodeo tonto hallar primitivas para derivarlas a continuación.
Z 1
4 Determinemos el x ∈ [1, 3] que hace máximo el valor de E
Ej. Sea E(x) = e−t dt .
5−2x y probemos que este valor máximo es mayor que 32 .
4
E 0 (x) = 0 + 2 e−(5−2x) > 0 ⇒ E es estrictamente creciente en todo R, con lo que el valor
R 1 −t 4
máximo se toma en x = 3 y será E(3) = −1 e dt (la primitiva no es calculable).
4 −t 4 dt = 2e > 23 .
R1 R1
t ∈ [−1, 1] ⇒ t 4 ≤ 1 ⇒ e−t ≥ e−1 ⇒ −1 e ≥ −1
−1 e dt

Z x
dt
Ej. Hallemos los x que hacen máximo y mínimo el valor de R(x) = en el intervalo [0,4] :
x/2 6+t 3
Integrando continuo y 2x , x derivables ⇒ R derivable si t ≥ 0 .
1/2 3(8−x3 ) f
R0 (x) = 1
6+x3
− 6+(x/2)3 = (6+x3 )(48+x3 )
= 0 si x = 2 . R
R2
Crece en [0, 2] y decrece en [2, 4] ⇒ máximo R(2) = 1 f.
El mínimo, que estará en uno de los extremos, será R(0) = 0 ,
pues R(4) = 24 > 0 , ya que el integrando es positivo.
R

R √x
Ej. Estudiemos en qué x ∈ [0, 2π] alcanza sus valores extremos la función: S(x) = 0 sent 2 dt .
√
1 Su derivada es S0 (x) = sen [( x)2 ] 2√
1
x
− 0 = sen
√ x , ∀x > 0 .
2 x
√π ≈ 0.8
Máximo y mínimo de S existen por ser continua. Los candidatos son
los extremos y los puntos en que S0 = 0 , es decir, x = 0 , x = π y
sen(x2) x = 2π . Con el signo de S0 se ve que S crece antes de π y decrece
√2π ≈ 2.5
–1 después, luego en ese punto se alcanza el máximo.
[No era necesario hallar S0 para decirlo: estaba claro viendo la gráfica de f (x) = sen x2 , pues
hasta x2 = π añadimos áreas positivas y a partir de entonces quitamos áreas bajo del eje x ].
Precisar cuál de los mínimos locales es el absoluto exige saber cuál de estos números es menor:
R √2π
S(0) = 0 ó S(2π) = 0 sen (t 2 ) dt
La gráfica de f sugiere que S(2π) > 0 , pero no podemos dar el valor exacto, por no tener f
primitiva elemental (veremos cómo aproximar numéricamente las integrales en la sección 5.5).

96
5.3. Cálculo de primitivas
Ya vimos en la anterior sección cómo calcular primitivas inmediatas, consecuencias directas de
las fórmulas de derivación (o casi inmediatas, teniendo cuidado con las constantes que pudieran
faltar). Nos dedicamos ahora a ver cómo hallar primitivas algo más complicadas. No existen
muchas más técnicas que las que veremos en esta sección. Que nos quede claro que la mayoría
de las primitivas no son calculables.
De la linealidad de la derivada se deduce inmediatamente para las primitivas que:
R R R R R
[ f (x) + g(x)]dx = f (x)dx + g(x)dx , c f (x)dx = c f (x)dx .
R √ x
 R√ R R x 8 3/2 7 x
Ej. 4 x + 6 + 5 sen x − 7 dx = 4 x + 6 dx + 5 sen x dx − 7 dx = 3 [x+6] − 5 cos x − log 7
(insistimos en que no lo escribiremos nosotros, pero que no olvidaremos que podemos añadir +K ).

Es falso que la integral de un producto sea el producto de las integrales por no serlo la derivada,
pero de la fórmula del producto ( f g)0 = f g0 + f 0 g obtenemos:
Integración por partes. Sean f 0 y g0 continuas (para que existan las primitivas). Entonces:
b R
f (x)g0 (x)dx = f (x)g(x) − f 0 (x)g(x)dx ; ab f (x)g0 (x)dx = f (x)g(x) a − ab f 0 (x)g(x)dx
R R R

Esto reduce el problema a calcular otra primitiva, que será más sencilla si f 0 y g lo son (o si
una de ellas lo es y la otra no es más complicada que la anterior).
Con la notación d f ≡ f 0 (x)dx , la integración por partes se escribe
R R
udv = uv − vdu .

u = x , dv = sen x dx
R h i R
Ej. x sen x dx = = −x cos x − (− cos x) dx = −x cos x + sen x .
→ du = dx , v = − cos x

xe−x dx = [ u = x , dv = e−x dx → du = dx , v = −e−x ] = −x e−x + e−x dx = −(x+1) e−x .

R R
Ej.
[las primitivas de sen x y e−x no son peores que ellas, pero la derivada del x sí es más sencilla].
Otras funciones que mejoran al derivarlas son los logaritmos (las potencias de x no se complican):
R√ √
x log |x|dx = u = log |x| , dv = x dx = 23 x3/2 log |x| − 23 x3/2 dx 3/2 2 log |x| − 4 .
  R  
Ej. x =x 3 9

Algunas veces conviene tomar g0 = 1 (es decir, dv = dx ):

Ej. log x dx = u = log x , dv = dx → du = dx
R   R
x , v = x = x log x − dx = x log x − x .
Z
x
= x arctan x − 21 log (1+x2 ) .
R
Ej. arctan x dx = [ u = arctan x , dv = dx ] = x arctan x − 1+x2

Otras veces hay que repetir la integración por partes:

R 2 x
x e dx = x2 ex − 2 x ex dx = x2 ex − 2xex + 2 ex dx = [x2 − 2x + 2] ex .
R R
Ej.
u↑ dv↑ u↑ dv↑
log x dx log x dx log x dx
= 21 [log x]2 [se podía haber hecho a ojo].
R R R
Ej. Otro truco: x = log x log x − x → x

Combinando las dos últimas ideas:

Ej. I = sen x ex dx = sen x ex − cos x ex dx = ex [sen x − cos x] − I ⇒ I = 12 ex [sen x − cos x] .
R R

u↑ dv↑ u↑ dv↑
R R
dx dx
→ du = dx, v = − 1x dx
 
Ej. Curiosidad: x = u = x, dv = x2
= −1 + x ¿ ⇒ ? 0 = −1 !!
[no olvidemos que hay una K arbitraria aunque no la escribamos].

97
 P(x)
Z
Primitivas de funciones racionales: Q(x) dx , con P y Q polinomios .
P
Si el gr P ≥ gr Q , dividimos: Q = C + QR con el resto R de grado menor que Q .
Vimos que
Q se puede escribir como producto de polinomios del tipo (x−a)m [raíces reales] y
n
x2 +cx+d [complejas], siendo m y n la multiplicidad de las raíces [ m = 1 si son simples].
[El problema es que (como vimos en 3.3), salvo para Q especialmente sencillos, realizar
esta descomposición es, en la práctica, imposible por ser imposible hallar sus raíces].
R
Se prueba que Q se puede escribir como suma de constantes por funciones del tipo:
1 1 x
, y , con 1 ≤ j ≤ m , 1 ≤ k ≤ n (llamadas fracciones simples).
(x−a) j (x2 +cx+d)
k
(x2 +cx+d)
k

Para ‘descomponer en fracciones simples’ QR (hallar la constante que acompaña a cada fracción)
basta resolver un sistema lineal de ecuaciones. Y así, el problema de integrar P/Q se reduce,
una vez factorizado Q , al de integrar el polinomio C y funciones como las últimas.
4 3 +5x2 −11x+4 R(x)
Ej. I = 4xx5−6x
Z Z
4 3 2
−x +x −3x +2x
dx = Q(x) dx (ya es 4 < 5 ) . Empezamos factorizando:

Q(x) = x (x−1)2 (x2 +x+2) [suerte hemos tenido] y descomponemos en fracciones simples:
R(x) A(x4 −x3 +x2 −3x+2)+B(x4 +x2 −2x)+C(x3 +x2 +2x)+(Dx+E)(x3 −2x2 +x)
= Ax + x−1
B
+ C
+ xDx+E
2 +x+2 =
Q(x) (x−1)2 x(x−1)2 (x2 +x+2)
Bm B Dn x+En n D x+E
Si hubiera (x−1)m escribiríamos x−1
 2

1
+ · · · + (x−1)m ; si (x +x−2) , x2+x+2 + · · · + 2
1 1
(x +x+2)
n .

Igualando coeficientes de x4 , x3 , x2 , x y la constante de ambos términos se obtiene el sistema:

A + B + D = 4 , −A +C − 2D + E = −6 , A + B +C + D − 2E = 5 , −3A − 2B + 2C + E = −11 , 2A = 4
Resolviéndolo: A = 2 , B = 1 , C = −1 , D = 1 , E = −1
R 2dx
R dx
R dx
R (x−1)dx
⇒ I= x + x−1 − (x−1)2
+ x2 +x+2
Z
dx
Las (x−a)m
son casi inmediatas. Más trabajo dan las otras. Primero se busca un logaritmo:
1 (2x−2)dx 1 (2x+1)dx
− 32 dx
R R R
2 x2 +x+2
= 2 x2 +x+2 x2 +x+2
h i
2
2
Y luego un arco tangente completando el cuadrado: x2 +x+2 = (x+ 12 ) + 47 = 7
4
2x+1
√
7
+1 .
√
(x−1)dx 2/ 7 dx
Z
= 12 log(x2 +x+2) − 3 √2
R
Por tanto: x2 +x+2 2 7 √ 2 →
([2x+1]/ 7 ) + 1
1
+ 12 log (x2 +x+2) − √3 arctan 2x+1


I = 2 log |x| + log |x−1| + x−1 √ .
7 7

x4 −5x2 +x+8 x+4

R R R
Ej. I = x3 +x2 −4x−4
dx = (x−1) dx + (x+1)(x+2)(x−2) dx [de nuevo las raíces eran sencillas].

x+4 A B C A(x+2)(x−2)+B(x+1)(x−2)+C(x+1)(x+2)
(x+1)(x+2)(x−2) = x+1 + x+2 + x−2 = (x+1)(x+2)(x−2)

Cuando haya tantas raíces reales, mejor que igualar coeficientes se hace x = a para cada raíz a :
1 1
x = −1 → −3A = 3 , A = −1 ; x = −2 → 4B = 2 , B = ; x = 2 → 12C = 6 , C = 2 2
2 −4|
h i
→ I = 21 x2 − x − log |x+1| + 12 log |x+2| + 12 log |x−2| = 21 x2 − 2x + log |x
|x+1| 2

Z
dx

Para hallar las primitivas de las fracciones simples más complicadas como (x2 +x+2)
n

se utilizarían fórmulas de reducción como la propuesta en problemas .

98
Cambios de variable:
Supongamos que buscamos una primitiva de f (g(x))g0 (x)dx (con f y g0 continuas). Si F es
R

una primitiva de f , por la regla de la cadena: (F ◦g)0 (x) = F 0 (g(x))g0 (x) = f (g(x))g0 (x) . Así
pues, F ◦g es la primitiva buscada. Basta pues integrar la f y evaluar el resultado en g(x) . Si
lo que queremos es la integral definida entre a y b su valor es F(g(b)) − F(g(a)) . Por tanto:
R g(b)
f (g(x)) g0 (x) dx = f (u) du u=g(x) ; ab f (g(x)) g0 (x) dx = g(a) f (u) du
R R  R

En la práctica se suele usar la notación de diferenciales: se escribe u = g(x) , du = g0 (x)dx

y si hay límites de integración es fácil recordar que: x = a → u = g(a) , x = b → u = g(b) .
En algunos casos la g0 (x) aparece explícitamente y es claro el cambio que hay que hacer:
1
sen3 2x cos 2x dx = [ u = sen 2x, du = 2 cos 2x dx ] = u du = 81 u4 = 81 sen4 2x .
R R 3
Ej. 2
[en casos tan sencillos no será necesario escribir la sustitución, es fácil ver a ojo que
sen4 2x es casi la primitiva; derivándola mentalmente se ve que falta el 1/8 ].

u = log x , du = dx
 
R 5 dx R log 5
Ej. e x log x = x = 1 du u = log | log 5| − 0 = log (log 5) .
x = e → u = 1 , x = 5 → u = log 5
[podíamos haber calculado la primitiva olvidando límites de
integración y sustituir al final, una vez deshecho el cambio].
R x√  R√
Ej. e ex − 1 dx = u = ex , du = ex dx = u − 1 du = 23 [u−1]3/2 = 23 [ex −1]3/2 .


u = ch x, du = sh x dx
R h i R
Ej. sh3 x e− ch x dx = 2 2 = [u2 −1]e−u du = [partes]
sh x = ch x − 1
= −u2 e−u + 2 ue−u du + e−u = [partes] = [1−2u−u2 ]e−u + 2 e−u du = −[1+2u+u2 ]e−u
R R

= − 1 + 2 ch x + ch2 x e− ch x .
 

u = sh2 x
h i
o partes directamente: sh3 xe− ch x dx = = − sh2 x + 2 ch x sh xe− ch x dx
R R
dv = shx e− ch x dx
= [partes] = − sh2 x e− ch x − 2 ch x e− ch x + 2 sh x e− ch x dx = −[2 + 2 ch x + sh2 x] e− ch x .
R

Pero en la mayoría de los casos no es tan evidente el cambio ni hay una clara du . La forma
del integrando puede sugerir hacer algún cambio u = g(x) . Para obtener entonces la f (u) se
despeja la x en función de u , se calcula el dx y se sustituyen x y dx en la integral (sin
olvidar el cambio de límites de integración, si los hay):
R9 √ √
cos x dx = u = x , x = u2 , dx = 2udu , x = 4 → u = 2 , x = 9 → u = 3 = 2 23 u cos u du =
  R
Ej. 4
R3
[partes] = 2[u sen u]32 − 2 2 sen u du = 2[u sen u]32 + 2[cos u]32 = 2[3 sen 3 − 2 sen 2 + cos 3 − cos 2] .
R√ R √u−1 √
ex −1 dx = [ u = ex , x = log u , dx = du u−1 = z , u=z2 +1 , du=2zdz

Ej. u ]= u du =
R 2 R √ √
= 2 z2z+1 dz = 2 1 − z21+1 dz = 2z − 2 arctan z = 2 ex −1 − 2 arctan ex −1

√
[con un poco más de vista podríamos haber hecho directamente z = ex −1 acabando antes].

[En los cambios de variable las funciones f y g0 deben ser continuas. Si g0 aparece explícitamente
es fácil ver que es así. Pero si no aparece se nos podría olvidar y cometer errores].

La práctica sugiere qué y cuándo sustituir. Para tipos concretos de funciones (trigonométricas,
con radicales...) hay cambios típicos que se sabe que dan buen resultado (los veremos a conti-
nuación) y que suelen llevar a la integración de funciones racionales de las que hemos visto.

99
Primitivas de funciones trigonométricas
[Aparecen muy a menudo, ya que están muy ligadas a la integración en polares].
R
Para integrar  x, cos x) dx , con R función racional en sen x y cos x , existe siempre un
R(sen
x

cambio u=tan 2 que la lleva a una racional, pero veamos antes una serie de casos más fáciles.
R
senm x cosn x dx :
k
sen2k+1 x = sen x (1 − cos2 x) y se hace u = cos x .
Si m ó n son impares: k
cos2k+1 x = cos x (1 − sen2 x) y se hace u = sen x .
Si m y n pares, se escriben en función del ángulo doble: sen2 x = 1−cos
2
2x
, cos2 x = 1+cos
2
2x
.

sen2 x cos3 x dx = (1−sen2 x) sen2 x cos xdx = [u = sen x] = (u2 −u4 ) du = 31 sen3 x − 15 sen5 x .
R R R
Ej.
o a ojo

Ej. cos 2x sen x dx = cos2 x sen x−sen3 x dx = 2 cos2 x sen x−sen x dx = cos x − 23 cos3 x .
R R  R 

R1
O bien sen a cos b = 12 [sen(a+b)+sen(a−b)] : 1 1
2 [sen 3x−sen x] dx = 2 cos x − 6 cos 3x .
↑
[Este tipo de igualdades se usan para integrar productos de senos y cosenos de distinto ángulo].
1R 1R 1R 1R 3x sen 2x
4 2 + sen324x .
R
Ej. cos x dx = 4 (1+cos 2x) dx = 4 dx + 2 cos 2x dx + 8 (1+cos 4x) dx = 8 + 4
R
La integral general R(sen x, cos x) dx se convierte en cociente de polinomios haciendo:
u = cos x , si R es impar en sen x [ es decir, si R(− sen x, cos x) = −R(sen x, cos x) ].
u = sen x , si R es impar en cos x [ es decir, si R(sen x, − cos x) = −R(sen x, cos x) ].
1 du
 2 
u = tan x cos x = 1+u 2 , dx = 1+u2 , si R(− sen x, − cos x) = R(sen x, cos x) .
x 2u 1−u2 2du
 
u = tan 2 sen x = 1+u2 , cos x = 1+u2 , dx = 1+u 2 , para cualquier R [último recurso].

Z Z
dx sen x dx R du 1R du 1 R du 1 u−1 1 cos x−1
Ej. sen x = 1−cos 2 x = [ u = cos x ] = u2 −1
= 2 u−1 − 2 u+1 = 2 log | u+1 | = 2 log | cos x+1 | .
2du/[1+u2 ]
Z Z
dx
= [ u = tan 2x ] =
R du x
O de otra forma: sen x 2u[1+u2 ]
= u = log | tan 2 | .
[Ha salido tan fácil por casualidad; las dos expresiones de la primitiva deben
coincidir salvo K arbitraria (con pocas cuentas se ve que son iguales)].
2
[1+u2 ] du
Z Z Z
dx
= cos4dx 1
R du R 2
Ej. cos3 x sen x x tan x
= [ u = tan x ] = u[1+u2 ]
= u + udu = log | tan x| + 2 tan x .
u=cos x
Z Z Z Z 
1/2 1/2
dx
= cos3sen x dx du
− 1u − u13 du

Más largo: cos3 x sen x x[1−cos2 x]
= u3 [u+1][u−1]
= · · · = u+1 + u−1
= 21 log |1−u2 | − log u + 2u12 = log | sen x| − log | cos x| + 2 cos
1
2x .

Peor todavía sería hacer u = sen x (también es impar en coseno) ó u = tan 2x ;



(1+u2 )3
Z 
por ejemplo, con el último cambio queda la complicada primitiva u(1−u2 )3
du .
Z π Z 0
dx du du
 
Ej. 2 = u = tan x , dx = 1+u2
= 2+u2
=0 1
__1
___
0 1+cos x 0
1/2 1+cos2x
[resultado evidentemente falso: el integrando es siempre positivo
!
y la integral debía ser un número positivo. No olvidemos que en
los cambios de variable las funciones f y g0 deben ser continuas. El cambio hecho (clásico, como
hemos dicho, para este tipo de integrales) es válido sólo hasta π2 ; sí es cierto que
Z π/2 Z ∞ Z ∞ √ i∞ √ √
dx du √1 1/ 2 du 1 u π 2 Rπ π 2
2
1+cos x
= 2+u 2 = 2
√ 2 = √ arctan √
2 2
= 4 → 0 = 2 .
0 0 0 1+[u/ 2 ] 0
π
pues el integrando es simétrico respecto a x = . Al ∞ que nos ha aparecido (que como siempre
2
representará un límite) le daremos más seriedad cuando estudiemos las integrales impropias].

100
Primitivas de irracionales
(las más simples; R función racional de x y de la raíz que se indica).
R √
√
R x, n ax + b dx se convierte en racional haciendo u = n ax + b .

u = [1+x]1/4 , x = u4 −1
 
R 9 5
x[1+x]1/4 dx = = 4(u8 −u4 )du = 4u9 − 4u5 = 49 [1+x]9/4 − 45 [1+x]5/4 .
R
Ej. 3
dx = 4u du
También se puede hacer por partes:
x[1 + x]1/4 dx = 54 x[1 + x]5/4 − 45 [1 + x]5/4 dx = 45 x[1+x]5/4 − 16 9/4 .
R R
45 [1+x]
√
u = x − 4 , x = u2 +4
Z 5  
= 01 u22u = 01 2(u−2+2) du
= 01 u−2
dx R du R R 2du R 1 4du
Ej. √ = −4u+4 (u−2)2
+ 0 (u−2)2
4 x−4 x−4 dx = 2u du
1 4 1

= 2 log |u−2|0 − u−2 0
= 2(1−log 2) .

R √

R x, a2 − x2 dx se convierte en trigonométrica haciendo x = a sen u .
R√
4 − x2 dx = [ x = 2 sen u , dx = 2 cos u du ] = 4 cos2 u du = 2u + sen 2u
R
Ej.
√ √
= 2u + 2 sen u 1− sen2 u = 2 arc sen 2x + 2x 4 − x2 .
R √
√
R x, x2 +a dx se convierte en racional haciendo u = x + x2 +a ,
1
puesto que (u−x)2 = u2 −2xu+x2 = x2 +a → x = u2 − 2u
a a


→ dx = 2 + 2u2 du .
√

El cambio u = x2 +a no sirve de nada pues vuelven a aparecer raíces al despejar la x .
Z √ 2 2
 
= u = x + x2 +1 , x = u 2u−1 , dx = 1+u
 R 2du
√dx u−1  √x

Ej. du = = log | | = log .

2u2 2
u −1 u+1 
x x2 +1 2 1+ x +1
Z √
Ej. √x dx = x2 + 1 [¡a ojo! , antes de ponerse a calcular a lo loco, miremos si es inmediata].
2 x +1
 √ 
Las primitivas con raíces ax2 +bx+c se reducen a las últimas completando cuadrados .
√ √
Recordamos que si las raíces son más complicadas como x3 + a ó 3 x2 + a , las integrales, son, en

general, no calculables. Esto no quiere decir que alguna, en particular, lo sea:

Z 7 R √ 3 √
√x dx = [t = x4 ] = 14 √tdt = [ u = t+1 ] = 12 (u2 − 1)du = u6 − u2 = 16 [x4 −2] x4 + 1 .
R
Ej. 4 t+1
x +1
Z
Pero no se podría hallar la primitiva de √dx4
 
.
x +1

Otro tipo de primitivas que se convierten en racionales mediante cambios de variable son:
R
R(ex ) dx , siendo R función racional de ex .

Haciendo u = ex se convierte en la racional R(u) du

R
u , pues dx = u .

ex dx A(1+u2 )+Bu2 +Cu = 1

h
dx du
Z Z Z Z  i
= Au + Bu+C

Ej. = = du =
1+e2x ex (1+e2x ) 2
u(1+u ) 1+u2
C = 0 , A = 1 , B = −A
Z Z
= du
u − u du
1+u2
= log u− 12 log(1+u2 ) = x − 12 log (1+e2x ) .
[Aunque (parecido) también se podía hacer con u = e2x → 1 R du
2 u(1+u) = 12 log u− 12 log(1+u) ↑

101
Más ejemplos variados de cálculo de integrales:
√
Z 1
Partes: dv = dx 1
h i
2x2 +6−6
1 Z π 3
log 3+x2 dx = = x log(3+x2 ) 0 −


Ej. 2 3+x 2 dx = 2 log 2 − 2 + 3 ,
0 u = log(3+x ) 0
Z 1
6 dx
√ Z 1 dx/√3 √ x 1

pues 3+x 2 =2 3 √ 2 = 2 3 arctan √ 0 .
3
0 0 1+(x/ 3 )

Z 4 √
√
Ej. log 1+ x dx . El aspecto del integrando sugiere hacer el cambio: t = x , dx = 2t dt .
0
Los nuevos límites de integración son x = 0 → t = 0 , x = 4 → t = 2 y la integral queda:
2 R 2 dt  2 2
2 02 t log(1+t) dt = t 2 log(1+t) 0 − 02 tt+1 = (t −1) log(1+t) − 12 (t−1)2 0 = 3 log 3 .
R
↑ &
t 2 −1+1 1
partes = t+1 = t − 1 + t+1
R
Ej. cos(log x) dx . Apetece, para simplificar el integrando, hacer log x = t , dx = et dt , que lleva a
R
et cost dt = et cost + et sent dt = et [cost+sent] − et cost dt
R R
R
⇒ et cost dt = 12 et [cost+sent] = 12 x cos(log x)+sen(log x) .
 

Z
Ej. √1 + x dx . En la página anterior tenemos el cambio que la lleva a racional: x = 3 sen u .
9−x2 Z √
→ 1+ 3 sen u x 2
3 cos u 3 cos u du = u − 3 cos u = arc sen 3 − 9−x .

− √−x dx2 .
Z Z
dx/3
Pero hubiéramos acabado antes observando que era ‘casi inmediata’: √
1−(x/3)2 9−x

ex =u partes
Z
e2x cos(ex ) dx = u cos u du = u sen u − sen u du = u sen u + cos u = ex sen ex + cos ex .
R R
Ej.

dx 2 dx 2 ex
Z Z Z
Ej. ch x = ex + e−x = dx = 2 arctan ex . [A ojo, no hay que hacer ningún cambio].
1 + e2x
Z 1
Ej. sh3 x dx = 0 , sin nada que calcular. El integrando es impar (como sh x ).
−1 R
ch2 x −1 sh x = 31 ch3 x −ch x .


Aunque se podría hallar la primitiva:
Z 1
dx
Ej.
−2 1+|x| . El |x| tiene dos expresiones en el intervalo, con lo que hay que dividir la integral:
R 0 dx R 1 dx 0 1
−2 1−x + 0 1+x = − log |1−x| −2
+ log |1+x| 0
= log 3 + log 2 .
Z 1
dx 1
Ej.
−2 1+x . Esta parece más inofensiva y podría pensarse que es: log |1+x| −2 = log 2 .
Pero esto es absolutamente falso. La integral se ha definido para funciones acotadas
y existía para funciones continuas (o a trozos) y los teoremas fundamentales eran
también para ese tipo de funciones. Nuestro integrando tiene una asíntota en x=−1
(dentro del intervalo) y la integral no existe.
[Ni siquiera existirá como una integral impropia de las que vamos a definir en la próxima sección].

102
5.4. Integrales impropias
La integral la hemos definido para funciones f acotadas en intervalos [a, b] finitos. Extendemos
la definición, primero para intervalos de integración no acotados [a, ∞) ó (−∞, b] . Como siempre
que aparece un ∞ aparecerá un límite en la definición:
Rb Rb
Supongamos que a f existe para todo b≥a . Si existe el lı́m f se le llama integral
R∞ R∞ b→∞ a
impropia de f en [a, ∞) , se representa por a f ó a f (x)dx y la integral impropia
Def.
se dice convergente. Si a∞ f no es convergente, se dice divergente.
R
Rb
f = lı́m ab f .
 R 
Análogamente se define −∞
a→ −∞

[La integral entre a y b existe ∀b , como sabemos, si por ejemplo f es continua (o continua
a trozos) en [a, ∞) ; como para cada b la integral es un número, tenemos una función de b y
tiene sentido hablar de su límite cuando b → ∞ ; este límite (el valor de la integral impropia)
será otro número si la integral converge].
R b dx b
= lı́m [− 1x ]1 = lı́m [1 − 1b ] = 1 .
R ∞ dx
Ej. 1 x2 = lı́m
b→∞ 1 x
2
b→∞ b→∞ deja un área
1 1/x infinita debajo
[la integral es convergente y su valor es 1 ].
R ∞ dx 1/x2
Ej. 1 x = lı́m [log x]b1 = lı́m [log b] diverge. 1
b→∞ b→∞ área=1

R ∞ dx  1

En general, 1 xs diverge si s ≤ 1 y converge si s > 1 hacia s−1 .
[es inmediato comprobarlo; para los mismos s converge la a∞ dx
R
Rb Rb xs ∀a
Ra
> 0 , pues
a e 1 son dos funciones de b que sólo difieren en la constante 1 ]
1 ax b 1
hacia − 1a y diverge si a ≥ 0 .
R ∞ ax ab
 
0 e dx = a lı́m [e ]0 = a lı́m [e −1] converge si a < 0
b→∞ b→∞
b
[Se suele abreviar [eax ]∞
0 en lugar de lı́m [eax ]0
b→∞
; pero no olvidemos que es un límite].

Aunque no sepamos calcular la primitiva podremos, en bastantes ocasiones, determinar si es o no

convergente (como ocurría con las series; incluso teníamos un criterio integral que relacionaba
unas y otras; los criterios de convergencia son muy parecidos).
Rb

Criterios para funciones positivas los damos para la a∞ ; son análogos para −∞
R
.
En todos suponemos que las funciones que aparecen son integrables en [a, b] ∀b .
Teorema:
R∞ R∞ R∞ R∞
Si 0 ≤ f (x) ≤ g(x) para x ≥ a , a g converge ⇒ a f converge, e a f≤ a g.
Rb Rb R∞
g 0 ≤ F(b) = a f≤ a g≤ a g ∀b ≥ a ⇒
f
F(b) creciente y acotada superiormente ⇒ F(b) tiene límite si b → ∞
a b (la última ⇒ se prueba como en las sucesiones).
!
[El teorema dice también que a∞ f divergente ⇒ a∞ g divergente,
R R
"a
g
F(b)
desde luego; pero como siempre, en este tipo de criterios, de que
a la pequeña converja o de que la gorda diverja, no se sigue nada; e
insistimos en que es para funciones positivas: si una f cualquiera es menor que otra de integral
convergente, no tiene que converger su integral, ya que podría irse a −∞ ].

103
Las comparaciones con ≤ son siempre más complicadas que las hechas por paso al límite:
Teorema:
f (x)
Si f y g son positivas para x ≥ a y x→∞
lı́m g(x) = c finito, entonces:
R∞ R∞
Si c > 0 , a g convergente ⇔ a f convergente.
c = 0 , a∞ g convergente ⇒ a∞
R R R∞ R∞
Si f convergente [es decir, a f diverge ⇒ a g diverge].
f (x)
Si c > 0 , para x ≥ M es 2c ≤ g(x) ≤ 3c c 3c
2 ⇒ 0 ≤ 2 g(x) ≤ f (x) ≤ 2 g(x) y podemos
aplicar el teorema anterior. Si c R= 0 , para x ≥ M Res 0 ≤ f (x) ≤ g(x) y de nuevo el
teorema. Además está claro que M∞ f converge ⇔ a∞ f converge.
Si el integrando f no es positivo, como en las series, conviene considerar el | f | :

R∞ R∞ f se dice absolutamente
Teorema: a | f | convergente ⇒ a f convergente 
integrable en [a, ∞) .
R   R   R
0 ≤ f +| f | ≤ 2| f | ⇒ a∞ f +| f | convergente ⇒ a∞ f = a∞ f +| f | − a∞ | f | convergente.
R

[log x]2 [log x]2

Z ∞
1 R ∞ dx
Ej. x dx diverge, pues si x ≥ 3 es x ≥ x e 3 x diverge.
3

Por paso al límite debemos utilizar la parte con c=0 ya que log x no se parece a ningún xs :
1/x R ∞ dx R ∞ [log x]2
→ 0 e 3 x divergente ⇒ 3 x dx diverge (mayor que divergente).
[log x]2 /x x→∞
R ∞ [log x]2
dx = 13 [log x]3
∞
También nos bastaba la definición: 3 x 3
→∞.
Z ∞ √
xdx x 1 x/ x5 −x+1
√ . Cuando x → ∞ , √
 
Ej. ∼ x3/2
es decir, 1/x3/2
→ 1 .
0 x5 −x+1 x5 −x+1 x→∞
R∞ 1
Como 1 3/2 converge, la dada también (no sabemos a qué número).
x

R ∞ −x2 −x2 2
dx (sin primitiva elemental) converge, pues ee−x = ex−x → 0 e 0∞ e−x dx converge.
R
Ej. 0 e x→∞
2
O bien, por desigualdades: si x ≥ 1 es e−x ≤ e−x y de aquí:
R ∞ −x R∞ R ∞ −x2
converge (⇔ 0∞ converge).
R
1 e dx converge (⇔ 0 converge) ⇒ 1 e
R ∞ −x2 1√
 
con técnicas de integrales dobles se puede ver que 0 e dx = 2 π .
sen(1/x)
sen 1x dx ∼ lı́m 1/x = lı́m+ sent
R∞ R ∞ dx  
Ej. 1 1 x divergente pues x→∞ t = 1 ⇒ la dada diverge.
t→0
Z ∞
sen x
 sen x  1 R ∞ dx 1


Ej. 3 dx es convergente pues  1+x 3 ≤ 1+x3 e 0 1+x3 converge ∼ cuando x → ∞ .

0 1+x x3

R ∞ sen x
√ dx no podemos concluir nada, ya que ∞ √1 diverge.
R
Ej. Aplicando la misma idea a 1 x 1 x
∞
√ dx = π + 2π + · · · ≡
R ∞ sen x
–
R R
Pero 0 x 0 π ∑ ak , donde 1/"x
k=1
a1
R kπ | sen
√ x|
R kπ dx
√ p 
|ak | = (k−1)π x
≤ (k−1)π
√
x
=2 kπ − [k−1]π . ! a2 2! 3! 4!
–
–1/"x
La serie es alternada, decreciente y con ak → 0 , con
lo que por Leibniz converge (y por tanto la integral). sen(x2 )
––
De aquí deducimos que –
"2!
"!
R ∞ sent
sen x2 dx = [ t = x2 ] =
R∞
√ dt también converge.
0 0 t
(¡a pesar de que f (x) no tiende a 0 si x → ∞ ! [esto no es como en las series]).

104
La segunda extensión de la integral es para f no acotada en un extremo del intervalo:

Supongamos que tb f existe para todo t ∈(a, b] . Se define ab+ f = lı́m+

R R Rb
t f si el límite
t→a
Def. existe y en ese caso la integral impropia se dice convergente.
R −
Análogamente: ab f = lı́m− at f .
 R 
t→b

(En vez de a+ y b− suele escribir a y b ; no olvidemos que la integral es impropia).

[No se pide que f esté acotada en (a, b] , ni siquiera que esté
definida en el punto a ; para que f sea integrable en [t, b] ,
debe, desde luego, estar acotada en cada intervalo de esa forma;
por ejemplo, si f es continua en (a, b] se tiene, para todo t ,
garantizada la existencia de la integral de f en [t, b] , aunque el a t b
límite puede no existir y divergir la integral impropia].
R 1 dx R1
Ej. lı́m+ t dx
0+ x2 = t→0 x2
= lı́m+ [ 1t − 1] no existe (la integral impropia diverge).
t→0
R 1 dx √
0+
√
x
= lı́m+ [2 − 2 t ] = 2 , converge (y su valor es 2 ).
t→0

Rb dx R a− dx
En general, se ve fácil que a+ [x−a]s e c [a−x]s convergen si s < 1 y divergen si s ≥ 1 .

[x − a]s tiene sentido para x < a si s = 13 , 72 , ... ; si s = 21 ó s = π la función no está definida .

Para este otro tipo de impropias existen criterios de convergencia totalmente análogos a los
vistos para las del primer tipo. Resumiendo (las de a+ ) y sin demostraciones:
Teorema:
Rb Rb Rb Rb
Si 0 ≤ f ≤ g en (a, b] , a+ g convergente ⇒ a+ f convergente e a+ f≤ a+ g .
f (x)
Sean f , g ≥ 0 en (a, b] y sea finito el lı́m+ g(x) = c , entonces:
x→a
Rb Rb Rb Rb
si c > 0 , a+ g converge ⇔ a+ f converge; si c = 0 , a+ g converge ⇒ a+ f converge.
Rb Rb
a+ | f | convergente ⇒ a+ f convergente.
R 1 cos2 x cos2 x 1 R1 1 cos2 x/x3/4

Ej. 0+ x3/4 dx converge, pues 0 ≤ x3/4
≤ x3/4
e 0+ x3/4 converge o porque 1/x 3/4 → 1 .
x→0+
R7 dx R 7 dx
Ej. 2+ x3 −8 diverge, pues se parece cerca de x = 2 a 2+ x−2 divergente:
1/[x3 −8] 1
1/[x−2] = → 1
x2 +2x+4 x→2 12
(o usando L’Hôpital).

R3 dx 1/ sen x R 3 dx
Ej. 0+ sen x . Cerca de 0 el sen x ∼ x : 1/x x→0→ 1 . Como 0+ x diverge, la dada diverge.

R ∞ sen x
Ej. La √ dx de antes, no plantea problemas en x = 0 , a pesar de anularse su denominador,
0+ x √
pues se parece cerca de 0 a x que no sólo converge, es continua.
R1 2 (log√x)2
Ej. 0+ (log x) dx es convergente, pues 1/ x
= (x1/4 log x)2 → 0 (lo sabemos desde 4.5),
x→0+
con lo que la nuestra es más pequeña que una convergente. Y podemos hallar su valor:
R1
(log x)2 dx = x(log x)2 − 2 log x dx = x(log x)2 −2x log x+2x ⇒ 2
R R
0+ (log x) dx = 2.

105
 R∞ R ∞ R b−
Hay otras integrales que reúnen más de un tipo de impropiedad: −∞ , a+ , a+ , ...
Cada integral de estas se dice convergente si, dividido el intervalo en subintervalos tales que en
cada uno de ellos haya una única impropiedad, todas las integrales resultantes convergen. Por
ejemplo, si f es continua en todo R:
R∞ R0 R∞ R∞ R0 R∞
−∞ f converge ⇔ −∞ f e 0 f convergen y su valor es −∞ f= −∞ f+ 0 f
b R ∞ R
[esta integral no se define como lı́m −b f que podría existir a pesar de ser −∞ f divergente;
b→∞
a ese límite de las integrales calculadas en intervalos simétricos [−b, b] , si existe, se le llama
valor principal de Cauchy de la integral y aparece en matemáticas más avanzadas].
R∞ R∞
Ej. −∞ sen xdx diverge, pues 0 sen xdx = lı́m [1 − cos b]
b→∞
R0
no existe (y tampoco existe −∞ sen xdx ). b senx
–b 0
[Sí existe el valor principal de Cauchy:
R∞ Rb
VP −∞ sen xdx = lı́m −b sen xdx = 0 ( sen x es impar)].
b→∞

R ∞ arctan x R ∞ R ∞ dx  arctan x/x+x2 

Ej. 0+ x+x2 . 1 es convergente pues comporta como 1 x2 convergente 1/x2
→ π2 .
x→∞
arctan x 1 R1
Cerca de 0 : x+x2
∼ 1+x con límite finito ⇒ 0+ converge.
R∞
Como convergen las dos, 0+ converge.
R ∞ dx R 1 R∞
Ej. 0+ xs = 0+ + 1 diverge ∀s : 1/x2
Si s > 1 converge la de ∞ , pero diverge la de 0+ ,
_
si s < 1 ocurre al revés y si s = 1 divergen ambas. 1 1/"x
R 1 dx 1/x
Ej. −1 x2 no es una integral normal y ni siquiera existe 1/x2
R0 R1 1
como impropia, pues no convergen ni −1 ni 0 .
área área
Tampoco existe el VP de Cauchy de la impropia, ! !
–1 1
f + b1 f ] .
Rb R

definido en estos casos por lı́m [ −1
b→0

R2
√x dx √x ∼ √x , pero 2∞ diverge pues ∼ 1x .
R∞
R

Ej. 1+ . 1+ converge pues √
x4 −1 x−1 x3 +x2 +x+1 1+ 2 x−1 ∞
Por tanto, la inicial diverge (insistimos en que deben converger las dos para ser convergente).
R ∞ 1−e−x2 1−e−x
2
dx converge, pues lo hacen 01+ (tiene límite en x=0 ) e 1∞ 1
R R

Ej. 0+ x2
es 0 ≤ x2
≤ x2
.
R ∞ x−1 −t
e dt . En x = 0 converge si y sólo si x > 0 (pues se parece a 0∞ t x−1 dt ).
R
Ej. Γ(x) = 0 t
t x−1 e−t t x+1 R ∞ dt R∞
En ∞ converge siempre: t −2
→ 0 ∀x
et t→∞ e 1 t2 converge. La 0 inicial converge ∀x > 0 .
La Γ(x) (función gamma) definida por esta impropia generaliza el factorial:
R ∞ x −t partes x −t ∞ R ∞ x−1 −t
0 t e dt = −t e ]0 + x 0 t
Γ(x + 1) = e = x Γ(x) ⇒ si n ∈ N ,
Γ(n+1) = nΓ(n) = n(n−1)Γ(n−1) = n(n−1) · · · 2 · 1 · Γ(1) = n! , pues Γ(1) = 0∞ e−t dt = 1 .
R

+ ±
R ∞ 1−cos x
Ej. 0+ x3 log x dx . Plantea problemas en 0 , 1 , ∞ . Para converger, deben hacerlo las cuatro.
+
R 1/2 dx 1/2
Analizamos todas. En 0 : ∼ 0+ x log x
= log(| log x|) t
−→ −∞ (diverge).
t→0 R +

1±
R dx R dx
En : ∼ log x ∼ x−1 (divergen ambas). En ∞ : ≤ dx
x3
(converge).

106
5.5. Integración aproximada
Como sabemos, funciones integrables pueden no tener primitivas elementales o exigir un cálculo
muy largo. Pero en muchas ocasiones, sólo se necesita el valor aproximado de una integral
definida (en otras, simplemente, cotas de dicha integral). Las Un y Ln de 5.1 (y algún teorema
con desigualdades visto en ella) nos daban ya alguna (mala) estimación, pero será más preciso
utilizar series de Taylor o utilizar las fórmulas sencillas (sobre todo para los ordenadores) de
los trapecios o de Simpson que veremos al final de esta sección.
Integración de series de Taylor.
Estas series se podían derivar término a término (en el intervalo de convergencia). Veamos que
también se pueden integrar término a término en ese intervalo (de nuevo como si se tratasen de
’polinomios infinitos’). Esto será consecuencia de los siguientes resultados:
Sea { fn } sucesión de funciones continuas que converge uniformemente
Teorema: Rb Rb
hacia f en [a, b] . Entonces a f = lı́m a fn .
n→∞
ε
Sea ε > 0 . Existe un N tal que si n ≥ N entonces | f (x)− fn (x)| < b−a para todo x ∈ [a, b] .
R b Rb  Rb Rb ε
Si n ≥ N ,  f (x)dx − fn (x)dx ≤ | f (x) − fn (x)|dx <
a a a dx = ε . a b−a
Este resultado es falso si la sucesión de funciones converge sólo puntualmente (el límite de las
integrales puede ser distinto de la integral del límite) como para la siguiente { fn }:
 2
 2n x , 0 ≤ x ≤ 1/2n
4
Ej. fn (x) = 2n − 2n2 x , 1/2n ≤ x ≤ 1/n .
0 , 1/n ≤ x ≤ 1

3
La gráfica de cada fn es un triángulo isósceles de altura n
sobre el intervalo [0, 1n ] y vale 0 en el resto de [0, 1] ; el área
2
encerrada por cada fn es 21 para todo n . El límite puntual de
las fn es f (x) = 0 para todo x ∈ [0, 1] ya que para cada x , a
1
partir de un N todas las fn (x) = 0 y fn (0) = 0 ∀n . Está claro
que { fn } no converge uniformemente y que se tiene:
R1 R1
0= 0 fn 6= lı́m f = 21 . 1/4 1/3 1/2 1
n→∞ 0

Como consecuencia inmediata de lo anterior, tenemos que:

Teorema:
∞ Rb ∞ Rb
Si ∑ fn converge uniformemente hacia f en [a, b] entonces a f=∑ a fn
n=1 n=1

∞ ∞ ∞
Ej. Como f (x) = ∑ senn2nx converge uniformemente en todo R, es
Rπ R π sen nx 2
0 f =∑ 0 n2 = ∑ [2n−1]3 .
n=1 n=1 n=1

Y en el caso particular de las series de potencias concluimos:

∞
Si f (x) = ∑ an xn para |x| < R ⇒
n=0
Teorema: Rx ∞ Rx ∞
an n+1
0 f (t)dt = ∑ 0 ant n dt = ∑ n+1 x = a0 x + a21 x2 + a32 x3 + · · · si |x| < R .
n=0 n=0

Pues en [−x, x] sabemos que la serie converge uniformemente.

[Fuera de (−R, R) la serie no convergerá y no servirá para aproximar niguna integral].
El conjunto de primitivas de f será, desde luego: f (x) dx = C + a0 x + a21 x2 + a32 x3 + · · · .
 R 

107
 R1
Ej. Calculemos aproximadamente 0 senx2 dx (función sin primitiva elemental). Tenemos que:
Rx R x 2 1 6 1 10 1 14 1 3 1 7 1 1
0 sent 2 dt = 11 15
0 [t − 3! t + 5! t − 7! t + · · · ]dt = 3 x − 42 x + 1320 x − 75600 x + · · · ∀x
R1 2 1 1 1 1
→ 0 sent dt = 3 − 42 + 1320 − 75600 + · · ·
y podemos aproximar la integral con las sumas parciales de esta serie alternada decreciente:
R1
0 ≈ 31 − 42
1
≈ 0.3095 con error menor que 1
1320 ≈ 0.0007 < 10−3
R1
0 ≈ 31 − 42
1 1
+ 1320 ≈ 0.310281 con error menor que 1
75600 ≈ 0.000013 ∼ 10−5
R1
0
1
≈ 31 − 42 1
+ 1320 1
− 75600 ≈ 0.310268158 con error menor que 1
9!·19 ≈ 0.000000145 ∼ 10−7
1
La misma serie de potencias nos da la integral para cualquier otro x . Por ejemplo, si x = 2 :
R 1/2 1 1 1 1 (converge mucho más rápidamente, pues
0 sent 2 dt =
− 5376 + 2703360
24 − 2477260800 + · · · cerca de x = 0 se parece más el desarrollo).
√
También vemos que si x = 2π ( ≈ 2.51 ) la integral es positiva (como sospechábamos en 5.2):
R √2π 2/3  2 4 4π 6
sent 2 dt = [2π]3 1− 2π7 + 2π

0 55 − 1575 +· · · ≈ 5.24 [1 – 2.82 + 3.54 – 2.44 + 1.06 – 0.31 + · · · ]
Las sumas parciales de la serie entre corchetes son: 1, – 1.82, 1.72, – 0.72, 0.34, 0.09, . . . (todo va
más lento ahora). Como es alternada decreciente (a partir de tercer término) su suma está entre
dos sumas parciales consecutivas, con lo que la integral es > 0. [Para dar su valor con un error
< 10−2 se ve que hay que sumar 8 términos (dos más) y se obtiene 0.43 ].
Como disponemos de su desarrollo de Taylor, aparte de las anteriores aproximaciones, podemos
realizar otras operaciones en la que aparezca la integral, como, por ejemplo, calcular algún límite
indeterminado:
Rx 1 8 1 8
3x 0sent 2 dt−x4 [x4 − 14 x +··· ]−x4 8
14 x +o(x ) 1
lı́m 8 = lı́m 1 = 8 8 = − 14
x→0 arctan x x→0 x8 − 3 x24 x +o(x )
Rx
3 sent 2 dt+3x sen x2 −4x3 2 2 2 2
lı́m [1+x16 ] 0 = lı́m 6 sen x +6x cos x6 −12x = · · · .


Por L’H más largo:
x→0 arctan 8x7 x→0 arctan 56x

R1 2
Ej. Encontremos cotas racionales de I = 0 g si g(x) = x2 e−x (de primitiva no calculable).
R1 R1
Las cotas más sencillas, pues claramente 0 ≤ g(x) ≤ 1 , son 0 = 0 0≤I≤ 0 1 = 1.
Podemos mejorar la cota superior hallando el máximo de g en [0, 1] :
2 R 1 −1 e>2.7
g0 (x) = 2x(1−x2 )e−x ⇒ máximo si x = 1 y g(1) = e−1 ⇒ I ≤ 0 e ≤ e−1 < 10 27 .
Si comparamos en [0, 1] con diversas funciones integrables:
1
g(x) ≤ x2 ⇒ I ≤ 13 x3 0 = 13 (mejor que la anterior)
2 2 1 e<2.8
g(x) ≤ xe−x ⇒ I ≤ − 21 e−x 0 = 12 [1 − e−1 ] < 12 [1 − 28 10 9
] = 28 (aún menor)
3 3 1
g(x) ≤ x2 e−x ⇒ I ≤ − 31 e−x 0 = 31 [1 − e−1 ] < 13 [1 − 28
10 3

] = 14 (más pequeña aún)
1
g(x) ≥ x2 e−x ⇒ I ≥ 01 x2 e−x dx = − [x2 + 2x + 2]e−x 0 = 2 − 5e−1 > 2 − 27 50 4
R 
= 27
partes
Pero si queremos obtener cotas con la precisión que necesitemos, lo mejor es usar Taylor:
I = 01 x2 − x4 + 21 x6 − 16 x8 + · · · dx = 31 − 15 + 14
1 1
R  
− 54 + · · · ∀x
1
⇒ 3 − 15 = 2
15 < 13 − 15 + 14
1 1
− 54 = 176
945 < · · · < I < · · · < 13 − 51 + 14
1
= 43
210 < 13 .
2 176 4
La cota inferior 15 es peor que la obtenida comparando, pero 945 > 27 ya la mejora.
43 43 3
Y la superior 210 es más pequeña que la menor de las anteriores: 210 < 14 .
[Con un ordenador se consigue mucha precisión ( I ≈ 0.189472 ), nosotros hemos conseguido
176 43
sólo deducir que 945 ≈ 0.186 < I < 210 ≈ 0.205 ; pero nos costaría poco sumar más términos].

108
 2x R1
Ej. Si h(x) = 8−x 2 , hallemos racionales que aproximen I = h con error menor que 10−2 .
0
1
Parece inútil aproximarla si podemos fácilmente dar el valor exacto: I = − log|8−x2 | 0 = log 87 .
El problema es que, sin calculadora, no sabemos el valor de ese logaritmo. Pero por Taylor:
log 1+ 17 = 71 − 2·49
1 1
− · · · serie de Leibniz → I ≈ 13 1 −2 .


+ 3·243 98 , con error < 729 < 10
Podríamos también desarrollar primero el integrando y luego integrar la serie:
∞  2 n ∞
2x x3 x5
8−x2
= 4x 1−x12 /8 = 4x ∑ x8 = 4x + 32 + 256 + · · · → I = 18 + 128
1 1
+ 1536 +··· = ∑ 1
n8n .
n=0 n=1
El problema de esta serie (que, desde luego, debe sumar lo mismo) es que no es alternada, lo
que hace menos fácil y mecánico estimar los errores.
∞ ∞ n
Sumando dos términos I ≈ 128 17
, el error cometido es ∑ n81n < 3·81 3 ∑ [ 18 ] = 3·7·8
1 −2 .
2 < 10
n=3 n=3

Rx
Ej. Dibujar la gráfica de r(x) = xx−1
4 +1 y hallar, si existen, los x en los que la función R(x) = 0 r ,

con x ∈ [0, ∞) , alcanza sus extremos.

0 3x4 −4x3 −1 4x2 [3x5 −5x4 −5x+3]
r ∈C∞ (R) . r(x) → 0 . r(x) > >
< 0 si x < 1 . r (x) = − 4 2 , r00 (x) = − 4 3 .
x→±∞ [x +1] [x +1]

P = 3x4 − 4x3 − 1 tiene 1 raíz positiva x+ [+ − −] y 1 negativa x− [+ + −].

P(−1) = 6 , P(0) = −1 , P(1) = −2 , P(2) = 15 ⇒ x− ∈ [−1, 0] y x+ ∈ [1, 2]
⇒ r decrece hasta x− , crece hasta x+ y decrece a partir de entonces.
x = −1 inflexión; en x = 0 no hay (no cambia de signo r00 ); los otros puntos de inflexión los
darían las raíces (ya no hay más enteras y negativas sólo la −1 ) de 3x4 −8x3 +8x2 −8x+3
se pueden hallar haciendo z = x + 1x .


–1 x–
1 x+ 2
3 1
Valores: r(−2) = − 17 , r(2) = 17 ,
r(−1) = −1, r(0) = −1 (Rolle confirma x− );
–1
r0 (−1) = − 32 , r0 (0) = 1 (otra vez x− ),...
A la vista de la gráfica de r : R decrece si 0 ≤ x ≤ 1 [añadimos áreas negativas] y luego crece
[lo mismo se deduce del signo ya analizado de R0 (x) = r(x) ]. El mínimo se da, pues, si x = 1 .
Aproximemos el valor de R(1) desarrollando por Taylor el integrando, que si |x| < 1 es:
?
[x−1][1−x4 +x8 −· · · ] = −1 + x + x4 − x5 − x8 + · · · → I1 = −1 + 12 + 51 − 16 − 91 + 10
1 1
+ 13 +···
En principio, quizás no podamos integrar hasta 1 , pero parece ir bien, pues la serie converge:
S1 = −1 < S5 ≈ − 0.578 < · · · < I1 < · · · < S7 ≈ − 0.401 < S3 = − 0.3
Podría no haber máximo de R en [0, ∞) . Si 1∞ fuese divergente (que no lo es, pues r(x) ∼ x−3
R

en el ∞ ), la R tendería a ∞ ; si fuese convergente y tendiese a un valor I2 > |I1 | , R tendería

hacia I1 + I2 > 0 (valor que no alcanzaría); y si converge hacia un I2 < |I1 | entonces el máximo
se alcanza en x = 0 (y vale R(0) = 0 ). Veamos que esto último es lo que sucede realmente:
r ≥ 0 en [1, ∞) . El criterio de comparación por desigualdades da cotas fáciles de la impropia:
∞
I2 ≡ 1∞ r < 1∞ x4x+1 = arctan x2 1 = 12 [ π2 − π4 ] = π8 < 0.4 , o bien:
R R 
∞
I2 < 1∞ x−1 = 3x13 − 2x12 1 = 12 − 13 = 16 < 0.17 , bastante mejor cota.
R 
x4
R1 R∞
R(x) → f+ f = I1 +I2 < 0 , según las cotas halladas ⇒ el máximo es R(0) .
x→∞ 0 1

[Con esfuerzo podemos hallar la primitiva R (el denominador lo factorizamos en 1.5):

Z x √ 2
√ √  √ √
t dt √
=· · ·= arctan x2 − 2
log x2 +√2x+1 − 2

4 arctan(x 2 +1)+arctan(x 2 −1)
R(x)= √
0 [t 2 +2t+1][t 2 − 2t+1] 8 x − 2x+1

y el valor exacto de ambas integrales. Con calculadora obtenemos: I1 ≈ −0.474 , I2 ≈ 0.149 ].

109
Fórmulas de los trapecios y de Simpson.
Para aplicar cualquiera de estos dos métodos no necesitamos la expresión analítica de f ; nos bastan
algunos de sus valores [situación que experimentalmente se presenta a menudo].
Trapecios:
Dividimos [a, b] en n partes iguales de anchura b−a
=h . T
n
R a+[k+1]h f(a+kh) f
Como aproximación de a+kh f tomamos el área del f
h f(a+[k+1]h)
trapecio T de la figura: 2 [ f (a+kh) + f (a+[k+1]h)] .
h
Entonces ab f será aproximadamente igual a la suma
R
h
a a+h b
de las áreas del los n trapecios: a+kh a+(k+1)h
Rb h h h
a f≈ 2 [ f (a) + f (a+h)] + 2 [ f (a+h) + f (a+2h)] + · · · + 2 [ f (a + [n−1]h) + f (a+nh)] ,
Rb h
  1 2 2 2 2 1
a f≈ 2 f (a) + 2 f (a+h) + 2 f (a+2h) + · · · + 2 f (a+[n−1]h) + f (b) a b

Simpson:
Una aproximación mejor se tendrá si, dividido [a, b] en un número par n = 2m de partes iguales de
longitud h = b−a b−a
n = 2m , en vez de sustituir cada trozo de f por una recta, la sustituimos por la parábola
que interpola la gráfica de f en tres puntos consecutivos:
x0 = a+kh , x1 = a+[k+1]h = x0 +h y x2 = a+[k+2]h = x0 +2h , f

es decir, por el polinomio: Q2 (x) = A0 + A1 (x − x0 ) + A2 (x − x0 )(x − x1 ) , Q2

con: A0 = f (x0 ) , A1 = h1 [ f (x1 )− f (x0 )] , A2 = 2h12 [ f (x2 )−2 f (x1 )+ f (x0 )] .
h h
Integrando Q2 se tiene tras algunos cálculos: x0 x1 x2
R a+[k+2]h R x0 +2h
Q2 (x)dx = 2hA0 + 2h2 A1 + 32 h3 A1 = h
 
a+kh f≈ x0 3 f (x0 ) + 4 f (x0 +h) + f (x0 +2h)
1 4 2 4 2 2 4 1
Y sumando las m integrales anteriores obtenemos: a b
Rb h
 
a f ≈ 3 f (a) + 4 f (a+h) + 2 f (a+2h) + 4 f (a+3h) + 2 f (a+4h) + · · · + 4 f (a+[n−1]h) + f (b)

Si se quiere utilizar con seriedad un método numérico se debe hablar del error cometido. Demos algún
resultado sin demostración. La estimación por trapecios es exacta si f (x) es una recta, función con
f 00 = 0 . No es de extrañar que el error dependa de f 00 . Puede probarse que:

Si | f 00 (x)| ≤ M2 para x ∈ [a, b] entonces: |error| ≤ 1

12 (b−a)M2 h
2

Se prueba que Simpson es exacto si f (x) = a + bx + cx2 + dx3 y que:

1
Si | f (4) (x)| ≤ M4 para x ∈ [a, b] entonces: |error| ≤ 180 (b−a)M4 h
4

Se ve que ambos métodos mejoran, como era esperable, cuando h → 0 , más rápidamente Simpson ya que
h4 tiende más fuertemente a 0 que h2 . Como las cuentas a realizar en ambos casos son casi las mismas,
será mejor acudir a Simpson si tenemos que aproximar una integral (hay métodos mucho mejores, pero
también más complicados).
R 1 4dx
Ej. Poco práctico, para comparar y ver si funcionan los métodos. Aproximemos 0 1+x2 (= π) :

Trapecios: h = 12 , n = 2 : 01 ≈ 44 1 + 2 54 + 21 = 31
R  
10 = 3.1

h = 14 , n = 4 : 01 ≈ 48 1 + 2 17
16
+ 2 54 + 2 16 4
R  
25 + 5 ≈ 3.1312

1 R1 4
1 + 4 54 + 12 = 47
 
Simpson: h= 2 ,n=2: 0 ≈ 6 15 ≈ 3.13
1 R1 4 16
+ 2 45 + 4 25
16
+ 54 ≈ 3.14157
 
h= 4 , n = 4 : 0 ≈ 12 1 + 4 17

110
Comparemos ahora los números que nos da trapecios y Simpson con los obtenidos por otros caminos en
los ejemplos anteriores de esta sección:
R1
Ej. Calculemos aproximadamente 0 senx2 dx (ya estimada utilizando Taylor):
1 R1 R1
h= 2 T. 0 ≈ 14 [0 + 2 sen 14 + sen 1] ≈ 0.334 S. 0 ≈ 16 [0 + 4 sen 41 + sen 1] ≈ 0.305
1 R1
h= 4 T. 0 ≈ 18 [0 + 2 sen 16
1
+ 2 sen 14 + 2 sen 16
9
+ sen 1] ≈ 0.316
R1 1 1
S. 0 ≈ 12 [0 + 4 sen 16 + 2 sen 41 + 4 sen 16
9
+ sen 1] ≈ 0.3099
1 R1 1 1 1
h= 6 T. 0 ≈ 12 [0 + 2 sen 36 + 2 sen 9 + · · · + sen 1] ≈ 0.313
R1 1 1 1
S. 0 ≈ 18 [0 + 4 sen 36 + 2 sen 9 + · · · + sen 1] ≈ 0.310205

1 R1 R1
h= 100 T. 0 ≈ 0.3105 S. 0 ≈ 0.3102683009
1 R1 R1
h= 1000 T. 0 ≈ 0.31026839 S. 0 ≈ 0.3102683017
[estos últimos valores exigen, desde luego, o una enorme paciencia o un ordenador].
Como f 00 (x) = 2 cos x2 − 4x2 sen x2 , f (4) (x) = (16x4 − 12) sen x2 − 48x2 cos x2 , en [0, 1] es
| f 00 | ≤ 6 , | f (4) | ≤ 4|4x4 − 3| + |48x2 | ≤ 60 → |error T| ≤ 12 h2 ; |error S| ≤ 12 h4
R√
[Para aproximar la integral de 5.2, 0 2π sent 2 dt , Simpson con n = 2 y n = 4 da,
respectivamente, 1.67 (la gráfica se parece muy poco a una parábola) y 0.42 ].

R 1 2 −x2
Ej. Para otra integral aproximada con Taylor I = 0 x e dx , Simpson da muy buenos resultados:

n = 2 → I ≈ 16 [0 + e−1/4 + e−1 ] ≈ 0.191

1 −1/16
n=4 → I≈ 1
12 [0 + 4 e + 12 e−1/4 + 49 e−9/16 + e−1 ] ≈ 0.18951 ].
[Lo largo de Simpson es acotar el error (tampoco sabemos con Taylor si sale serie no alterna-
da)].
R1 2x
Ej. Hallemos también con estos métodos algún racional que aproxime I = 0 h , h(x) = 8−x2
.

Probablemente Simpson daría un error < 10−2 con ya con h = 21 , pero necesitaríamos acotar la
h(4) ,
lo que es largo. Probemos con Trapecios que hay que derivar menos:
2
8+x2 ]
h0 = 2 2 , h00 = 4 x[24+x 00 100 100 2
3 → en [0, 1] es |h | ≤ 73 → |error| ≤ 12·343 h →
[8−x2 ] [8−x2 ]
1
→ I ≈ h(0) + 2h( 12 ) + h(1) = 41 0 + 31 8
+ 27 = 434
59
con error < 10−2 .
   
basta tomar h= 2

Ej. Por último, aproximemos con Simpson I1 = 01 r , para la r(x) = xx−1

R
4 +1 .

Tomemos h = 12 : I1 ≡ 01 r ≈ 16 − 1 − 4·8 49
R  
17 + 0 = − 102 ≈ – 0.480 (sin cota del error).
Número coherente con las cotas que calculamos con Taylor:
S5 ≈ − 0.578 < · · · < I1 < · · · < S7 ≈ − 0.401 ,
Y bastante similar al ‘exacto’ obtenido con la primitiva y utilizando la calculadora: I1 ≈ −0.474 .

111
5.6. Aplicaciones
Áreas planas
Ya vimos que la integral no es exactamente un área, sino una suma de áreas con signo. Por tanto,
para hallar el área encerrada entre el eje y=0 y la gráfica de una función f habrá que sumar las
integrales de f en los intervalos en que esté por encima del eje y restar las integrales cuando f
quede por debajo. Esto es equivalente a integrar | f | . Así:
Ej. Hallar el área de la región encerrada entre el eje horizontal y la gráfica de f (x) = x3 − x .
R1 R0 R1 R1 R1 |f|
(x − x3 )dx = 1 f
Área = −1 | f | = −1 f− 0 f = −2 0 f =2 0 2
f impar
 R1 
la −1 f (que es 0 por ser f impar) no representa el área sombreada . –1 1

Ej. Determinar el área de la región acotada por los ejes y la gráfica de tan(x–1)
h(x) = tan (x−1) [la de tan x trasladada una unidad a la derecha].
R1 R1 0
Área = 0 |h| = − 0 tan(x−1)dx = log | cos (x−1)| ]10 = − log(cos 1) > 0 1–!/2 1 1+!/2

(porque cos 1 < 1 ; las áreas deben salir siempre positivas).

Más en general, el área comprendida entre las gráficas de f y g en el intervalo [a, b]

Rb
viene dada por a | f − g| .

Ej. Determinar el área de la región encerrada entre las gráficas de f (x) = x2 − 2x y g(x) = x .
Las gráficas se cortan si x2 − 2x = x ⇔ x = 0 (y = 0) ó x = 3 (y = 3).
3
En [0, 3] la gráfica de g está por encima de la de f , por tanto:
R3 R3 –
g
A= [g − f ] = (3x − x2 )dx = 92 . x=1–! –
0 0 x=1+!
O bien de otra forma (más complicada en este caso, pero mejor en 3
√ –1 f
otros), integrando respecto a y : y = x2 − 2x ↔ x = 1 ± 1+y →
R0  √ √ √
1 + 1+y − (1 − 1+y) dy + 03 1 + 1+y − y dy
 R  
A = −1
2 3
0 h i
= 43 (1+y)3/2 −1 + y + 23 (1+y)3/2 − y2 = 92 .

0

Ej. Hallar (sin calculadora) con un error menor que 0.04 el valor aproximado del área de la región
acotada por la gráfica de h(x) = arctan x2 y la recta y = π4 .
h par , h → π
2 , h0 (x) = 2x
1+x4
(crece si x > 0 ) → !/2
|x|→∞
R1 h
corta y = π4 si x = ±1 → Área = π
2 −2 0 arctan(x2 )dx . !/4

Hallar la primitiva es posible pero largo:

–1 0 1
2x2 dx
R R
arctan(x2 )dx = x arctan(x2 ) − 1+x4
= ···
[y los log y arctan del resultado no podríamos evaluarlos sin calculadora].
Trapecios o Simpson dan valores desconocidos del arctan (y sería largo acotar el error). Mejor
es integrar el desarrollo de Taylor [se puede llegar hasta x=1 ] y usar la serie alternada que sale:
2 01 x2 − 13 x6 + 51 x10 − 17 x14 + · · · dx = 23 − 21
2 2
− · · · ≈ 47 , con error < 55
2 2
R  
+ 55 < 50 = 0.04 .
Por tanto, Área ≈ 1.571 − 0.571 = 1.00 con error < 0.04 (con ordenador: Área ≈ 0.974991 ).

112
Áreas en coordenadas polares.
Un punto P del plano de coordenadas cartesianas (x, y) se
puede describir además por un par de coordenadas polares y
y r
(r, θ ) siendo r la distancia de P al origen y θ el ángulo en "+arctan–x y
x !=arctan– x
radianes comprendido entre el semieje de las x positivas y x
la semirrecta que une el origen con P . Unas coordenadas se y
pueden obtener de otras utilizando que:
x = r cos θ , y = r sen θ ↔ r = x2 + y2 , θ = arctan xy [+π] , si θ ∈ (− π2 , π2 ) θ ∈ ( π2 , 3π
p  
2 ) .

! =# Para hallar el área de una región R como la del dibujo, acotada por
r=f(!) las semirrectas θ =α , θ =β y la curva r= f (θ ) , con f (θ )≥0 , di-
vidamos el ángulo β−α en n partes iguales (de longitud ∆θ ). Co-
mo el área de un sector circular de radio r y ángulo θ es r2 θ /2 , si
!=" mk y Mk son el mínimo y el máximo de f (θ ) en cada sectorcillo,
r=M1 se tiene que el área de cada uno de ellos está comprendida entre
r=m1 m2k ∆θ y Mk2 ∆θ → ∑ m2k ∆θ ≤ área de R ≤ ∑ Mk2 ∆θ . Como estas
son las sumas inferior y superior de f 2 en [α, β ] deducimos que:
2
Área de R = 12 αβ f (θ ) dθ .
R 
3

Ej. Hallar el área de la región acotada por la curva r = 3 + cos θ .

1 R 2π 2
R π 19 1 19π
A= 2 0 [3 + cos θ ] dθ = 0 [ 2 + 6 cos θ + 2 cos 2θ ]dθ = 2 .
–2 0 1 4
R no es el círculo de centro (1, 0) y radio 3 , con
√ área 9π y expresión
en polares: r2 −2r cos θ −8 =0 → r = cos θ + cos2p θ ; la curva dada
en cartesianas es muy complicada: x2 + y2 − x = 3 x2 + y2 .

De forma similar a las del área en polares se prueban las otras fórmulas de la sección:
Longitud de la gráfica R bq longitud=L
f(x)
L= a 1+[ f 0 (x)]2 dx .
de f en el intervalo [a, b] :
(Lo natural es probar la fórmula estudiando las ‘integrales de línea’).
a b

1
Ej. Hallar la longitud del tramo de parábola y = x2 que une los puntos (0, 0) y (1, 1).
y=x2 √
R √ h √ i R √ 2 ]2 √
0 1 L = 01 1+4x2 dx = u = 2x+ 1+4x2 = 81 12+ 5 [1+u u3
du = 2
5
+ log(2+ 5)
4 ≈ 1.48

Ej. Probar que la longitud L del tramo de y = x3 que une (0, 0) y ( 12 , 18 ) cumple 21 < L < 16
9
.
1/2 √
y0 = 3x2 → L = 0 1 + 9x4 dx (primitiva no calculable).
R
y=x3
5/4 f
f (x) ≡ [1+9x4 ]1/2 = 1+ 92 x4 − 81 8 4 √1
8 x + · · · si |9x | < 1 ⇔ |x| < 3 1/8 1
9 5 1/2
→ L = [x + 10 x − 98 x9 + · · · ]0 = 12 + 320
9 9
− 4096 +··· , 1/2
[válido por estar [0, 1/2] dentro del intervalo de convergencia]. La
1/2
serie de L es decreciente y alternada a partir de segundo término
1 169 9
→ 2 = S1 < S3 < · · · < L < · · · < S2 = 320 < 16 .
De otra forma (la integral puede describir el área limitada por f en [0, 1/2] ):
1 9
área rectángulo = 2 < L área trapecio = 16 [es fácil ver que f es convexa en ese intervalo].
[La acotación L > 1/2 era clara geométricamente antes de hacer ninguna cuenta].

113
Volúmenes (sencillos):
El instrumento natural para calcular volúmenes son las integra-
les múltiples del cálculo en varias variables, pero para hallar al- A(x)
gunos bastan integrales de funciones de una variable.
Volumen de un sólido que se extiende desde x = a hasta x = b ,
Rb a x b
conocida el área A(x) de cada sección plana: V = a A(x)dx

Ej. Un sólido tiene base circular de radio 2. Cada sección producida por un
x plano perpendicular a un diámetro fijo es un triángulo equilátero. Calcular
–2 2 2
el volumen del sólido.
x 2 √ √ √
A(x)= área triángulo de base 2 4−x2 = 3(4−x2 ), V =2 02 A(x)dx = 323 3
R

–2 f(x)
En particular, volumen de un sólido de revolución engendrado
al girar en torno al eje x la región comprendida entre la gráfica
a b
de f ( f (x) ≥ 0 ) y el eje x en [a, b] . El área de cada sección
[círculo de radio f (x) ] es
[ f (x)]2 dx .
Rb
A(x) = π[ f (x)]2 . Por tanto, V = π a

Ej. El volumen obtenido al girar la región determinada por

g(x) = 1x y el eje x en el intervalo [1, b] , b > 1 es
Rb 1 2 1
1 Vb = π 1 [ x ] dx = π[1 − b ] . Observemos que Vb −→ π :
b→∞
b
el volumen del sólido infinito es finito (impropia convergente).

1 Rb
Valor medio de una función en un intervalo; se define: M = b−a a f (x) dx

(si f ≥ 0 , es la altura de un rectángulo de base b − a M

y área igual a la limitada por la gráfica de f ) a b

2!/" Ej. Hallar el valor medio de f (x) = A sen ωx en el semiperiodo [0, ωπ ] :

!/" R ω/π
M = ωπ 0 A sen ωx dx = 2A el valor medio en [0, 2π


π ω ] es 0 .

Trabajo de una fuerza variable: un punto se mueve en el eje x sometido a una fuerza f (x)
función sólo de x . El trabajo realizado por f para mover el punto desde a hasta b es
Rb
T= a f (x) dx
0 b
Ej. El trabajo preciso para estirar un muelle una longitud b desde su
posición de equilibrio es 0b cxdx = 12 cb2 .
R

Sea una varilla de densidad lineal variable ρ(x) que ocupa desde a hasta b . Su masa m , su
centro de gravedad x∗ y su momento de inercia I respecto a 0 son:
Rb Rb Rb 2
m= a ρ(x)dx , x∗ = a xρ(x)dx , I = a x ρ(x)dx

Rb
Distancia recorrida en el intervalo de tiempo [a, b] por un móvil de velocidad v(t) : D = a v(t)dt .

114