Problemas de Integración Múltiple
Problemas de Integración Múltiple
PRÓLOGO
La presente publicación contiene una colección de problemas, resueltos en
su mayor parte, sobre integración múltiple, e integrales impropias.
Muchos de los problemas resueltos han sido propuestos en diversos exámenes
de la asignatura Ampliación de Cálculo de segundo curso (plan 76) de esta
Escuela de Ingenieros Industriales, …gurando en el enunciado el curso académico
al que corresponden.
Las soluciones se presentan en detalle y a muchas de ellas se han añadido
explicaciones adicionales, con el …n de completar el entendimiento del problema
y sus relaciones con la teoría correspondiente.
Uno de los objetivos principales que hemos pretendido con la redacción de
este manual es mostrar a los alumnos cómo debe presentarse la solución de un
problema. La exposición de cualquier trabajo, y más aún la presentación de
la solución de un problema de matemáticas, debe hacerse de manera clara y
ordenada, manteniendo a lo largo del escrito la coherencia de todo lo que se
escribe, que es uno de los índices más signi…cativos del grado de comprensión
de la materia objeto del problema.
Desde cualquier punto de vista, es preferible no escribir nada si no se entiende
nada de un enunciado, que presentar un ejercicio con el que se pretende dar la
sensación, mediante unas páginas caóticamente emborronadas de fórmulas y
frases incoherentes, de que algo se sabe. Con el ánimo de reforzar el valor de
esta opinión, recogemos aquí las palabras de Serge Lang, distinguido matemático
y profesor de reconocido prestigio:
No se suele reconocer que algunas di…cultades de la enseñanza de las matemáti-
cas son análogas a las de la enseñanza de una lengua extranjera. Me parece
esencial que los estudiantes tengan que escribir sus ejercicios de matemáticas
en frases completas y coherentes. Una gran parte de sus di…cultades con la
matemática surge de su caótico uso de los símbolos y fórmulas matemáticas,
aislados de toda proposición signi…cativa y todo cuanti…cador apropiado. Así
pues, todos los trabajos que se presenten deben ser legibles y estar bien ordena-
dos. No deben sugerir que una mosca recién salida del tintero se haya arrastrado
por el papel.
No se nos oculta el hecho de que, para un alumno principiante, la lectura de
los problemas de la colección que aquí se presenta puede provocar el desánimo
y la sensación de impotencia ante su resolución. Las matemáticas son difíciles
y no existe un camino cómodo para su aprendizaje. Continuamente, cuando
se explica la resolución de un problema matemático, plantean los estudiantes
la cuestión de cómo se les puede ocurrir a ellos la idea precisa que conduce al
método adecuado. Su sentimiento es que el único recurso es la pura imitación
de la solución de un problema análogo, sin ninguna con…anza en lo que su
propia inteligencia les pueda sugerir como camino fecundo. Esta es una actitud
equivocada, ya que el esfuerzo personal de meditación, re‡exión y estudio que
debe hacer el alumno, no puede, a nuestro juicio, suplirse de ninguna manera.
Por ello aconsejamos que el estudiante intente con perseverancia, antes de leer
la solución que aquí se propone, elaborar su propia solución, en la seguridad
de que el tiempo que a ello dedique siempre estará bien empleado, aún cuando
no consiga completar del todo su tarea. Solo después de este trabajo previo,
debe afrontarse la lectura de la solución propuesta, para que sirva como guía en
las partes que no se hayan podido resolver o como comprobación de resultados.
2
Debe también tener claro el alumno que las soluciones que aquí se proponen no
son las únicas, ni necesariamente las mejores y que es posible que él encuentre
por sí mismo otras que le resulten más satisfactorias.
Para …nalizar, queremos agradecer por anticipado cualquier indicación que
se nos haga sobre el contenido de esta publicación, que sirva para su mejora en
posibles futuras ediciones.
3
INTEGRAL MÚLTIPLE
TEORÍA GENERAL
1.1. Aplicando la de…nición, demostrar que la función de…nida por f (x; y) = xy,
es integrable en el cuadrado [0; 1] [0; 1] y calcular el valor de su integral sobre
ese cuadrado.
8x 2 Q : f = f + f ; f + (x) 0 ; f (x) 0:
INTEGRACIÓN REITERADA
1.4. Un teorema de derivación bajo el signo de integración.
Sea f : [a; b] [c; d] ! R una función real con derivadas parciales continuas.
Se de…ne la función g : [a; b] ! R mediante:
Z d
g(x) = f (x; y) dy:
c
1.5. Sea la funcion ' real de variable real de…nida sobre R f 1; 1g por:
Z
'(a) = log(1 2a cos x + a2 )dx
0
para y 2 2; 2 :
1.14. Calcular la integral de la función f (x; y) = xy, sobre la región del primer
cuadrante limitada por las parábolas y = x2 y x = y 2 .
T = (x; y) 2 R2 ; x 0; y 0; x + y 1 :
T = (x; y; z) 2 R3 ; x 0; y 0; z 0; x + y + z 1 :
1.17. Sea T el triángulo del plano limitado por la recta x + y = 1 y los dos ejes
de coordenadas y sea f : [0; 1] ! R una función integrable en [0; 1]. Comprobar,
para todo par p y q de enteros no negativos, la validez de la fórmula
Z Z Z 1
p!q!
xp y q f (x + y) dxdy = tp+q+1 f (t) dt:
T (p + q + 1)! 0
Obtener una fórmula semejante para integrales triples.
6
1.19. Calcular el valor medio de las distancias entre dos puntos elegidos ar-
bitrariamente de un segmento.
1.20. Calcular el valor medio de las áreas de todos los triángulos inscritos
en una circunferencia.
DESIGUALDAD DE CAUCHY-SCHWARZ
1.24. Sean f; g : A ! R dos funciones integrables sobre el conjunto acotado
A medible-Jordan. Demostrar la desigualdad de Cauchy-Schwarz:
Z 2 Z Z
fg f2 g2
A A A
1.29. Calcular
Z Z Z
cos(ax + by + cz) dx dydz ; a; b; c 2 R:
x2 +y 2 +z 2 1
8
Sea A el cuadrado [0; 1] [0; 1]. Denotemos por L(A) al paralelogramo imagen
de este cuadrado por la transformación lineal L. Demostrar que el área de este
paralelogramo es j det Lj.
Nota.–Este ejercicio es la comprobación, en un caso particular sencillo, del
siguiente resultado general.
Sea L : Rn ! Rn una aplicación lineal y A un subconjunto de Rn aco-
tado y medible–Jordan. Denotemos por L(A) al conjunto imagen de A bajo la
aplicación L. Se puede demostrar que siempre se veri…ca:
siendo
M = f(x; y) 2 R2 ; x2 + y 2 1; x 0; y 0g:
siendo
D = f(x; y) 2 R2 ; x2 + y 2 1; x 0; y 0g:
x+y+z =1
siendo
D = f(x; y; z) 2 R3 ; x2 + y 2 + z 2 R2 ; z 0g:
y se pide:
a) Expresando la suma x + y + z en términos de u, v y w, determinar la
imagen por la aplicación T del cubo unidad Q de R3 , dado por:
Q = f(u; v; w) 2 R3 ; 0 u 1; 0 v 1; 0 w 1g
(ETSII 1993-94).
0 a r ; 0 b r ; 0 c r
y r = xi + yj + zk.
(ETSII 1996-97).
Se pide:
a) Expresar In en términos de In 2 : Obtener In en términos de n:
b) Para cada entero positivo n; calcular el producto In In 1 :
c) Como aplicación, obtener la fórmula de Wallis:
2
(n!) 22n p
lim p =
n !1 2n! n
de…nida por:
x2 + y 2 a2 ; 0 z x2 y 2 + a2 ; (a > 0):
4x2 + (y 6)2 = 4 ; z = 0:
x2 + y 2 = ax ; x2 + y 2 = bx ; x2 + y 2 = cy ; x2 + y 2 = dy
x2 + y 2 = 1 ; (x 1)2 + (y 1)2 = 1:
13
CENTROS DE GRAVEDAD.
1.69. Supóngase que un sólido M de masa m y centro de gravedad G está
formado por dos piezas M1 y M2 de masas m1 y m2 y centros de gravedad G1 ,
G2 respectivamente.
Demostrar que
1
G = (m1 G1 + m2 G2 ):
m
Determinar el centro de gravedad del sólido homogéneo de…nido por las condi-
ciones:
x2 + y 2 + z 2 4a2 ; (x a)2 + y 2 + z 2 a2 :
1.71. Calcular las coordenadas del centroide de la porción del plano interior a
la elipse
x2 y2
+ =1
a2 b2
y situada en el primer cuadrante.
15
S = f(x; y; z) 2 R3 ; x2 + y 2 + z 2 2 ; x2 + y 2 zg:
1.75. Determinar las coordenadas del centro de gravedad del sólido limitado
por el paraboloide z = x2 + y 2 y el cono z 2 = x2 + y 2 , cuya densidad en cada
punto es proporcional al cuadrado de la distancia de dicho punto al eje OZ.
1.76. Determinar la posición del centroide del sólido homogéneo de…nido por
x2 + y 2 + z 2 2 ; x2 + y 2 1 ; z 0:
1.77. Determinar la posición del centroide del sólido homogéneo de…nido por
x2 + y 2 z 4:
1.78. Se considera el sólido homogéneo limitado por las super…cies cuyas ecua-
ciones son:
z = x2 + y 2 ; z = 2ax + 2y + 2:
Determinar el valor de a para que la distancia de su centro de gravedad al eje
OZ sea 2.
(ETSII 1987-88).
1.82. En una esfera maciza de radio a existe una distribución de masa cuya
densidad en cada punto es proporcional a la distancia de dicho punto a uno
…jo de la super…cie de la esfera. Se pide determinar la posición del centro de
gravedad de este sólido.
(ETSII 1990-91).
1.84. Determinar el centro de masa del sólido limitado por dos esferas con-
céntricas de radios 0 < a < b y un plano que pasa por el centro de ambas,
sabiendo que la densidad en cada punto es proporcional a la distancia de dicho
punto al centro.
(ETSII 1991-92).
1.85. Se considera el cono formado por las semirrectas que parten del punto
(0; 0; 10) y que se apoyan en la circunferencia situada en el plano 3x 4z = 0,
uno de cuyos diámetros tiene por extremos los puntos (0; 0; 0) y (8; 0; 6). Hallar
el centro de gravedad del sólido homogéneo limitado por dicho cono y el plano
z = 0.
(ETSII 1996-97).
17
p
1.86. Una circunferencia de radio a gira sobre el exterior de otra de radio a
q
p
donde es irreducible. Se pide:
q
a) Encontrar unas ecuaciones paramétricas de la curva descrita por un punto
de la circunferencia de radio a (epicicloide).
b) En el caso p = q = 1 situando adecuadamente el origen, obtener unas ecua-
ciones de dicha curva en coordenadas polares (cardiode).
c) Calcular el centro de gravedad de la super…cie encerrada por la curva del
apartado b.
a2 x2 + y 2 + z 2 b2 ; z 0:
x2 + y 2 + z 2 = a2 ; z 0:
(ETSII 1997-98)
MOMENTOS DE INERCIA.
1.88. Teorema de Steiner. El momento de inercia de un sólido M respecto de
una recta r es igual al momento de inercia de M respecto de una recta paralela
a r y que pasa por su centro de gravedad, más el producto de la masa de M por
el cuadrado de la distancia que separa las dos rectas.
1.89. En una esfera de radio R existe una distribución de masa con densi-
dad en cada punto proporcional a la distancia de ese punto al centro de la
esfera, con constante de proporcionalidad k. Se considera el sólido formado por
la parte de la esfera interior al cono de vértice el centro de la esfera y semiángulo
en el vértice igual a . Se pide:
a) Calcular la masa de ese sólido.
b) Determinar el valor de para que el centro de gravedad del sólido esté
a distancia R=2 del centro de la esfera.
c) Calcular el momento de inercia del sólido respecto de su eje de simetría.
(ETSII 1988-89).
1.90. En un hexágono regular cuyo lado tiene una longitud l hay una dis-
tribución de masa. La densidad en cada punto es inversamente proporcional
a la distancia de ese punto al centro del hexágono, siendo k la constante de
proporcionalidad. Se pide:
a) Calcular la masa total del hexágono.
b) Calcular el momento de inercia del hexágono respecto de su centro.
c) Determinar la posición del centro de masa de la mitad del hexágono,
suponiéndolo dividido por una de sus diagonales.
(ETSII 1989-90).
18
1.94. Sobre el triángulo isósceles que tiene por vértices los puntos de coor-
denadas ( a; 0), (a; 0) y (0; h) existe una distribución super…cial de masa de
densidad constante igual a . Se pide:
(a) Determinar los momentos de inercia de esta distribución respecto de cada
uno de los dos ejes de coordenadas.
(b) Determinar el momento de inercia de la distribución respecto de una
recta cualquiera que pasa por el origen, expresándolo en términos del ángulo
polar de esta recta.
(c) Calcular el área que encierra la curva cuya ecuación en polares es
1
=p
I( )
1.95. Calcular el momento de inercia respecto del eje de giro del toro sólido
homogéneo engendrado por un círculo de radio a al girar alrededor de un eje
situado en su mismo plano y que dista una distancia d > a de su centro.
(ETSII 1995-96).
SÓLIDOS DE REVOLUCIÓN
1.98. Sea el dominio rayado en la …gura adjunta, donde los lados curvos
son arcos de circunferencia de radios 1 y 2. Determinar la posición del centroide
de . Calcular el volumen engendrado por al girar alrededor de OX y al girar
alrededor de OY .
y=x
1.5
Y 1
0.5
Ω
0
0 0.5 1 1.5 2
X
Calcular el volumen del sólido limitado por la región como base y limitado
superiormente por la super…cie esférica centrada en el origen y de radio 2.
1.104. a) Encontrar una fórmula para calcular el centroide del sólido obtenido
al girar un recinto plano alrededor de un eje no secante y situado en su mismo
plano.
b) Idem para calcular cualquier momento sobre el sólido es decir
Z Z Z
xl y m z n dxdydz
V
= f(x; y) 2 R2 ; xy 4 ; jx yj 3; x 0; y 0g:
= 2a(1 + cos )
x = y = z:
22
en el supuesto de que la integral anterior sea convergente para cada valor del
parámetro real x.
Este problema tiene por objeto determinar la convolución de dos campanas
de Gauss.
a) Calcular para cada valor real del parámetro p la integral
Z +1
2
I(p) = e t +pt dt
1
b) Dados los números reales y , considérense las dos funciones ' y de…nidas
a partir de f y g mediante
Calcular I( ) para todos los valores del parámetro real para los que la in-
tegral es convergente y especi…car dichos valores de . Estudiar la continuidad
de I( ) en = 0.
INTEGRALES EULERIANAS
2.6. Funciones Gamma y Beta de Euler. Las funciones Gamma y Beta
de Euler se de…nen mediante las integrales
Z +1 Z 1
x 1 t
(x) = t e dt ; (x; y) = tx 1 (1 t)y 1 dt
0 0
a) Demostrar que (x) está de…nida para x > 0. Demostrar (x; y) está
de…nida para x > 0, y > 0.
b) Demostrar que (p + 1) = p (p), p > 0. Demostrar que (n) = (n 1)!,
n = 1; 2; : : :.
c) Esbozar razonadamente la grá…ca de (x).
d) Demostrar que
Z +1
tx 1
(x; y) = dt
0 (1 + t)x+y
x1 + x2 + + xn 1 ; x1 0 ; x2 0 ; : : : ; xn 0
(ETSII 1987-88).
y calcularla si es convergente.
(ETSII 1987-88).
(ETSII 1991-92).
(ETSII 1993-94).
2.18. Se trata de encontrar una función f (t) continua, real, de variable real,
de…nida para t > 0, tal que satisface la siguiente ecuación integral:
Z x
f (t) p
p dt = x x ( )
0 x t
Se pide:
a) Calcular Z y p 1
G(y) = x x p dx
0 y x
b) Dibujar el recinto de integración e invertir el orden de integración en la
siguiente integral reiterada:
Z y Z x
f (t) 1
F (y) = p dt p dx
0 0 x t y x
c) Calcular la integral simple:
Z y
1
p dx ; t y
t (x t) (y x)
d) Transformar la ecuacion ( ) en otra ecuación que relacione F (y) con G(y).
Obtener la función f (x).
(ETSII 1998-99).
(ETSII 1999-2000).
(p) (1 p) = ; 0<p<1
senp
y despues integrar entre 0 y 1=2:
4) Para cada > 0, calcular
Z +1
F( ) = log( (x))dx
Indicación: Calcular F 0 ( ):
(ETSII 1999-2000).
a) Demostrar que:
1
F 0( ) = F( )
2
b) Obtener F ( ).
c) Como aplicación, calcular la integral:
Z +1
2 sen 2ax
e x dx
0 x
R +1
c) Probar que 0 g(a; x) dx es uniformemente convergente para a 2 [ ; +1)
con > 0:
d) Obtener la derivada F 0 (a) para a > 0:
e) Obtener F (a) para a > 0:
f) Probar que F (a) es uniformemente convergente para a 0:
g) Calcular
Z +1
sen x
F (0) = dx :
0 x
Z +1
sen x
2.24. Sabiendo que dx = , calcular:
0 x 2
Z +1 Z +1
senxy sen2 x
1) dx ; 2) 2
dx
Z0 +1 x Z0 +1 x 4
1 cos x sen x
3) 2
dx ; 4) dx
Z0 +1 x Z 0 x2
+1
sen4 x sen3 x
5) 4
dx ; 6) dx
0 x 0 x3
Z +1
sen x
2.25. Sabiendo que dx = , calcular:
0 x 2
Z +1 Z +1
cos ax sen ax
1) 2
dx ; 2) dx
0 1+x 0 x(1 + x2 )
2.26. Estudiar la convergencia y calcular cuando sea posible cada una de las
siguentes integrales impropias
Z +1 Z +1
1) e ax cos bxdx ; 2) e ax senbxdx
0 0
Z +1 Z +1
3) xe ax cos bxdx ; 4) xe ax senbxdx
0 0
Z +1 Z +1
5) x2 e ax cos bxdx ; 6) x2 e ax senbxdx
0 0
2.27. Este problema tiene por objeto desarrollar una demostración de la cono-
cida fórmula Z +1
2 p
e t dt = ( )
1
2.28. Se considera la función real de dos variables reales dada por la representación
integral: Z +1
arctan(xt) arctan(yt)
f (x; y) = dt
0 t2
y se pide:
a) Determinar el dominio de la función, es decir, el conjunto de parejas de
números reales (x; y) para las cuales es convergente la integral anterior.
b) Admitiendo la posibilidad de derivar dentro de la integral cuando xy 6= 0,
@2f
calcular la expresión explícita de la función .
@x@y
c) Determinar explícitamente el valor de f (x; y) para cada punto (x; y) del
dominio de f .
(ETSII 1990-91).
2.30. 1) Calcular
Z +1
dt
dt:
0 (x2 + t2 )
2) Calcular
Z +1
dt
dt ; n = 1; 2; :::::::
0 (x2 + t2 )n
Indicación: Derivar la integral respecto de un parámetro y razonar por induc-
ción.
3) Calcular
Z +1
dt
2 n
dt
0 1 + tn
4) Como aplicación de lo anterior, calcular
1 3 5 (2n 3) p
lim n
n!1 2 4 6 (2n 2)
(ETSII 1997-98).
30
2.33. Calcular Z Z
2 2
xye (x +y ) dx dy
M
donde
M = (x; y) 2 R2 ; 0 x 2y
2.35. a) Determinar los valores del parámetro real para los que existe la
integral triple impropia
Z Z Z
dxdydz
3 1 + (x2 + y2 + z2 )
R
2.36. Estudiar la existencia de la integral según los valores del parámetro real
a 0. En los casos en que sea posible, calcularla.
Z Z Z
dxdydz
I(a) = 2 + y 2 + (z
2 2
x +y +z 12 x a)2
(ETSII 1990-91).
2.37. Estudiar la existencia de la integral según los valores del parámetro real
a 0. En los casos en que sea posible, calcularla.
Z Z Z
dxdydz
I(a) = p
2 2 2
x +y +z 1 x + y 2 + (z a)2
2
Calcular la integral
Z Z Z
dxdydz
q
x2 +y 2 +z 2 1 2 2 2
(x a1 ) + (y a2 ) + (z a3 )
Interpretar físicamente.
(ETSII 1991-92).
32
siendo
M = (x; y) 2 R2 ; 3x2 + 3y 2 + 2xy 1
Especi…car los valores de para los cuales I( ) es convergente.
(ETSII 1997-98).
donde
M = x 2 Rn : (xt Qx) 1
siendo Q una matriz real cuadrada n n, simétrica y de…nida positiva.
2.42. Sea f una función real de variable real y de clase C 1 (R), tal que
Z +1
f (u)
du es convergente
1 u
PRIMERA PARTE
SEGUNDA PARTE
INTEGRALES MÚLTIPLES
1.1. Para cada entero positivo n consideramos la partición Pn del cuadrado
[0; 1] [0; 1] en los n2 subcuadrados que resultan al dividir cada intervalo [0; 1]
en n subintervalos iguales. La …gura adjunta muestra esta partición cuando se
toma n = 5.
1
Y
0.8
p'ij
j 0.6
R(i,j)
j-1 0.4
p''
ij
0.2
0
0 0.2 0.4 0.6 0.8 1
i-1 i
X
Utilizando la notación anterior, el subcuadrado Rij tiene como vértices los pun-
tos de coordenadas:
i 1 j 1 i j 1 i 1 j i j
; ; ; ; ; ; ; :
n n n n n n n n
ij (i 1)(j 1)
Mij = ; mij = :
n2 n2
Para esta partición la correspondiente suma superior de Darboux es:
n X
X n Xn X n n n
ij 1 1 X X
Sf (Pn ) = Mij A (Rij ) = = i j
i=1 j=1 i=1 j=1
n2 n2 n4 i=1 j=1
2
1 n(n + 1) 1 1 1
= = + + :
n4 2 4 2n 4n2
De igual modo se calcula la suma inferior de Darboux para esa partición que
resulta ser:
1 1 1
If (Pn ) = + 2:
4 2n 4n
Observemos que:
1
Sf (Pn ) If (Pn ) =
n
35
1 1 1 1
inf Sf (Pn ) = lim + + = :
n n!1 4 2n 4n2 4
Dado que estas sumas superiores son sólo una parte de las sumas superiores
posibles se cumple que:
1
inf Sf (P ) inf Sf (Pn ) =
P n 4
Por otra parte no es difícil comprobar que la sucesión de sumas inferiores If (Pn )
es monótona creciente. Para ello basta observar que la función:
1 1 1
'(t) = +
4 2t 4t2
tiene por derivada:
1 1 t 1
'0 (t) ==
2t2 2t3 2t3
que es positiva para t 1, de manera que '(t) es creciente en el intervalo [1; 1).
Por consiguiente se cumple:
1 1 1 1
sup If (Pn ) = lim + 2 = :
n n!1 4 2n 4n 4
Por idéntica razón que en el caso de las sumas superiores, se cumple ahora que:
1
sup If (P ) sup If (Pn ) = :
P n 4
La cadena de desigualdades:
1 1
sup If (P ) inf Sf (P )
4 P P 4
fuerza a que los dos términos intermedios sean iguales a 1=4, es decir que:
ZZ
1
xy dxdy = :
[0;1] [0;1] 4
Por lo tanto:
Si jaj < 1 :
'(a) = cte = '(0) = 0
Si jaj > 1; se tiene que:
'(a) = 2 log jaj + cte = log a2 + cte
Con el …n de determinar la constante, obsérvese que si jaj > 1 entonces j1=aj < 1
y por lo tanto:
Z
2 1
0 = '(1=a) = log(1 cos x + 2 )dx
0 a a
Z
= log(a2 2a cos x + 1) log a2 dx
0
= '(a) log a2
Se deduce …nalmente que:
'(a) = log a2
Por lo tanto:
Z y
1
G(y) = G(0) + ytgy + log 1 + tg2 y dy
0 2
Z y
= y log(cos y) + log(cos y) + log(cos y)dy
0
= y log(cos y)
x2 y2
+ = 1:
9 4
Invirtiendo el orden de integración se obtiene
Z Z Z x=3 Z y=2p1 (x2 =9) !
y y
dxdy = dy dx
D x+3 x=0 y=0 x+3
p
Z 3 y=2 1 (x2 =9) Z
dx y2 2 3
= = (3 x) dx = 1:
0 x+3 2 y=0 9 0
x+y =1 ; x+y = 1 ; x
y = 1 ; x y = 1:
Z Z Z 0 Z 1+x Z 1 Z 1 x
2 2 2
(x+y) dxdy = (x + y) dy dx+ (x + y)2 dy dx = :
M 1 1 x 0 x 1 3
Nota.- El aspecto geométrico del recinto de integración sugiere efectuar el giro
de ejes coordenados:
" #
u p1 p1 x
= p1 2 2
v p1 y
2 2
1.11.
Z Z p Z 1 Z 1 p
jy x2 jdxdy = dx jy x2 jdy
[0;1] [0;1] 0 0
Z 1 Z x2 p Z 1 Z 1 p
= dx x2 ydy + dx y x2 dy
0 0 0 x2
Z 1 Z 1
2 3 2 3 1
= x dx + 1 x2 2
dx = + :
0 3 0 3 6 8
1.12. Si f (x; y) es una función continua, las dos posibles maneras de emplear
la integración reiterada sobre el triángulo
conducen a la igualdad
Z a Z x Z a Z a
f (x; y) dy dx = f (x; y) dx dy:
0 0 0 y
1
1.13. Respuesta: e .
e
1.14. Por simetrías de la función y del recinto con respecto a la diagonal
principal Z 1 Z x Z 1
1
I=2 dx xydy = x3 x5 dx = :
0 x2 0 12
Nótese que para p = q = 0 resulta el valor 1=2 que coincide con el área del
triángulo.
42
t=x+y x=x
Z Z Z 1 Z t
xp y q f (x + y) dxdy = xp (t x)q f (t)dx dt:
T 0 0
El cambio x = tu en la integral simple proporciona:
Z 1 Z 1 Z 1 Z 1
tp up tq (1 u)q f (t)tdu dt = tp+q+1 f (t)dt up (1 u)q du :
0 0 0 0
f ( 1; 2) = f ( 2; 1)
sen
g( ) =
1 cos
es una función decreciente en [0; 2 ] ; será g( 1 ) g( 2 ), lo que permite concluir
que en tal región, f ( 1 ; 2 ) 0: Asi pues:
Z 2 Z 1
r2
M = d 1 (sen ( 1 2) sen 1 + sen 2 ) d 2
4 2 0 0
Z 2
r2
= (2 2 cos 1 1 sen 1 ) d 1
4 2 0
Z 2
r2 r2 3r2
= 2
(2 1 sen 1 ) d 1 = 2
6 = :
4 0 4 2
1.22. Respuesta: :
96
1.23. Respuesta: :
35
1.24. La demostración es, salvo en un detalle, idéntica formalmente a la
demostración de la desigualdad de Schwarz para n–uplas de números reales.
Consideremos la función real de variable real de…nida mediante:
Z Z Z Z
I( ) = ( f + g)2 = 2 f2 + 2 fg + g2 :
A A A A
Z
Salvo cuando f 2 = 0 la función I( ) es un trinomio de segundo grado
A
que además es no negativo, ya que surge de integrar una función no negativa.
Supongamos por el momento que
Z
f 2 > 0:
A
45
y sigue siendo una función no negativa para todo real. Puesto que se trata
de un polinomio en de grado menor o igual que uno, esto solo es posible si el
coe…ciente del término en se anula. Es decir,
Z
fg = 0
A
f= min
kL2
C :
T
Por otra parte si se elige una velocidad constante, que debe ser la velocidad
media L=T , el consumo es:
Z T 2
L kL2
C=k dt =
0 T T
y por lo tanto, la velocidad constante proporciona el consumo mínimo. Ob-
sérvese que si únicamente se admiten regímenes de velocidades que varíen con
continuidad, esta ley de velocidades es la única que proporciona el mínimo con-
sumo.
1.28.
Y y=x
y=ax a>1
yx=b
y=ax 0<a<1
Haciendo el cambio
y
u=
; v = xy
x
transformamos el recinto D en el rectángulo [a; 1] [0; b], de forma inyectiva en
el interior del rectángulo ya que:
r
v p
x= ; y = uv:
u
b ja 1j b2 jlog aj
I1 = ; I2 = :
2 4
Obsérvese que al cambiar a por 1=a; la segunda integral no se modi…ca.
Esto se debe a la simetría del recinto con respecto a la diagonal principal y a la
simetría de la función con respecto a la misma diagonal.
48
1.29. Denotemos por a al vector (a; b; c). Efectuaremos un giro de los ejes
cartesianos que haga coincidir uno de los ejes con la recta soporte del vector a.
Esto signi…ca que hemos cambiado las variables (x; y; z) por las nuevas variables
(u; v; w) en la forma
2 3 2 3
x u
4 y 5 = P 4 v 5 con u = ax + by + cz :
kak
z w
La matriz P es una matriz ortogonal y por lo tanto, el valor absoluto del deter-
minante jacobiano del cambio de variables es igual a jdet P j = 1. Nótese además
que el recinto de integración no ha variado en esta transformación geométrica,
por lo que la integral que deseamos calcular es:
Z Z Z Z Z Z
I= cos(ax+by+cz) dxdydz = cos (ukak) dudvdw
x2 +y 2 +z 2 1 u2 +v 2 +w2 1
1.30. Nótese que con un tolerable abuso de notación, denotamos con el mismo
símbolo a la aplicación lineal y a su matriz asociada. Sea A el cuadrado
[0; 1] [0; 1]. Según sabemos, la imagen L(A) de este cuadrado por la transfor-
mación lineal L es el paralelogramo determinado por los vectores columna de la
matriz de la transformación, es decir los vectores
u = ai + cj v = bi + dj
La …gura adjunta representa el cuadrado y su imagen.
(a+b,c+d)
Y V
1 (b,d)
(a,c)
O 1 X O U
49
Por consiguiente:
m(L(A)) = jdet Lj
Como el área de A es obviamente 1, se observa la validez de la proposición en
este caso especial.
a4
1.41. Respuesta: sh1:
8
1.42. a) Dado que los valores de las variables u, v, w se situan en el intervalo
[0; 1], es claro que los valores de las variables x, y, z son números no negativos.
Dicho en otros términos, la imagen del cubo unidad está situada en el octante
positivo del espacio R3 . Por otra parte se tiene:
x+y+z = u + (1 u)v + (1 u)(1 v)w
= [1 (1 u)] + (1 u)v + (1 u)(1 v)w
= 1 (1 u)(1 v) + (1 u)(1 v)w
= 1 (1 u)(1 v)(1 w):
La misma consideración anterior sobre las restricciones de los valores de u, v, w
permite asegurar que el producto (1 u)(1 v)(1 w) es positivo y menor que
1, con lo cual:
0 x + y + z 1:
La imagen del cubo unidad está por consiguiente contenida en el tetraedro
S = f(x; y; z) 2 R3 ; x 0; y 0; z 0 ; x+y+z 1g:
Nótese además que la desigualdad anterior es estricta a menos que el punto
(u; v; w) sea un punto del interior del cubo, lo cual quiere decir que los puntos
del interior del cubo se transforman en puntos del interior del tetraedro.
La …gura adjunta representa grá…camente el cubo unidad Q y el tetraedro
S.
D
A' Z
W
B B'
C'
C
V U Y X
51
Esta inclusión permite a…rmar que ningún punto del interior del tetraedro es
un punto frontera del conjunto T (Q), y por consiguiente que todos los puntos
del interior del tetraedro son imagen de algún punto del cubo unidad por la
transformación T . Por consiguiente la imagen del cubo unidad Q es exactamente
el tetraedro S.
b) Es sencillo calcular el jacobiano de la transformación. Se tiene:
2 3
1 0 0
JT (u; v; w) = det 4 v 1 u 0 5
(1 v)w (1 u)w (1 u)(1 v)
= (1 u)2 (1 v):
Para un punto (x; y; z) del interior del tetraedro S, los valores de u, v, w quedan
unívocamente determinados por los de x, y, z, de manera que puede asegurarse
que la transformación es inyectiva en el interior del cubo Q. Ya hemos visto en
52
1.43. Efectuamos un giro de ejes, de…nido por una matriz constante P or-
togonal, que haga coincidir uno de los nuevos ejes con la recta soporte del
vector …jo dado a = (a; b), tal y como se indica en la …gura.
U
(a,b)
O X
53
Obsérvese que jdet Aj es el volumen del paralelepípedo cuyas aristas son los vec-
tores a, b, c.
es decir
n 1
In = In 2:
n
Por lo tanto
n 1 n 1n 3
In = In 2 = In 4:
n n n 2
Si n es par
n 1n 3 31 n! n!
In = ::::: I0 = 2 = :
n n 2 42 (n!!) n 2
2! 2n
Si n es impar
n 1 2
n 1n 3 42 (n 1)!! 2 ! 2n
In = ::::: I1 = 2= :
n n 2 53 n! n!
b)
n 1n 2 n 2
In In 1 = In 2 In 3 = In 2 In 3 :
n n 1 n
Si n es par
Z Z
n 2n 4 2 2 2 2 2
In In 1 = ::::: I2 I1 = sen 'd' sen'd' = 2= :
n n 2 4 n 0 0 n2 n
55
Si n es impar
Z Z
n 2n 4 1 1 1 2
In In 1 = ::::: I1 I0 = sen'd' d' = 2 = :
n n 2 3 n 0 0 n n
c) Es evidente que,
0 In+1 In n2N
y por lo tanto fIn g es una sucesión decreciente y acotada inferiormente, lo que
asegura que es convergente. Del resultado
2
In In 1 =
n
se obtiene
lim In = 0
n !1
y como
n
2 = nIn+1 In nIn2 nIn In 1 =2
n+1
cuando n tiende a in…nito obtenemos
lim nIn2 = 2
n !1
d) Si n es par
n
2 2 2 2 2
In In 1 In 2 In 3 :::::I2 I1 = ::::: = n :
n n 2 4 2 2!
Si n es impar
n+1
n 1
2 2 2 2 2
2 !2n
In In 1 In 2 In 3 :::::I1 I0 = ::::: = :
n n 2 3 1 n!
xn
cos '1 =
relaciones que determinan '1 de manera única si el sen '1 es positivo, es decir
en (0; ) : Este razonamiento es aplicable a '2 ; ::::; 'n 2 : Para 'n 1 se tiene:
x1
sen'n 1 =p
x21 + x22
x2
cos 'n 1 =p
x21 + x22
que determina 'n 1 de manera única en (0; 2 ).
b)
sen'1 cos '1
J 2 = =
cos '1 sen'1
sen'1 ::sen'n 1 cos '1 ::sen'n 1 ::: sen'1::: cos 'n 2 sen'n 1 sen'1 ::: cos 'n 1
sen'1 :: cos 'n 1 cos '1 :: cos 'n 1 ::: sen'1 ::: cos 'n 2 cos 'n 1 sen'1 :::sen'n 1
J n = ::::: ::::: ::: ::::: ::::: :
sen'1 cos '2 cos '1 cos '2 ::: 0 0
cos '1 sen'1 ::: 0 0
sen'1 ::sen2 'n 1 cos '1 ::sen2 'n 1 ::: sen'1 :: cos 'n 2
2 sen 'n 1 sen'1 ::sen'n 1 cos 'n 1
sen'1 :: cos2 'n 1 cos '1 :: cos2 'n 1 ::: sen'1 :: cos 'n 2
2 cos 'n 1 sen'1 ::sen'n 1 cos 'n 1
::::: ::::: ::: ::::: :::::
sen'1 cos '2 cos '1 cos '2 ::: 0 0
cos '1 sen'1 ::: 0 0
:
sen'n 1 cos 'n 1
sen'1 ::sen2 'n 1 cos '1 ::sen2 'n 1 ::: sen'1 :: cos 'n 2 sen2 'n 1 sen'1 ::sen'n 1 cos 'n 1
sen'1 ::sen'n 2 cos '1 ::sen'n 2 ::: sen'1 ::: cos 'n 2 0
::::: ::::: ::: ::::: :::::
sen'1 cos '2 cos '1 cos '2 ::: 0 0
cos '1 sen'1 ::: 0 0
:
sen'n 1 cos 'n 1
57
( 1)n 1
sen'1 ::::sen'n 2
sen'1 :::sen'n 2 cos '1 :::sen'n 2 ::: sen'1 ::: cos 'n 3 sen'n 2 sen'1 ::: cos 'n 2
sen'1 ::: cos 'n 2 cos '1 ::: cos 'n 2 ::: sen'1 ::: cos 'n 3 cos 'n 2 sen'1 :::sen'n 2
::::: ::::: ::: ::::: :::::
sen'1 cos '2 cos '1 cos '2 ::: 0 0
cos '1 sen'1 ::: 0 0
Si n es impar,
n 1 n 1
n 3
2rn 2
2 !2n 2 2 ( n 2 1 )!2n n
V (En ) = = r :
n n 2! n!
d) Consideremos los recintos
x1 + x2 + :::::::: + xn < 1
Dn
xk > 0 k = 1; 2; ::::; n
Tn (0; 1) 0; :::::: 0;
2 2
n
Cn (0; 1)
58
y las transformaciones:
(x1 ; x2 ; ::::::::; xn ) = n (u1 ; u2 ; ::::::::; un ) = u21 ; u22 ; ::::::::; u2n
1 = x1 + x2 + :::::::: + xn 1 + xn
1 2 = x1 + x2 + :::::::: + xn 1
::::::: ::::::::
1 2 :::::::: n 1 = x1 + x2
1 2 :::::::: n 1 n = x1
e) El jacobiano de n vendrá dado por:
1
J n = [J n (u1 ; u2 ; ::::::::; un )] J n ; '1 ; ::::::::; 'n 1 J n ( 1; 2 ; ::::::::; n ) =
[J n (u1 ; u2 ; ::::::::; un )] J n ; '1 ; ::::::::; 'n 1
= =
J n ; '1 ; ::::::::; 'n 1
2n u1 u2 ::::::::un (cos n2 + sen n2 ) n 1
senn 2 '1 ::::::sen'n 2
= =
2n sen'1 ::::::sen'n 2 sen'n 1 cos '1 :::::: cos 'n 1
n
senn 1 '1 :::sen2 'n 2 sen'n 1 cos '1 ::: cos 'n 1 n 1 senn 2 '1 :::sen'n 2
= =
(cos n2 + sen n2 ) sen'1 ::::::sen'n 2 sen'n 1 cos '1 :::::: cos 'n 1
n n
= (cos + sen ) 2n 2 sen2n 4 '1 ::::::sen2 'n 2 =
2 2
n n
= (cos + sen ) n1 1 n2 2 n3 3 :::: 2n 2 n 1
2 2
f) El volumen del recinto Dn es:
Z Z Z Z
n 1 n 2 n 3
::: dx1 dx2 ::::dxn = ::: 1 2 3 :::: 2n 2 n 1 d 1 :::::d n =
Dn 1 (D
n)
Z 1 Z 1 Z 1 Z 1
n 1 n 2
= 1 d 1 2 d 2 :::: n 1d n 1 d n =
0 0 0 0
1 1 1 1
= ::::: 1 = :
nn 1 2 n!
59
1.49. En la …gura hemos elegido un sistema de referencia cuyos ejes son las
tres aristas del prisma, estando la base cuyos lados miden a y b en el plano
z = 0.
a
b
Este gr…co corresponde a los valores:
a=3; b=2; p=1; q=2; r=8
60
z = px + qy + r
1.50. Supongamos que los conjuntos A y B quedan descritos por los valores de
z comprendidos en el intervalo [z0 ; z1 ]. Entonces:
Z Z Z Z z1 Z Z
vol (A) = dxdydz = dz dxdy
A z0 Az
Z z1 Z z1
= m(Az ) dz = m(Bz ) dz
z0 z0
Z z1 Z Z Z Z Z
= dz dxdy = dxdydz = vol (B):
z0 Bz B
a4
1.52. Respuesta: :
2
1.53. La super…cie x2 + y 2 z 2 = a2 es un hiperboloide cuyo eje de revolución
es el eje OZ. El volumen vendrá dado por:
Z Z Z Z Z Z Z Z Z
V = dxdydz = dxdydz + dxdydz
S C E
8
1.54. Respuesta: :
6
1.55. La ecuación de la super…cie del cono viene dada en forma paramétrica
por:
2 3 2 3 2 3 2 3 2 3
x 0 e
x e
x cos t
4 y 5 = 4 0 5 + 4 ye 5 = 4 ye 5 = 4 (6 + 2sent) 5
z 6 ze 6 6 (1 ) 6 (1 )
5
Y
3 Dxy
2
0
0 2 4 6 8
X
1 x2 + y 2 x2 + y 2
:u= v=
x y
63
1.57. Respuesta: 1:
2
1.58.- La segunda condición se puede escribir como:
La curva que limita esta región es una parábola de eje oblicuo. Efectuando el
giro: " #
u p1 p1 x
= p12 p1
2
v 2 2
y
2
W
1.5
0.5
0
-2 -1 0 1 2
-0.5
U
-1
-1.5
-2
64
Z p
2
2 u2 4
A= 2 p du = :
0 2 2 3
1.59. Teniendo en cuenta que las dos super…cies que limitan el sólido tienen
simetría de revolución con respecto al eje OZ la masa vendrá dada por
Z Z Z Z Z Z p 2 2 3 x y
M= (x2 + y 2 + z 2 ) dxdydz = dxdy (x2 + y 2 + z 2 )dz
D Dxy (x2 +y 2 )=2
Por lo tanto,
Dxy = (x; y) 2 R2 : (x2 + y 2 ) 2 :
Efectuando un cambio a coordenadas polares:
Z Z Z p 3 2
2
M= p d d ( + z 2 )dz
[0; 2] [0;2 ] 2 =2
y …nalmente
Z Z p
3 2
1 5 1 2 3 1 7
M = p 3 + (3 )2 d d
[0; 2] [0;2 ] 2 3 24
Z 2 p
1 p 1 3 1 2 1 3 108 3 77
= 2 u 3 u + (3 u) 2 u u du = :
0 2 6 4 48 30
en donde ( )
2
2 2 2 R
Dxy = (x; y) 2 R : (x + y ) :
2
Efectuando un cambio a coordenadas polares:
Z Z p
p 4 3 3 R3
VT = 2 R2 2 d d = R
[0; 2 ] [0;2 ]
R 3 2
y …nalmente p
3 R3
V = :
2
65
con
Dxy = (x; y) 2 R2 : (x2 + y 2 ) 2x ; (x2 + y 2 ) 2y :
X
Dxy
en donde
x2 y2 z2
E = f(x; y; z) 2 R3 ; 2
+ 2 + 2 1g
a b c
es el elipsoide y E + representa la parte de dicho elipsoide contenida en el octante
positivo x 0, y 0, z 0.
Denotemos Ix a la integral
Z Z Z Z Z r
2 2 x2 y2
Ix = 8k x dxdydz = 8kc x 1 2
dxdy:
E+ S+ a b2
Hemos utilizado la notación
x2 y2
S + = f(x; y) 2 R2 ; 2
+ 2 1 ; x 0 ; y 0g:
a b
Si en la integral doble anterior efectuamos el cambio de variables
x = ar cos ; y = brsen
que transforma el recinto S + en el rectángulo f( ; r) ; 0 =2 ; 0 r 1g,
resulta:
Z =2 ! Z
1 p 4 k 3
3 2
Ix = 8ka bc cos d r3 1 r2 dr = a bc:
0 0 15
Por razones de simetría, es evidente que
Z Z Z Z Z Z
4 k 3 4 k
Iy = 8k y 2 dxdydz = ab c ; Iz = 8k z 2 dxdydz = abc3 :
E + 15 E+ 15
De los cálculos anteriores se obtiene …nalmente la masa del elipsoide
4 k
M= abc(a2 + b2 + c2 ):
15
La condición M < V proporciona …nalmente los valores admisibles de la constante
k:
5
0<k< 2 :
a + b2 + c2
O
Y
X
67
x21 x2n
= (x1 ; : : : ; xn ) 2 Rn ; + + 1
a21 a2n
x21 x2
2 + 22 = 1
a1 a2
68
Con el …n de calcular el área que encierra la elipse, seccionamos dicha elipse con
la recta x2 = t. La longitud del segmento interceptado es:
s
t2
2a1 1
a22
El área de la elipse viene dada por tanto por la integral:
Z a2 s Z 1p
t2
V2 = 2a1 1 dt = 4a1 a2 1 t2 dt:
a2 a22 0
X3
X2
X1
69
Cn 1 a1 : : : an 1 1
a2n
70
Vn = 2In Cn 1 a1 : : : an
u = sen ; dv = cosn 2
sen d
resulta la igualdad:
1
In = In 2 In
n 1
que proporciona la relación de recurrencia:
n 1
In = In 2:
n
Para valores pares de n se tiene entonces que:
2k 1 2k 3 3 1
I2k = I0
2k 2k 2 4 2
y puesto que
Z =2
I0 = d =
0 2
se cumple que
(2k 1)!!
I2k = :
(2k)!! 2
Análogamente cuando n es impar se tiene:
2k 2k 2 4 2
I2k+1 = I1
2k + 1 2k 1 3 1
71
y puesto que
Z =2
I1 = cos d = 1
0
se cumple que
(2k)!!
I2k+1 = :
(2k + 1)!!
Observemos que se tiene:
(2k)!! (2k 1)!! 1
I2k+1 I2k = =
(2k + 1)!! (2k)!! 2 2k + 1 2
y también
(2k 1)!! (2k 2)!! 1
I2k I2k 1 = =
(2k)!! 2 (2k 1)!! 2k 2
Las igualdades anteriores permiten a…rmar que para todo entero positivo n se
cumple que:
In In 1 = :
2n
A partir de la observación anterior, la recurrencia entre las constantes Cn se
puede escribir del modo siguiente:
2
Cn = 2In Cn 1 = 2In (2In 1 Cn 2 ) = 4In In 1 Cn 2 = Cn 2
n
Para n par se tiene entonces:
k
2 2 2
C2k = C2 =
2k 2k 2 4 k!
Análogamente para n impar se tiene:
2 2 2 2k+1 k
C2k+1 = C3 =
2k + 1 2k 1 3 (2k + 1)!!
que son las expresiones del volumen de la hiperesfera unidad, en términos de su
dimensión.
Recordemos ahora que, mediante la aplicación de la integración reiterada,
hemos demostrado que el volumen del hiperelipsoide viene dado por:
Vn = a1 : : : an Cn
xt Qx 1
X Y
z
-2
2
0
4
2
v 6
4
8 u
6
10
8
x = xg ( ) + u cos
y = yg ( ) + usen
z = zg ( ) + v
el Jacobiano será:
cos 0 x0g ( ) usen
J = sen 0 yg0 ( ) + u cos = x0g ( )sen yg0 ( ) cos u
0 1 zg0 ( )
74
nos queda:
J = s( ) u
el módulo del Jacobiano será
jJ j = s( ) + u
siempre que
s( ) + u 0
lo cual se asegura si el sólido esta contenido en la región x 0, y 0:
El volumen ahora vendra dado por:
Z Z Z Z Z Z Z Z Z
V = dxdydx = jJ j dudvd = (s( ) + u) dudvd
V 1 (V ) 1 (V )
Z 2
Z Z Z 2
Z Z
= s( )d dudv + d ududv
0 0
y nos queda
Z 2
Z 2
Z Z Z 2
V = s( )A( )d + d ududv = s( )A( )d
0 0 0
A (R; d) = 2 (VS VC )
75
donde VS es el área del sector circular y VC es el área del triangulo. El área del
sector circular vendrá dada por:
h d
VS = R2 = R2 arccos = R2 arccos
2 R 2R
…nalmente
p ! p
2 4R2 d2 d 4R2 d2
A (R; d) = 2R arctg
d 2
V (R; d) = 2 (VS VC )
…nalmente
1 2 4R2 d2 d
V (R; d) = 2 R (2R d) = 16R3 12dR2 + d3
3 3 4 2 12
Z Z Z
n 1
VSn = senn 2
'n 1 :::sen'2 :d d'n 1 ::::d'1
[0;R] [0; ] [0; ] :::: [0;2 ]
Z
n 1
= RVEn 1
(R) senn 2
'd'
n 0
Si n es impar
0n 1 n 1
1
1 @X 2
n k
X
2
n k (n 2k)'i A
senn ' = n ( 1) e(n 2k)'i
( 1) e
(2i) k n k
k=0 k=0
n 1
1 X
2
n k e
(n 2k)'i
e (n 2k)'i
= n 1 ( 1)
(2i) k 2i
k=0
n 1 n 1
( 1) 2 X
2
n k
= ( 1) sen (n 2k) '
2n 1 k
k=0
Z n 1 n
X
2
1
k
n ( 1) 2
n ( 1)
sen 'd' = [1 cos (n 2k) ]
0 2n 1 k n 2k
k=0
77
para n par
0 n n
1
1
1 @ n n X
2
n k
X1
2
n k (n 2k)'i A
senn ' = n n ( 1) +
2
( 1) e(n 2k)'i
+ ( 1) e
(2i) 2 k n k
k=0 k=0
0 n
1
1
2 @ n n X
2
n k e
(n 2k)'i
+ e (n 2k)'i
= A
n ( 1) + ( 1)
2
n
(2i) 2 k 2
k=0
0 n
1
n 1
n X
2
( 1) 2 @ n n k
= n ( 1)
2
( 1) cos (n 2k) 'A
2n 1 2 k
k=0
0 1
Z n n
X1 n ( 1)k
2
( 1) 2 n n
senn 'd' = n 1 @ n ( 1) 2 + sen (n 2k) A
0 2 2 k n 2k
k=0
4
1.71. Respuesta: Cg = (a; b):
3
4 3
1.72. Respuesta: Cg = (0; r) ; Cg = (0; 0; r):
3 32
1.73. Elegimos como sistema de referencia un sistema de referencia cartesiano
en el que hacemos coincidir el origen con el centro de la esfera y el eje OZ con
el eje del cono.
El volumen del sólido en cuestión es:
Z Z Z Z Z Z
2
V = dxdydz = sen ' d d d'
M [0;R] [0;2 ] [0; ]
2 3
= R (1 cos )
3
habiendo efectuado un cambio a coordenadas esféricas en la última integral.
Por razones de simetría es evidente que el centro de gravedad de este sólido está
situado sobre el eje OZ y por lo tanto,
xG = 0 ; yG = 0:
Puesto que
Z Z Z Z Z Z
3
z dxdydz = sen ' cos ' d d d' = R4 sen2 ;
M [0;R] [0;2 ] [0; ] 4
resulta …nalmente que:
RRR
z dxdydz 3R
zG = R R R M = cos2 :
M
dxdydzV 4 2
79
1.74. El sólido esta limitado por una esfera y un paraboloide que se cor-
tan en
z + z 2 = 2 ) z = 1:
En la resolución llamaremos h a este valor para obtener expresiones generales
de los volumenes referentes a un paraboloide y a un casquete esférico.
3.5
Z
3
2.5
1.5
0.5 h
0
-2 -1 0 1 2
Puesto que
Z Z Z Z h Z r
2 2 7
zdxdydz = z dz + z(r2 z 2 )dz = h3 + r2 h2 =
S 0 h 3 4 12
80
Z
a
Ty
Y
X
Su volumen es
Z Z Z Z Z Z !
a=2
V = dxdydz = 4 dxdz dy
S 0 Ty
Este ejercicio puede también resolverse considerando las secciones del sólido S
con planos paralelos a XOY . Tales secciones son elipses que pueden represen-
tarse en términos de la coordenada z:
x2 y2
+ = 1:
(a2 =4)(1 (z=b))2 a2 =4
Y
C
η
ξ
A O a B X
jb + a abj ab b a
p = p :
2
a +b 2 2
a +b 2
x = cos ; y = sen
Por razones de simetría es evidente que las coordenadas del centro de gravedad
de este sólido son G = (0; yG ; 0). Con el …n de calcular la coordenada yG , comen-
zaremos calculando su volumen, puesto que se trata de un sólido homogéneo.
87
x = cos ; y = sen
teniendo en cuenta que en el intervalo [0; =2] de variación del ángulo se tiene
que cos 0:
Z =2 Z 2a Z =2
32 16a3
V =4 d t2 dt = (1 cos3 ) d = (3 4):
0 2a cos a3 0 9
Finalmente:
Iy 24 a
yG = = :
V 5 3 4
M= k(b4 a4 ):
2
La simetría de la distribución de masa y de la forma geométrica del sólido per-
miten a…rmar que las coordenadas de su centro de gravedad son G = (0; 0; zG ).
Con el …n de determinar la coordenada zG hemos de calcular la integral:
Z Z Z p
Iz = kz x2 + y 2 + z 2 dxdydz
Z Z SZ
4
= k sen ' cos ' d d'd = k(b5 a5 ):
[a;b] [0; =2] [0;2 ] 5
Finalmente:
Iz 2 b5 a5
zG = = :
M 5 (b4 a4 )
o lo que es lo mismo:
2 2
x 4 y
+ = 1
4 5
3x
z =
4
y en forma paramétrica:
8 9
< x = 4 + 4 cos t =
y = 5sent ; t 2 [0; 2 ] :
: ;
z = 3 + 3 cos t
89
Así pues, el cono descrito en el enunciado viene dado por las ecuaciones paramétri-
cas: 2 3 2 3 2 3
x 0 4 + 4 cos t
4 y 5=4 0 5+ 4 5sent 5
z 10 7 + 3 cos t
eliminando t y en las ecuaciones anteriores obtenemos la ecuación cartesiana
de dicho cono:
2 2
x + z 10 8y
+ p = 1:
z 10 2 10 (z 10)
El volumen vendrá dado por:
Z Z Z Z 10 Z Z
V = dxdydz = dz dxdy
0 Exy
Z 10
p p
2 10 2 250 10
= (z 10) dz = :
0 8 3
La coordenada z del centro de gravedad será:
Z Z Z Z 10 Z Z
3 3
zg = p zdxdydz = p zdz dxdy =
250 10 250 10 0 Exy
Z 10 p
3 2 10 2 5
= p z (z 10) dz = :
250 10 0 8 2
La coordenada x del centro de gravedad será:
Z Z Z Z 10 Z Z
3 3
xg = p xdxdydz = p dz xdxdy =
250 10 250 10 0 Exy
Z 10 p
3 2 10 2 15
= p (z 10) (z 10) dz = :
250 10 0 8 2
Proponemos a continuación otro método de obtención del centro de gravedad.
Suponiendo que el corte de la super…cie del cono con el plano z = 0 es una elipse,
(suposición que es claramente cierta y está demostrada en la anterior resolución).
podemos calcular el eje mayor mediante:
a 8
= =) a = 20
10 10 6
Claramente el centro de la elipse estará en (10; 0; 0) y al eje menor lo llamaremos
b: El sólido está descrito por las ecuaciones:
2 3 2 3 2 3 8 9
x 0 10 + 10 cos t < 0 1 =
4 y 5=4 0 5+ 4 b sent 5 0 t 2
: ;
z 10 10 0 1
Utilizando como nuevas variables ; t; nos queda:
ZZZ Z 1 Z 1 Z 2
2 100b
V = dxdydz = 100b d d dt = :
0 0 0 3
90
1.86. a)
Y
D E
y(t) pt/q
t
B
C
O x(t) A X
Y D
t
F C
t
t
t t
G O X
91
x( ) = 2a cos (1 + cos )
y( ) = 2asen (1 + cos )
donde evidentemente es el ángulo polar.
Calculando la distancia al origen
2
( ) = x2 ( ) + y 2 ( ) = 4a2 (1 + cos )2 ;
( ) = 2a(1 + cos ):
c) RR RR 2
xdxdy cos d d 10 a3 5
xg = R R = RR = = a:
dxdy d d 6 a2 3
1.87. Las simetrías del sólido respecto a los planos XZ e Y Z aseguran que
las coordenadas yg y xg del centro de gravedad son nulas. El volumen del sólido
es la diferencia entre el volumen de la semiesfera exterior y el de la interior, es
decir
2
V = b3 a3
3
utilizando coordenadas esféricas para calcular zg
Z Z Z
1 3
zg = cos 'sen'd d'd
V [a;b] [0; 2 ] [0;2 ]
3 b4 a4
= b4 a4 =
4V 8 b3 a3
Para calcular la coordenada zeg del centro de gravedad de la super…cie, hace-
mos el límite de la expresión anterior cuando b tiende a a
3 b4 a4 1 b3 a
zeg = lim+ = lim =
b!a 8 b3 a3 b!a+ 2 b2 2
92
d Z
r
X
G
Y
La integral Z Z Z
y (x; y; z) dxdydz
M
se anula puesto que representa el producto de la masa de M por la coordenada
y de su centro de gravedad, que es precisamente el origen de coordenadas.
de la …gura.
Y
O
M= kR4 (1 cos ):
2
b) Por consideraciones de simetría es obvio que el centro de gravedad del
sólido está situado en su eje de simetría, ya que la densidad de masa está
distribuida con simetría respecto de ese eje. Por consiguiente, sólo es necesario
calcular la cota zG del centro de gravedad. Esa cota coincide con la distancia
al centro de la esfera.
Por de…nición se tiene que:
Iz
zG =
M
donde Iz es el momento estático de primer orden de la distribución de masa
respecto del eje OZ y M es la masa del sólido, que acabamos de calcular. En
coordenadas esféricas se tiene que:
Z Z Z Z Z Z
4 2
Iz = z dm = k sen ' cos ' d d'd = kR5 sen2 :
5
94
La identidad:
sen3 ' = (1 cos2 ')sen '
permite realizar con facilidad el cálculo:
Z
1
sen3 ' d' = (cos3 3 cos + 2):
0 3
Las otras dos integrales son triviales. Resulta así la expresión siguiente para el
momento de inercia pedido:
O
95
La masa de cada uno de los triángulos que componen el hexágono viene dada
por la integral Z Z
k dxdy
m=2 p
T x2 + y 2
siendo T el triángulo rectángulo rayado en la …gura. La descripción de ese
triángulo en coordenadas polares es:
h
0 ; 0
6 cos
p
en donde hemos utilizado la notación h = (l 3)=2. Por consiguiente, haciendo
un cambio a polares en la integral anterior, se obtiene la masa de cada uno de
los triángulos equiláteros:
Z =6 Z h= cos ! Z =6 p
d 3
m = 2k d d = 2kh = kl log 3:
0 0 0 cos 2
2 1 1
i= + log 3 kh3
3 3 4
4
I = 6i = + log 3 kh3 :
3
Dividiendo por la masa del triángulo obtenemos la distancia al centro del hexá-
gono del centro de masa de cada triángulo, que resulta ser:
h
d= :
2
Podemos ahora calcular la posición del centro de masa del semihexágono propuesto
en el enunciado. Para ello basta determinar la posición del centro de masa de tres
masas puntuales iguales a la masa de cada triángulo situadas respectivamente
en los centros de masa de los tres triángulos que constituyen el semihexágono.
La …gura adjunta ilustra esta situación.
X
97
resulta:
Z Z Z Z
(x2 + y 2 ) dxdy = 12 (x2 + y 2 ) dxdy
Sz Tz
Z p Z p ! p
lz ( 3=2) ( 3=3)x
2 2 5 3 4
= 12 dx (x + y ) dy = lz :
0 0 96
h
r
Y
O
Sustituyendo:
Z h
1 5 2 2 3
I1 = (R2 z 2 )2 dz = h R h + R4 h ;
0 5 3
Z R
1 5 2 2 3 8 5
I2 = (R2 z 2 )2 dz = h + R h R4 h + R :
h 5 3 15
La condición I1 = I2 proporciona para h la ecuación:
Esta es una función continua y derivable en R; f (0) = 4R5 < 0, f (R) = 4R5 >
0; lo que asegura la existencia de al menos una raíz en [0; R]: Puesto que
Z Z
X Y
d a
La densidad del toro viene dada por la función (x; y; z) = kjzj, en donde k > 0
es la constante de proporcionalidad.
Tal y como se ha situado el sistema de referencia, el momento de inercia
pedido ha de calcularse respecto del eje OZ, así que está de…nido por la integral:
Z Z Z Z Z Z
I= (x2 + y 2 )kjzj dxdydz = 2k (x2 + y 2 )z dxdydz
T T+
Z Z Z Z
2
Iy = jxj dxdy = 2 x2 dxdy
T T+
Z Z ! Z
a h(1 x=a) a
2 x a3 h
= x dy dx = 2 h x2 1 dx = :
0 0 0 a 6
En nuestro caso:
Z Z
1 2 1 12 =2
d
A= d = :
2 0 ( ah=6)h2 cos2 + a2 sen 2 ah 0 h2 cos2 + a2 sen 2
En esta última integral hacemos el cambio de variable tg = u, obteniéndose
…nalmente: Z +1
12 du 6
A= = :
ah 0 a2 u2 + h2 a2 h2
1.96. Teniendo en cuenta que la distancia del punto (e x; ye; ze) a un eje con vec-
tor director (x; y; z) viene dado por la norma del producto vectorial de ambos
vectores, resulta:
Z Z Z 2
k(x; y; z) (e x; ye; ze)k
I(x;y;z) = 2 de
xdey de
z
S k(x; y; z)k
RRR 2
S
k(yez ze y ; ze
x xe z ; xey ye x)k de xdey de
z
= 2
k(x; y; z)k
RRR 2 2 2 2
S
y ze + z ye 2yze y ze + z 2 x
e2 + x2 ze2 2xze xze + x2 ye2 + y 2 x
e2 xye de
2xye xde
y de
z
= 2
k(x; y; z)k
x2 Ix + y 2 Iy + z 2 Iz 2yzPyz 2xzPxz 2xyPxy
= 2
k(x; y; z)k
102
2
I(x;y;z) k(x; y; z)k = x2 Ix + y 2 Iy + z 2 Iz 2yzPyz 2xzPxz 2xyPxy = 1
2 1
k(x; y; z)k =
I(x;y;z)
1
k(x; y; z)k = p
I(x;y;z)
es decir la distancia al origen del punto (x; y; z) es justamente el inverso de la
raiz cuadrada del momento de inercia con respecto a un eje cuyo vector director
es el (x; y; z)
1.97.
Z Z Z Z Z b Z Z
2 2 2 z3 2b3
Ix = y +z dxdydz = y z+ dxdy = 2by 2 + dxdy =
S 3 b 3
Z Z
2b3 2 3
3 2
2 b r br4 1
= r + 2b sen2 d d = + = br2 4b2 + 3r2
3 [0;r]x[0;2 ] 3 2 6
Z Z Z
1
Iy = x2 + z 2 dxdydz = Ix = br2 4b2 + 3r2
S 6
Z Z Z Z Z
2 2 r4
Iz = x +y dxdydz = 2b x2 + y 2 dxdy = 2b2 = br4
S 4
Z Z Z
Pxy = xydxdydz = 0
S
Z Z Z
Pxz = xzdxdydz = 0
S
Z Z Z
Pyz = yzdxdydz = 0:
S
El elipsoide de inercia será:
1
br2 (4b2 + 3r2 ) x2 + y 2 + br4 z 2 = 1:
6
El punto de intersección de la recta x = y = z con el elipsoide es:
s
1
x =y =z = 1 2 2 2
6 br (4b + 9r )
y en coordenadas polares:
RR 2
RR 2
[1;2] [0; 4 ]
cos d d [1;2] [0; 4 ] sen d d
xg = 3 ; yg = 3 ;
4 4
es decir p p
14 2 14 2 2
xg = ; yg = :
9 9
Aplicando el teorema de Guldin:
p p
7 2 7 2 2
Vy = 2 xg A( ) = Vx = 2 yg A( ) = :
3 3
El volumen limitado por la esfera y con base en es:
Z Z Z p4 x2 y 2 Z Z p
V = dxdydz = 4 x2 y 2 dxdy
0
Z Z p
p 3
= 4 2 d d = :
[1;2] [0; 4 ] 4
1.99.
2
1
T
S
0
-1 0 1 2 3
X
-1
-2
104
1.100. La curva
x2 + y 2 2xy 2x + 6y = 0
se puede escribir en la forma
2
(x y) 2x + 6y = 0
u 1 1 1 x
=p :
v 2 1 1 y
En las nuevas coordenadas (u; v), el eje OX pasa a ser la recta u = v y la curva
dada tiene la ecuación p p
u2 2 2u + 2v = 0:
3
V
Y 1
0
-1 0 1 2 3
U
-1
Z p
2 Z p
2 2u
p u
2 p
1 2 2
vg = du vdv =
A 0 u 5A
en donde A es el área de S:
La coordenada x del centro de gravedad será:
vg + ug 11
xg = p =
2 30A
8t1 ; t2 2 [0; =2] ; t1 6= t2 ; se veri…ca que, (x(t1 ); y(t1 )) 6= (x(t2 ); y(t2 )):
1.5
X(t)
1
Y
X
0.5 Y(t)
0
0 0.5 1 1.5
t
Obsérvese que ambas funciones alcanzan su pmáximo para el pmismo valor del
parámetro, t0 2 [0; =2], tal que sen t0 = 2=3, cos t0 = 1= 3. Los valores
máximos de estas funciones son:
r
4 2 2
xM = x(t0 ) = ; yM = y(t0 ) = p :
3 3 3 3
106
0.8
0.6
Y
A
0.4
0.2
B
O 0.5 1
X
Z =2 Z =2
2 3 4
V = 2 ( )sen d = a3 sen 6 2 sen d
3 0 3 0
Z =2 Z
28 a3 7 6 28 a3 1 6
= sen cos d = u (1 u2 )3 du
3 0 3 0
212 a3
= :
3:7:11:13
Por consiguiente la abscisa del centroide del sólido viene dada por
Rb 2
xf (x) dx
xG = Ra b
a
f 2 (x) dx
Como aplicación de esta fórmula determinaremos la posición del centroide de
una semiesfera. La función f viene en este caso dada
p
f (x) = r2 x2 ; x 2 [0; r]
siendo r el radio de la esfera. Las integrales resultantes son inmediatas:
Z r Z
2 2 1 4 r 2 2
x r x dx = r r x2 dx = r3
0 4 0 3
108
1.105.
Y
4
2
Ω
1
0
0 1 2 3 4 5
X
Z Z Z 1 Z x+3 Z 4 Z 4=x Z 4 Z x 3
ydxdy = dx ydy + dx ydy dx ydy =
0 0 1 0 3 0
Z Z 4 Z
1 1 8 1 4
= (x + 3)2 dx + dx (x 3)2 dx =
2 0 1 x2 2 3
37 1
= +6 = 12:
6 6
Z Z Z 1 Z x+3 Z 4 Z 4=x Z 4 Z x 3
xydxdy = xdx ydy + xdx ydy xdx ydy =
0 0 1 0 3 0
Z 1 Z 4 Z 4
1 8 1
= x(x + 3)2 dx + dx x(x 3)2 dx =
2 0 1 x 2 3
Z 1
1 11
= 3x2 + 9x dx + 16 log 2 = + 16 log 2:
2 0 4
Por lo tanto:
11 4
xg = + log 2:
48 2
110
X
r
así que la ecuación del arco de circunferencia que gira alrededor del eje OX es:
p
3 p R R
y= R + R2 x2 ; x 2 ; :
2 2 2
El momento de inercia del sólido S respecto del eje OX es
Z Z Z Z R=2 Z Z
I= (y 2 + z 2 ) dxdydz = (y 2 + z 2 ) dydz dx
S R=2 Sx
resulta
a2 b=6 2a
xG = = :
Area (S) 3( 2)
Por razones de simetría, se tiene que:
2b
yG = :
3( 2)
V1 = Area (S) 2 xG :
Por otra parte, considérese el sólido engendrado al girar la parte de elipse con-
tenida en el primer cuadrante, alrededor de OY . Este sólido es la mitad del
elipsoide de semiejes a, a, b y por lo tanto, su volumen es
14 2
V = a b:
23
Finalmente, consideremos el cono de revolución obtenido al girar alrededor de
OY el triángulo que determina la cuerda C con los ejes coordenados en el primer
cuadrante. Su volumen es:
1 2
V2 = a b:
3
Puesto que V = V1 + V2 , resulta:
2 2 1 2 ab
a b= a b + 2 xG ( 2):
3 3 4
Despejando se obtiene xG = 2a=3( 2).
1.108. Necesitamos calcular los momentos de inercia con respecto a los ejes
112
z
-2
-4 -4
-2
u 0 8
6
2
4
4 2 v
0
Z Z Z Z Z Z Z
1 1 3 1
Ix = y 2 + z 2 dxdydz = 2 ; vu2 dvdu+2 ; v 3 dvdu
S 2 2 2 2
Z Z Z Z Z
2 2
Iy = x +z dxdydz = 4 (3=2; 1=2) v 3 dvdu
S
Z Z Z
Iz = x2 + y 2 dxdydz = Ix
S
Z Z Z
Pxy = xydxdydz = 0
S
Z Z Z
Pxz = xzdxdydz = 0
S
Z Z Z
Pyz = yzdxdydz = 0:
S
Puesto que
Z Z Z Z 2a(1+cos )
vu2 dvdu = d 4
sen cos2 d
0 0
Z
32a5 5
= (1 + cos ) sen cos2 d
5 0
5 Z
32a5 X 5
= cosk+2 sen d
5 0 k
k=0
5 5
X
32a 5 2
=
5 k k+3
k=0;2
4096a5
=
105
113
Z Z Z Z 2a(1+cos )
3 4
v dvdu = d sen3 d
0 0
Z
32a5 5
= (1 + cos ) sen 1 cos2 d
5 0
5 Z
5 X
32a 5 4096a5
= cosk sen d
5 0 k 105
k=0
5
5 X
32a 5 2 4096a5
=
5 k k+1 105
k=0;2
1024a5
=
35
sustituyendo resulta
1024a5
Iy = 2
35
4096a5 1024a5 11264 5
Ix = Iz = 2 + = a :
105 35 105
El elipsoide de inercia será:
11264 5 2 2048 5 2
a x + z2 + a y = 1:
105 35
El punto de intersección de la recta x = y = z con el elipsoide es:
r
105
x =y =z =
28672 a5
la distancia al origen de dicho punto es:
r
315
d=
28672 a5
y por lo tanto aplicando el ejercicio anterior, se tiene que el momento de inercia
con respecto a la recta x = y = z es
28672 5
Ir = a
315
114
INTEGRALES IMPROPIAS
2.1.1) Puesto que la función
sen x
f (x) =
xp
no está acotada en un entorno de cero para algunos valores de p; aislaremos las
dos posibles singularidades considerando los dos sumandos:
Z +1 Z 1 Z +1
sen x sen x sen x
p
dx = p
dx + dx
0 x 0 x 1 xp
Teniendo en cuenta que
sen x 0 si p < 1
lim =
x!0+ xp 1 si p = 1
y que
sen x
xp sen x
lim+ 1 = lim+ = 1;
x!0 p 1 x!0 x
x
aplicando el criterio de comparación se obtiene que la integral que …gura en el
primer sumando no es impropia si p 1; es convergente si 1 < p < 2 y es
divergente si p 2:
Analicemos la segunda integral. Si p 0; se veri…ca la negación del criterio
general de convergencia
hx i
0
(9" = 1) (8x0 > 1) 9x1 = + 1 2 ; x2 = x1 +
2
Z x2 Z x1 + Z x1 + Z
sen x sen x
dx = dx sen xdx = sen xdx = 2 1 = "
x1 xp x1 xp x1 0
Z A Z ] 1 Z (k+1)
[ AX ] 1 R (k+1)
[ AX
jsen xj jsen xj jsen xj k
jsen xj dx
dx dx+ dx p
1 xp 1 xp k xp (k + 1) p
k=1 k=1
[X
] A
1
2
= p p
(k + 1)
k=1
p 0 0<p 1 1<p<2 p 2
Z 1
sen x
dx No Impropia No Impropia Conv. Abs. Diverge
Z +1
0 xp
sen x
p
dx No Converge Converge Conv. Abs. Conv. Abs.
Z1 +1 x
sen x
dx No Converge Converge Conv. Abs. Diverge
0 xp
2) La función
cos x
g(x) =
xp
no está acotada en un entorno de cero para algunos valores de p: Por lo tanto
consideramos:
Z +1 Z 1 Z +1
cos x cos x cos x
p
dx = p
dx + dx
0 x 0 x 1 xp
Teniendo en cuenta que:
cos x 0 si p < 0
lim =
x!0+ xp 1 si p = 0
1
Z +1 Z X 1 Z (k+1)
sen x sen x sen x
dx = dx+ dx
1 xp 1 xp k=1 k xp
Z X 1
sen x
dx+ ( 1)k Ak
1 xp k=1
y que
cos x
xp
lim 1 = lim cos x= 1
x!0+ x!0+
xp
se deduce del criterio de comparación que el primer sumando no es una integral
impropia si p 0; es una integral impropia convergente si 0 < p < 1 y divergente
si p 1:
Analicemos la segunda integral. Si p 0; se veri…ca de manera análoga
a la integal con el seno la negación del criterio general de convergencia y por
consiguiente la integral
Z +1
cos x
dx es divergente.
1 xp
jcos xj 1
xp xp
se deduce inmediatamente la convergencia absoluta del segundo sumando si
p > 1 consideremos ahora 0 6 p 1
Z Z 1 Z 1
R (k+1)+ =2
+1
jcos xj =2
jcos xj X (k+1) + =2
jcos xj X k + =2
jcos xj dx
dx = dx+ dx p
1 xp 1 xp k + =2 xp (k + 1) p
k=0 k=0
1
X 2
= p p
(k + 1)
k=0
p 0 0<p<1 p=1 p 1
Z 1
cos x
dx No Impropia Conv. Abs Diverge Diverge
Z 0 xp
+1
cos x
p
dx No Converge Converge Converge Conv. Abs.
Z1 +1 x
cos x
dx No Converge Converge Diverge Diverge
0 xp
p
3) Haciendo el cambio x = u
Z +1 Z +1
sen u
sen x2 dx = 1 du
0 0 2u 2
117
Z +1 Z +1
2 cos u
cos x dx =
1 du
0 2u 2 0
De los apartados (a) y (b) se deduce que las integrales de Fresnel convergen
pero no absolutamente.
[X
2A
] Z (k+1)
2
k
e 2 log(cos x)2 dx
k
k=0 2
[X] 2A
Z 2
k
= e 2 ( 2) log(cos x) dx:
k=0 0
4 1 1
= p log 2 :
3 2 2
Asi pues la integral I es convergente por lo tanto
Z A [X
2A
]
k
x
e log(cos x)2 dx Ie 2 :
0 k=0
Haciendo A ! +1 :
Z 1 1
X
x k 1
e log(cos x)2 dx Ie 2 =I
0 k=0
1 e 2
2.3. Sea Z A
xq sen 2 x
I(A) = xp e dx
0
118
] 1Z
[ AX (k+1) Z A [A] Z
X (k+1)
q 2 q 2
x sen x xq sen2 x
xp e x sen x
dx I(A) = xp e dx xp e dx
k=0 k 0 k=0 k
Además,
] 1Z
[ AX (k+1) [A] Z
X (k+1)
p (k+1)q q
sen 2 x p p (k )q sen2 x
(k ) e dx I(A) (k + 1) e dx
k=0 k k=0 k
] 1Z
[ AX [A] Z
X
p (k+1)q q
sen 2 x p+1 p p (k )q sen2 x
(k ) e dx I(A) + (k + 1) e dx
k=0 0 k=1 0
y=x
2.5
y=2x/ π
1.5
Y
1
y=senx
0.5
0
0 1 2 3
X
resulta
] 1 Z
[ AX =2 [A] Z
X =2 2
p (k+1)q q
x2 p+1 p p (k )q 4x2
2 (k ) e dx I(A) + 2 (k + 1) e dx
k=0 0 k=1 0
2
y teniendo en cuenta que e t es positiva, pueden acotarse las integrales susti-
tuyendo el intervalo de integración por otro menor o mayor respectivamente:
[ AX
] 1 p Z =2 [A]
X p p+1 Z 1
2 (k ) t2 p+1 (k + 1) t2
q=2
e dt I(A) + q=2
e dt:
(k + 1) q=2 0 (k ) 0
k=0 k=1
k=1
b)
Z +1 Z +1
(' )(x) = '(t) (x t) dt = f (t )g(x t ) dt:
1 1
Haciendo el cambio =t
Z +1
(' )(x) = f ( )g(x ) d = (f g)(x ):
1
y2 Z +1
2
+ r y
e 2b2 b2 ( 2a12 + 2b12 ) d
(f g)(y) = q e
1 1 1
2a2 + 2b2
y2 y2 p y2
e 2b2 p 4b4 ( 2a12 + 2b12 ) = pba 2 e 2(a2 +b2 )
= q e :
1
+ 1 a2 + b2
2a2 2b2
Finalmente resulta:
p (x )2
2 2(a2 +b2 )
(' )(x) = (f g)(x ) = ab p e :
a + b2
2
p 1
(1+ x)(1 x) 1 1
lim = lim p =p
x!1 p1 x!1 1+ x 1+
1 x
8 R1 9
Z 1 < p1 p1 dx
>0 =
dx 0 + 1
x x2
I( ) = p = R1
0 1+( 1) x x2 : p1 p dx
<0 ;
0 1 1
x+x2
8 R1 9
>
< p1 q dx
>0 >=
0 +1 2 1 2
( ) (x )
= R1 2 2
:
>
: p1 q dx
<0 >
;
0 1 2 +1 2
(x 2 ) ( 2 )
resulta:
Z 0 arccos 1
1 sen 1+
I( ) = p d = p si > 0;
arccos( 11+ ) jsen j
Z
1 0 +1
2 sh
+1 arg ch 11+ arg ch 1
I( ) = p +1
d = 2
+1
p = p 1+ si < 0:
arg ch ( 11+ ) 2 sh 2
Estudiemos la continuidad de I( ) en = 0:
q 1 2
1 1 p
arccos 1+ 1 ( 11+ )
2 (1+ )2
(1+ )
lim I( ) = lim p = lim 1 = lim 1 =2
!0+ !0+ !0+ p !0+ p
2 2
q 1 2
1 1p
arg ch 1+ ( 11+ )
2
1 (1+ )
2
(1+ )
lim I( ) = lim p = lim+ = lim+ =2
!0 !0 !0 p1 !0 p1
2 2
Z 1
dx p
I(0) = p = 2 1 x j10 = 2:
0 1 x
tx 1
e t
tx+1
lim 2
= lim = 0;
t!+1 t t!+1 et
Ahora bien,
Z 1 Z 0
x 1 y 1
2 (x; y) = t (1 t) dt = (1 u)x 1 y 1
u du = 1 (y; x):
1 1
2 2
122
R1
La integral 1 (x; y) tiene el mismo caracter que 02 tx 1 dt; :así que converge si
x > 0 cualquiera que sea y: Por lo tanto 2 (x; y) converge si y > 0 cualquiera
que sea el x: Concluimos que la función Beta converge si x > 0 e y > 0:
b)
Z +1 Z +1
p t p t +1
(p + 1) = t e dt = t e j0 + ptp 1
e t
dt = p (p)
0 0
(n) = (n 1) (n 1) = (n 1) (n 2) (n 2) = :::::::: =
Z +1
t
= (n 1)! (1) = (n 1)! e dt = (n 1)!
0
c)
(2) = (1) = 1
Puesto que la función es continua en x > 03 , se veri…ca que
También,
1 1
lim (x) = lim ( ) = lim n ( + 1) = lim n (1) = +1
x!+0 n!1 n n!1 n n!1
50
40
1
10
GAMMA(X)
GAMMA(X)
30
20
10
0
10
0
0 2 4 6 8 10 0 1 2 3 4 5 6
X X
e) 1)Haciendo el cambio u = 1 t
Z 1 Z 0
(x; y) = tx 1
(1 t)y 1
dt = (1 u)x 1 y 1
u du = (y; x)
0 1
4) Para p; q 1
Z +1 Z +1 Z A Z A
p 1 x q 1 y
(p) (q) = x e dx y e dy = lim dx xp 1 q 1
y e (x+y)
dy
0 0 A!+1 0 0
Z Z
= lim xp 1 q 1
y e (x+y)
dx dy = lim I [S(A)]
A!+1 [0;A] [0;A] A!+1
(A) = (x; y) 2 R2 : x 0; y 0; x + y A
resulta
Z A Z 2A
(p; q) lim v p+q 1
e v
dv lim I [S(A)] (p; q) lim v p+q 1
e v
dv
A!+1 0 A!+1 A!+1 0
(p) (q)
(p; q) = :
(p + q)
se deduce que
(p; q) = (p + 1; q) + (p; q + 1)
125
y aplicando (2)
q
(p; q) = (p + 1; q) + (p; q + 1) = (p + 1; q) + (p + 1; q) =
p
p+q p+q p+qq+p+1
= (p + 1; q) = (q; p + 1) = (q + 1; p + 1)
p p p q
p + q q + p + 1 (p + 1) (q + 1) p+qq+p+1 p (p)q (q)
= =
p q (p + q + 2) p q (p + q + 1) (p + q) (p + q)
(p) (q)
= :
(p + q)
40
35
30
25
20
15 V
10
0 0
0
2 2
Y X
4 4
2.7. Z 1
(p) (p) p 1
(p; p) = = tp 1
(1 t) dt:
(2p) 0
1+u
Haciendo el cambio t = 2 ; resulta:
Z 1 Z 1
(p) (p) p 1 p 1 du 1 p 1
= (1 + u) (1 u) = 1 u2 du:
(2p) 1 22p 1 22p 1 1
Operando
(p + 21 ) (p) 22p 1
1
(2p) = 22p 1
= p (p + ) (p)
( 12 ) 2
n
cubo [0; 1] en el recinto Tn
x1 = 1 2 ::::: n 1 n
x2 = 1 2 ::::: n 1 (1 n)
:::::::::::::::::::::::::::
xn 1 = 1 2 (1 3)
xn = 1 (1 2 )
sustituyendo en la integral:
Z Z n
!pn+1 1
X
p1 1 p2 1 pn 1
::: x1 x2 : : : xn 1 xk dx1 dx2 : : : dxn
Tn k=1
Z Z
p1 1 p2 1 pn 1 pn+1 1 n 1
= :: ( 1 :: n) 1 :: n 1 (1 n) :: [ 1 (1 2 )] (1 1) 1 :: n 1 d 1 ::d n
n
[0;1]
Z 1 Z 1 Z 1
p1 +::+pn 1 pn+1 1 p1 +::+pn 1 1 pn 1 p1 1 p2 1
= 1 (1 1) d 1 2 (1 2) d 2 ::: n (1 n) d n
0 0
0 10
n Z
Y 1 Pn+1 k
n
Y n+1
Xk
( j=1 pj ) 1 pn+2 k 1 @ A
= k (1 k) d k = pj; pn+2 k
k=1 0 k=1 j=1
Haciendo el cambio
1
x1 = a1 u1p1
1
x2 = a2 u2p2
a1 a2 :::an p11 1 1
1 1
1
:::::: cuyo J' = u u2p2 :::unpn
1
pn 1
p1 p2 :::pn 1
xn 1 = an 1 un 1
1
pn
xn = an un
resulta:
Z Z
n a1 a2 :::an p11 1 1
1 1
1
VE(n) = 2 :: u u2p2 :::unpn du1 du2 ::::dun
Tn p1 p2 :::pn 1
1 1 1
::: (1)
n a1 a2 :::an p1 p2 pn
= 2 :
p1 p2 :::pn 1
+ ::: + 1
+1
p1 pn
127
2.9.Haciendo el cambio x = e t
Z 1 n Z 0 Z 1
1
xm log dx = e mt n
t e t
dt = e (m+1)t n
t dt
0 x 1 0
xk e x log k
xk :
Puesto que Z +1
xk dx = +1
0
para estos valores de k, resulta que la función no es integrable. Es fácil com-
probar que cuando k > 1, la integral es convergente. En tal caso, si se utiliza el
cambio de variable t = x log k, se obtiene que:
Z +1 k Z +1
x 1 (k + 1)
x
dx = k+1
tk e t dt = :
0 k (log k) 0 (log k)k+1
que implica que la función integrando está acotada en las proximidades del cero.
La función exponencial nunca toma valores negativos, por lo que se cumple que:
y2
1 e 1
:
y2 y2
Por lo tanto, el integrando es una función positiva mayorada por una función
integrable en el intervalo [1; +1[. De todo lo anterior se deduce que la integral
pedida es convergente.
Para efectuar su cálculo, efectuaremos una integración por partes eligiendo
2
u = 1 e y , dv = (1=y 2 )dy, con lo que se tiene:
Z " #+1 Z
+1 y2 y2 +1 p
1 e 1 e y2
dy = +2 e dy = :
0 y2 y 0
0
cuyo jacobiano es
p+q
a pq
J= u(1=p) 1 (1=q) 1
v
pq
resultando:
p+q
pq +1
Z Z
4a
V = (1 u v)u(1=p) 1 (1=q) 1
v dudv
pq u+v 1 ; u;v 0
( p+q
p+q p+q
pq +1 pq +1
4a p+q 1 1 4a pq ) (2) (1=p) (1=q)
= ;2 ; = p+q
pq pq p q pq ( pq + 2) ( p+q
pq )
p+q
4a pq +1 (1=p) (1=q)
= :
pq ( p+q
pq + 2)
tx 1
(1 e t ) +x 1
lim = lim t (1 e t ) = 1 si = x
t!0+ 1=(t ) t!0+
la integral Z 1
tx 1
(1 e t ) dt
0
tiene el mismo carácter, en lo que se re…ere a la convergencia, que la integral
Z 1
1
x
dt
0 t
haciendo el cambio t = nu
Z 1
n
In (x) = nx 1 u x 1
(1 u) n du = nx (x; n + 1) =
0
x
n (x) (n + 1) nx (x)n!
= = =
(x + n + 1) (x + n) (x + n 1) ::::::x (x)
nx n! nx n!
= = Q
n :
(x + n) (x + n 1) ::::::x
(x + k)
k=0
a)
Z 1
1x 1! 1
I1 (x) = tx 1
(1 t) dt = =
0 Q
1 (x + 1)x
(x + k)
k=0
Z 2 2
t 2x 2! 2x 2
I2 (x) = tx 1
1 dt = =
0 2 Q
2 (x + 2)(x + 1)x
(x + k)
k=0
Z 3 3
t 3x 3! 3x 6
I3 (x) = tx 1
1 dt = =
0 3 Q
3 (x + 3)(x + 2)(x + 1)x
(x + k)
k=0
c) Sabiendo que
x
1 lim x log(1 1
) 1
= lim ex log(1 ) = ex!+1
1 x
lim 1 x =
x!+1 x x!+1 e
se obtiene:
2 !( nt ) 3t 2 !( nt ) 3t
n
t 1 1
lim 1 = lim 4 1 n
5 = 4 lim 1 n
5 =e t
n!1 n n!1
t
n!1
t
y=1-x
2
0
-2 -1 0 1 2
X
-1
Y
-2
resultando:
t2 t t t
)e n
(1 1 e n
n2 n
y teniendo en cuenta que todos los terminos son positivos,
n n
t2 t t
1 e 1 e t:
n2 n
8
(1+h)n , 1+nh
0
-1 -0.5 0 0.5 1
h
-2
-4
-6
-8
n n
t2 t t2 t t
1 n e 1 e 1 e t:
n2 n2 n
Multiplicando por tx 1 e integrando:
Z n Z n n Z n
x 1 t2 t t
t 1 e dt tx 1
1 dt tx 1
e t dt:
0 n 0 n 0
1 2
e 1 ; 06 61
g( ) = e 1+ ! g0 ( ) = e + 1 ! g 00 ( ) = e
g es convexa alcanza el mínimo en = 0 y g(0) = 0 por consiguiente g( ) 0 para todo :
El anterior argumento de la convexidad nos permite también demostrar la desigualdad de
Bernouilli
133
Z n Z n n
1 t
(x) lim tx+1 e t dt lim tx 1
1 dt (x)
n!1 n 0 n!1 0 n
y …nalmente: Z n n
t
(x) lim tx 1
1 dt (x):
n!1 0 n
d a 1 u du
I = = u 2
0 a(1 + )c 0 a(1 + 1
c
u ) (1 u)
Z 1
1 b
1 c 2 b
a +1
1 b b b b
= ua (1 u) du = ;c con >0 c> :
a 0 a a a a a
Finalmente
1 b b
I= ;c ; a > 0; b > 0; ac > b:
a a a
2) En la integral (3) hacemos b = m + 1; a = n; c = 1 resultando
Z +1 m
x dx 1 m+1 m+1
= ;1 = n > m + 1 > 0:
0 1+x n n n n nsen m+1n
y por lo tanto,
Z =2
1 +1 +1
tg x dx = ;1 = +1
0 2 2 2 2sen 2
y
x
y
t
Z y Z x Z y Z y
f (t) 1 1 1
F (y) = p dt p dx = f (t) p p dx dt
0 0 x t y x 0 t x t y x
c) Tranformemos la integral anterior en una beta haciendo un cambio lineal que
transforme el intervalo [t; y] en el [0; 1] ; es decir x = t + u(y t)
Z y Z 1
1 (y t)du 1 1
p dx = p = ; = :
t (x t) (y x) 0
2
u(y t) (1 u) 2 2
2.19. 1) La integral
Z 1
(x 1) xk
I(k) = dx
0 log x
aparentemente es impropia en 0 y en 1 pero si analizamos detenidamente la
función subintegral en un entorno de 1 tenemos
(x 1) xk (x 1) 1
lim = lim = lim =1
x!1 log x x!1 log x x!1 1=x
(x 1) xk xk 0 si k > 0
lim+ = lim+ = k 1
x!0 log x x!0 log x limx!0+ kx1=x = limx!0+ kxk = 1 si k < 0
lo que asegura que la integral solo es impropia para valores negativos de k. Para
estos valores, observando la acotación
Z Z 1=2 Z 1=2
1 1=2 xk (x 1) xk xk
dx 5 5 dx
2 0 log x 0 log x 0 log x
se deduce que la integral pedida tiene el mismo carácter que la integral
Z 1=2 Z log 2 ku Z +1
xk e +1 si k + 1 6 0
dx = e u du = u 1 e (k+1)u du =
0 log x +1 u log 2 convergente si k > 1
y Z +1
1 1
2
du = :
log 2 u log 2
Resumiendo
k>0 la integral no es impropia
0>k> 1 la integral es impropia y convergente
1>k la integral diverge.
para valores de k en el intervalo [b; +1) para todo b > 1; se tiene que:
Z 1 Z 1 Z 1
(x 1) xk log x 1 1
I 0 (k) = dx = (x 1) xk dx = xk+1 xk dx =
0 log x 0 0 k+2 k+1
de donde
k+2
I(k) = log + cte
k+1
136
y puesto que
k+2
lim I(k) = 0 ) cte = 0 ) I(k) = log :
k!+1 k+1
c) Aunque no se pide, vamos a demostrar que
lim I(k) = 0:
k!+1
0 x 2 [0; 1)
lim fk (x) =
k!1 1 x = 1:
(x 1) xk
M bk ! 0
log x
137
Puesto que
1 x
5 0;
N 0 (x) =
x2
resulta que N (x) es decreciente y puesto que
lim N (x) = log x 1 + 1=x = 0;
x!1
podemos asegurar que N (x) > 0 si x 2 [0; 1] : Por lo tanto, g 0 (x) > 0 si x 2 [0; 1] ;
así que g(x) es creciente, lo que conduce a:
x 1
0 5 g(x) = 5 1:
log x
Entonces
Z 1 Z 1
1) xk(x xk+1 1
0 5 I(k) = dx 5 xk dx = = ; k+1>0
0 log x 0 k + 1 k + 1
Z 1
(x 1) xk 1
0 5 lim I(k) = lim dx 5 lim = 0:
k!+1 k!+1 0 log x k!+1 k + 1
2 log (1 + a) + c ; a>0
I(a) =
2 log (1 a) + c ; a<0
y como I es impar:
2 log (1 + a) ; a>0
I(a) = = sgn(a) log (1 + jaj)
2 log (1 a) ; a<0 2
2)
Z Z =2 Z =2
2
=2 cos x x
J = log(senx)dx = x log(senx) j0 x dx = dx
0 0 senx 0 tgx
Z +1
arctgu du
= = I(1) = log 2:
0 u 1 + u2 2
3) Si 0 < x < 1;
Por lo tanto,
Z 1=2 Z 1=2 Z 1=2
1
log (x)dx + log (1 x)dx = log log sen ( x) dx
0 0 2 0
Z 1 Z =2
1 1 1
log (x)dx = log log (senu) du = log :
0 2 0 2 2
4)
( + 1)
F 0 ( ) = log ( + 1) log ( ) = log = log si >0
( )
Z
F( ) F (1) = log xdx si >0
1
Z 2 Z 1 Z 1 Z 1 Z 1
F (1) = log (x)dx = log (u+1)du = log [u (u)] du = log u du+ log (u) du = F (0) 1
1 0 0 0 0
Z Z 1 Z Z
1
F ( ) = F (1)+ log x dx = F (0)+ log u du+ log x dx = log + log x dx
1 0 1 2 2 0
Finalmente
1
F( ) =
log + log si >0
2 2
Se puede observar que la función es continua por la derecha en 0.
esto quiere decir que F (x) es creciente y por consiguiente inyectiva, lo cual nos
indica que la solución de
F (x) =
si existe, es única. Calculemos los extremos de F (x) :
Z +1 h i+1 Z +1
1=t2 1=t2 2 1=t2
sup F (x) = 1 e dt = t 1 e + e dt =
x2R 0 0 0 t2
Z 0
s2
p
= 0+ 2e ds =
+1
como F es impar, p
inf F (x) =
x2R
x2 x2
e cos( x) e
y puesto que Z p
+1
x2
e dx = ;
0 2
F ( ) admite una mayorante convergente, lo que implica que converge uniforme-
mente para todo : Además,
x2
@ e cos( x)
x2
= xe sen( x)
@
x2
es continua en todo R2 y mayorada por xe y teniendo en cuenta que
Z +1 2
x2 e x +1 1
xe dx = j0 = ;
0 2 2
la integral Z +1
x2
xe sen( x)dx
0
140
y puesto que Z p
+1
x2
F (0) = e dx = ;
0 2
se obtiene p
2
F( ) = : e 4
2
c) Considérese la integral impropia dependiente del parámetro :
Z +1
2 sen x
G( ) = e x dx:
0 x
Esta integral es convergente y además,
@ e x2 sen x
x x2
=e cos( x)
@
es continua en todo R2 y la integral
Z +1
x2
e cos( x)dx
0
Integrando
Z 2a Z 2a p
2
p Z 2a 2
G(2a) = G(0) + F ( )d = 0 + e 4 d = e 4 d
0 0 2 2 0
2.23. a)
Z +1
ax sen x
F (a) = e dx (a 0)
0 x
141
Puesto que Z +1 x
e 1
e x
dx = j+1
0 = >0
0
R +1
0
e ax sen x dx converge uniformemente en [ ; 1) ; >0
d) De los apartados anteriores se obtiene:
Z +1 Z +1
F 0 (a) = g(a; x)dx = e ax sen xdx
0 0
Z +1
= e ax cos x j+1
0 + ae ax cos xdx
0
Z +1
= 1 + ae ax sen x j+1
0 + a2 e ax sen xdx
0
2 0
= 1 a F (a):
Despejando:
1
F 0 (a) = a>0
1 + a2
e) Integrando en la expresión anterior:
F (a) = C arctga:
Puesto que
Z +1 Z +1
ax sen x ax sen x
jF (a)j = e dx e dx
0 x 0 x
Z +1 ax
e 1
e ax
dx = j+1
0 = ; a>0
0 a a
resulta que
lim F (a) = 0 = C )C= :
a!+1 2 2
142
Finalmente
F (a) =
arctga; a > 0:
2
f) Aplicando un conocido teorema del valor medio7 , para 0 < x1 < x2 :
Z x2 ax1 Z
ax sen x e
sup e dx = sup sen x dx
a>0 x1 x a>0 x1 x1
ax1
e 2
sup 2 = :
a>0 x1 x1
Por lo tanto
Z x2
ax sen x 2
lim sup e dx = lim = 0;
x1 !+1 a>0 x1 x x1 !+1 x1
R +1
lo que implica que 0 e ax senx x dx converge uniformemente en el intervalo
[0; 1) :
g) Por la convergencia uniforme, F es continua para a 2 [0; 1) y entonces
Z +1
sen x
lim+ F (a) = = F (0) = dx:
a!0 2 0 x
Por lo tanto
G0 (y) = y jyj + c ; G0 (0) = 0 ) c = 0
y de aquí
Z +1
sen4 xy 3
G(y) = dx = jyj + c ; G(0) = 0 ) c = 0:
0 x4 3
Finalmente Z +1
sen4 x
G(1) = dx = :
0 x4 3
6) Sea
Z +1
sen3 xy
H(y) = dx:
0 x3
Derivando:
Z +1 Z +1
3sen2 xy cos xy 3sen xysen2xy
H 0 (y) = 2
dx = dx
0 x 0 2x2
Z
3 +1 cos xy cos 3xy
= dx:
4 0 x2
Z +1
00 3 sen xy + 3sen xy 3 3
H (y) = dx = 2 = ; y > 0:
4 0 x 4 2 4
De aquí se obtiene
Z +1
3sen2 xy cos xy 3
H 0 (y) = 2
dx = jyj +c ; H 0 (0) = 0 ) c = 0:
0 x 4
Por lo tanto
Z +1
sen3 xy 3
H(y) = dx = y jyj +c ; H(0) = 0 ) c = 0
0 x3 8
Finalmente Z +1
sen3 x 3
H(1) = dx = :
0 x3 8
Todas las integrales de este ejercicio pueden resolverse integrando adecuada-
mente por partes.
2.25.Sea Z +1
sen ax
F (a) = dx:
0 x(1 + x2 )
Derivando respecto del parámetro se obtiene:
Z +1
0 cos ax
F (a) = dx
0 1 + x2
ya que se veri…can las condiciones siguientes, que aseguran la validez de la
derivación
cos ax
2) es continua en todo R2 ;
1 + x2
Z +1 Z +1
cos ax dx
3) La integral 2
dx esta mayorada por = ; así que
0 1 + x 0 1 + x2 2
converge uniformemente.
Z +1
xsenax
Derivando nuevamente, y teniendo en cuenta que dx con-
0 1 + x2
verge uniformemente en [ ; 1) , > 0 :
Z +1 Z +1
00 xsenax x2 senax
F (a) = dx = dx
0 1 + x2 0 x (1 + x2 )
Z +1 Z +1
senax senax
= dx + dx =
0 x 0 x (1 + x2 )
= + F (a) ; a > 0:
2
Así pues, F (a) satisface la ecuación diferencial
F (a) = c1 ea + c2 e a
+ :
2
Utilizando la continuidad de F (a) y F 0 (a); se obtiene
F (0) = 0 ) c1 = 0
Z +1
0 dx
F (0) = = ) c2 =
0 1 + x2 2 2
y por lo tanto
a
F (a) = 1 e ; a>0
2
F 0 (a) = e a
; a>0
2
Por las simetrías de las funciones
Z +1
sen ax jaj
F (a) = dx = sgn(a) 1 e
0 x(1 + x2 ) 2
Z +1
0 cos ax jaj
F (a) = dx = e :
0 1 + x2 2
xn e ax
cos bx xn e ax
xn e x
;
junto con Z Z
1 1
n x yn e y
(n + 1)
x e dx = n+1
dy = n+1
0 0
aseguran que la integral J admite una mayorante de Weierstrass convergente y
puesto que la función subintegral es continua en todo R2 ; se puede derivar con
respecto a los parámetros todas las veces que se desee en [ ; 1) R cualquiera
que sea > 0:
Z +1 Z
1 a +1 ax
1) I (a; b) = e ax cos bxdx = e ax senbx j+1 0 + e senbxdx
0 b b 0
Z
a 1 ax a +1 ax
= e cos bx j+10 + e cos bxdx
b b b 0
a a 2
= 2
I (a; b)
b b
Z +1
a
I (a; b) = e ax cos bxdx = 2 :
0 a + b2
Z +1 Z
b +1 ax b
2) e ax senbxdx = e cos bxdx = 2 :
0 a 0 a + b2
Z +1 Z +1 ax Z +1
ax @e
senbx @ ax
3) xe cos bx = dx = e senbxdx
0 0 @b @b 0
@ b a2 b2
= = 2:
@b a2 + b2 (a2 + b2 )
Z +1 Z +1 ax Z +1
ax @e
cos bx @ ax
4) xe senbxdx = dx = e cos bxdx
0 0 @b @b 0
@ a 2ab
= = 2:
@b a2 + b2 (a2 + b2 )
Z +1 Z +1 ax Z +1
@xesenbx @
5) x2 e ax
cos bxdx = dx = xe ax
senbxdx
0 0 @b @b 0
@ 2ab 2a a2 2b2
= = 3 :
@b (a2 + b2 )2 (a2 + b2 )
Z +1 Z +1 ax Z +1
@xesenbx @
6) x2 e ax
senbxdx = dx = xe ax
senbxdx
0 0 @a @a 0
2 2
@ 2ab 2b 2a b
= 2 = 3 :
@a (a2 + b2 ) (a2 + b2 )
147
f (x) + k = [g(x)]2
así que basta con determinar f en el primer cuadrante. En lo que sigue, supon-
dremos que x 0, y 0.
Puesto que
arctg (xt)arctg (yt) (xt)(yt)
lim = lim = xy
t!0+ t2 t!0+ t2
la integral que de…ne a la función f presenta una singularidad únicamente en
+1.
Por otra parte,
arctg (xt)arctg (yt) ( =2)2
(t ! +1)
t2 t2
148
1 x2 1 y2 1
= ; (x2 6= y 2 ):
(1 + x2 t2 )(1 + y 2 t2 ) x2 y2 1 + x2 t2 x2 y 2 1 + y 2 t2
Z +1 Z
@2f dt 2 du
(x; x) = 2 =
@x@y 0 (1 + x2 t2 ) 0 x (1 + tg 2 u)
Z 2
2 cos u
= du =
0 x 2x
@2f
(x; y) =
@x@y 2 (x + y)
c) Integrando respecto de la variable x en la fórmula que proporciona la
derivada segunda @ 2 f =@x@y, se tiene:
Z
@f dx
= + k(y) = log(x + y) + k(y):
@y 2 x+y 2
e integrando ahora respecto de la variable y, resulta:
Z
f (x; y) = log(x + y) dy + g(x) + h(y) = (x + y) log(x + y) + g(x) + h(y)
2 2
Aplicando el criterio de la mayorante de Weierstrass se prueba la convergencia
uniforme en cualquier compacto de R2 de la integral que de…ne a f , lo que
permite asegurar la continuidad de dicha función,.Por lo tanto,
h i
0 = f (x; 0) = lim f (x; y) = lim (x + y) log(x + y) + g(x) + h(y)
y!0 y!0 2
Análogamente
h(y) = y log y g(0):
2
Finalmente
g(0) = h(0):
2.30. 1) Integrando
Z +1 Z +1
dt 1=xdt 1 t
I(x) = = h i= arctg j+1 = :
0 x + t2
2
0 x 1+ t 2 x x 0 2x
x
I20 (x) 3 1 3
I3 (x) = = =
2 2x 4 4x x4 4 4 x5
8 Observese la continuidad de la función subintegral así como la de sus derivadas con re-
I30 (x) 3 5 3 5 1
I4 (x) = = =
2 3x 4 4 6x x6 2 2 4 6 x7
que sugiere la fórmula
3 5 2n 3 1
In (x) =
22 4 6 2n 2 x2n 1
d) Operando formalmente
Z +1 Z +1 p
dt t2 3 5 2n 3 p
lim 2 n = e dt = = lim n
n!1 0 1 + tn 0 2 n!1 22 4 6 2n 2
3 5 2n 3 p 1
lim n= p
n!1 2 4 6 2n 2
Cn = (x; y) 2 R2 : x2 + y 2 n2 :
Entonces
Z Z Z Z Z Z
2 2
(x2 +y 2 ) 2
e (x +y ) dx dy = lim e dx dy = lim e d d
R2 n!1 Cn n!1 [0;n] [0;2 ]
2
e
= 2 j+1
0 = :
2
Por otra parte, la sucesión de conjuntos siguientes también recubren monóton-
amente todo R2
y por lo tanto
Z Z Z Z Z +1 Z +1
2 2
(x2 +y 2 ) x2 y2
e (x +y ) dx dy = lim e dx dy = e dx e dy = :
R2 n!1 Cn 1 1
Se deduce que Z Z p
+1 +1
x2 1 x2
e dx = e dx =
0 2 1 2
Cosideremos ahora la sucesión de conjuntos que recubre monótonamente todo
Rn
Km = X 2 Rn : X t QX m :
151
Entonces
Z Z Z Z
X T QX X T QX
e dX = lim e dX
Rn m!1 Km
Dm = fY 2 Rn : jyk j m k = 1:::::ng :
2.32.
Z Z Z Z
2 2
lim sen(x + y ) dx dy = lim senx2 cos y 2 + seny 2 cos x2 dx dy
r!+1 Qr r!+1 Qr
Z +1 Z +1
= 2 senx2 dx cos y 2 dy
1 1
resultado …nito, ya que las integrales de Fresnel son convergentes. Por otra parte
Z Z Z Z
2 2
lim sen(x + y ) dx dy = lim sen 2 d d = lim cos r2 :
r!+1 Cr r!+1 Qr r!+1
2.33. Sea
Mr = (x; y) 2 R2 ; 0 x 2y ; x2 + y 2 r2 :
Es sabido que toda matriz real y simétrica tiene todos sus valores propios reales
y es diagonalizable por semejanza. Si la matriz es de…nida positiva, esos valores
propios son números reales positivos. Además se puede siempre encontrar una
base ortonormalizada de vectores propios, lo que signi…ca que la matriz que
establece la relación de semejanza puede elegirse ortogonal. En términos más
precisos, si son 1 ; : : : ; n los valores propios de la matriz A y si introducimos
la matriz diagonal
D = diag f 1 ; : : : ; n g
existe una matriz ortogonal P (es decir, tal que P P T = I) que relaciona las
matrices A y D mediante
D = P T AP:
Utilizando este resultado, hagamos en la integral I el cambio de variables
de…nido por x = P y. En esta notación se han identi…cado x e y con las
matrices columna cuyos elementos son las componentes de los vectores x e y.
Con este cambio de variables se tiene que:
n
X
T T T T 2
x Ax = y P AP y = y Dy = i yi :
i=1
Recordando que Z +1 p
t2
e dt =
1
un sencillo cambio de variable conduce a:
Z +1 r
2
e i yi dyi = :
1 i
b) Calculemos la entropia de
0 1
X
g(x1 ; : : : ; xn ) = exp @ aij xi xj A :
1 i;j n
ya que g es siempre positiva y menor que uno lo cual hace que la función subinte-
gral de la entropia sea siempre positiva. Utilizando el mismo cambio de variables
del apartado anterior, que transforma la matriz A en la matriz diagonal D, re-
sulta:
Z Z
Hg = g log g = exp xT Ax log exp xT Ax dx1 : : : dxn
Rn Rn
0 1 !
Z Xn Xn
= exp @ j yj
2 A i yi
2
dy1 : : : dyn
Rn j=1 i=1
0 1
n
X Z n
X
= i exp @ j yj
2A
yi 2 dy1 : : : dyn :
i=1 Rn j=1
y2
Integrando por partes, con u = y, dv = ye dy, se calcula la integral
Z +1 +1 Z +1 r
y2 1 y2 1 2 1
y2 e dy = ye + e y dy = :
1 2 1 2 1 2
2.35.–a) Para una función positiva, su integral impropia sobre todo R3 puede
de…nirse como el valor límite de las integrales sobre esferas centradas en el ori-
gen, supuesto que ese límite sea …nito. En términos más precisos, si f : R3 ! R
es una función continua y positiva en todo R3 , se de…ne:
Z Z Z Z Z Z
f (x; y; z) dxdydz = lim f (x; y; z) dxdydz
R3 r!+1 Er
0 r ; 0 2 ; =2 ' + =2:
(x) (1 x) =
sen x
resulta:
J = :
2 sen 23
Recordando la relación entre la integral impropia simple que acabamos de cal-
cular y la integral triple impropia, se llega a la fórmula:
Z Z Z
dxdydz 2 2
=
2 2
R3 1 + (x + y + z )
2 sen 23
que completa la solución del primer apartado del problema.
b) La forma cuadrática
2 30 1
5 2 0 x
5x2 + 4xy + 2y 2 + 3z 2 = (x; y; z) 4 2 2 0 5 @ y A
0 0 3 z
es de…nida positiva, puesto que los autovalores de su matriz asociada A son
1 = 1, 2 = 3, 3 = 6.
Existe por lo tanto una transformación lineal ortogonal T de R3 tal que si
(x; y; z) = T (u; v; w), se tiene la igualdad
5x2 + 4xy + 2y 2 + 3z 2 = u2 + 6v 2 + 3w2 :
Aplicando este cambio de variables a la integral del segundo apartado del enun-
ciado y teniendo en cuenta que el valor absoluto del jacobiano de una transfor-
mación lineal ortogonal es 1, resulta:
Z Z Z Z Z Z
dxdydz dudvdw
2 + 4xy + 2y 2 + 3z 2 )
= 2 + 6v 2 + 3w 2 )
:
R 3 1 + (5x R 3 1 + (u
156
2n
k=1
Z Z Z Z Qn 1
dx1 : : : dxn 2n rn k=1 uk du1 : : : dun
2
=p lim ::: n 2
:
Rn 1 + q(x1 ; : : : ; xn ) det A r!1 Tn 2 1 + r (u1 + + un )
Haciendo el cambio de Dirichlet
Z Z Z 1 n n2 1 n Z
dx1 : : : dxn 1 r 1 d 1 Y 1 n+12 k 1 1
= p lim k (1 k)
2
d k
Rn 1 + q(x1 ; : : : ; xn ) det A r!1 0 2n 1 + r2 1 k=2 0
Z r2 n 1 Y n
1 2 d n+1 k 1
= p lim ;
det A r!1 0 1 + k=2
2 2
Z 1 n n
1 1 v 2 1 dv 2
= p n
det A 0 1+v 2
n n
1 1 n n 2 1 1 2
= p ;1 n =p n n
det A 2 2 2 det A sen 2 2
El integrando es una función positiva, lo que nos permite elegir para su cálculo
el recinto
se tiene si a 6= 0
Z Z Z
dxdydz
I(a) =
x2 +y 2 +z 2 1 x2 + y 2 + (z a)2
Z Z Z 2
sen '
= lim 2 sen2 '
d d'd
"!0+ I" [0; ] [0;2 ] + ( cos ' a)2
Z Z
2 sen '
= 2 lim 2
d' d
"!0+ I" 0 + a2 2a cos '
Z
2
= lim [log ja + j log ja j] d :
a "!0+ I"
Si 0 < a 1
Z Z Z
dxdydz
I(a) =
x2 +y 2 +z 2 1 x2 + y 2 + (z a)2
Z 1 Z a Z 1
= 2 log (a + ) d lim 2 log (a )d lim 2 log ( a) d
a 0 "!0+ 0 "!0+ a+
Z 1 Z 1 2
2 1 d
2 log (a + ) d = log (a + ) 0
=
0 0 (a + )
Z 1
2 1 a2
= log (a + ) 0
a+ d
0 a+
" #1
2
2 ( a) 2
= log (a + ) a log (a + )
2
0
" #1
2
2 ( a)
= a2 log (a + )
2
0
2 2 1 2a
= 1 a log (a + 1) + a log a :
2
158
Z a Z a 2
2 a d
lim 2 log (a )d = lim log (a ) 0
+ =
"!0+ 0 "!0+ 0 (a )
Z a e
2 a a2
= lim log (a ) 0
+ a+ d
"!0+ 0 a
" #a
2
2 ( + a) 2
= lim log (a ) a log (a )
"!0+ 2
0
" #a
2
2 ( + a)
= lim a2 log (a )
"!0+ 2
0
" #
2
2 2 (2a ) a2
= lim 2a log + a log a
"!0+ 2
3a2
= a2 log a :
2
Z Z Z
dxdydz
I(a) =
x2 +y 2 +z 2 1 x2 + y 2 + (z a)2
1+a 1 a2 1+a
= 1 a2 log + 2a = log +2
a 1 a a 1 a
Si a > 1
Z Z Z
dxdydz
I(a) = 2 2 + (z
x2 +y 2 +z 2 1 x + y a)2
Z 1 Z 1
= 2 log (a + ) d 2 log (a )d
a 0 0
8" #1 " #91
< ( a)
2
( + a) =
2
2
= a2 log (a + ) 2
a2 log (a )
a: 2 2 ;
0 0
1 2
2 a+ 1 a a+1
= a2 log +2 a = log +2 :
a a 0 a a 1
159
El integrando es una función positiva, lo que nos permite elegir para su cálculo
el recinto
f(x; y; z) ; x2 + y 2 + z 2 (a ")2 g [ f(x; y; z) ; (a + ")2 x2 + y 2 + z 2 1g
Efectuando un cambio a coordenadas esféricas y utilizando la notación
I" = [0; a "] [ [a + "; 1]
se tiene si a 6= 0
Z Z Z
dxdydz
I(a) = p
x2 +y 2 +z 2 1 x2 + y 2 + (z a)2
Z Z Z 2
sen '
= lim p d d'd
"!0+ 2 sen2 ' + ( cos ' a)2
I" [0; ] [0;2 ]
Z Z !
2 sen '
= 2 lim p d' d
"!0+ 2
+ a2 2a cos '
I" 0
Z q q
2 2
= lim (a + ) (a ) d
a "!0+ I"
Z
= lim (ja + j ja j) d :
a "!0+ I"
Si 0 < a 1
Z Z Z
dxdydz
I(a) = p
x2 +y 2 +z 2 1 x2
+ y 2 + (z a)2
Z a Z 1
2
= lim 2 d + lim 2a d
a "!0+ 0 "!0+ a+
2 3 1 2
= a + a a3 = a + :
a 3 3
Si a > 1
Z Z Z Z 1
dxdydz 2
I(a) = p = 2 2d = :
x2 +y 2 +z 2 1 x2 + y2 + (z a)2 a 0 3a
160
p
Sustituyendo a por a21 + a22 + a23 obtendremos la segunda integral que nos
pide.
La integral anterior es el potencial electrico creado en el punto a por el in-
terior de la esfera unidad uniformemente cargada. Obsérvese que el valor del
potencial en el interior de la esfera no es constante. Está claro que la esfera
no puede ser un conductor ya que si lo fuera las cargas estarían uniformemente
distribuidas en la super…cie de la esfera para minimizar la energía electrostatica.
a b 1 0
= PT P ; j det P j = 1:
b c 0 2
x = u2= ; y = v 2=
cuyo jacobiano es
4
J= u(2= ) 1 (2= ) 1
v
u = r cos ; v = rsen
resultando:
Z Z
4 r(2= )+(2= ) 1
I = lim (cos )(2= ) 1
(sen )(2= ) 1
d dr:
R!+1 [0; =2] [0;R] 1 + r2
>0 ; >0
4 0
xt Qx = xt P t DP x = x
et De
x , donde D = :
0 2
Haciendo el cambio indicado nos queda:
Z Z p
1 (4e x21 + 2ex22 )
In ( ) = 2 2
1de
x1 de
x2 ;
Mn (4e
x1 + 2ex2 )
donde
1
Mn = e 2 R2 ;
x x21 + 2e
(4e x22 ) 1 :
n
El recinto Mn es el comprendido entre dos elipses así que hacemos un cambio
de variables a coordenadas elípticas
1
e1 =
x 2 cos t 1
p1 sent
; jJ'j = p
e2 =
x 2 2 2
La integral queda:
Z Z p
1 2
In ( ) = 2
p d dt ; donde
Sn 2 2
1
Sn = ( ; t) 2 R2 ; 2
1;0 t 2 :
n
Finalmente:
Z 1
p Z 2 Z 1
p
1 2 1 2
In ( ) = 2
p d dt = p 2
d
1
n
2 2 0 2 1
n
163
2
Haciendo el cambio =u 2 d = du; se obtiene
Z 1 p Z 1p
1 2 1 u
In ( ) = p 2
d = p du;
2 n1 2 2 n1 u
y por lo tanto
Z 1 1 1
I( ) = lim In ( ) = p (1 u) 2 u du = p ;1 :
n!1 2 2 0 2 2 2
Es evidente que el límite existe solamente cuando la beta converge es decir para
< 1:
Se puede simpli…car el cálculo observando que la función subintegral es
constante en las elípses:
3x2 + 3y 2 + 2xy = ;
así que el área de los conjuntos
g 1
M n = x 2 R2 : (3x2 + 3y 2 + 2xy)
n
es:
1
An ( ) = p ( ) =) A0n ( ) = p :
2 2 n 2 2
Entonces:
Z Z p Z 1 p
1 (3x2 + 3y 2 + 2xy) 1
In ( ) = dxdy = p d :
Mn (3x2 + 3y 2 + 2xy) 1
n
2 2
2.41. Siguiendo los pasos del ejercicio anterior. Seleccionemos una sucesión
de conjuntos Mr que recubre de modo monótono al conjunto M :
1
Mr = x 2 Rn ; (xt Qx) 1 :
r
Calculemos: Z Z p
1 (xt Qx)
Ir ( ) = :: dx1 ::dxn
Mr (xt Qx)
Como Q es diagonalizable simultaneamente por congruencia y semejanza existe
un cambio de variable y = P x que es una isometría (ya que P 1 = P t y por
consiguiente jdet P j = 1 ) tal que :
xt Qx = xt P t DP x = y t Dy;
en donde D es la forma diagonal de Q. Haciendo el cambio indicado nos queda:
Z Z p
1 (y t Dy) 1
Ir ( ) = :: t Dy)
1dy1 :::dyn ; donde Mr = y 2 Rn ; (y t Dy) 1
Mr (y r
p
Siguiendo el apartado c) del ejercicio 2.35. Haciendo el cambio u = Dy
q
Z Z p Z Z 2
1 (y t Dy) 1 1 kuk
Ir ( ) = :: dy 1 ::dy n = p : : : 2 du1 : : : dun ;
Mr (y t Dy) det Q Pr kuk
164
donde
1 2
Pr = u 2 Rn ;
(ut u) = kuk 1 :
r
Teniendo en cuenta las simetrías del integrando y haciendo el cambio
n
1 1 Y 1
uk = vk2 ; k = 1::::n ; J' = vk 2 ;
2n
k=1
p Q
n 1
n Z Z 1 (v 1 + + vn ) vk 2
2 1 k=1
Ir ( ) = p lim ::: dv1 : : : dvn ;
det Q r!1 Tr 2n (v 1 + + vn )
donde
1
Tr = u 2 Rn ; (v 1 + + vn ) 1 :
r
Haciendo el cambio de Dirichlet
Z 1 p n
1 n Z
Y 1
1 1 1 1
2
d 1
n+1 k
1 1
Ir ( ) = p lim k
2
(1 k)
2
d k
det Q r!1 1=r 1 k=2 1=r
n
Y
1 n 1 n+1 k 1
= p ; ;
det Q 2 2 2 2
k=2
n 1 n n n+1
1 2 2
2 1 2
2
= p n+1 n =p n+1 n
det Q 2 2
det Q 2 2
Esta claro que el límite existe solamente cuando la beta converge es decir para
< n2
Z +1 Z +1 Z R Z b
f (ax) f (bx)
2) dx = F (x)dx = lim dx f 0 (xy)dy:
0 x 0 R!+1 0 a