2.1.5 EcuacidiferencelementEDRainvil3ed
2.1.5 EcuacidiferencelementEDRainvil3ed
2.1.5 EcuacidiferencelementEDRainvil3ed
com/
't
l
r .
.. " , .
,>
I-¡~
,l'
t
¡. Ecuaciones
diferenciales
I
I
1,
I
elementales
'1
Earl D. Rainville
,,
"
·'1
I
http://carlos2524.jimdo.com/
,
Indice general
CAPÍTULO 2
Ecuaciones de primer orden y primer grado 29
5. Soluciones generales de ecuaciones diferenciales ordina-
rias - 6. Separación de variables - 7. Sobre la forma de las
... soluciones - 8. La notación EXP U - 9. Funciones homogé-
neas - 10. Ecuaciones con coeficientes homogéneos - 11.
Ecuaciones exactas - 12. Métodos de"solución - 13. La ecua-
ción lineal de primer orden.
CAPÍTULO 3
Aplicaciones elementales 59
14. Velocidades de escape de la Tierra - 15. Ley de enfria-
miento de Newton - 16. Conversión química simple.
CAPÍTULO 4
Tópicos adicionales sobre ecuaciones de primer orden y
primer grado 69
17. Factores integrantes obtenidos por inspección - 18. De-
terminación de factores integrantes - 19. Sustitución sugerida
por la ecuación dada - 20. Ecuación de Bernoulli - 21. Coe-
ficientes lineales en las dos variables - 22. Soluciones que
involucran integrales no elementales.
7
http://carlos2524.jimdo.com/
8 Indice general
CAPÍTULO 5
Trayectorias ortogonales 95
23. Trayectorias 'ortogonales; coordenadas rectangulares - 24.
Trayectorias ortogonales; coordenadas polares - 25. Potencial
eléctrico - 26. Temperaturas de estado estable - 27. Flujo de
un fluido bidimensional en el estado estable.
CAPÍTULO 6
Funciones hiperbólicas 109
28. Definición de_las funciones hiperbólicas - 29. Fórmulas
básicas de trigonometría hiperbólica.
CAPÍTULO 7
Ecuaciones düerenciales lineales 115
30. La ecuación lineal general - 31. Independencia lineal -
32. El wronskiano - 33. Solución general de una ecuación
homogénea - 34. Solución general de una ecuación no homo-
génea - 35. Operadores diferenciales - 36. Leyes fundamen-
tales de operación - 37. Algunas propiedades de los operadores
diferenciales - 38. La enésima derivada de un producto.
CAPÍTULO 8
Ecuaciones lineales con coeficientes constantes 133
39. Introducción - 40. La ecuación auxiliar; raíces distintas
- 41. La ecuación auxiliar; raíces repetidas - 42. Definición
de EXP z para z imaginaria - 43. La ecuación auxiliar; raí-
ces imaginarias.
CAPÍTULO 9
Ecuaciones no homogéneas: coeficientes indeterminados 147
44. Construcción de una ecuación homogénea a partir de una
solución especificada - 45. Solución de una ecuación no ho-
mogénea - 46. Método de coeficientes indeterminados - 47
Solución por inspe~ción.
CAPÍTULO 10
Operadores düerenciales inversos 165
48. El cambio de la exponencial - 49. El operador l/f(D) -
50. Evaluación de [l/f(D)]ellZ - 51. Evaluación de
(D2 + a2)-:.1 sen ax y (D2 + a 2fl cos ax - 52. Evaluación
de [l//(D)]x'" - 53. Observaciones adicionales sobre el mé-
todo operacional.
http://carlos2524.jimdo.com/
Indice genera l 9
CAPÍTULO 11
La transformada de Laplace 183
54. El concepto de la transformada - 55. Definición de la
transformada de Laplace - 56. La transformada de funciones
elementales - 57. Funciones seccionalmente continuas - 58.
Funciones de orden exponencial - 59. Funciones de clase A
- 60. Transformadas de derivadas - 61. Derivadas de trans-
formadas - 62. La función gama - 63. Funciones periódicas.
CAPÍTuLO 12
Transformadas inversas 207
64. Definición de una transformada inversa - 65. Función
escalón - 66. Un teorema de convolución - 67. · Fracciones
parciales - 68. Problemas elementales de valores a la fronte-
ra - 69. Ecuaciones integrales especiales.
CAPÍTuLO 13
Aplicaciones 237
70. Vibración de una cuerda - 71. Vibraciones no amorti-
guadas - 72. Resonancia - 73 . Vibraciones amortiguadas -
74. El péndulo simple - 75. Vigas.
CAPÍTULO 14
Sistemas de ecuaciones 261
76. El método de la transformada de Laplace - 77. El método
del operador diferencial.
CAPÍTULO 15
Circuitos eléctricos y redes 269
78. Circuitos - 79. Redes simples.
CkPÍTULO 16
. ariación de parámetros y otros métodos 279
80. Variación de parámetros - 81. Solución de y" + y = f(x)
- 82. Ecuación lineal general de segundo orden - 83. Reduc-
ción de orden usando factores del operador - 84. Tercer
método: cambio de variable dependiente.
CAPÍTULO 17
La existencia de soluciones 293
85. Observaciones preliminares - 86. Naturaleza de las solu-
ciones de una ecuación particular - 87. Teorema de exis-
tencia.
http://carlos2524.jimdo.com/
10 Indice general
CAPÍTULO 18
Ecuaciones de orden uno y ':.llayor grado 299
88. Factorización del miembro izquierdo - 89. Soluciones
singulares - 90. Ecuación con discriminante e - 91. Ecuación
con discriminante p - 92. Eliminación de la variable depen-
diente - 93. Ecuación de Clairaut - 94. Eliminación de la
variable independiente.
CAPÍTULO 19
Ecuaciones especiales de orden dos 321
95. Variable dependiente perdida - 96. Variable indepen-
diente perdida - 97. Catenaria.
CAPÍTULO 20
El método de series de potencias 329
98. Ecuaciones lineales y series de potenciaos - 99. Ecuaciones
de mayor orden y grado - 100. Convergencia de las series de
potencias - 101. Puntos ordinarios y puntos singulares - 102.
Validez de las soluciones cerca de un punto ordinario - 103.
Soluciones cerca de un punto ordinario.
CAPÍTULO 21
Soluciones cerca de puntos regulares singulares 345
104. Puntos regulares singulares - 105. Ecuación indicial -
106. Forma y validez de las soluciones cerca de un punto re-
gular singular - 107. Ecuación indicia! con diferencia no
entera de raíces - 108. Diferenciación de un producto de fun-
ciones - 109. Ecuación indicial con raíces iguales - 110.
Ecuación indicial cuya diferencia de raíces es un entero po-
sitivo, caso no logarítmico - 111. Ecuación indicial cuya
diferencia de raíces es un entero positivo: caso logarítmico
- 112. Sumario.
CAPÍTuLo 22
Ecuaciones de tipo hipergeométrico 379
113. Ecuaciones que se tratarán en este capítulo - 114. Fun-
ción factorial - 115. Ecuación hipergeométrica - 116. Poli-
nomios de Laguerre - 117. Ecuación de Bessel con índice no
entero - 118. Ecuación de Bessel con índice entero - 119.
Polinomios de Hermite - 120. Polinomios de Legendre -
121. Ecuación · confluente hipergeométrica.
http://carlos2524.jimdo.com/
fndice general 11
CAPÍTULO 23
Tópicos adicionales sobre soluciones en series de potencias 391
122. Relaciones de recurrencia de muchos términos; puntos
ordinarios - 123. Relaciones de recurrencia de muchos térmi-
nos; punto regular singular - 124. Limitación de los errores
en los cálculos - 125. Solución para x grandes - 126. Ecua-
ción con un punto irregular singular - 127. Papel de las series
divergentes - 128. Solución de ciertas relaciones lineales de
recurrenCla - 129. Ecuaciones no homogéneas.
CAPÍTULO 24
Métodos numéricos 423
130. Observaciones generales - 131. Método del incremento
- 132. Método de aproximación sucesiva - 133. Mejora del
método procedente - 134. Tratamiento puramente gráfico -
135. Series de potencias.
CAPÍTULO 25
Ecuaciones diferenciales parciales 433
136. Notas sobre las ecuaciones diferenciales parciales - 137.
Ecuaciones diferenciales parciales en matemáticas aplicadas
- 138. Métodos de separación de variables - 139. Problema
sobre la conducción del calor en una plancha.
CAPÍTULO 26
Conjuntos ortogonales 449
140. Ortogonalidad - 141. Conjuntos simples de polinomios
- 142. Polinomios ortogonales - 143. Ceros de polinomios
ortogonales - 144. Ortogonalidad de los polinomios de Le-
gendre - 145. Otros conjuntos ortogonales.
CAPÍTULO 27
Series de Fourier 457
146. Ortogonalidad de un conjunto de senos y cosenos - 147.
Series de Fourier: Un teorema de desarrollo - 148. Ejem-
plos numéricos de series de F:ourier - 149. Series senoidales
de Fourier - 150. Series cosenoidales de Fourier - 151. Aná-
lisis numérico de Fourier - 152. Electrocardiografía - 153.
Artificios para obtener mayor rapidez en la convergencia.
http://carlos2524.jimdo.com/
12 Indice general
CAPÍTULO 28
Problemas de valores a la frontera 489
154. Ecuación térmica unidimensional - 155. Verificación
experimental de la validez de la ecuación térmica - 156. Tem-
peraturas superficiales que varían con el tiempo - 157. Con-
ducción térmica en una esfera - 158. Ecuación simple de
onda - 159. Ecuación de Laplace en dos dimensiones.
CAPÍTULO 29
Propiedades adicionales de la transformada de Laplace 511
160. Series de potencias y transformadas inversas - 161. La
función error - 162. Funciones de Bessel- 163. Ecuaciones
diferenciales con coeficientes variables.
CAPÍTULO 30
Ecuaciones düerenciales parciales; métodos de
transformación 527
164. Problemas de valores a la frontera - 165. La ecuación
de onda - 166. Difusión en un sólido semiinfinito - 167. Va-
riables canónicas - 168. Difusión en una losa de anchura fi-
nita --169. Difusión en un sólido en el que una de sus cuartas
partes tiende a infinito.
CAPÍTULO 1
Definiciones, eliminación de
constantes arbitrarias
(4) ou = h2
ot
2U
C
+ 02U\
ox2 oy2)'
d 2° dO 1
(5) L dt! + R d; + e i = Ew cos wt,
a2v a2 V
(6) ax2 + ay2 = 0,
(7) (d
-
dX2
2
Wy -xy-+w=O
dw
dx '
d3x dx
(8) dy3 + x dy - 4xy = 0,
(9) ~~ + 7 (~~y - 8y = 0,
d 2y d 2x_
( 10) dt2 + dt 2 - x,
of af
(11 ) x ay + y ax = nlo
Cuando una ecuación contiene una o más derivadas con respecto a
una variable particular, a esa variable se le llama variable indepen-
diente. Una variable es llamada dependiente si una derivada de esa
variable aparece en la ecuación diferencial.
En la ecuación
d2 ° dO 1
(5) L dt! + R d~ + e i = Ew cos wt
(6)
§2] Definici"ones 15
o
dx
(X2 + r) dy -2xy = 0,
EJERCICIO ORAL
2. DEFINICIONES
Una relación que define una función explicita se llama una solu-
ción de una ecuación diferencial si satisface la ecuación en el sentido
o
( 1)
_
d 21y 2
o_o 0+ d _ 6y = 4e2"
dX2 dx
+ 2e 2"
o
- 6e 2" = °
'
EJERCICIOS ORALES
1. ~: + k x = o. 2
d2 y
~! + Ri = E.
d2 x
9. x dt 2 - y dt 2 = el 10. L
( 1)
-y e- 2Z es:t
(4) -y I
- 2e-"2Z 3e sz = O.
-y" 4e- 2Z gesz
Ya que e-"2Z y eS'" no pueden ser cero, la ecuación (4) puede escribirse
de otra manera, sin los factores e-2Z y eSz, como
y 1 1
y' -2 3 = O
y" 4 9
de la cual, la ecuación diferencial
y" - y' - 6y = °
se obtiene inmediatamente.
Este último método tiene la ventaja de hacernos ver fácilmente que
la eliminación de las constantes el, e2, . . . ,en de una relación de la
forma
21
EJERCICIOS
5. cy2 =
X2 + y. SOL. 2xy dx - (y + 2X2) dy = O.
6. x = A sen (wt + (3); w un parámetro que no debe ser elimi:1ado.
2
d x
SOL. -d ' +wx = 2
O.
t-
2
d x
7. x = C1 cos wt + C2 sen wt; w un parámetro. SOL . dt2 + w x = O.
2
8. y = cx + c + 1.
2
SOL. Y = xy' + (y') + 1. 2
4. FAMILIAS DE CURVAS
U na relación que contiene un pará;:netro, así como una o ambas
coordenadas de un punto en un plano, representa una fanúlia de cUIVas,
una curva correspondiente a cada valor del parámetro. Por ejemplo, la
ecuación
http://carlos2524.jimdo.com/
(1)
o
(2) X2 + y2 - 2c (x + y) = O,
puede interpretarse como la ecuación de una familia de circunferencias,
cada una teniendo ~u centro sobre la recta y = x y pasando por el
origen. La figura 1 muestra algunos elementos o miembros de esta
familia.
y
FIGURA 1
encontramos que
(x + y)(2x dx + 2ydy) - ( X2 + y2) (dx + dy) = O.
(x + y )2
Por tanto
(3) (X2 + 2xy - y2) dx - (X2 - 2xy + y2) dy = O
http://carlos2524.jimdo.com/
.'
"
§4) '... ~'.
•
i(.5·t'·~
r
"
:t'V
,v
~1~\'·~
Famil ias de curvas
\ \V '" \ ~ •
por hf e a.8i'Ó\(h)~
Obse~~~iUos que la ecuación (3) asocia a cada punto (x, y) en el
plano, ~~a pendiente definida.
dy _ X2 + 2xy _ y2
(4)
dx X2 - 2xy _ y2'
excepto donde el denominador en el miembro derecho de (4) es nulo.
Un examen del método que hemos usado para ir de (4) a (1) muestra
que en cualquier punto (x, y) la ecuación (4) da la pendiente de la
curva de la familia (1) que pasa por el punto en cuestión.
Es conveniente ahora hacer una comprobación. Cuando el deno-
minador en el miembro derecho de la ecuación (4) es nulo, la curva
pasa por el punto que debe tener una tangente vertical. De
X2 - 2xy - y2 = O
vemos que
(5) y= (v2 - l)x
o
( 6) y= - (ti + 1 )x.
En la figura 2, las líneas rectas (5) Y (6) aparecen junto con la
familia (1). Se ve que las rectas ( 5) Y (6) cortan a los miembros
y=H+.j2)x
FIGURA 2
http://carlos2524.jimdo.com/
----------~~~~~~-----------x
FIGURA 3
x + yy' = b.
y'
Entonces
y'[l + yy" + (y') 2] - y" ( x + yy') = O
(y') 2
la ecuación diferencial deseada es
xy" - (y')3 - y' = O FIGURA 4
EJERCIC10S
En cada uno de los ejercicios obténgase la ecuación diferencial de la
familia de curvas planas descritas y bosquéjense algunos miembros repre-
sentativos de la familia.
1. Rectas que pasan por el origen. SOL. y dx - x dy = O.
2. Rectas que pasan por el punto fijo (h, k). La h y k no deben
eliminarse. SOL. (y - k) dx - (x - h ) dy = O.
3. Rectas con pendiente y la intercepción con el eje y, iguales.
SOL. y dx - (x + 1) dy = O.
4 .. Rectas con la pendiente y la intercepción con el eje x iguales.
SOL. (y') 2 = xy' _ y.
5. Rectas con la suma algebraica de las intercepcior:es iguales a k.
SOL. (xy' - y) (y' - 1) + ky' = O.
~. \
l' http://carlos2524.jimdo.com/
CAPÍTULO 2
Ecuaciones de
priIller orden y priIller grado
6. SEPARACIóN DE VARIABLES
2 dx _ y dy -
x y+3-'
°
o
(4) x2 dx - [ 1 - y +3 3] dy = O.
ahora podemos
http://carlos2524.jimdo.com/
y llegamos a ( 7)
El problema de cambiar una forma de la solución en otra, es un
problema que surge frecuentemente cuando dos o más personas resuel-
ven la misma ecuación diferencial y desean comprobar los resultados
obtenidos. A menos de que la solución obtenida sea destinada a un uso
específico, hay poca razón para preferir una forma u otra, excepto
por consideraciones de compactación, simetría u otras cualidades esté-
ticas. Esta preferencia es esencialmente de carácter individual. La
sección 7 contiene explicaciones más amplias sobre la forma de las
soluciones.
x + y = -1
1 - xy ,
o
(11 ) xy - x - y - 1 = O.
\
\
\
\ , .....
'--
'\
IO,-ll~
FIGURA 5
http://carlos2524.jimdo.com/
Una diferencia entre las soluciones (10) Y (11) puede verse obser-
vando que una máquina computadora, dada la ecuación diferencial (8)
Y buscando una solución que pase por el punto (O, - 1) , dibujará
solamente la rama ' izquierda de la curva de la figura 5. La barrera
( asíntota) en x = 1 hace que la máquina no compute la otra rama de la
hipérbola (11).
Consideremos la ecuación
( 1) 2x(y + 1) dx - y dy = O.
Separando las variables en la ecuación (1) obtenemos
podemos ver que la solución (3) es útil para el rango y > - 1, mientras
que la solución (4) es útil para y < -1. En fin, dos caminos son
posibles para escribir la solución combinando las ventajas de (3) con
las de (4).
Ya que
_ C2 dy _ dy
dln[c"l (y + 1)] - C2 (y + 1) - Y + l'
la solución general de la ecuación (1) puede escribirse en la fonna
(5)
Ajustando la solución (5) a nuestro propósito inmediato, ponemos
x = 0, y = - 2 Y encontramos que
o= - 2 - In ( - C2) ' .
o ..
X2 = y +2- In [- (y + 1 )].
En la ecuación (5), C2 se hará negativa cuando una condición a la
frontera requiera que y < - 1 Y positiva cuando se necesite que y > - 1.
De la ecuación (5) no es difícil obtener la relación
(6) (y + l)e Z
' = Cs eY
como otra fonna de la solución de la ecuación diferencial (1).
En este libro dejaremos usualmente las soluciones en forma similar a
( 5) o (6). Dondequiera que una respuesta tal como la ecuación (3)
se encuentre está implícito que el que usa esta ecuación (3) deberá ajus-
tarla para su propósito práctico. El resultado deberá verificarse directa-
mente en la ecuación diferencial de tal fonna que la validez de la
transición de una fonna de solución a otra, no intervenga en la de-
mostración.
8. LA NOTACIÓN EXP U
La función e" se usa en nuestro trabajo frecuentemente, y . algunas
veces el exponente u es complicado. Un simple ejemplo tal como é 1r y2
puede ser indeseable en la impresión tipográfica, porque los exponentes
http://carlos2524.jimdo.com/
EJERCIClOS
En los ejercicios 1-26 obténgase la solución general.
1. (4 + x)y' = y3. SOL. 2y2 In[e(4 + x)] =
-1.
2. exp (y2) dx + x 2y dy O. = SOL. x exp (_y2) + 2 =ex.
3. cos x cos y dx + sen x sen y dy = O. SOL. sen x = ecos y.
~~~=b~ ~ ~=~
http://carlos2524.jimdo.com/
8. CZ(y - 1) dx + 2(eZ + 4) dy = O.
SOL. (y - 1)2 (e" + 4) = c2.
9. dr = e(r sen O dO - cos O dr). SOL. r(l + ecosO) = c.
10. (xy - x) dx + (xy.+ y) dy = O.
SOL. (y - 1) exp (x .+ y) = c(x + 1).
11. (y +1) dx = 2xy dy. SOD. e2Y = cx(y + 1)2.
12. x 2 dx+y(x - 1) dy=O.
13. (xy.+ x) dx = (X2y2 X2 + + y2 + 1)dy,.
SOL. In (X2 + 1) = y2 - 2y + 41n [c( y + 1)].
14. X cos 2 y dx + tan y dy = O. SOL. X2 + tan 2 y = c2. •
15. x yy' = eY.
2
SOL. x(y + 1) = ( 1 + cx)eY.
16. tan y dy = sen x dx.
2 3 SOL. cos X - 3 cos x = 3 (tan y - y + e).
3
17. y' = cos2 x cos y. SOL. 4ln (sec y + tan y) = '2x + sen 2x + c.
· 18. y' = y sec X. SOL. Y = c(sec x + tan x).
19. dx = t (1 + t 2) sec 2 x dt. SOL. 2x + sen 2x = e + (1 + t 2) 2.
20. (e 2Z + 4) y' = y. SOL. ya(1 + 4e- 2") = c2 •
21. a dJ3 + f3 da + aJ3(3 da + dJ3 ) = O. SOJ.,. caJ3 = exp (-3a - J3).
22. (1 + lnx) dx + (1 + In y) dy = O. sOL. xlnx + ylny = c.
23. x dx - ...¡ a 2 - X2 dy = O.
SOL. Y - e = - ...¡ a2 - x 2, la mitad inferior de la
circunferencia x2 + (y - e ) 2 = a2 •
24. x dx + va 2
- X2 dy = O.
SOL. Y - e = ...¡ a 2 - x 2 , la mitad superior de la
circunferencia X2 + (y - e) 2 = a2 •
y+c
25. ' a 2 dx = x ...¡ X2 - a 2 dy. SOL. x = asee--
a
26. y In x In y dx + dy = O. SOL . x In x + In In y = x + c.
En los ejercicios 27-33, obténgase la solución particular que satigaga las
condiciones a la frontera indicadas.
dr
27. dt= -4rt ; cuando t = 0, r = ro. SOL. r = rQ exp ( - 2t2).
28. 2xyy' = 1 + y2; cuando x = 2, Y = 3. SOL . y2 = 5x - 1.
29. xyy' = 1 + y2; cuando x = 2, Y = 3. SOL. 5X2 - 2y2 = 2.
30. y' = x exp (y - X2) ; cuando x = 0, y = O.
SOL. 2e-1I = 1 + exp ( _X2).
31. xy2 dx + eZdy = O ; cuandox~ oo,y~f.
SOL. Y = e"j(2CZ - x - 1).
32. (2a2 - r 2 ) dr =
r a sen OdO; cuando O = 0, r = a.
dv
33. v - = g; cuando x = XQ, V = VQ. SOL. v2 ~ VQ2 = 2g(x - Xo).
dx
http://carlos2524.jimdo.com/
37
9. FUNCIONES HOMOGÉNEAS
Los polinomios en los cuales todos los términos son del mismo grado,
tales como
x? - 3xy + 4y2,
(1) X3 + ya,
x y + 7y
4
5, _
(2)
EJERCICIOS ORALES
7. exp G). 3y
8. tan-o
x
X
(X2 _ y2)i
15. Y2 tan-. 16. (X2 _ y2) t'
y
a
17. a - 4b'
+ 4b 18. In:: .
y
19. x In x - y In y. 20. xlnx - xIny.
(2) dy
dx
+ g(2:)
x
= °
Esto sugiere la introducción de una nueva variable v al poner y = vx.
Entonces (2) será
dv
, (3) x - + v + g (v) = 0,
dx
en la cual las variables son separables. Podemos obtener la solución de
(3) por el método de la sección 6, insertando y/x por v, y así llegamos
a la solución de (1). Hemos mostrado que la sustitución y = vx trans-
formará la ecuación (1) en una ecuación en v y x en la que las variables
son separables.
El método anterior hubiera tenido el mismo éxito de haber usado
x = xy para obtener de (1) una ecuación en y y v. Ver el ejemplo b)
más adelante.
(X2 -
2
x v + 2 2
X V ) dx - x 2v(vdx + xdv) = 0,
de la cual el factor X2 puede suprimirse inmediatamente. Hecho esto
debemos resolver
http://carlos2524.jimdo.com/
.,
(1 - V +v 2
) dx - V(V dx + X dx) = O,
O
(1 - v) dx - xv dx = O.
De aquí, separando las variables obtenemos
dx + v dv = O.
x v - 1
Entonces de
dX+[l +_l_Jdv =O
x v- 1
la solución será
In x + v + In (v - 1) = In e,
o
x(v - l)e V = c.
(y - x)exp(~) = c.
EJEMPLO b): Resolver la ecuación
(5)
Otra vez, los coeficientes en la ecuación son homogéneos y de grado
dos, Podríamos usar y = vx pero la simplicidad relativa del término
dx en (5) sugiere que pongamos
x = vy.
Entonces dx = V + y dv, y la ecuación (5) se reemplaza por
dy
y4(~ + 1) = C;
esto es,
(7)
Ya que el miembro izquierdo de la solución (7) no puede ser negativo,
podemos, por motivos de simetría, cambiar la constante arbitraria por
Cl \ escribiendo
y2(2x2 + i) = C.1 4 •
Para el estudiante que va a resolver la ecuación (5), vale la pena
usar y = vx. Este método nos lleva directamente a la ecuación
(v 3 + 2v) dx + x( v + 1) dv =
2
O.
Frecuentemente en ecuaciones con coeficientes homogéneos, no tiene
importancia el uso de y = vx o x = vy.
EJERCICIOS
En los ejercicios 1-19 obténgase la solución general.
: 1. (x - 2y) dx + (2x + y) dy = O.
SOL. In (X2 + y2) + 4 Aretan (y/x) = e.
2. 2(2x 2 + y2) dx - xydy = O. SOL. x 4 = e2 (4x2 + y2).
3. xydx - (X2 + 3y2) dy =O. SOL. X2 = 6y2In (y/e).
4. x 2y' = 4X2 + 7xy + 2y2. SOL. X2(y + 2x) = e(y + x).
5. 3xydx + (X2 + y2) dy = O.
6. (x - y) (4x + y) dx + x(5x - y) dy = O.
SOL. _x_(_y_+-x.)-=--= -c (y - 2xJ.
7. (5u - u) du + (3u - 7u) du =O. SOL. (3u + U)2 = e(u - u).
8. (x 2 + 2xy - 4y2) dx - (X2 - 8xy - 4y2) dy = O.
SOL. X2 + 4y2 = e (x + y) .
9. (x 2 + y2) dx - xydy = O. SOL. y2 = 2x2In (x/e).
10. v 2 dx + x(v - 4x) du = O. SOL. xu 2' = e(u ~ 2x).
11. (2x + y)2dx = xydy. SOL. x (x + y) = e exp (y/x).
3
(A)
43
x-
aF + y-
aF = kF.
ax ay
y
oN 02F
ox - oxoy·
Otra vez del cálculo sabemos que
aF2
02F
oyox - oxoy'
siempre que esas derivadas parciales sean continuas. Por tanto, si (1)
es una ecuación exacta, entonces
oM oN
(4) ay =a;.
De este modo, para que (1) sea exacta es necesario que (4) se satisfaga.
Mostraremos ahora que si la condición (4) se satisface, entonces
( 1) es una ecuación exacta. Sea 9' (x, y) una función tal que
09' = M
OX •
( 8)
Y
oF
(9) -
oy = N = x3 + 2y.
Queremos determinar la F de la ecuación (8). Integrando ambos
miembros de ( 8) con respecto a x, manteniendo y constante se obtiene
( 10)
donde la constante arbitraria usual en la integración indefinida, es ahora
necesariamente una función T (y), hasta aquí desconocida. Para deter-
minar T (y), usamos el hecho de que la función F de la ecuación (10)
debe también satisfacer la ecuación (9).
En consecuencia
x 3 + T'(y) = x 3 + 2YJ
T'(y) = 2y.
En·la obtención de T(y) no se necesita la constante arbitraria, ya que
una se introduce en el miembro derecho de la solución F = c. Entonces
T(y) = i,
y de (10)
F = X3 y - 3X2 + y2.
Finalmente la solución de la ecuación (7) es
X3 y - 3X2 + y2 = C.
Debido a que (13) es más simple que (12), y por motivos de variedad
principiamos la determinación de F de la ecuación (13) . .
Inmediatamente de (13)
F == _ix2y2 - 2xy + Q(x),
donde Q(x ) se determinará de (12). La última da
_xy2 - 2y + Q'(x) = 2X3 - xy2 - 2y + 3,
Q'(x) = 2X 3 + 3.
Por tanto
Q(x} = i X4 + 3x,
y la solución deseada de (11) es
- ix2y2 - 2xy + i X + 3x =
4 ie,
o
x4 - x 2i - 4xy + 6x = c.
EJEMPLO e): Resolver la ecuación del ejemplo a) por inspección.
Supóngase que hemos examinado
(7) 3x(xy - 2) dx + (x 3 + 2y) dy = O
y hemos encontrado que es exacta. Entonces, podemos escribir la solu-
ción general por inspección; esto es, por una cuidadosa observación del
miembro izquierdo de (7), encontraremos una función de la cual (7)
es la diferencial total.
Primero, el término 3x2y dx sugiere la diferencial de x 3y. De aquí
buscamos en (7) el término necesario compañero de x 3 dy y agrupamos
los dos términos. Cualquier término tal como - 6x dx, que contiene:
solamente una variable, es una diferencial exacta tal como se ha puesto.
Por tanto, acomodaremos la ecuación (7) en la forma
(3x2y dx + x dy)3
- 6x dx + 2y dy = O,
de la cual se sigue que
X3y - 3X2 +i = c.
EJEMPLO d): Resolver la ecuación
(11) (2x 3 - xy2 - 2y + 3) dx - (X2y + 2x) dy = O
del ejemplo' b~ por inspección.
El agrupamiento
http://carlos2524.jimdo.com/
EJERCICIOS
Examínense cada una de las siguientes ecuaciones para saber si son exac-
tas y resuélvase la ecuación. Las ecuaciones que no son exactas pueden,
por supuesto, resolverse por los métodos discutidos en las secciones prece-
dentes. .
1. (x + y) dx + (x ,- y) dy =O. SOL. .+
X2 2xy - y2 =c.
2. (6x + y2) dx + y(2x -3y) dy = O. SOL. 3X2 + Xy2 - y3 ~ I c.
3. (2xy - 3X2) dx + (X2 + y) dy = O. SOL. X2y - X3 + 'b2 =c.
4. (y2 - 2xy + 6x) dx - (X2 - 2xy + 2) dy = O.
SOL. Xy2 - X2y + 3X2 - 2y = c.
5. (2xy + y) dx + (x 2 - x) dy O. = SOL. Y- =CX(X - 1)-3.
6. (x - 2y) dx + 2(y - x) dy =
O. SOL. X2 + 2y2 = 4xy + c.
7. Haga el ejercicio 6 por otro método.
8. (2x - 3y) dx + (2y - 3x) dy = O. SOL. X2 + y2 = 3xy + c.
9. Hacer el ejercicio 8 por otro método.
10. v(2uv 2 - 3) du + (3U 2V 2 - 3u + 4v) dv = O.
V(U V - 3u + 2v ) = c.
2 2
SOL.
11. (cOS 2y - 3X2y2) dx +
(COS 2y - 2x sen 2y - 2x3y ) dy = O.
SOL. i sen 2y + x cos 2y - X3y2 = c.
12. (1 + y2) dx + (x 2y + y) dy = O.
SOL. 2 Arctan x + In (1 + y2) = c.
13. (1 + y2 + xi) dx + (X2y + y + 2xy) dy = O.
SOL. 2x + y2(1 + x)2 = c.
14. (w 3 + wz 2 - z ) dw + (Z3 + w 2z - w) dz = O.
SOL. (w + z2)2 = 4wz + c.
15. (2xy - tan y) dx + (X2 - x sec 2 y) dy = O.
SOL. x2y - x tan y = o.
16. (cosxcosy - cotx) dx - senxsenydy O. =
, SOL. senxcosy = In (csenx) .
17. (r + senO - cos O) dr + r( sen O + cosO) dO = O.
SOL. r + 2r (sen O - cos O) = c.
2
yexp () Pdx)
La solución de (8) es
(9)
y debe ser comprobado particularmente, ya que la verificación de este
resultado es muy fácil.
De (9) obtenemos (8 ) por diferenciación. Entonces, la ecuación
diferencial .original se sigue de (8 ) por un simple ajuste. Así (9) es
una solución de la ecuación original.
( 10) dx + (~ -
Y
1) x dy = - '¿ dy.
y'
Ahora
J(~ - 1) dy = 3 In y - y + el,
así que un factor integrante de la ecuación ( 10 ) es
exp (3lny - y) = exp(~ lny ) • e-V = [exp (lny)]3. e-Y = ie-v.
~plicando este factor integrante a la ecuación (10) llegamos a la
ecuación exacta
y3e-Y dx + y2 (3 - y) e-Yx dy = - 2y2e-Y dy,
de la cual obtenemos
xy3e-V = - 2 S y2 e-Y dy
= 2ie-v + 4ye-V + 4e-v + c.
De este modo pasamos a escribir la solución
xy3 = 2y2 + 4y + 4 + c.eY.
EJERCIClOS
En los ejer6cios 1-25, encuéntrese la solución general.
1. (x 5 + 3y) dx - x dy = O. SOL. + cx
_2y = x 5 3
•
2. 2(2xy + 4y - 3) dx + (x + 2)2dJ' = o.
SOL. y = 2(x + 2)-1 + c(x + 2) -4.
http://carlos2524.jimdo.com/
6. u dx + (1 - 3u) X du =
3u duo SOL. XU = ceS\< - u - !.
7. y' = X - 4xy. Resolver por dos métodos.
SOL. 4y = 1 + c exp ( -2x2).
8. y' = csc X + Y cot X . SOL. Y = c sen cos x.
X -
9. y' = cse X - Y cot X. SOL. y sen X = X + ·c.
10. (2xy + X2 + x 4 ) dx - (1 + X2) dy = O.
SOL. Y = (1 + X2) (c + X -,. Arctan x) .
11. (y - cos 2 x) dx + cos x dy = O.
SOL. y(secx + tan x) = c + x - cosX.
12. y' = x - 2y cot 2x. SOL. 4y sen 2x = c + sen 2x - 2x cos 2x.
13. (y - x + xycotx) dx + x.dy = O.
SOL. xy sen x = c + sen x - x cos x .
14. dy -my = clem:e, donde Cl y m son consta."ltes.
dx
SOL. Y = (C1X + C2) em:e.
15. dy _ m2Y = cl.emJ:e, donde C1, 11'!0¡" m2 son constantes y mI =1= m2.
dx
Cl
SOL. Y = cs em,,,, + C2e"""', donde C3 r= _----'0._
m1- m2
16. v dx + (2x + 1 - vx) dv = O. SOL. xv 2 = V + 1 + ceV •
17. X(X2 + l)y' + 2y = (X2 + 1)3.
SOL'. x2y = i{x2 + 1)3 + C(X2 + 1).
18. 2y (y2 - x) dy = dx. SOL. X = y2 - 1 + cexp (_y2).
19. (1 + xy) dx - ( 1 + X2) dy = O. SOL. Y = x + c(1 + x2)i.
20. 2y dx = (X2 - 1) (dx - dy) .
SOL. (x - l)y = (x + l)[c + x - 2Jn (x + 1)].
21. dx - (1 + 2x tan y) dy = =
O. SOL. 2x cos2 Y Y + c + sen y cos y.
22. (1 + cos x) y' = sen x (sen x + sen x cos x - y) .
, SOL. Y = (1 + cosx) (c + x - sen x) .
23. y' = 1 + 3y tan x. SOL. 3y cosa x =
c + 3 sen x - sen s x.
24. (x + a)y' = bx - ny; a, b, n son constantes con n =1= O, n =1= -1.
SOL. n(n + l)y = b(nx - a) + c(x + a) -".
25. Re~uélvase la ecuación del ejercicio 24 para los casos n = O Y
n = -1.
SOL. Si n = O, Y = bx + c - ab In (x + a).
Si n =- 1, Y =
ab + c (x + a) + b (x + a) In (x + a).
26. E:tl la forma estándar dy + Py dx =
Q dx póngase y vw, de este =
modo obténgase
w(dv + Pv dx) + v dw = Q dx.
Entonces, escogiendo primero v tal que
dv + Pv dx = O
http://carlos2524.jimdo.com/
EJERCICIOS COMPLEMENTARIOS
11. ~~ = sec 2
x sec3 y. SOL. 3 tan x + e = 3 sen y - sens y.
http://carlos2524.jimdo.com/
I
http://carlos2524.jimdo.com/
56
I
I
I
----~--------~~--------~---x
FIGURA 6
\
http://carlos2524.jimdo.com/
,
CAPITULO 3
Aplicaciones ele:rnentales
dv k
a - -
---¡¡-ro
La aceleración es negativa porque la velocidad es decreciente. Por lo
tanto, la constante k es negativa. Cuando r = R, entonces a= -g.
la aceleración de la gravedad en la superficie de la Tierra.
Por tanto,
k
- g = R2'
de lo que se obtiene
dv dr dv dv
a=-=--=v-
dt dt dr dr'
de tal forma que la ecuación diferencial para la velocidad se ve que
tiene la forma
v dv = _ gR~
(1 )
dr r2 •
Se aplica el método de separación de variables a la ecuación (1) Y se
obtiene inmediatamente la solución
2 R2
v 2 =-g- + c.
r
Supóngase que la partícula deja la superficie de la Tierra con una
velocidad va. Entonces v = va cuando r = R, de estas condiciones es
fácil determinar el valor de la constante C, que es
e= v0
2
- 2g R.
Entonces, una partícula proYectada en una dirección radial hacia
afuera de la superficie de la Tierra, con una velocidad inicial Va, viajará
con una velocidad v dada por la ecuación
2 R2
(2) v
2
= ~+ v 0
2
- 2gR.
r
Es de considerable interés determinar cuándo la partícula escapará
a la atncción gravitacional de la Tierra. En la superficie de la Tierra,
http://carlos2524.jimdo.com/
62
(2) x = XOe-kt
Vamos ahora a estipular otra condición más, que nos sirve para
determinar k. Supóngase que se sabe que al terminar el primer medio
minuto, esto es en t = 30 seg., dos tercios de la cantidad original Xo se
han convertido. Vamos a determinar cuánta sustancia no convertida
queda en t = 60 seg.
Cuando dos tercios de la sustancia se han convertido, un tercio
permanece sin convertirse. Por lo tanto, x =
!xo cuando t 30. La =
ecuación ( 2) nos da ahora la relación
!xo = xoe- 30k
~ In 3. Entonces si t está
de la cual se encuentra fácilmente que k vale 30
medida en segundos, la cantidad de sustancia no convertida está dada
por la ecuación
(3) x = Xo exp ( - ~t In 3) .
30
En t = 60
x = xo exp ( - 2 In 3) = Xo ( 3 ) = 2.xo.
-2
9
EJERCICIOS
68
FIGURA 8
CAPÍTULO 4
Tópicos adicionales sobre
ecuaciones de
pri:rner orden y primer grado
69
http://carlos2524.jimdo.com/
(9)
( 10)
~d(~) - d (xy) = O
y y x2y2
de la cual la solución es
t(~y + x~ =~.
Finalmente, podemos escribir la solución deseada de la ecuación
(7) como
xa + 2y = cxi.
EJEMPLO e): Resolver la ecuación
3x2y dx + (y4 - x 3 ) dy = O.
Ya que dos términos en los coeficientes de dx y dy son de tercer
grado y el otro coeficiente no es de tercer grado, reagrupamos los
ténninos para obtener
http://carlos2524.jimdo.com/
d (~) = v du - u dv
V v2
Por tanto, dividimos cada término de nuestra ecuación entre y2 y
obtenemos
d(X3) - x 3 dy
:JI~....'..:..:......!...-2---':~...L. + y2 dy = 0,
y
o •
d (;3) + y2 dy = O.
EJERCICIOS
aM au aN au
u-+M-=u-+N - ,
ay oy ax ax
o
(3 ) u ( aM _ aN) = N ou _ M au .
ay ax ox ay
Además, tomando el argumento anterior a la inversa, puede verse
que si u satisface la ecuación (3), entonces u es un factor integrante
de la ecuación ( 1). Hemos "reducido" el problema de resolver la ecua-
ción diferencial ordinaria (1) al problema de obtener una solución
particular de la ecuación diferencial parcial (3).
No se ha ganado mucho ya que aún no hemos desarrollado métodos
para atacar una ecuación como la (3). Por tanto, volvemos atrás, al
problema de resolver ecuaciones diferencial ordinarias y restringimos
a u a que sea una función de una sola variable.
au au
Primero, sea u una función de x solamente. Entonces ;- = O Y -
ux ay
du D
se transf orma en ax' e aqm ' (3) d
se re uce a
u (OM _ aN) = N du
ay ax dx'
o
(5) -1 (aM
---aN)' =f(x)
N ay ox '
entonces el factor integrante deseado es u = exp (S f (x) dx).
Por un argumento similar, suponiendo que u es una función de y,
sólo llegamos a la conclusión de que si
(6) -1 (aM
---aN) =g(y),
M oy ax
entonces un factor integrante para la ecuación (1) es u = exp ( - S
g(y) dy).
Nuestros dos resultados están expresados en las siguientes reglas.
,I
;
/
http://carlos2524.jimdo.com/
se obtiene de inmediato.
http://carlos2524.jimdo.com/
1 (oM oN) _ x + y+ 1 _ 1
M Ty - ox - - y(x + y +l ) - - y '
Por tanto exp (In y ) = y es el factor integrante deseado para (9).
Usando el factor integrante, escribimos '
(xy2 + y3 + y2 ) dx + (X2y + 3xi + 2xy) dy = O,
o
(xy2 dx + x2y dy) + (i dx + 3xy2dy) + (i dx + 2xy dy) = O.
Se encuentra entonces que la solución de (9) es '.
ix2y2 + xy3 + X'y 2 = ic,
o
xi(x + 2y + 2) = c.
1. (OM _ oN) = x + 2y - 1 = 1.
N oy OX x + 2y - 1
En consecuencia eft es un factor integrante para (10) . Entonces
(xye'" + y e"') dx +
2
(xe'" + 2ye'" - eft ) dy = O
es una ecuación exacta. Agrupando los términos de la siguiente manera:
http://carlos2524.jimdo.com/
EJERCICIOS
1. (X2 + y2 + 1) dx + x (x - 2y) dy = O.
SOL.
I • X2 _ y2 + xy - 1 = cx.
2. 2y(x 2 - y + x) dx + (X2 - 2y) dy = O. SOL. Y(X2 - y ) = ce-2".
3. y (2x - y + 1) dx + x(3x - 4y + 3) dy = O.
SOL. xy3( X - y + 1) = c.
4. y (4x + y ) dx - 2 (x 2 - y ) dy = O. SOL. 2X2 + xy + 2ylny = cy.
5. (xy + 1) dx + x (x + 4y - 2) dy = O.
SOL. xy + In x + 2y2 - 2y = c.
6. (2y2 + 3xy - 2y + 6x) dx + x (x + 2y - 1) dy = O. '
S OL. X2(y2 + xy - y + 2x) = c.
/ 7. y(y + 2x - 2) dx - 2(x + y) dy = O. SOL. y( 2x + y) = ce".
8. y 2dx + (3xy + y2 - 1) dy = O. SOL. y2(y2 + 4xy - 2) = c.
9. 2y(x + y + 2) dx + .(y2 - X2 - 4x - 1) dy = O.
SOL. X2 + 2xy + y2 + 4x + 1 = cy.
10. 2(2y2 + 5xy - 2y + 4) dx + x( 2x + 2y - 1) dy = O.
SOL. X4(y2 + 2xy - y + 2) = c.
11. 3(X2 + y2) dx + X(X2 + 3y2 + 6y) dy = O.
SOL. X(X2 + 3y2) = ce-I/.
12. y (8x - 9y) dx + 2x( x - 3y) dy = O. SOL. x3y(2x. - 3y) = c. ...-
13. Hacer el ejercicio 12 por otro método.
14. y(2x 2 - xy+ 1) dx + (x - y) dy O. =
y(2x - y) = c exp (_X2).
SOl;.
15. El teorema de Euler (ejercicio 34, página 42) sobre funciones homo-
géneas, establece que si F es una función homogénea de grado k en
x y 'YJ entonces
aF aF
x ax + y ay = kF.
En la ecuación
(2) (1 + 3x sen y) dx - X2 cos y dy = 0,
la presencia de sen y. y su diferencial cos y dy, Y el hecho de que y
aparece en la ecuación solamente de esta manera, nos lleva a poner
sen y = w y obtener la ecuación diferencial en w y x. (Véase nuestro
ejemplo b.)
( v - 1) dv + (v + 2) dy = O.
Ahora las variables pueden separarse. De la ecuación en la forma
v - 1
- - 2 dv
v+ + dy = 0,
obtenemos
y entonces
. v - 3ln (u + 2) + y + e = O.
Pero v = x + 2y, así nuestro resultado final es
x + 3y + e = 3 In (x + 2y + 2).
EJEMPLO b): Resolver la ecuación
(2) (1 + 3xseny) dx- x 2
cosydy = O.
. Haciendo sen y = w. Entonces cos y dy = dw y (2) se convierte en
(1 + 3xw) dx - X2 dw = 0,
una ecuación lineal en w . De la fonna estándar
http://carlos2524.jimdo.com/
3 dx
dw - - wdx=-
X Jf"
se ve que un factor integrante es
exp ( - 3 In x) = x-3
La aplicación del factor integrante produce una ecuación exacta
x- 3 dw - 3x--4w dx = x-5 dx,
de la cual obtenemos
o
4xw = cx4 - 1.
En consecuencia (2) tiene la solución
4x sen y = cx 4
- 1.
de la cual
(1 - n) y-n dy = dz.
De este modo la ecuación en x y z es
dz + (1 - n) pz dx = (1 - n) Q dx,
una ecuación lineal en la forma estándar. De aquí que cualquier ecua-
ción de Bemoulli puede resolverse con la ayuda del cambio anterior
de la variable dependiente (a menos que n = 1, en la que la sustitución
no es necesaria).
http://carlos2524.jimdo.com/
EJERCICIOS
dy -
5. dx = sen (x + y). SOL. + e = tan (x + y) - sen (x + y) .
X
Considérese la ecuación
http://carlos2524.jimdo.com/
o
(vz - V)3 = c(vz + v).
Ahora vz = u, así nuestra solución aparece como
( u - V)3 = c(u + v) .
http://carlos2524.jimdo.com/
EJERCICIOS
2. (x - 4y - 9) dx + (4x +y- 2) dy = O. x _ 1
S OL. In [( x - 1)2 + (y + 2)2] - 8 Arctan-- y+2
= c.
3. (2 x - y ) dx + (4x + y - 6 ) dy=0 .
SOL. (x +y-3)s=c ( 2x+y-4)2.
4. (x - 4y - 3) dx - (x =
6y - 5) dy = O.
SOL . (x - 2y - 1)2 = c (x - 3y - 2).
5. ( 2x + 3y - 5) dx + (3x - y - 2) dy = O. Resolverla por dos mé~
todos.
6. 2dx + ( 2x - y +
3) dy = O. Usar un cambio de variable.
SOL. y + e = - In ( 2x - y + 4) .
7. R esuélvase la ecuación del ejercicio 6 usando el hecho de que la
ecuación es lineal en x.
8. (x - y + 2) dx + 3 dy O. =
SOL. X + e 31n (x - y + 5) . =
9. R esuelva el ejercicio 8 por otro método.
10. (x + y - 1) dx + ( 2x + 2y + 1) dy = O.
SOL. X + 2y + e = 3 In (x + y + 2) .
11. (3x + 2y + 7) dx + ( 2x - y) dy = O. Resolverla por dos métodos.
12. (x - 2) dx + 4 (x + y - 1) dy = O.
S OL. 2 (y + 1) = - (x + 2y) In [c (x + 2y)]'
13. (x - 3y + 2) dx + 3 (x + 3y - 4) dy = O.
In [(x - 1) 2 + 9 (y - 1)2] -2 Arctan 3 ty -=- 11 ) = c.
I 14. (6x - 3y
S OL.
+ 2) dx -
( 2x - y - 1) dy = O.
SOL. 3x - y +
e = 5 In ( 2x - y + 4 ).
15. (9x - 4y + 4) dx -
( 2x - y 1) dy = O. +
SOL. y - 1 = 3 (y - 3x - 1) In [e (3x - y + 1)].
16. (x +3y - 4) dx (x + +
4y - 5) dy = O.
S OL. Y - 1 = (x +
2y - 3) In [e (x 2y - + 3 r
17. (x +2y - 1) dx - (2x +
y -5 ) dy = O.
S OL . (x - Y - 4)3 = c (x y - + 2) .
18. (x -1 ) dx - (3x - 2y - 5 ) dy =0.
SOL. ( 2y - x 3)2 = e (y - x+ + 2) .
En les ejercicios 19-2 2 obténgase la solución par ticular indicada.
19. ( 2x - 3y + 4) dx + 3 (x - 1) dy O; cuando x =
3,y 2. = =
x - 1
SOL. 3 (y - 2) = -2 (x - 1) In - 2- '
20. La ecuación del ejercicio 19, pero con la condición: cuando x = -1,
1- x
y = 2. SOL . 3 (y - 2) = -2 (x - 1) In - - .
2
21. (x +y- 4 ) dx - (3x - y - 4 ) dy = O; cuando x = 4, y = 1.
x - y
SOL. 2 (x + 2y -
6) = . 3 (x - y) In - - ·
3
22. La ecuación del ejercicio 21, pero con la condición : cuando x = 3,
y = 7. SOL. Y - 5x + 8 = 2 (y - x) In y - x .
4
http://carlos2524.jimdo.com/
(3) 2
fer x = _1
v 7i
f'" exp ( -[32) df3
o
* Véase E. D . Rainville, Unified Calculus and Analytic C eom etry. Nueva York:
The M acmillan. 1961 , páginas 605-607.
http://carlos2524.jimdo.com/
~
xy = e +
J "' sen w
o
x2n
--dw,
w
o y = cx-
1
+ ~o ( - 1) n (2n + 1) (2n + 1) !
3. y' = 1 - 4x 3 y.
4. (y cos 2
X - X sen x) dx
SOL. Y = exp ( _ x 4 ) [e
+ sen x cos x dy = O.
+ J: exp (w) df3 ]
MISCELÁNEA DE EJERCICIOS
CAPÍTULO 5
Trayectorias ortogonales
97
-------+~--~--+44_------x
FIGURA 9
EJERUCIOS
8 8
FIGURA 10 FIGURA 11
http://carlos2524.jimdo.com/
100
I /'
I /
__,I____-+I____~~------+_----+_--x
\
"
FIGURA 12
o '
sen 8 dr - r( 1 + cos 8) d8 = O.
Separando las variables, usamos la forma
dr _ ( 1 + cos 8) d8 = O
r sen 8 '
o
dr _ ese 8 d8 - cos 8 d8 = O.
r sen 8
De aquí
1n r - In (csc 8 - cat 8) -In sen 8 = In b,
así, las trayectorias ortogonales de la familia original de cardioide:.
serán
r = b sen 8 ( esc 8 - cot 8 ) ,
o
r = b ( 1 - cas 8) .
EJERClCIOS
(2)
métodos que hemos estudiado aquí para determinar las líneas de flujo,
cuando se conocen las curvas equipotenciales.
FIGURA 13
http://carlos2524.jimdo.com/
de aquí se concluye que las líneas de flujo son círculos con centros sobre
el eje y y que pasan a través de los puntos ( 1,0 ) y ( -1 ,0 ). Ver figura 13.
T = TI + T 2 -a TI Y
Arctan-.
X
FIGURA 14
(1 ) 1/1 (x, y) = c
107
o 8=0
FIGURA 15
(1 - ~) cos e dr - (r + De sen dO = O,
o
(4) (r2 - 1) cos Odr - r( r2 + 1) sen OdO = O.
Entonces la ecuación diferencial de las líneas de corriente es
/
(5) (r 2 + 1) senO dr + r(r2 - 1) cosO dO = O.
De (5) se sigue que
r + 1 d + cos O de = o
r( r - 1) r sen e '
o
(r - 1) dr + cos ede = O,
2
2r -
r( r2 - 1) sen e
de tal forma que necesitamos solamente la ecuación
(6) :r dr _ dr + cos e de = O.
r--l r sen e
La solución de (6) es
In (r2 - 1) - In r + In sen e = In k J
que puede escribirse como
(7) (r 2 - 1) sen e = kr.
La ecuación ( 7) representa la deseada familia de líneas de corriente ,
algunas de estas líneas están mostradas en la figura 15.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
CAPÍTULO 6
Funciones hiperbólicas
(1 ) senhx=
e--r
2 '
y el coseno hiperbólico de x, simbolizado * cosh x, está defini-
do por
e· - e4
(2) coshx = .
2
El uso de sÍmbJlos y nombres tan similares a los de trigo-
nometría puede parecer imprudente. Alguna jUf>tificación
aparecerá en la siguiente sección. donde se muestra que las
así
(1 ) cosh 2 X - senh2 x = 1,
una identidad similar a la identidad bien conocida cos2 x + sen2 x = 1
en trigonometría circular. Véanse también los ejercicios 2 y 3.
Directamente de la definición encontramos que
y = senhu
es equivalente a
esto es,
d du
(2) - senhu = cosh u-o
dx dx
http://carlos2524.jimdo.com/
o
http://carlos2524.jimdo.com/
112 Funciones hiperbólicas [Cap. 6
x = a cos t; y = a sen t
son ecuaciones p ar ámetricas de la circunferencia X2 + y2 = a 2•
En 10s textos de trigonometría elemental se da generalmente una re-
presentación lineal de las seis funciones trigonométricas en una figura
basada en una circunferencia de radio unitario. Se puede obtener una
representación lineal similar';':' para las seis funciones hiperbólicas.
Las fórmulas usadas para definir el seno y el coseno hiperbólicos
también tienen análogos en el estudio de las funciones trigonométricas.
La definición dada en la página 141 para la función exponencial con
exponentes puramente imaginarios puede emplearse para deducir las
fórmulas pertinentes. Los resultados se muestran en el ejercicio 24.
EJERCICIOS
6. cosh 2 y = .} (cosh 2y + 1) .
+
7. cosh 2y = cosh 2 y senh 2 y
= 2 cosh 2 y - 1
= 2 senh 2 y +1.
8. senh 2y = 2 senh y cosh y.
9. senh (x + y) = senh x cosh y + cosh x senh y;
senh (x - y) = senh x cosh y - cosh x senh y.
10. cosh (x + y) = cosh x cosh y + senh x senh y;
cosh (x - y) = cosh x cosh y - senh x senh y.
11. tanh (x ) tanh x +tanh y .
+y 1 + tanh x tanh y'
tanh x - tanh y
tanh (x - y)
1 - tanh x tanh y'
d du
12. die cosh u = senh u dx'
d th h2 du
16 . dx co U = - csc U dx'
17. Obténgase la ecuación diferencial formada por eliminación de las
constantes el y e2 de y = el senh (kx) e2 cosh (kx) . +
d'
SOL. -...:1 - k2y = o.
d X2
18. Resuélvase la ecuación
(ysenhx - y + 2x) dx + (coshx - x) dy = O.
SOL. y cosh x - xy + X2 = c.
'19. Resuélvase la ecuación xy dx +
(x'Z + 1) (1 - Y tanh y) dy = O.
SOL. y2 (X2 + 1) = e cosh 2 y.
20. Resuélvase la ecuación y' senh x = 1 - 2y cosh x.
SOL. y senh 2 x = e + cosh x.
21. Encuéntrense las trayectorias ortogonales a las curvas sen x cosh y = c.
SOL. COS x senh y = k.
22. Defínase el seno hiperbólico i:1Verso de x, escribiendo y = senh-1 x,
por la ecuación.
I x = senhy
o
(A) x = t(e V - e-V
Póngase (A) en la forma
e2V - 2xev - 1 = O
y de este modo conclúyase que
y = senh- l x = In (x + V + 1) X2
CAPÍTULO 7
Ecuaciones diferenciales lineales
32. EL WRONSKIANO
Con las definiciones de la seCClon 31 en mente, obtendremos ahora
una condición suficiente para que n funciones sean linealmente indepen-
dientes en un intervalo a ::; x ::; b.
Supongamos que cada una de las funciones f,(X), f2(X),· .. ,fn(X)
es diferenciable al menos ( n - 1) veces en el intervalo a ::; x ::; b. En-
tonces de la ecuación
( 1)
por diferenciación sucesiva se sigue que
eJ; + e2f; + .. . + enf: = 0,
eJ:' + e2f~ + ... + e"f~ = 0',
EJERClCIOS
1. Obténgase el wronskiano de las funciones
1, x, x 2 , • • • , Xn-l para n > 1.
SOL. W = 0!1!2! ... (n - 1)!.
2. Demuéstrese que las funciones eX, e 2X, e 3X, son linealmente indepen-
dientes. SOL. W = 2e 6X =1=- O.
3. Demuéstrese que las funciones eX, cos X, sen X son linealmente inde-
pendientes. SOL. W = 2ex =1=- O.
4. Por determinación de las constantes Cl, C2' Ca, C4, no todas cero, y
tales que cdl + cz/2 + c3f3 + c4f4 = O idénticamente, demuéstrese
que las funciones
http://carlos2524.jimdo.com/
§33] Solución generol de una ecuación homogénea 119
(2)
Sea la relación
(4 ) y = f(x)
Considérese la solución
(7 )
de la cual
d"y d"-'y
(2) Ay = ao- +
dxn
a,--+
dX'O-l
.
Los coeficientes Ck, a,·, . . . , an en el operador A pueden ser funciones de
x, pero en este libro solamente usaremos los operadores con coeficientes
constantes.
Se dice que dos operadores A y B son iguales si, y sólo si, producen el
mismo resultado cuando cada uno de ellos actúa sobre la función y. Esto
es, A = B si, y sólo si, Ay = By para todas las funciones y que posean
las derivadas necesarias para las operaciones involucradas.
El producto AB de dos operadores A y B está definido como el opera-
dor que produce el mismo resultado que se obtendría usando el operador
B seguido del operador A. Así A By = A (By) . El producto de dos ope-
radores diferenciales siempre existe y es un operador diferencial.
Para operadores con co.eficientes constantes, pero no generalmente
para los de coeficientes variables, es verdad que AB = BA.
A (By) = (D + 2) (3 :~ - y)
2
= 3 d y _ dy + 6 dy _ 2
dX2 dx dx y
2
= 3 d y + 5 dy - 2
dX2 dx y
= (3D 2 + 5D - 2)y.
De aquí AB = (D + 2) (3D - 1) = 3D + 5D -
2
2.
Ahora considérese BA. Actuando sobre y, el operador BA produce
Por tanto
BA = 3D 2 + 5D - 2 = AB.
EJEMPLO b): Sea e= xD + 2, Y H = D - 1. Entonces
d2 d d
= x~ - x~+ 2~ - 2y
dX2 dx dx
2
dy dy
= x-+
dX2
(2 - x)- - 2y,
dx
Así
eH = XD2 + (2 - x)D - 2.
Por otro lado
esto es,
He = XD2 + (3 - x)D - 2.
Es digno de mencionarse que cuando tenemos dos operadores e y H
(uno de ellos con coeficientes variables), su producto depende del orden
de los factores. Sobre este tema, véanse también los ejercicios 17-22 de la
siguiente sección.
La suma de dos operadores diferenciales se obtiene expresando cada
uno de ellos en la forma
A = 3D2 - D +x - 2
y
entonces
A +B = (3 + X2) D + 3D + x + 5.
2
http://carlos2524.jimdo.com/
EJERCICIOS
así
(3) f (D) e1WlJ = eP"'f( m).
Si m es una raíz de la ecuación f (m) = O, entonces en vista de la
ecuación (3),
f (D)em" = O.
(2D 2 + 5D - 12) Y = O.
Ahora considérese el efedo del operador (D - m) sobre el producto
de em: y una potencia de X. Ya que
(D - m) (x"e m,,) = k~lem" + mx"em" - mx"e""',
obtenemos (D - m) (x"e m,,) = kx"-le m".
Entonces (D - m) 2 (x"em,,) = k(D - m) (~lem.,)
= k(k - 1 )~2em.,.
Repitiendo la operación, llegamos a
(D - m)"(x"e mz ) = k (k - 1) . .. 2 . lxOe fMJ
= k!e=.
http://carlos2524.jimdo.com/
128 Ecuaciones diferenciales l ineales [Cap 7
o
(3)
Otro de estos pasos nos lleva al resultado
(4) y( 3 ) = U( 3 )V ( 0) + 3U(2)V(1) + 3U(1)V(2) + U(0)V(3l,
11. Obténgase y(4) de y = x sen 2x. SOL. y(4) = 16x sen 2x - 32 cos 2x.
12. Obténgase y(4) de y = x cos 2x.
13. Obténgase y(4) de y = x 2 e"'.
14. Obténgase y(4) de y = x 2 e-X.
15. Obténgase y(4) de y = X2 sen x.
SOL. y(4) = X2 sen x - 8x cos x - 12 sen x.
16. Obténgase y(4) de y = X2 cos X.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
CAPÍTULO 8
Ecuaciones lineales
con coeficientes constantes
39. INTRODUCCIÓN
(1 )
dny
ao-
dx"
dn-ly
+ al - -+
dXn-l
.. .
dy
+ an-l-
dx
+ any = °
'
puede representarse en la forma
(2) f(D)y = 0,
133 .
http://carlos2524.jimdo.com/
134 Ecuaciones 'l ineales con cQéficientes constantes [Cap. 8
La ecuación auxiliar es
http://carlos2524.jimdo.com/
3= 2Cl - 2cz.
De las ecuaciones simultáneas para Cl y C2 concluimos que Cl = t
y c~ = - !. Por tanto
x = ¡(eH - e- zt ),
EJERCICIOS
(6)
Además, si f (D) contiene el factor (D - b) "', entonces la ecuación
(1 ) f(D)y = O
puede escribirse
(7) g(D) (D - b)"y = O
donde g(D) contiene todos los factores de f(D) excepto (D - b)n.
Entonces cualquier solución de
(4)
es también solución de (7) Y por tanto de (1).
Ahora estamos en la posibilidad de escribir la solución de la ecuación
( 1) siempre que la ecuación auxiliar tenga raíces reales solamente. Cada
raíz de la ecuación auxiliar es distinta de todas las otras raíces, o
pertenece al conjunto de las raíces ig.uales. Correspondiendo a una raíz
m. distinta de todas las otras, hay la solución
(8) yi = Ci exp (m ix),
y correspondiendo a las n raíces iguales ml, mil, . . . ,mt>, cada una
igual a b, hay la solución
(9)
La suma de tales soluciones (8) y (9) produce el número propio de
soluciones, un número igual al orden de la ecuación diferencial, debido
o que hay una solución por cada raíz de la ecuación auxiliar. Las solu-
ciones, obtenidas de este modo pueden probarse que son linealmente
independientes.
o
y = Cle" + (C2 + C3X + c.x2)e 2".
EJEMPLO b): Resolver la ecuación
4 3 2
d y + 2 d y + d y = O.
dx 4 dx 3 dX2
La ecuación auxiliar es
m4 + 2m + m3 2
= 0,
con raíces m = 0,0, - 1, - 1. De aquí, la solución deseada es
y = Cl + C2X + c3e-" + c1xe-".
EJERCICIOS
En los ejercicios 1-20 e:lcuéntrese la solución general.
1. (D2 - 6D + 9)y = O. SOL. y = (Cl + C2X) eS".
2. (D2 + 4D + 4)y = O.
3. (4D3 + 4D2 + D )y = O. SOL. y = Cl + (c 2 + cax ) exp ( - ix).
4. (D3 - 8D2 + 16D )y = O.
5. (D4 + 6D3 + 9D2) Y O. = SOL. y =
Cl + C2X + (Ca + C4X) c 3".
6. (DS - 3D2 + 4)y = O.
7. (4D3 - 3D + l )y = O. SOL. y = cle-" + (C2 + C3X) exp (ix ) .
8. (D4 - 3D3 - 6D2 + 28D - 24 )y = O.
9. (4D4 - 4D3 - 23D2 + 12D + 36)y = O.
SOL. y = (c l + c2x)e2" + (c s + c.x) exp ( - ~x).
10. (D3 + 6D2 + 12D + 8)y = O.
11. (DS - DS) y = O. SOL. y =
Cl + C2X + cax2 + C4e" + cse-";
O y = Cl + C2X + C3X2 + C6 cosh x + C7 senh x.
12. (DS - 16D3)y = O.
13. (4D4 + 4Da - 3D2 - 2D + l)y = O.
14. (54D4 - 27D3 - 9D2 + 7D - l)y = O.
15. (D4 - 2Da - 3D2 + 4D + 4)y = O.
SOL. y = (Cl + C2X) e-" + (cs + c.x) e2".
16. (4D4 + 4D3 - 3D2 - 2D + l)y = O.
17. (4DS + 4D4 - 9D3 - 11D2 + D + 3)y = O.
18. (DS - 15D3 + lOD2 + 60D - 72)y = O.
19. (D4 + 2D3 - 6D2 - 16D - 8)y = O.
SOL. y = (Cl +C2X) e-2"+ca exp [( 1 + V3) x] + c. exp [( 1 - V3) xJ.
20. (D4 - 2D3 - 5D2 + 2)y = O.
En los ejercicios 21-26 encuéntrese la solución particular indicada.
21: (D2 + 4D + 4) Y = O, cuando x = O, Y = 1 Y y' = - 1.
SOL. Y = (1 + x-)c 2
".
http://carlos2524.jimdo.com/
_.; 'x"
(2) ea: _ '" ---¡.
,,=0 n.
Si ahora tentativamente ponemos x = if3 en (2) como una definición de
eifl, obtenemos
ifl _ if3 i2 f32 i S f32 i 4p ... i"f3n
(3) e - 1+1!+21+31+41+ +nl+
http://carlos2524.jimdo.com/
§42) Definición del exp z para z imaginaria 141
(4)
o
(4)
e'(3 = ~
'X)
+ t. ~ea -':-::-~--'-:c-:--:--
"=0 ( 2k)! k=o (2k + 1)!'
Pero las series del miembro derecho de (5) son precisamente aquellas
para cos f3 y sen f3 tal como se desarrollan en cálculo. De aquí, llegamos
al resultado tentativo
(6) e i (3 = cos f3 + i sen /3.
El lector deberá darse cuenta de que las manipulaciones anteriores
no tienen significado por sí mismas en esta etapa (suponiendo que las
series infinitas con términos complejos no son una parte del contenido
de las matemáticas elementales). Lo que se debe entender es esto: las
manipulaciones formales anteriores han sugerido la definición (6) .
Combinando (6) con (1) podemos sacar una definición razonable de
exp (a + if3), específicamente
(7) exp (a + if3) = ea ( cos f3 +i sen f3 ), a y f3 reales.
R{;emplazando f3 por ( - f3) en (7) obtenemos un resultado que será de
valor en la siguiente sección,
exp (a - if3) = ea (cos f3 - i sen f3) .
Es interesante e importante observar que cpn la definición (7), la
función ett para z compleja retiene muchas de las propiedades que posee
http://carlos2524.jimdo.com/
142 Ecuaciones lineales con coeficientes constantes [Cap. 8
la función efIJ con x real. Tales asuntos son a menudo estudiad;k; en deta-
lle en libros sobre variables complejas.-K Aquí necesitamos conocer en
o
particular que si
y = exp (a + ib) x,
con a, b y x reales, entonces
(D - a - ib)y = O.
El resultado deseado se obtiene de inmediato por diferenciación, con
respecto a x, de la función
y = ea.z( cos bx + i sen bx).
EJERCICIOS
21. -
~
~
+ 2b -~
& + k x = 0, k > b > 0, cuando t = 0, x =
2
°&= y-
~
Ve}.
Vk
~--
MISCELÁNEA DE EJERCIC.IOS
CAPÍTULO 9
Ecuaciones
no hOIllogéneas:
Coeficientes indeterIllinados
. [ (m - 1)2 + 4)2 = 0,
o
m4 - 4m 3 + 14m 2
- 20m + 25 = O.
De aquí que la ecuación deseada es
(D 4 - 4D s + 14D 2 - 20D + 25)y = O.
Obsérvese que en todos estos problemas una solución correcta (pero
indeseable) puede obtenerse insertando raíces adicionales de la ecuación
auxiliar.
EJERCICIOS ORALES
l' '»- m = 0, 0, 1
de la ecuación auxiliar. La función complementaria es
http://carlos2524.jimdo.com/
(3)
Como la solución general de (1) es
y = yc + YP
donde yc está dada en (3) Y Yv es una solución particular cualesquiera
de (1), todo lo que resta por hacer es encontrar una solución p"articu-
lar de (1).
El miembro derecho de (1),
(4) R(x) = 3e'" + sen x
es una solución particular de una ecuación diferencial lineal homogénea
cuya ecuación auxiliar tiene las raíces
(5) m' = 1,Cti:i.
Por tanto la relación (4) es una solución particular de la ecuación
(6) (D - 1) (D 2 + 1) R = O.
Deseamos convertir (1) en una ecuación diferencial lineal homogénea
con coeficientes constantes, debido a que conocemos cómo resolver
tales ecuaciones. Pero, por ( 6), el operador (D - 1) (D 2 + 1) anulará
el miembro derecho de (1). Por tanto, aplicamos ese operador a ambos
miembros de la ecuación ( 1) y obtenemos
(7) (D - 1) (D 2 + 1) D 2
(D - 1) Y = O.
Cualquier solución de ( 1) debe ser una solución particular de (7) .
La solución general de ( 7) puede escribir inmediatamente de las raíces
de su ecuación auxiliar, esas raíces son m = O, O, 1 de (2) y los valores
m' = 1, -+- i de (5). De este modo la solución general de (7) es
(8) Y= Cl + C2X + c3e'" + c4xe'" + Cs cos x + C6 sen x.
Pero la solución general deseada de (1) es
(9) y = yc + yp,
donde yc = Cl + C2X + C3e'",
las Cl, C2, C3 son constantes arbitrarias como en (8). Así, debe existir una
solución particular de (1) que contenga al menos los términos sobrantes
de (8). Usando diferentes le~ras como coeficientes para enfatizar que no
Son arbitratias, concluimos que (1) tiene una solución particular
http://carlos2524.jimdo.com/
YP = A + Bx + C cos 2x + E sen 2x
en la cual Cl y C2 son constantes arbitrarias, mientras que A, B, C y E
serán determinadas numéricamente, tales que YP satisfaga la ecuación ( 7 ) .
Como
Dyp = B - 2C sen 2x + 2E cos 2x
y D2yp = - 4C cos 2x - 4E sen 2x,
la sustitución directa de yp en (7) produce
(8) - 4C cos 2x - 4E sen 2x +B- 2C sen 2x + 2E cos 2x - 2A
-2Bx - 2C cos 2x - 2Esen 2x = 2x - 40 cos 2x.
Pero (8) es una identidad en x, así que debemos igualar los coeficientes
de cada uno de los conjuntos de funciones linealmente independientes
cos 2x, sen 2x, x, 1 que aparecen en la identidad. De este modo se
slgue que
http://carlos2524.jimdo.com/
-6C + 2E =
40,
-6E - 2C = O,
- 2B = 2,
B - 2A = O.
de la cual A = -t y B = O.
La solución general de (9) es
y = Cl COS X + cz sen x - tx cos x.
Ahora
y~ = A( - xe-4 + e-4) + 2Be 2X,
y~ = A (xe-4 - 2e-4 ) + 4Be 2X ,
y;" = A( -xe-4 + 3e-4) + 8Be 2X
Entonces
y~' - y~ = 2Ae-4 + 6Be 2X
,
-1 =
° = C1
C2 -
+ C2 + C3 + t,
Cs + 3,
2= C2 + Cs - 2,
de las cuales C1 = -~,
2
C2 = 0, Cs = 4. Por tanto el resultado final es
Un detalle importante que algunas veces pasan por alto los estudian-
tes, es que la solución general, la y de ( 11), es la que debe satisfacer
las condiciones a la frontera.
EJERCICIOS
4. (D2 + 3D + 2)y = 1 + 3x + X2 .
.t5. (D 2 + 9)y = 5e" - 162x.
SOL . y = Cl COS 3x + C2 sen 3x + fe'" - 18x.
6. (D2 + 9)y = 5e" - 162x 2.
SOL. y = Cl COS 3x + C2 sen 3x + ie" - 18x 2 + 4.
7. y" - 3y' - 4y = 30e". SOL. y = cle"}) .+ C2 e-x - 5e"'.
8. y" - 3y' - 4y = 30e ". 4
SOL. y = (Cl + 6x) e4" + C2 e-x.
9. (D2 - 4)y = e " + 2.
2
SOL . Y = cl e- " + (C2 + !x)e
2 2X - i·
10. (D 2 - D - 2)y = 6x + 6e-x.
SOL . y = Cle-X + C2e2X - 3x + 3/2 - 2xe- x.
11. y" - 4y' + 3y = 20 cos X.
SOL . Y = cle x + C2eax + 2 cos x - 4 sen )O.
12. y" - 4y' + 3y = 2 cos x + .4 sen x. SOL. y = cle x + C2esx + cos x.
13. y" + 2y' + y = 7 + 75 sen 2X1.
SOL . Y = e-"(Cl + C2X) + 7 - 12 cos 2x - 9 sen 2x.
14. (D2 + 4D + 5yY = 50x + 13e aa1 .
SOL. y = e- 2X (cl COSX + C2Se:lx) + 10x - 8 + ie 3 ",.
15. (D2 + l)y = COSX. SOL. Y = ClCOSX + C2sen x + ixsenx.
16. (D2 - 4D + 4 )y = e2~. SOL. y = e2X (Cl + C2X + }x2 ).
17. (D2 - 1) y = e-X (2 sen x + 4 cos x) .
18. (D2 - l )y = 8xe x. SOL . y = cle-X + e"'(c2 - 2x + 2X2).
19. (Da - D )y =X. SOL. Y =Cl + C2e'" + C3r - !X2.
20. (DS - D2 + D - l)y = 4senx.
SOL. y = Cle'" + (C,2 + x) cos x + (C3 - x) sen x.
21. (D3 + D 2 - '4D - 4 )y = 3e-X - 4x - 6.
SOL . Y = cle '" + C2e-
2 2X + (C3 - x)r + x + l
22. (D4 - l )y = 7x • 2
y=
y y'
0, y'
3.
- 4
y y" = 10.
=
44. y" + 2y' + 5y 8e--"', cuando x 0, y =
y y'= ° = 8.
En los ejercicios 45-48, obténgase de la solución particular i:1dicada
el valor de y y el valor de y' en x = 2. .
45. y" + 2y' + y = x; en x = 0, y = - 3, Y en x = 1, y = -1.
SOlJ. En x = 2, y = e- y y' = 1. 2
por inspección.
Por ejemplo, si R (x) es una constante Ro y si b" =1= 0,
Ro
(2) YP = -
b"
es una solución de
(3) (boD" + b1D,,-1 + ... + bn) y = R{); bn =1= 0, Ro constante,
debido a que todas las derivadas de YP son cero, tal que
(boDn + b1Dn-l + ... + bn)yp = bnRo/b" = Ro.
°
Supongamos que b" = en la ecuación (3). Sea Dky la derivada de
más bajo orden que aparece en la ecuación diferencial. Entonces, la
ecuación diferencial puede escribirse como
(4 ) ( boD" + ... + bn-¡J)k ) y = R {) ; bn-k =1= 0, Ro constante.
Ahora, Dkx" = k! es constante, así todas las derivadas de mayor orden
de x" son cero. Esto nos conduce a que (4) tiene una solución
Rox"
(5) YP=-klb
. n-k '
entonces
(DS + 4D 3
)
7x 3
2{ = °+ 7.6
4 . 24 = 7.
es una solución.
Finalmente, observemos que si yl es una solución de
http://carlos2524.jimdo.com/
f62 Ecuaciones no homogéneas: coeficientes indeterminados [Cap. 9
y y2 es una solución de
entonces
es una solución de I
f (D)y = Rl ( X) + R2(X).
Se sigue fácilmente que la tarea de obtener una solución particular de
• f(D)y = R (x)
puede dividirse en partes tratando de separar los términos de R (x)
independientemente, si ello es conveniente. Veremos los ejemplos siguien-
tes. Esta es la base del "método de superposición" que desempeña un
papel importante en las matemáticas aplicadas.
y que
ya = 2\ sen 4x
satisface
(D2 - 9)ya = - sen 4x.
De aquí que
yp = - ~e'" - ;x + 2\ sen4x
es una solución de la ecuación ( 13 ).
EJERCICIOS
1. Demuéstrese que si b =1= a, entonces
(D2 + a2)y = sen bx
tiene la solució:J. particular y = (a 2 - b 2 )-1 sen bx.
2. Demuéstrese que la ecuación
(D2 + a 2 ) y = sen ax
. no tiene solución de la forma y = A sen ax, con A constante. En-
cuéntrese una solución particular de la ecuación.
X
SOL. Y = - 2a cos ax.
CAPÍTULO 10
Operadores diferenciales inversos
Por tanto
o
yP = (2x 3 - 3X2 + ~x)e3!JJ
y la solución general de la ecuación original (3) es
y = Cle'" cos 2x + cze'" sen 2x + (2x 3
- 3X2 + :x) e3'"
2
a
(11 )
U na solución particular de (11) se obtiene fácilmente:
e-"Yl = ¡¡Ox s,
así que
(12)
La ecuación (10) se trata como en el capítulo 9. Hacemos
y2 = Ax + B.
Entonces Dy2 = A, Y de (10) es fácil encontrar que A = 7, B = 12.
De este modo una solución particular de (10) es
(13) y2 = 7x + 12.
Usando (12) Y (13), y las raíces de la ecuación auxiliar de (8 ), la
solución general de (8 ) puede ahora escribirse como:
y= (Cl + C2X) e" + ~xse" + 7x + 1~.
EJERCICIOS
- ,.
49. EL OPERADOR 111(D)
Buscando una solución particular de
( 1) f(D)y = R(x),
es natural escribir
1
(2) Y = f(D) R ( x)
Y
(2)
(3 ) 1 az _ e I() =1= O
I(D)e - / (a)' a •
Ahora de (1) se sigue que
uz UZ
f(D) e = f(a)e = e"'''.
f (a) f(a)
http://carlos2524.jimdo.com/
senax= - - - -
2i
entonces
=--
x eai'" + e-4i'"
2a 2
http://carlos2524.jimdo.com/
§51] Eva luación de ID2 + a 2)- l sen ax y ID2 + a 2)- l cas ax 173
tal que
1 x
(2) =-------,.sen ax = - -cos ax.
D2 + a 2
2a
La comprobación de (2) se deja como ejercicio.
Otro resultado útil,
1 x
(3 ) =-- -2 cos ax = - sen ax,
D2 + a 2a
se puede obtener del mismo modo.
EJERCICIOS ORALES
EJERCICIOS
f<n) (a) *
donde n es el entero no negativo más pequeño para el cual
O. Véase C. A. Hutchinson, "A.-lOther Note on Linear
Operators", Amer. Math. Monthly, 46, 1939, página 161.
1 "
f(D) x
se evalúa correctamente si reemplazamos l/f (D) por su desarrollo for-
mal en series de potencias ascendentes de D, llevando al cabo después
las operaciones indicadas. Por una evaluación entendemos que si yp es el
resultado de la evaluación, entonces f(D)yp = x".
La extensión de este método para la evaluación de [1/f (D) ] P(x) ,
. donde p ( x) es un polinomio en x, deberá ser evidente.
Dos comentarios se imponen aquí. Primero, si f (D) tiene un factor
Dk, k > O, entonces las potencias negativas de D aparecerán en el des-
arrollo. Estas son intérpretadas como integraciones iteradas según ilus-
tramos abajo.
Ahora
entonces
YP = D2 +
1 X2 =
1
(1 - D2) X2 = X2 - 2
.
en la cual todos los términos que involucran a D"'X2 para m > 2 han
sido eliminados debido a que cada uno es cero.
Entonces
1 5
2 4 x5
YI' = D2 _ 1 . :0 = (- 1- D - D - D
6
- • ) 60'
de la que
X5 XS
YI' = - 60 - 3" - 2x.
= ( - D-s - JFl - D _ D S ~ . . . ) x 2.
como antes.
+
YI' = D2 b_
2 (2x - x
3
).
El carácter formal del proceso que hemos usado hace imperativa la veri-
ficación de nuestro resultado final por sustitución directa en la ecuación
diferencial original.
Se deja la verificación al estudiante.
La brevedad de la solución anterior se debe a la omisión de casi todo
el trabajo, por ejemplo, la división larga y la reducción de (3) a (4).
Antes de ser demasiado entusiasta de este método, el estudiante deberá
resolver la ecuación (2) por el método del capífulo 9, así como efectuar
. los pasos omitidos en el desarrollo anterior. Cuando se haga esto se pueden
comparar los dos procesos.
Parece haber un tipo de ecuación para el cual este método es mucho
más ventajoso que los vistos previamente. Véase el siguiente ejemplo.
* Para aquellos que la vida haya sido especialmente dura, es admisible que
15!/7! = 259 459 200.
http://carlos2524.jimdo.com/
178 Operadores diferenciales inversos [Cap. 10
(2 )
(3)
en e! que Al, A2, ... ,A" son constantes. Se puede demostrar * que en
[l/f(D)]R (x), e! operador 11t(D) puede ser reemplazado por e! IIÚem-
bro derecho en (3). Entonces, la fórmula (2) o cualquier otra herra-
mienta matemática adecuada puede usarse para calcular cada térrrúno
de! tipo Ak( D - mktlR (x).
Aun si f(D) involucra factores repetidos, se puede efectuar un cálcu-
lo similar al realizado en el caso anterior. Entonces, el desarrollo en
fracciones parciales requiere la evaluación de términos de la forma
(D - m)-nR (x) con n> 1. Haciendo
yp = (D - mtnR (x) ,
entonces
(D - m)"yp = R (x) ,
en la cual el cambio de la exponencial trabaja perfectamente. Así, de
e-tll4: (D - m)nyp = e-"''''R (x )
se sigue que
D" ( e-"''''yp) = e-m",R ( X ) ,
de tal forma que el valor de yp puede obtenerse con n integraciones.
En un paso de la deducción de la fórmula para el cambio de la ex-
ponencial, demostramos (sección 48) que para cualquier y
f(D)y = eG.2:f(D + a)(e-a",y).
Escojamos
_ az 1 V( )
y - e f(D + a) x.
Entonces
f(D)y = eG.2:f (D + a) [,(D ~ a) V(x) ]
= eazV(x).
* E. L . Ince, Ordinary Differential Equations. Londres: Longmans-Green, 1927,
páginas 140-141.
http://carlos2524.jimdo.com/
Por tanto
1
- - eU"V(x) = y
f(D) ,
o
1 1
(4) _ _ [eU"V(x)] = ea.: V(x).
f(D) f(D + a)
La relación (4) puede emplearse para hacer el cambio de la exponen-
cial automáticamente al resolver una ecuación tal como
(5) (D 2 - 2D + 5)y = e'" cos 2x.
Aquí
- X 1 2
y" - e (D + 1) 2 _ 2 (D '+ 1) + 5 cos x
o
1
_
Y'P - e D 2
x
+ 4 cos 2x.
Entonces, empleando (3), de la sección 51, obtenemos
yp = ixex sen 2x.
EJERCICIOS
lO!
15. (D6 - l)y = x lO • SOL. Y = _ x lO - 4! x 4
16. (D3 + 3D2 - 4)y = 16x3 + 20 2
•
SOL. y = -4x3 _ 5X2 _ 18x _ 2/ .
MISCELÁNEA DE EJERCICIOS
°
34. La ecuación del ejercicio 33 con las condiciones de que cuando
x = 0, y = y y' = O. SOL. 4y r c- 2X2 - 2x + 1 - e-2".
En los ejercicios 35-37, obténgase de la solución particular indicada
el valor de y y el valor de y' en x = 1. .
35. y" + 2y' + y = x + 2, cua:J.do x = 0, y = 1 Y y' = O.
SOL. Enx = l,y:;::: (e - l) j eyy' = (e - l)je.
36. y" + 2y' + Y = x + 2, cuando x = 0, y = O Y y' = O.
En x = 1, y= '(e-l) je y y' = 1.
37.
SOL.
(D2 + l)y = 3, cuando x = 7rj2, Y =
SOL. En x = 1, Y
°=
y y' = O.
0.4756 Y y' = -1.6209.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
CAPÍTULO 11
La transformada de Laplace
(1) T{F(t)} =J 00
-00
K(s, t)F( (t) dt = I(s).
183
http://carlos2524.jimdo.com/
184 la transformada de laplace [Cap. 11
( 1) L{F (t)} = r o
e-stF ( t) dt.
L{eTct } =
,
1""
o
e-(B-k)t 'dt
+ e.
Como
f az
e senmx x -
d - ·eai: ( a sen mx - m cos mx)
a +m
2 2
Para 5 positiva, e-8t ~ O cuando t ~ oo. Más aún, sen kt y cos kt están
limitadas cuando t ~ oo. Entonces (3) da
_ 1(0 - k)
L{senkt} - 0- 2 k2 ,
5 +
o
le.
(4) L{sen kt} = 5> O.
5
2
+k 2 ,
El resultado
5
(5 ) L{cos kt} = 2 k2 , 5> O,
5 +
puede obtenerse de manera similar de la fórmula elemental
d - e=(a cos mx + m sen. mx)
(6)
f az
e cos mx x - 2 + 2
a m + c.
EJEMPLO e): Obtener L{f>} para n entero positivo.
Por definición
L{t"} = ro
e-stt" dt.
nt"-l dt - -1 e-st
s
De este modo obtenemos
S> O,
o
(8) L{t"} = 2: L{t"-l}, S> O.
S
L{tn-1} = n - 1 L{t"-2}.
S
así
(9 )
( 10) L{t n} = ~
S,,+l ,
s > O. \
L{H(t)} = ( [stH(t) dt
= ro
e-stt dt + J'" e- st 5 dt.
4
http://carlos2524.jimdo.com/
188 La transformada de Laplace [Cap. 11
De este modo
4e-48 e- 48 1 5e- ol8
L{H (t)} = - -5 - -52+ 0+ -
52
- 0 +5-
EJERCICIOS
= O, t> 7r.
2(1 - e-1TS )
SOL.
52 +
--'---::----;-
4--'--' s > O.
( 1) 1rOOo e- 8 t F ( t) dt
existe para algún rango de valores de 5. Por tanto, buscamos una clase
amplia de funciones para las cuales podamos probar que para todas
ellas la integral (1) existe.
Uno de nuestros intereses reconocidos en la transformada de Laplace
es su utilidad como herramienta en la solución de problemas que tienen
que ver con aplicaciones más o menos elementales, particularmente en
problemas de valores a la frontera en ecuaciones diferenciales. Por
tanto, no vacilamos en restringir nuestro estudio a funciones F (t) que
sean continuas o aun diferenciables, excepto posiblemente en un conjun-
to discreto de puntos, en el rango semi-infinito t ~ O.
http://carlos2524.jimdo.com/
190 La transformada de Laplace [Cap. 11
, /
I
I que en cada uno de los subintervalos:
I
I
I 1. F ( t) es continua en el intervalo
N
I
I
I
I
I
I
I
I
I
I
I
I
I
abierto c < t < d,
2. F ( t ) tiende a un límite cuando t
I I I I
I I I , se aproxima a cada extremo del interva-
I
O 2 3 5 6
lo desde adentro del intervalo; esto es,
FIGURA 17 Lim F ( t) Y Lim F ( t) existen.
t-+c+ t~d-
Sabemos que
(2) J'' ' e -ct dt
o
converge para c > O. Esto de·spierta nuestro interés por funciones P (t)
que están, para t ( t ~ lo), esencialmente lirrútadas por alguna expo-
nencial ebt , tales que el integrando en ( 1) se comporte como el integrando
en ( 2) para s suficientemente grandes.
(5) r o
e-stP (t) dt = r o
e-stP(t) dt + r ~
e-stP(t) dt.
(7)
existe para s > b.
L{P (t)} = r e-stP(t) dt = f(s)
. t
t -+", e
3
Lun;-¡¡¡= bbt
t-+", e
2
L'1m 3t = L'llT!.b 6t = L'1m 6 =
t-+",
2bt
e t-+",
b3bt
e
°
.
L{r i } = 2 r o
exp ( -sf32 )d{3, s> O.
r
o
exp ( - i) dy = 1- y-:;'
Por tanto.
( 1) L{rll } = 2s~' i y-;'
= GY, s> O,
(1 ) L{F'(t)} = ro
e-8t F'(t) dt
F' (t) dt
_se- 8 t dt F (t)
( 1) f(s) = r
o
e-8t F (t ) dt
se sigue que
( 2) f'(s) = r o
( -t )e-8t F(t) dt.
( 6)
De la fóm1Ula conocida
S
s
2
+ k~ = L{cos kt}
(7)
http://carlos2524.jimdo.com/
200 la transformada de laplace [Cap. 11
S2 ~P = L{¡senkt}
para obtener
S2 +(S2p ++ k2)2
k2 - S2 {1
= L ksenkt - tcoskt ,
}
de.la cual se sigue que
(1)
r
cuando {3 ~ oo. Así
(4) 1'(x + 1) = x e- f3{3"-ld{3 = x1'(x).
= n(n - 1) (n - 2) .. ·2·1·1'(1).
= n!1'(I).
http://carlos2524.jimdo.com/
L{r-i} = r(¡) .
s
ro)!.De aquí
Pero ya sabemos que L{t1} = (--;
c (
- s(3) d(3
, = 1 - exp ( -se)
Jo exp S
http://carlos2524.jimdo.com/
203
¡P (t, e)
o C 2C 3C 4C
FIGURA 19
en consecuencia, de ( 11 )
(15)
L{
Q(
te) }
J
= !s [ 1 + exp2 ( - se)
-IJ =1-. s
l-exp ( - se).
1 + exp ( - se)
Q (t,c)
o C 2C 3C 4C
-1
FIGURA 20
http://carlos2524.jimdo.com/
204 la transformada de laplace [Cap. 11
EJERClCIOS
1 (ñ)~
1. Demostrar que L{ t! } = 25 ~ , 5 > O.
15 ~ ,
2. Demostrar que L {ti' } = 8f! (ñ)! 5 > o.
3. Use la ecuación (4), página 197, para deducir L{senkt}.
4. Use la ecuación (4), página 197, para deducir L{coskt}.
5. Compruebe las tra:lsformadas conocidas de sen kt y cos kt una con-
tra otra, empleando el teorema 9, página 197.
6. Si n es un entero positivo, obtener L{tnek t } de la transformada cono-
cida L{ek t } empleando el teorema 13, página 199.
n!
SOL. (5 _ k)n+l' 5> k.
2k (3s - k 2)
2
7. Encontrar L{t 2 sen kt} . SOL. (52 + P )3 , 5> o.
2S (s2 - 3k 2)
8. Encontrar L{t 2 cos kt}. SOL. (S2 + k2) 3 ' 5> O.
9. Para la función
F(t) = t + 1, o ~ t ~ 2,
= 3, t> 2,
grafique F (t) Y F' (t ) . Encuentre L{ F (t) }, Encuentre L{ F' (t)} por
dos métodos.
SOL. L{F' (t)} = ,1(1 - e- 28 .) , 5> O.
10. Para la función
H (t ) = t + 1, o ~ t ~ 2,
= 6, t> 2,
encuentre lo mismo que en el ejercicio 9.
11. Definida una función de onda triangular T (t, e) por
T(t, e) = t, o ~ t ~ e,
= 2e - t, e < t < 2e;
T (t + 2e, e) = T (t, e).
Describa T (t, e) y encuentre su transformada de Laplace.
1 es
SOL. -; tanh - .
s- 2
http://carlos2524.jimdo.com/
§63) Funciones periódicos 205
= O, -r. <t<-
2r.;
w w
F(t + ~) = F (t).
w 1
SOL . S2+ w2. (S7r) .
l-exp - -
w
14. Encontrar L{F(t)} donde F(t) = t para 0< t < w y F(t + w) =
F( t).
SOL. .!.. - ~ exp ( -sw) = 1. + ~ (1 _ coth~).
52 S 1 - exp ( - sw ) S2 25 2
15. Probar que si L{F(t ) } = f(s) y si F(t) It es de clase A.
CAPÍTULO 12
Transformadas inversas
207
http://carlos2524.jimdo.com/
208 Transformadas ¡n-, ¡ersas [Cap. 12
Entonces
(6) L-
1
t2~ F} = cos kt.
= r e-st[e"tp(t)] dt.
t2: =
Como sabemos que L-l k2} senkt, procedemos como SIgue:
Escribimos
L -1 {
52 + 65+ +1 25} -_ L (5 + 3)+21+ 16} .
5 -1 { 5
Entonces
L-
1
t2 : ~ ~ 2s} = e-stL- t: ; ~6} 1
= e- 3t
(cos 4t - i sen 4t).
EJERCICIOS
1
e 3t sen t.
2. S2 _ 6s + 10' SOL.
s
e-t(cos2t-isen2t).
3. S2 + 2s + s· SOL.
25 - 1
e~2t (2 cos st - sen st) .
8. + 4s + 29'
$2
SOL.
3s + 1
-1 J 9. (s + 1) . 4
SOL. e-t G i ts) .
t2 -
L-l{
(s + a1)2 + b2}=..!.b e-atsenbt.
14. Demuéstrese que
L-l{f(as)) = ~ F G).
16. Para a> 0, demuéstrese que de L -l{f(s)} = F(t) se sigue que
y y
/'
I
I
I
I
I
I
o 01 e
FIGURA 22 FIGURA 23
L{y ( t)} = r ca
t 2e-8t dt + r
2
6e- st dt.
S2
S8
}.
0:;;t<1
en el cual a(t - 1) = O Y a(t - 3) = O. Encontramos
( 10) g (t) =3, O:;;t<1.
EJERCICIOS
t >
71",
71". SOL.
S2
3(1 + tr''')
S2 +9 .
9. Obténgase y descríbase una transformada inversa de Laplace de
s s
F ( t ) = L -1 { (s :381) a}
evaluar F (2), F (5) , F (7 ).
SOL. F (2 ) = O, F (5) = 2e- 2 , F (7) = 8e-4.
http://carlos2524.jimdo.com/
216 Transformadas inversas [Cap. 12
(3 ) f (s) = r o
e-stF (t) dt.
(4) g(s) = ro
e- 8 !3G (f3) d,f3,
se sigue que
(7) L-l{e- S !3f(s ) } = F (t - (3 )a(t - (3 ) ,
en la que a es la función escalón discutida en la sección 65. La ecua-
ción ( 7) significa que
(8) e-SPI(s) = [ e-stF (t - (3 )a (t - (3) dt.
L-
1
{~} = 1,
S3 + 4s + 3s s 5 + 1 s + 3'
2
... Véase E. D. Rainville, U nified Calculus and Analytic Geometry. Nueva York :
The Macmillan, 1961, páginas 357-364.
http://carlos2524.jimdo.com/
220 Transformadas inversas [Ca p. 12
3
( 3) 5s - 6s - 3 = Al + A2+ A3+ ~ + Bz
S3(S + 1) 2 S S2 S3 S+ 1 (s + 1)"
Correspondiendo a un factor denominador (x - yr, debemos suponer,
en general, r fracciones parciales de la forma
~+ A2 • + ...+ Ar l"
16 + :i. + BIS + C I = B 2s + C2
S(S2 + 4) 2 S
. , S2 + 4 (S2 + 4) 2'
http://carlos2524.jimdo.com/
§67] Fracciones parciales 221
De la identidad
16 = A (S2 + 4 )2 + (BIS + e l) S(S2 + 4 ) + (B2S + e 2)S,
no es difícil encontrar los valores A = 1, B1 = -1, B2 = - 4, el = o,
e 2 = O.
Por tanto obtenemos
L -1 { 16
S(S2 + 4 ) 2
} = L- 1 {lS __+
s~
S_ _
4
4s
(S2 + 4 ) 2
}
= 1 - ces 2t - t sen 2t
Es posible obtener fórmulas para los desarrollos en fracciones parcia-
les de las fracciones racionales que serán tratadas en esta sección. Tales
fórmulas son muy útiles en la teoría y no particularmente ineficientes
en la práctica. Las téc~cas elementales anteriores, si se aplican inteli-
gentemente, son eficientes en problemas numéricos y son los únicos
métodos de fracciones parciales que serán prEsentados aquí.
EJERCICIOS
En los ejercicios 1-10, encué:1trese una transformada inversa de la f(s)
dada..
1
1. -.- -.
s- +
as
SOL.
s+2
2. S 2 _ 6s 8· + SOL.
1) + 5s - 4.
3 2s2 2+ e- 2t •
+ S2 - 2s
· s~
SOL. el _
:1 4 2s + 1 2
1 + e-t - 3te-t •
6'j · s (s +1 )2 · SOL .
4s + 4
5. S 2 (S _ 2) . SOL. 3e2t - 3 - 2t.
1
6. S3 ( S2 + 1) . SOL. .~p - 1 + cos t.
5s - 2
7. S2 (s + 2) (s - 1) .
SOL. t - 2 + el + e-2t •
del cual
(2)
http://carlos2524.jimdo.com/
diferentes.
Podemos escribir
4 (s + 1) + 6 3
y (s)= (S+1 )2 + (S+1 )4
4 6 3
= s+ 1 + (s+ 1)2+ (s+ 1)4
Empleando la transformada inversa, obtenemos
(8)
En este problema hace su aparición otra faceta del poder del método
de la transformada de Laplace. El hecho de que la función .p(t) en la
ecuación diferencial tenga derivadas discontinuas, hace difícil el uso
del método clásico de coeficientes indeterminados, en tanto que estas
discontinuidades no interfieren en lo absoluto con la simplicidad del
método de la transformada de Laplace.
Para tratar este problema, hagamos L{x(t)} = h(s). Necesitamos
obtener L { .p (t ) }. En términos de la función o: podemos escribir, de (1 0)
(11) .p(t)=4t - 4(t -1 )0:(t -1 ), t~O.
s2h(s) - s - 0+ 4h (s)
del cual
(12)
Ahora
4 1 1
(18)
+ -f
Bit
x(t) = A coskt + -senkt
k k o
F (t - 13 ) senkf3 df3,
o
(19) x(t) = A cos kt
Bit
+ k sen kt + k L F(f3) sen k, (t - 13) df3.
(22) w'(O) = B
y esperamos detenmnar B posteriormente usando la condición de que
w ( l ) = -l.
El problema transformado es
1
s2g(S) - s( - 3) - B + 2[sg (s) - ( -3)J + g(s)
del cual
_ -3 (s + 1) +B- 3 1
(23) g (s)- (s +1 )2 +s2(s +1 )"'
pero, usando el desarrollo en fracciones parciales tendremos
1 2 1 2 1
--,,-,-------,--.,. = - - + - + - - + - -----,.
S2(s + 1) 2 S S2 S + 1 (s + 1) 2 ,
tal que
(24)
de donde se obtiene
(2 5) w(x) = x - 2 - e-:C + (B - 2)xe-".
Tenemos aún que imponer la condición de que w ( 1) = -1. De (25),
con x = 1, obtenemos
- 1 = 1 - 2 - e- 1 + (B - 2) e-t,
de modo que B = 3.
Por tanto nuestro resultado final es .
(26) w (x ) = x - 2 - e-31 + xe-".
El problema en el ejemplo e) puede resolverse eficientemente por los
métodos' del capítulo 9. Véanse también 10s ejercicios del 21 al 42
siguientes.
EJERCICIOS
(2) F (t) = 4t - 3 r o
F({3) sen (t - (3) d{3.
(3)
o
MD = t + ~sen2t = F (t) = MI,
como se deseaba.
Es importante darse cuenta que la ecuación original
Ya que LUX2} = ly
Z2
L{x} = -; , podemos aplicar el operador Len
z
toda la ecuació~n (5) Y obtener
hez) = 1_ h(z)
Z3 Z2
de la cual
o
h z - 1 Z2 + 1 - Z2 _ 1 Z
( ) - Z (Z2 + 1) Z(Z2 + 1) - ~ - Z2 +1
http://carlos2524.jimdo.com/
§69] Ecuaciones integrales especiales 231
Entonces
g (x) = L-1{! __+z_}1
Z Z2
o
(7) g (x) = 1 - COSX.
= h 2
- [(X - y)(y - senY)J: - 1: (y - seny) dy
= ir - 0- [iy2+ cosyJ :
= ti2 - iX2 -:- cos x + 1 = 1 - cos x = MI.
EJERCICIOS
2. F (t) = + L
-t
1 F(,8) sen (t - ,8) d,8. SOL. F(t) = 1 + ·w.
3. F (t ) = t + J: F (t - ,8)e-¡3 d,8.
4. F (t ) = 4t 2
J: F(t - ,8)e- d,8.
- P
SOD. F(t) = -1 + 2t + 2t + e-
2 2t
•
(A )
L
t
F (t - 13 ) e-fJ df3 = t
j
http://carlos2524.jimdo.com/
235
TABLA DE TRANSFORMADAS
f (s ) = L{F (t ) } F (t)
f (s - a) eatF (t)
f(as + b) ~ exp ( - ~) F G)
_1 e- C8 , a(t - e) = O, °st<c
c>O
s = 1, t ~ e
1
- 1
s
1 tn
-
S1H-l n!
1 t~
x> -1
SX+l' I' (x + 1)
S1 (7rt) 1
-1- e-at
s+a
/
http://carlos2524.jimdo.com/
236
-
I (s) = L{F (t) } F (t)
1 {"e-'L t
-
(s + a)n+l n!
k
sen kt
S2 + k2
s
cos kt
S2 + k2
k
S2 _ k2 senh kt
s
S2 _ k2 cosh kt
2k 3
sen kt - Ict cos kt
(S2 + k 2) 2
2ks
t sen kt
(S2 + k2) 2
1 - e- t
In(l + +) --
t
In s + k 2 senh kt
s-k t
.-
k sen kt
Arctan - --~-
s t
http://carlos2524.jimdo.com/
CAPÍTULO 13
Aplicaciones
( 1) ( ) = xos + Vo + f (s)
u. S S2 + W l + 13 2'
o
( 5) u (s) = sXo + Vo + Aw
S2 + W (S2 + W) (S2 + ( 2)
El último término de .( 5) nos conduce a transformadas inversas dife-
rentes según que w = f3 o w =1= f3. El caso w = f3 nos lleva a resonancia.
Esto lo discutiremos en la sección siguiente.
Si w =1= f3 la ecuación (5) nos lleva a
(6 )
sXo + Vo
u (s) = S2 + W
Aw
+ w2 + f32
(1+ W - S2 +1)2 .
S2 w
esto es, el peso oscila entre los puntos 5 pulgadas abajo y arriba de E.
El periodo es trr seg.
EJEMPLO e): Una cuerda, con una constante de 0.75 libras por pie
está sobre una tabla grande y lisa (sin fricción). Un peso de 6 libras se
ata a la cuerda y está en reposo (velocidad cero) en la posición de equi-
- librio. Se aplica una fuerza de 1.5 libras al soporte a lo largo de la línea
de acción de la cuerda durante 4 segundos. Analizar el movimiento.
Debemos resolver el problema
( 11 ) 3
6
2
x" ( t) + ix ( t) = H ( t) ; x(O) = 0, x'(O) = 0,
8
S2U ( s) + 4u ( s) = - (1 - ¡,r-48),
S
o
_ 8( 1 - e-48)
u(s) - S(S2 + 4)
= 2 (1 __r+4
S
_s_)? (1 _ e-48).
La solución deseada es
(13) x(t) = 2(1 - cos2t) - 2[1 - cos2(t + 4)]a(t - 4)
Por ·supuesto, la solución (13) puede descomponerse en las dos rela-
Clones
(14 ) Para O ::; t ::; 4, x (t) = 2 (1 - cos 2t) ,
(15) Para t >4 x(t) = 2[cos 2(t - 4) - cos 2t],
si estas formas· parecen más simples de usarse.
La verificación de la· solución (13) o (14) Y (15) es directa. El
lector puede demostrar que
Lim x(t) = Lim x(t) = 2( 1 - cos 8) = 2.29
72. . RESONANCIA
En el problema de las vibraciones no amortiguadas de una cuerda,
uno de los ejemplos que encontramos fue
http://carlos2524.jimdo.com/
§72] Resonancia 245
Que (5) satisface todas las condiciones de (3) puede verificarse rápi-
damente.
En la solución (5 ) los términos proporcionales a cos [3t y sen Pt tienen
un límite, pero el término con {Jt cos f3t puede hacerse tan grande como se
quiera con una elección apropiada de t. Esta construcción de grandes
amplitudes en la vibración cuando w = P en (1) se llama resonancia.
EJERCICIOS
22. Una cuerda es alargada 1.5 pulgadas por un peso de 4 libras. Supón-
gase que el peso es llevado 3 pulgadas por debajo del ' punto de
equilibrio y entonces soltado repentinamente. Describir el movimiento
si hay, actuando sobre la cuerda, una fuerza de 8 sen 16t.
SOL. x = 1(1 - 8t) cos 16t + i sen 16t.
23. Para el movimiento del ejercicio 22, encontrar los cuatro primeros
tiempos en los que la velocidad es cero y determinar las posiciones
correspondientes. SOL. t = i , 7Tj16, 7Tj8, 37Tj16 (seg) y
x = +0.11, +0.14, -0.54, +0.93 (pies), respectivamente.
( 2) ( ) - xo (s + y) + vo + YXo + {(s)
u s - (s + y) 2 + f32 _ y2 (s + y ) 2 + (32 - y2 .
Trataremos ahora el caso cuando y > O; esto es, una fuerza retardataria
no desdeñable está presente. Sabemos que la forma de la transformada
inversa de u (s) depende de que (3 > y, (3 = y, 0(3< y.
Si f3 > y, (32 - y2 > O, así que hacemos
(3 )
y llegamos al resultado
(5) x( t) = e-yt[xocos8t + (vo + yxo )8-1sen8t~ + lh (t)
en el que 0/1 (t) es una transformada inversa del último término de (4).
La función 0/1 ( t) puede siempre escribirse en términos de F (t) por el
teorema de convolución, por ello conduce a complicaciones indeseables
para los,casos más comunes de F ( t). En la solución (5), la presencia del
http://carlos2524.jimdo.com/
§73] Vibraciones amortiguada s 24 9
}(
FI GU RA 29
http://carlos2524.jimdo.com/
250 Ap licaciones [Cap. 13
FIGURA 30
xt
()_
-
iL-1 { S - 4.4 }
(s + 0.8)2+ 63.36 .
= i exp (-0.8t)L-
1
+-6~·.~6}
{/
EJERCICIOS
pies por segundo, del peso móvil. Determínese la posición del peso
como una función del tiempo. SOL. x = _ie-4t sen 8t.
6. Una cuerda es alargada 6 pulgadas por un peso de tres libras. Este
peso está atado a la cuerda y principia su movimiento desde el
equilibrio CO::1 una velocidad hacia arriba de 12 pies/seg. La resis-
tencia del aire actúa como una fuerza retardataria igual en magnitud
a 0.03Ivl. Encontrar la ecuación de movimiento.
SOL. x = _1.5e-o. 16t sen 8t.
7. Una cuerda es alargada 6 pulgadas por un peso de dos libras. Hay
una fuerza amortiguadora presente de una magnitud igual a la mag-
nitud de la velocidad. Actuando sobre la cuerda hay una fuerza
motriz (2 sen 8t). Si, en t = O, el peso es soltado desde un punto
situado tres pulgadas por debajo del punto de equilibrio, encontrar
su posición para t > O. X = (i + 4t) e-
SOL. st - i cos 8t.
8. Una cuerda es alargada 10 pulgadas por un peso de 4 libras. El peso
principia su movimiento desde 6 pulgadas abajo del punto de
equilibrio con una velocidad hacia arriba de 8 pies/seg. Si el
movimiento se desarrolla en un medio que actúa como una fuerza
retardataria de magnitud llvl, describir el movimiento.
SOL. x = e- t [0.50 cos 6.1t - 1.23 sen 6.1t].
9. Para el ejemplo 8, encontrar los tiempos a los que la velocidad es
cero (los primeros tres), y la posición del peso en cada una de estas
ocasiones.
SOL. t 1 = 0.3 sec., Xl = -12 pulgadas; t:2 = 0.8 sec.,
X2 = +6 pulgadas; ta = 1.3 seg, X3 = -4 pulgadas.
~ dx
2
( A) b dx = w.
g dt 2 + dt
Resolver la ecuación (A) dadas las condiciones
(B) cuando t = O, x = O Y v = Vo.
Emplee a = bg/w.
http://carlos2524.jimdo.com/
254 Aplicaciones [Cap. 13
(1)
http://carlos2524.jimdo.com/
§75] Vigas 255
75. VIGAS
Considérese una viga de longitud 2c, como la mostrada en la figu-
ra 32. Denótese la distancia, a partir de uno de los extremos de la viga,
por x, y su deflexión por y. Si la viga está sujeta a una carga vertical
W ( x), la deflexión y deberá satisfacer la ecuación
d4
(1 ) El --.14 = W(x) para O < x < 2c
dx '
I
I
- - - . - - - - - - - -I r - - - - - - - -
I
C 2C
x_
FIGURA 32
( 8) u (s) = El ro
e-B,"y(x) dx.
= Wo (~ _
B
l + e -: ) •
e s S2 s-
Por tanto la ecuación diferencial (6) se transforma en
http://carlos2524.jimdo.com/
258 Apl icocí ones [Cap. 13
e s S2 S2 '
de la cual obtenemos
( 10)
EJERCICIOS
d8
(B) cuando t = O, 8 = 80 Y = dt =
ti) <Uo·
CAPÍTULO 14
Sistemas de ecuaciones
EJERCICIOS
En los ejercicios del 1 al 8, aplíquese el método de la transformada de
Laplace.
1. x" (t) - 3x'(t) - y'(t) + 2y(t) = 14t + 3,
x'(t) - 3x(t) + y'(t) = 1; x( O) = O, x'(O) =0, y(O) = 6.5.
SOL. x(t) = 2 - iet - ie 3t _ e- 2t ,
y(t) = 7t + 5 - e t + %e- 2t .
2. 2x'(t) + 2x(t ) + y'(t) - y(t) = 3t,
x'(t) + x(t) + y'(t) + y(t) = 1; x(O) = 1, y(O) = 3.
SOL. x(t) = t + 3e- t - 2e-3t ,
y(t)= 1 - t + 2e-3t
3. x'(t) - 2x(t) - y'(t) - y(t) = 6e 3t ,
2x'(t) - 3x(t) + y'(t) - 3y(t) = 6e 3t ;x(0) = 3,y(0) = O.
SOL. x(t) = (1 + 2t)e t + 2e~t
=
y(t) (1 - t)e t _ esto
4. x"(t) + 2x(t) - y'(t) = 2t + 5,
x'(t) - x(t) + y'(t) + y(t) = -2t - 1; x(O) = 3, x'(O) = O,
y(O) = -3. SOL. x(t) = t + 2 + e-2t + sen t,
y(t) = 1 - t - 3e- 2t - cost.
5. Las ecuaciones del ejemplo de la sección 76, página 261, con las
condiciones iniciales x(O) = O, x'(O) = O, y(O) = 1, y'(O) = O.
SOL. x(t) = -t - i sen t + ¡sen 2t,
y(t) = 1 + ~ cos t - ~ cos 2t.
6. Las ecuaciones del ejemplo de la sección 76, página 261, con las
condiciones iniciales x(O) = 9, x'(O)= =
2, y(O) 1, y'(O)= O.
SOL. X (t) = - t + 15 cos t - 5 sen t - 6 cos 2t + 4 sen 2t,
y(t) = 1 + 2 cos t + 6 sen t - 2 cos 2t - 3 sen 2t.
7. x"(t) + y'(t) - y(t) = O,
2x'(t) - x(t) + z'(t) - z(t) = O,
x' ( t) + 3x ( t ) + y' (t) - 4y ( t) + 3z ( t ) = O;
x(O) = O,x'(O )= 1,y(0) = =
O,z(O) o.
8. x." ( t ) - x ( t) + 5y' ( t ) = f3 ( t ),
y"(t) - 4y(t) - 2x'(t) = O,
en la que f3(t) = 6t, O S; t S; 2,
= 12, t> 2;
x(O) = O, x'(O) = O, y(O) = O, y'(O)= O.
http://carlos2524.jimdo.com/
266 Si stema s de ecuaci ones [Cap. 14
..-
http://carlos2524.jimdo.com/
CAPÍTULO 15
Circuitos eléctricos y redes
78. cmCUITOS
( 1) LI' ( t) + RI ( t) + C-1Q ( t) = E (t ) ,
en la cual
(2) 1 ( t) = Q' (t ) .
Las ecuaciones (1) y (2), con las condiciones iniciales
(3) ICO)- = 0, Q(O) = 0,
constituyen el problema a resolver.
Para retener el símbolo convencional L para el número de henrys de
la inductancia del circuito, denotaremos en esta sección por Lt el ope-
rador de Laplace para el cual se usa L en todas las otras secciones del
libro.
http://carlos2524.jimdo.com/
s (sL
1+ R) = R1(1 - ~
1)
s + RL-1 ,
tal que
i(s) = ~G - + kL-1)
s [1 - exp (-toS)].
Entonces
FIGURA 35
y tenemos, de (8)
( 11 )
il( S) !i Ra + S(L2 + La)
- S L zL3(S + al ) (S + az ) .
El miembro derecho de la ecuación ( 11) tiene un desarrollo en frac-
ciones parciales
http://carlos2524.jimdo.com/
(12)
tal que
Il(t) Ao + Alexp (-alt) + A2 exp (-a2t).
EJEMPLO b): Para la red mostrada en la figura 36, escribir las ecua~
ciones para la determinación de las corrientes 11, 1 2 e 13, y la carga Qa.
Suponer que cuando el interruptor está cerrado todas las _corrientes
y las cargas son cero. Obtener el problema transformado.
Usando las leyes de Kirchhoff escri-
bimos las ecuaciones
(13) 11=12+13,
dl2
( 14 ) Rl I 1 + L -;¡¡
2 = E sen wt,
(16) 13 = dQ3.
dt
Nuestro problema consiste de las cuatro ecuaciones, de la ( 13) a la
( 16) con las condiciones iniciales
(17) 12(0) = O, 13(0) = O, Qa(O) = O.
Sea L{Ik(t)} =
ik(s) , k =
1,2,3, Y L{Qa(t)} Entonces el = qa(S).
problema transformado es el de resolver el sistema algebraico
(18) i1 - iz - i3 = O,
(21) is = sqa.
EJERCICIOS
1. Para el circuito RL de la figura 34, página 271, encontrar la corrien-
te I si el elemento de corriente directa E, no se quita del circt!ito.
SOL. 1 = ER-1[1 - exp ( -RtL-l)].
http://carlos2524.jimdo.com/
En los ejercIcIOs del 8 al 11, úsese la figura 33, página 270, con
E (t) = E sen (¡Jt y con las siguientes notaciones empleadas para
simplificar las fórmulas
_ Ji. b2 _ 2 __ 1_ f.l2 __ 1_ _ 2
a - 2L' - a LG , 1" - LG a ,
FIGURA 38
FIGURA 39
http://carlos2524.jimdo.com/
278 Circuitos eléctricos y redes [Cap. 15
15. Para la red de la figura 39, escríbanse las ecuaciones para la determi-
'n ación de la carga Qs y las corrientes 11 , 12 , 1:,¡. Supóngase que todas
estas cantidades son cero al tiempo cero. Transfórmese el problema
a la forma algebraica.
FIGURA 40
16. Para la red de la figura 40, encuéntrense las ecuaciones para la deter-
minación de las corrientes. Supóngase que todas las corrientes son
cero al tiempo t = O. Encuéntrese i 1(s) =
L{11(t)} y analícese el ca-
rácter de 11 (t) sin hallar explícitamente la función.
http://carlos2524.jimdo.com/
CAPÍTULO 16
Variación de paráIlletros
y otros Illétodos
279
http://carlos2524.jimdo.com/
280 Variación de parámetros y otros métodos [Cap. 16
Por supuesto
yc = Cl COS X + C2 sen x.
Vamos a buscar una solución particular por el método de variación de
parámetros. Hagamos
(9) y = A cosx + B sen x)
de lo cual se obtiene
y' = -Asen x + B cos x + A' cos x + B' sen x.
En seguida hagamos
(10) A' cosx + B' sen x = 0,
tal que
y' = - A sen x + B cos x.
Entonces
(11) y" = - A cos x - B sen x - A' sen x + B' cos x.
Ahora eliminaTemos y combinando las ecuaciones (9) Y (11) con la.
ecuación original ( 8). Obtenemos entonces la relación
(12) -A' sen x + B' cosx = secx tanx.
De (12) Y (10) se elimina fácilmente A'. El resultado es
B' = tan x
así que
( 13) B = Insecx)
en la que la constante arbitraria ha sido descartada porque estamos
buscando sólo una solución particular para sumarla a nuestra función
complementaria yc) previament~ determinada.
De las ecuaciones (12) Y (10) se sigue también fácilmente que
A' = - sen x secx tan x)
o
A' = - tan 2 x.
Entonces
A = -Jtan2xdx = J(1 - sec 2 x) dx)
así que
( 14 ) A = x - tan x,
otra vez desdeñando la constante arbitraria.
http://carlos2524.jimdo.com/
282 Variación de parámetros y otros métodos [Cap. 16
imponemos la condición
(18) A'e" + B'e 21t = O.
Entonces
(19)
de lo que se sigue que
(20) y" = Ae" + 4Be + A'e" + 2B' e
21t 2
"
1
A'e'" =
A' =
Entonces
Similarmente
así que
o
B = - e-'" + In (1 + e-"') .
Entonces, de (17),
yp = e'" In (1 + e-4:) - e'" + re 2
'" In (1 + e-"').
El término ( - e"') en J!p puede absorberse en la función complementa-
ria. La solución general de la ecuación (16) es
y = Cae'" + c2e2'" + (e'" + e,2"') In (1 + e-"')
La ecuación (16) puede resolverse igualmente sin el método de
variación de parámetros, empleando adecuadamente el método de cam-
bio de exponencial de la sección 48.
y si elegimos
EJERCICIOS
En los ejercicios 1-20, úsese el método de variación de parámetros,
sección 80.
1. (D 2 + l )y = csc x cot X .
SOL. y = él COS X + C2 sen x - x sen x - cos x In sen x.
2. (D2 + l )y = cotx. SOL. y = Yc - senxln (cscx + cotx).
3. (D2 + l )y = sec x .
SOL. y = Cl COS X + C.2 sen x + x sen x + cos x In cos x.
4. (D2 + l )y = sec 2 x.
S O L. y = Yc - 1 + sen x In (sec x + tan x) .
5. (D 2 + l)y = sec3 x. SOl;. Y = Yc + t sec x.
6. (D2 + l)y = sec4 x.
SOL. y = Yc - ~ + ~ sec 2 x + fr sen x In (sec x + tan x) .
7. (D 2 + l )y = tan x. S OL . y = Yc - cos x ln (sec x + tanx ).
8. (D2 + l)y = tan 2 x. . ~"
S'ÜL. Y = YI - 1, donde Yl es la solución del ejercicio 4.
9. (D2 + l )y = sec x csc X .
SOL. y = Yc - cos x In (sec x + tan x) - sen x In (csc x + cot x) .
10. (D2 + l )y = sec 2 xcscx.
SOL. y = Yc - sen x In (csc 2x + cot 2x ).
11. (D2 - 2D + 1))> = e~Z(e:t + 1)-2. SOL. Y = Yc ,+ e"' ln ( 1 + eZ).
12. (D2 - 3D + 2 )y = e /(l + e2"').
2Z
SOL. Y = Yc + e'" Arctan (e-"') - ie 2'" In (1 + e- 2"').
13. (D2 - 3D + 2)y = cos (e-"') . SOL. y = Yc - e2"cos (e-") .
14. (D2 - l )y = 2(1 - e- )-!.2Z
SOL. y = Cle'" + C2e-Z - e'" Arcsen (e-Z) - ( 1 - e- 2Z ) •.
15. (D2 - 1) Y = e- " sen e-z.
2 SOL. y = Yc - sen e-" - e'" cos e-Z.
16. (D2 - 5D + 4)y = 6/( 1 + e- 2").
SOL. y = Yc + 2e" Arctan e--" - e2" + é" In (1 + e- 2") .
17. (D 2 - l)y = (1 + e-Z)-2.
SOL. y = Yc - 1 + xe-" + e-X In ( 1 + e-Z) .
18. (D - 1) (D - 3 )y = cos e- • z
(D + 1) (D - 1) Y =
e" +2 e_ , Z
hacemos
(5) (D-1 )y = w.
Entonces necesitamos resolver
2
(D + 1) w = e" + e-z '
que puede escribirse como
2 dx
(6) dw + w dx = "
e
_,,'
+e
'Un factor integrante de (6) es e". Entonces
2e" dx
e"dw + we"dx=---
e" + e-"
es exacta, por tanto
2e" dx = J 2e 2
"
e"w =
Je" + e-" e2 " + 1 dx'
de lo cual se tiene
e"w = In (e 2" + 1),
como de costumbre, la constante arbitraria será desdeñada ya que
sólo estamos buscando una solución particular.
En seguida insertamos la expresión para w en la ecuación (5), la cual
se transforma entonces en
( D - 1) Y = e-Z In (e 2" + 1) ,
dy - ydx = e-"In (e2 " + 1) dx.
Aquí un factor integrante es e-", de tal forma que escribimos
e-Z dy - ye--!.C dx = e- 2" In (e 2" + 1) dx,
e-"y - Je- 2" In (e 2" + 1) dx.
Se puede integrar por partes para reducir la anterior expresión a
la forma
(7) e-"y = - ie- 2 " In (e2" + 1) + J dx
e 2" +1
Entonces escribimos
2
e-"y = - ie- " In (e
2
" + 1) + J[1 - e2"e: 1J dx
http://carlos2524.jimdo.com/
§84] Tercer método, cambio de variable dependiente 289
EJERCICIOS
EJERCICIOS VARIOS
3. (D2 - 1) y = e "'(3 tan e" + e'" sec e"). SOL. y = Ye + e'" In sec e"'.
2 2
CAPÍTULO 17
La existencia de soluciones
----------~~~~~~--------~'x
FIGURA 41
tiene una solución que pasa por el punto (xo, yo) y que esa solución es
única cerca de (xo, yo) .
Es importante que la solución bajo discusión pueda realmente exhi-
birse. La función y(x) es el límite de una sucesión de funciones y,,(x);
n = 1, 2 . . " definida como sigue:
y = y(x) = Limy,,(x) ,
"....'"
donde
YI(X) = yo ·+ [f(t,yo) dt,
Y2(X) = yo + r"'o
f[t,YI(t)] dt,
k=l 2k In k x
Vn(X) = 1 + ~ -kl-'
y que
y = exp (2 In x) = X2.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
CAPÍTULO 18
Ecuaciones de orden uno
y Illayor grado
Si una ecuación
( 1) f( x, y, y') = O
es de grado mayor que uno, puede ser posible factorizar
el miembro izquierdo en partes, cada una de grado uno. El
problema de resolver (1) se reemplaza entonces por dos o
más problemas' de tipo más simple. Estos últimos deben poder
resolverse por los métodos de los capítulos 2 y 4.
Ya que y' deberá ser elevada a potencias tanto en el ejem-
plo como en los ejercicios, vamos a simplificar la nomencla-
tura y escribir de una manera común usando p por y' :
P= dy
dx·
EJEMPLo: Resolver la ecuación diferencial
(2) xyp2 + (x + y)P + 1 = O.
El miembro izquierdo de la ecuación (2) es fácilmente fac-
torizable.
299
http://carlos2524.jimdo.com/
300 Ecuacion es de orden uno y mayor grado [Cap. 18
se sigue que
ydy + dx = °
(5)
La ecuación (4) puede ser escrita como
(2) xyp2 + (x + y) P + 1 = O
tenía la solución general
(5 )
y
(6) y = - In (C2X).
19. Demuéstrese que de la familia (5) anterior, la única curva que pasa
por el punto (1, 1) es y = (3 - 2x) i Y que esta solución es válida
para x < ~.
20. Demuéstrese que de la familia (6) anterior, la única curva que pasa
por el punto (1, 1) es y = 1 - ln .x y que esta solución es válida
para 0< x. .
--------------~t_----~~~----- x
FIGURA 42
21. Demuéstrese que si una solución de (2) es puesta junto a las solu-
ciones que provienen de (5) Y (6), entonces las pendientes de las
curvas deben ser iguales donde las partes se unen. Demuéstrese enton-
ces que las partes deben unirse en un punto de la recta y = x .
22. Demuéstrese que la curva determinada por
y = (3 - 2x )" para x ~ 1,
Y = 1 - In x para 1 ~ x
(2) ydy = dx
Va 2 - i '
o
ydy
(3) . / 2 2 = dx,
va - y
o (si la división entre V a2 - y2 no puede efectuarse)
(4)
De (2) se SIgue que
( 5)
Y de (4) que
(7) y= a o y = -a.
Gráficamente, las soluciones (5) son semicircunferencias izquierdas
con radio a y centradas sobre el eje x; las soluciones (6) son semicir-
cunferencias derechas de radio a, centradas sobre el eje x. Podemos
combinar (5) Y (6) en
(8)
que es llamada la solución general de (1).
http://carlos2524.jimdo.com/
304 Ecuaciones de orden uno y mayor grado [Cap. i 8
--~----~--t--r--rt4---r-+-~~---i---X
FIGURA 43
(1) ¡(x, y, p) = O; P = dy
dx'
en la cual el miembro izquierdo es un polinomio en x, y y p. Puede no
ser posible factorizar el miembro izquierdo en factores que sean ellos
mismos polinomios en x, y y p. Se dice entonces que esta ecuación es
irreducible.
http://carlos2524.jimdo.com/
§90] Ecuación con discri minante e 305
(2) tp(x,y, e ) = O.
( 3)
(5) Otp = O
oc .
Ahora tenemos dos ecuaciones, tp = O Y :: = O, las cuales deberán ser
satisfechas por x, y y e. Estas dos ecuaciones deberán tomarse como
las ecuaciones paramétricas deseadas. Ellas contienen cualquier envol-
vente de la familia original de curvas, tp = O, que pueda existir. Afor-
tunadamente, no necesitamús expresar estas ecuaciones en la forma
x = x(e), y = y (e).
http://carlos2524.jimdo.com/
306 Ecuaciones de orden uno y mayor grado [Cap. 1 B
de los puntos (X, y) para los cuales (1 ) tiene al menos dos de sus
raíces iguales.
La condición de que la ecuación (1) tenga al menos dos raíces
iguales como una ecuación en p es que tanto I= O como ~~ = O. Estas
dos ecuaciones en las tres variables x, y y p son las ecuaciones paramé-
tricas de una curva en el plano xy donde p desempeña el papel de un
parámetro. La ecuación que resulta cuando se elimina p de las eoua-
ciones paramétricas I = O Y 01 = O se llama la ecuación con discrimi-
nante p. op
Si una envolvente de la solución general de I = O existe, estará con-
tenida en la ecuación con discriminante p. No se da aquí ninguna
prueba de esto último. * Para nosotros, la ecuación con discriminante
p es útil en dos formas. Cuando se obtiene una ecuación singular en el
transcurso de la solución de una ecuación, como ocurre en las secciones
de la 92 a la 94, la ecuación con discriminante p, nos sirve como mé-
todo de comprobación. Si ninguno de nuestros métodos de tratamiento
conduce a una solución general, entonces la ecuación con discriminante
p ofrece curvas que pueden ser soluciones particulares (incluyendo la
solución singular) de la ecuación diferencial. Entonces, la ecuación con
discriminante p puede ser probada para soluciones posibles de la ecua-
ción diferencial. Tales soluciones particulares no contribuyen, por su-
puesto, al proceso de encontrar la solución general.
La ecuación con discriminante p puede contener soluciones singulares,
soluciones que no son singulares y curvas que no son soluciones en todo
su dominio.
EJERCICIOS
1. Para la ecuación cuadrática
f = Ap2 + Bp + C = 0,
donde A, B, C son funciones de x y y, demuéstrese que el discriminante
p, obtenido al eliminar p de I = °y op01 = 0, es la ecuación conocida
B2 - 4AC = O.
2. Para la ecuación cúbica p3 + Ap + B = 0, demuéstrese que la ecua-
ción con discriminante p es 4As + 27B2 = o.
* Para detalles mayores sobre las soluciones singulares y los discriminantes, véase
E: L. Ince, Ordinary Difterential Equations. L ondres: Longmans-Green, 1927, pá-
gmas 82-92.
http://carlos2524.jimdo.com/
308 Ecuaciones d e orden uno y mayor grado [Cap . 18
(8) 2p = x~~.
Primero consideremos (8), la cual nos conduce a
así que
(9)
Entonces las ecuaciones (4) y (9), con p como parámetro, constituyen
una solución de (4), considerada esta ecuación como una ecuación
diferencial con p = dy/ dx.
En este ejemplo es fácil eliminar p de las ecuaciones (4) y (9). Efec-
tuando esta eliminación tenemos
x . C2 X4 - 3y . cx 2 + 9X2 = O.
Descartando por el momento la posibilidad de que x = 0, podemos
representar nuestra solución como
2
C X3 - 3cy +9= 0,
o
3cy = 2
C X3 + 9.
A,hora hagamos e = 3k para obtener
(10)
La ecuación 10, con k como una constante arbitraria se llama la
solución general de la ecuación diferencial (4).
http://carlos2524.jimdo.com/
x)dP
1 ( x+--
p=p--+
P p2 dx'
o
l=X(1
P
+~)dP.
p2 dx
Por tanto hemos llegado a la ecuación diferencial
http://carlos2524.jimdo.com/
312 Ecuaciones de orden uno y mayor grado [Cap. 18
dx = p2 + 1 dp
X P ,
de la cual
In x = tpz + In p + In c,
o
(14) x = cpexp (iPZ).
De la ecuación (13) obtenemos
y=x(PZ-l)/p.
Entonces, empleando (14), encontramos y en la forma
(15)
Las ecuaciones (14) Y (15) son ecuaciones paramétricas de la solu-
ción g.eneral de la ecuación diferencial original. Al bosquejar una curva
de la familia de soluciones, para cada valor de p insertado en ( 14) Y
(15), se obtienen un par de valores de x y y. Por tanto, para cada
p empleada obtenemos un punto sobre la curva deseada.
El parámetro p podría eliminarse resolviendo la cuadrática ( 13)
para p sustituyéndola en ( 14). El resultado podría ser más difícil de
bosquejar que la misma curva en forma paramétrica, ( 14) y (15).
El lector deberá usar la ecuación con discriminante p para demos-
trar que la ecuación ( 13) no tiene una solución singular.
Sucede frecuentemente, como se ve en las respuestas para algunos
de los ejercicios siguientes, que es inconveniente obtener y en términos de
p y c solamente. Puede ser mejor dejar la solución en la forma
(16 ) rp (x, p) = O
Y
(17) y = g ( x, P).
Entonces, en los cálculos, elegimos p y determinamos x de la ecuación
( 16), después de lo cual puede encontrarse y de la ecuación (17),
usando los valores tanto de x como de p.
(3) dp = O
dx
o
(4) x + f'(P) = 0.
La solución de la ecuación diferencial (3) es, por supuesto, p = c,
donde c es una constante arbitraria. Volviendo a la ecuación diferencial
( 1 ), podemos ahora escribir su solución general como
(5) y = cx + f (c) ,
y el resultado puede verificarse fácilmente por sustitución directa en la
ecuación diferencial ( 1). Nótese que (5) es la ecuación de una familia
de líneas rectas.
Considérese ahora la ecuación (4). Ya que f(P) y f'(P) son fun-
ciones conocidas de p, las ecuaciones (4) y ( 1)·, juntas constituyen un
conjunto de ecuaciones paramétricas que dan x y y en términos del
parámetro p. De hecho, de la ecuación (4) se sigue que
( 6) x= -f' (P),
la cual, combinada con la ecuación (1), nos da
(7) y = f (P) - pf'( P) .
Si f (P) no es una función lineal de p ni una constante, puede demos-
trarse (ejercicios 1 y 2 siguientes) que (6) y (7) son ecuaciones para-
métricas de una solución no lineal de la ecuación diferencial ( 1). Ya
que la solución general (5) representa una recta para cada valor de c,
la solución (6 ) y (7) no puede ser un caso especial de (5); es una
solución singular.
http://carlos2524.jimdo.com/
314 Ecuaciones de orden uno y mayor grado [Cap. 18
________________ ~~~~~~---x
FIGURA 44
X2p2 _ 2xyp + y2 - 1 - p 2 = O.
Entonces, es claro que la ecuación es de la forma
(12)
y que puede ser separada en dos ecuaciones, cada una de la forma de
Clairaut. Entonces, la solución general de (11) se obtiene reempla-
zando p en todas partes por una constante arbitraria c. Esto es,
( 13) (X2 - 1) c2 - 2xyc + y2 - 1 = O
es la solución general de (11). La solución (13) está compuesta de dos
familias de rectas,
(14)
Y
(15)
De la ecuación con discriminante p para (11) obtenemos inmediata-
mente la solución singular
(16)
EJERCICIOS
1. Sea ex un parámetro, pruébese que si f" (a) existe, entonces
(17) x =- f'(a) , y = fea) - af'(a)
es una solución de la ecuación diferencial y = px + f (p). Sugeren-
cia: Empléese dx y dy para obtener p en términos de y entoncesª
demuéstrese que y - px - f(p) se anula idénticamente.
2. Probar que si f" (a) =F 0, entonces (17) no es un caso especial de
la solución general y = ex + f (e). Sugerencia: Demuéstrese que la
pendiente de la gráfica de una solución depende de x, mientras que
la pendiente de la gráfica de la otra no depende de x.
En los ejercicios del 3 al 30, encontrar la solución general y tam-
bién la solución singular, si ella existe.
3. p2 + x 3p - 2x"2y = O. SOL. c2 + cx2 = 2y; sol. sing., 8y = _ x 4 •
4. p2 + 4x p - 12x y = O.
5 4
SOL. 12y = =_
c(c + 4x 3 ) ; sol. sing., 3y x6 •
5. 2Xp3 - 6yp2 + x 4 = =
O. SOL. 2c Sx 3 1 - 6c 2p; sol. sir{g., 2y = X2.
6. p2 - xp + y = O. SOL. y =
ex - c2; sol. sing., X2 = 4y.
7. =
Y px + kp2. SOL. Y =
ex + kc 2; sol. sing., X2 = -4ky.
8. X p2 + 3xp + 9y = O.
8
SOL. XS(y + c2 ) + e = =
O; sol. smg., 4x 6y 1.
http://carlos2524.jimdo.com/
316 Ecuaciones de o rden uno y mayor grado [Cap. 18
(1 ) f (x,y,P) =O;p=ix,
dx 1
dy -p
Esto no es, por supuesto, nada completamente nuevo. Ya que (1)
es una ecuación de primer orden, los papeles de las variables depen-
diente e independiente pueden intercambiarse mediante el uso de la
relación
~~ = (~;rl.
EJEMPLO: Resolver la ecuación diferencial
(2) yp2 _ Xp -t- 2p = O.
Primero, resolviendo para x se obtiene
(3) x = yp +-2y
p
,
diferenciando entonces ambos miembros de la ecuación (3) con res-
pecto a y, se obtiene
! + p + Y (1 -~)
P
dp = O.
p2 dy
Esto nos conduce a la ecuación
(4)
de la que
De aquí que
In y = - %nI ( 1 + P2) + 2 In p + In C J
así que
(5) y = cp<¿( 1 + p2)_3/2.
De la ecuación (3) se sigue que
x = p-ly (P2 + 2).
Entonces, empleando (5) se puede escribir
( 6)
http://carlos2524.jimdo.com/
318 Ecuaciones de orden uno y mayor grado [Cap. 18
EJERCICIOS
EJERCICIOS VARIOS
9. 4x sp2 + 12 x yp + 9 =
4 O.
SOL. x 3(2ey - 1) = e2; sol. sing., X3y2 = 1.
10. 4y2p3 - 2xp + y = O.
SOL. y2 = 2e(x - 2e 2) ; sol. sing., 8x 3 = 27y4.
11. p4 + xp - 3y = O. SOL. 5x = 4p3 + ep! y 15y = 9p4 + ep ~ .
12. Hágase el ejercicio 11 por otro método.
13. Xp2 + (k - x - y) p + y = O, la ecuación del ejercicio 5, página 25.
SOL. xe 2 + (k - x - y) e + y = O;
" sol. sing., (x - y) 2 - 2k(x + y) + k2 = O.
14. X p3 - 2xyp2 + y2p + 1 = O.
2
S OL. x 2 e3 - 2xye2 + y 2 e + 1 O; sol. sing., 27x
= = _4 y3.
15. 16xp2 + 8yp + y6 = O. SOL. r(e 2 x + 1) = 2e; soL sing., xy4 = 1.
16. Xp2 - (x 2 + I)P + x = O. SOL. X2 = 2(y - el )'y = In (e2x) .
17. p3 - 2xp - y = O. SOL. 8x = 3p2 + ep-~ y4y = p3 - epI".
18. Hágase el ejercicio 17 por otro método.
19. = =
9xy4p2 - 3ysp - 1 O. SOL . ey3 e2x - 1; sol. sing., y6 = -4x.
20. X2p2 - (2xy + l )p + y2 + 1 = O.
SOL. x 2e2 - (2xy + 1) e + y2 + 1 = O; sol. sing., 4X2 - 4xy - 1 = O.
21. X6p2 = 8(2y + xp) . SOL. e2x2 = 8 (2x2y - e ) ;sol. sing., x4 y= -1.
22. X2p2 = (x - y) 2 SOL. x(x - 2y) = el, Y = - x In (e·2x).
23. Xp3 - 2yp2 + 4X2 = O. Véase el ejercicio 10 anterior.
SOL. X2 = 4e(y - Be 2); sol. sing., Bya = 27x4 •
24. (p + 1)2(y - px) = 1.
SOL. (e + 1)2(y - ex) = l;sol.sing., 4(x + y)3 = 27x 2.
25. p 3 - p2 + xp - y = O.
SOL. y = ex + e3 - e2; sol. sing., con ecuaciones paramétricas,
x= 2a - 3a2y y = g2 - 2aS,.
26. Xp2 + y(1 - x)p - y2 = O. SOL . xy = el, x = In (e2Y)'
27. yp2 - (x + y) P + y = O.
SOL. py = e exp (p-l) Y px = Y(P2 - P + 1) ; sol. sing., y = x.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
CAPÍTULO 19
Ecuaciones especiales de orden dos
x ~: - p3 - p. = O.
La solución de (6) es
Y - C 2 -- -(c1 2
-+- - X2)~.,
o
(7)
La ecuación (7) es la solución general deseada de la ecuación dife-
rencial (3). Nótese que al dividir entre p al principio del proceso
hemos descartado las soluciones y = k(p = O), donde k es una constante.
http://carlos2524.jimdo.com/
Variable independ iente perdida 323
§96]
d P)_
(2) f ( y, p, P dy - O.
" _ P dp
y - dy
como antes, entonces la ecuación (3) se transforma en
(4)
en la cual las variables p y y son fácilmente separables.
De (4) se sigue que
http://carlos2524.jimdo.com/
324 Ecuaciones especiales de orden dos [Ca p . 19
pdp dy _ O
p2+1+-Y - ,
de la cual
fIn (p2 + 1) + In Y = In C1,
así que
(5 )
Resolvemos la ecuación (5) para p y encontramos que
( C1"2 _ y2P
p = + .
y •
Entonces
o
+y( C1 2 - y2)-i dy = dx.
Por tanto
-+- (C1' 2 - y2)!i = X - C2.
de donde obtenemos el resultado final
(x - · C2)".l + y2 = C1 2•
EJERCICIOS
2
3. x y" + (y' ) 2 - 2xy' = 0, cuando x = 2, y = 5 Yy' = 2.
SOL. X2 = 2 (y - 3 ). ./
4. yy" + (y' ) 2 = O. - SOL. Véase el ejercicio 20, página 26.
5. y2y" + (y')a = O. SOL. =
X C1Y - In (C2Y).
6. (y + 1) y" = (y') 2. . / SOL. Y+ 1 =
C2eC1"'.
7. 2ay" + (y')a = O. SOL. (y - c2)2 = 4a (x - C1) '
8. Hágase el ejercicio 7 por otro método.
9. xy" = y' + x 5, cuando x = 1, y = i y y' = 1.
SOL. 24y = x 6 + 9x2 + 2.
10. xy" + y' + x = O,cuandox = 2,y = -1 yy' = - l
SOL. y = - iX2 + In (ix) .
11.Y" = 2y(y') 3. . . _ / SOD . y3 = 3(C2 - X - C1Y)'
12. yy" + (y' ) 3 - (y')"2 O.= SOL. =
X Y - c1 In (C 2Y)'
13. y" + f3 2y = O. Compruébese el resultado resolviendo en dos formas
la ecuación.
14. yy" + (y') a = O. SOL. =
X C1 + y In (C2Y).
15.y"cosx = y'. SOL. y = C2+ c1 In (1 - senx).
16. y" = x (y' )2, cuando x = 2, y = i7T y y' = - i. OL. x = 2 cot y.
http://carlos2524.jimdo.com/
1+
2
SOL. Y = 1 + X2 - In - 2 x-) .
(
38. y" =2x + (X2 - y')2. SOL. 3y = x 30 + C2 - 31n (x + Cl) .
39. (y") 2- 2y" + (y') 2- 2xy' + .'\:'2 = 0, cuando x =O,y = !yy'= 1.
SOL. 2y = 1 + X2 + 2 sen x.
40. (y") 2 - xy" + y' = O.
SOL. Solución general: 2y = ClX 2 - 2C12X + C2;
familia de soluciones singulares: 12y = x 3 + k.
41. (y" ) 3 = . 12y' (xy" .- 2x' ) .
SOL. Solución general: y = Cl( X - Cl)3. C2; +
familia de soluciones singulares: 9y = x' + k.
42. 3yy'y" = (y') 3 - 1. SOL. 27cl(y + Cl )2 = 8(x + C2)3.
43. 4y (y' ) 2y" = (y' ) 4 + 3. sou. 256cl(y - Cl )3 = 243 (x - C2)4.
http://carlos2524.jimdo.com/
326
97. CATENARIA
Supongamos que se tiene un cable de peso uniformemente distri-
buido w (lb por pie) suspendido entre dos soportes A y B, como se
indica en la figura 45. El cable se combará y tendrá su punto más
bajo en V tal como lo indica la figura. Deseamos determinar la curva
formada por el cable suspendido. La curva se llama catenaria.
FIGURA 45
------+------------------------x
o
FIGURA 46
Supóngase que s representa la longitud (en pies) del cable medida des-
de V al punto variable P con coordenadas (x, y). Entonces la porción
del cable desde V hasta P está su jeta a las tres fuerzas mostradas en
la figura 46. Estas fuerzas son: a) la fuerza gravitacional ws (lb)
http://carlos2524.jimdo.com/
Catenaria 327
§97]
(3) tan O = dy
dx'
Podemos eliminar la tensión variable T 1 de las ecuaciones (1) Y
( 2) Y obtener
(4)
(5) tan O =~ .
a
(8)
http://carlos2524.jimdo.com/
CAPÍTULO 20
El método de series de potencias
y
r
Jo
Z
'f (y) dy = f
n={)
anx"+1,
n + 1
-R < X < R.
esto es, la serie es correctamente diferenciable e integrable en el sentido
de que la serie de las derivadas de los términos separados converge a la
derivada de la suma de la serie original, y de modo semejante para
la integración. Es importante que el intervalo de convergencia perma-
nezca inalterable. No estamos tratando aquí con comportamientos de
convergencia en los puntos extremos de ese intervalo.
Examinemos ahora más de cerca la razón por la cual una serie tiene
un intervalo de convergencia particular más bien que otro. Un ejemplo
elemental del cálculo es
ca
1
( 1)
1- x
~ x",
n=()
-1 < x < 1.
Es razonable sospechar que el mal comportamiento de la función
1/ ( 1 - x) en x = 1 es el que determina el hecho de que el intervalo
de convergencia termine en x = 1. Que
b este intervalo puede extenderse más allá
de este punto no ~ evidente. El punto
x = - 1 no se modificará en esta oca-
sión; el intervalo se termina en ese
extremo debido a la condición de que
debe ser simétrico alrededor de x = O.
--r---------r-------~r--a
Reemplacemos ahora x en ( 1) por
( - X2) para obtener
( 2)
-j
- l<x<1.
FIGURA 47 el intervalo de convergencia es otra
vez - 1 < x < 1, pero la función
f (x) = 1/( 1 + X2) se comporta adecuadamente para toda x real. N o
es claro qué es lo que hace que el intervalo de convergencia termine en
x = 1. El hecho es que necesitamos considerar a x como una variable
compleja, para entender qué es lo que hay que hacer aquí.
Usaremos el diagrama ordinario de Argand para números complejos.
Sea x = a + ib, con a y b reales, con i = Y=1, y asociado con el
punto ( a, b ) en el plano, el número a + ib. Ahora márquense sobre
http://carlos2524.jimdo.com/
§101] Puntos ordinarios y puntos singulares 333
EJERCICIOS
Para cada ecuación, enlistar los puntos singulares en el plano finito.
1. (X2 + 4)y" - 6xy' + 3y =O. SOL. X = 2i, -2i.
2. x(3 - x)y" - (3 - x)y' + 4xy = O. SOL. X = 0,3.
3. 4y" + 3xy' + 2y = O. SOD. Ninguno.
4. x(x - 1)2y" +3xy' + (x - l)y = O. SOL. X = O, 1.
5. x 2y" + xy' + (1 _ X2)y = O.
6. x4y" + Y = O.
7. (1 + X2)y" - 2xy' + 6y = O.
8. (X2 - 4x + 3)y" + xy - 4y = O.
9. x 2(1 - x) 3y" + (1 + 2x)y = O.
10. 6xy" + (1 - X2)y' + 2y = O.
http://carlos2524.jimdo.com/
§102] Val idez de los solucio"'es cerco de un punto ord inario 335
11. 4xy" + y = O.
12. 4y" + Y = O.
13. X2(X2 - 9)y" + 3xy' - y = O.
14. x2(1 + 4x 2)y" - 4xy' + y = O.
15. (2x + 1) (x - 3)y" - y' + ( 2x + l)y = O.
16. X3(X2 - 4PY' + 2 (x 2 - 4 )y' - xy = O.
17. X(X2 + l )2y" - Y = O.
18. (X2 + 6x + 8 )y" + 3y = O.
19. (4x + l )y" + 3xy' + y = O.
( 2)
que contiene dos constantes arbitrarias, por ejemplo ao y al, y que con-
verge dentro de una circunferencia con centro en x = O Y se extiende
hacia el punto singular (o puntos singulares) cercanos a x = O. Si la
ecuación diferencial no tiene puntos singulares en el plano finito, la solu-
ción ( 2) es válida para todo x finito. Queda, por supuesto, el asunto
de encontrar a n para n ~ 2. Esto constituye la mayor parte de nuestro
trabajo al resolver una ecuación particular.
Es necesario darse cuenta que el teorema citado establece que la
serie involucrada converge dentro de una cierta circunferencia. No
establece que la serie diverge fuera de esa circunferencia. En un ejem
plo particular puede suceder que la circunferencia de convergencia se
extienda más lejos que el mínimo dado por el teorema. En cualquier
caso, la circunferencia de convergencia pasa por un punto singular de
la ecuación - que puede no ser el punto singular más cercano.
(2)
o
00 00
(4 ) ¿ n (n -1 )anxn - 2 - ~ (n + l )( n + 4)anxn = O.
n=O n=O
/
http://carlos2524.jimdo.com/
338 El método de series de potencias [Cap. 20
. _ 4 . 6 . 8 . .. ( 2k + 2)
k ~ 1. a2a4aS ' •• a 2k - 2.4. 6 . . . (2k) a Oa2a4'" a.2k-2,
se simplifica de inmediato a
k ~ 1: a 2k = (k + 1) ao,
k ~ 1:
(2)
(11 )
y de (10) se obtiene
00 00 00
~ n (n - 1)anvn-2 + ~ (n - 4)anv'Hl = 0,
n=O n=O
as = - --ao
3.2
-4
a5 = - --a2 =
5 ·4
-2 .
°
-1
a6 = - --as
6·5
a1 = - -
7 .- °
6 a.- = ° as = () a5 = °
2
ag = - --a6
9·8
3
a10 = - - - a7
10·9
= ° all = °
3k - 7
ask = - 3k( 3k _ 1) aSk-3 a3/m = 0, k 2:: 2 aak+2 = 0, k 2:: 1.
(
14) Y = ao [1 + ~ ~(-=-,1)~/<['---'-(--::- . --::',---'(-'-3
-_4-'--.:)('------:-:1:-':-)-:-:;'-;:-;:2:- k _-_7--:-,)-!:,J(,---X_--:-:1-f)_3k]
:c _
.. =1 [3· 6 . 9 ... ( 3k) ] [2 . 5 . 8 ... (3 k - 1) ]
+ a1[(x - 1) + -!(x - 1)4J.
La ecuación diferencial original no tiene puntos singulares en el plano
finito, entonces la serie de (14) es convergente para toda x finita. En
los cálculos, 'Por supuesto, es de mayor utilidad en la vecindad del
punto x = 1.
El coeficiente de (x - 1) 3k es suficientemente complicado como para
garantizar que puede ser simplificado. En el producto 3 . 6 . 9 . . . (3k)
hay k factores, cada uno múltiplo de 3. Llegamos entonces al resultado
3 . 6 . 9 ... (3k) = 3/« 1 . 2 . 3 ... k) = 3/<k!.
Además todos, menos los dos primeros factores dentro del paréntesis
rectangular del numerador, aparecen también en el denominador. Con
un pequeño tanteo, probando con los términos k = 0, 1, 2, ya que los
factores que se cancelan no aparecen sino hasta k > 2, puede demos-
trarse que
'" 4( -l) /«x - 1)37<
(15) Y = ao ~o3/«3k _ 1) (3k -4)k! + al[ (x - 1) + i(x - 1)4J.
EJERCICIOS
co ( - 3) kX2k+l ]
+ al [ x + k~l 3 . .5 . 7 ... (2k + 1) ;
válida para toda x finita.
2. (1 + 4X2 )y" - By = O.
co ( -1 ) k+122kx2k+1
SOL. y = ao(1 + 4X2) + al ::0 4P _ 1
válida para Ixl < 1-
3. (1 +x 2
)y" - 4xy' + 6y = O.
SOL. Y = ao(1 - 3X2) + a¡(x - ix3);
válida para toda x finita.
4. (1 + x 2)y" + lOxy' + 20y = O.
SOL. y= ;0 k~O (-l)"k(k + 1) (2k+ 1) (2k + 3)x21c
+ ~1 k~ (- 1) k (k + 1) (k + 2) (2k + 3) x3 k+1;
válida para Ixl < 1.
5. (X2 + 4)y" + 2xy' - 12y = O.
00 3(-1)k(k+1)X2lc]
SOL. y = ao [ 1 + k~l 22k (2 k _ 1) (2k _ 3)
+ al(x + 3
l52X .); válida para Ixl < 2.
6. (X2 - 9 )y" + 3xy' - 3y = o.
[3 . 5 . 7 . . . (2k + 1) ]X 2k ]
SOL. y = ao [ 1 - ::1
co
(18)k(2k _ l )k!
+ a1X; válida para Ixl < 3.
7. y" + 2xy' + 5y = O.
_ [ ~ (-1 ) k[5 . 9 . 13 ... (4k + 1) ]X 21c ]
SOL. y - ao 1 + (:1 (2k) !
. [ ~ (-1)k[7. 11.15 ... (4k + 3)]X 2k+1]
+ al x + {:1 (2k + 1) ! ;
válida para toda x finita.
8. (X2 + 4)y" + 6xy' + 4y = o.
SOL. y ~ ao [ 1 + 22kel
co (-l)k(k + 1)x ]
2k
+ al [ x+ el 00 ( _ 1 ) k (2k + 3) X 2k+1]
3 . 22k
, .
; valIda para Ixl < 2.
9. 2y" + xy' - 4y = o.
SOL. + X2 + 1\ x 4 )
Y = ao(l
co 3 ( _ 1) kX2k+1
+ al eo 22k k! (2k - 3) (2k - 1) (2k + 1) ;
válida para toda x finita.
10. (1 + 2X2)y" - 5xy' + 3y O. =
_ [1 ~ 3 ( - 1) k[ ( - 1) . 3 . 7 . . . (4k - 5) ]x 2k ]
SOL. Y - ao + k~l 2kk!(2k _ 3) (2k _ 1)
+ al (x + ix 3
,) ; válida para Ixl < 1/ Y2.
http://carlos2524.jimdo.com/
§ 103] Soluciones cerca de un punto ordinario 343
+ el
00 (_1 )kX4k+1 ]
+ al [ x 22k k! .5.9. 13 . . . (4k + 1) ;
válida para toda x Jinita.
12. y" - 2 (x + 3 ) y' - 3y = O. Resuélvase alrededor de x = -3.
~ 3 . 7 . 11 ... (4k - 1) (x + 3) 2k]
" SOL. Y = ao [ 1 + {;1 (2k) ! J '
00 ( - 1) k (x - 2) 31<+1 ]
+ al [ (x - 2) + ~13kk![4. 7 . 10 .. . (3k + 1)] ;
válida para toda x finita.
14. (1 - 4x 2) y" + 6xy' - 4y = O.
SOL. Y = ao(1 + 2X2 ) + al [ x -
00
k! (4k2 _ 1) el
1. 5 . 9 . . . (4k - 3) X 2k+1]
;
válida para /x/ < t.
15. (1 + 2X2 )y" + 3xy' - 3y = O.
SOL. Y = alx + ao [ 1+ el
00 ( _ ) k+1 3 . 7 . 11 . . . (4k - 1) x 2k ]
2k( 2k _ l ) k! ;
válida para /x/ < 1/ V2.
16. y'" + x2y" + 5xy' + 3y = O.
00 ( -1 ) k X 3k ]
SOL. Y = ao [ 1 + k~l 2 . 5 .8 ... (3k - 1)
1) kx3l<+1]
00 ( _
+ al [ x + ~1
3kk!
+~ [ X2 + el
00 (_I)k x 3l<+2
4 . 7 . 10 ... (3k + 1) ;
]
SOL. .~ = - 15 2 k! k
-1 ) k (k 1) X2k+l 00 ( +
.. (2k
+ al k=o
1)
~
;
1 .3·5· +
válida para x finita.
18. Resuélvase la ecuación y" +
y = O tanto por series como por métodos
elementales y compárense los resultados.
http://carlos2524.jimdo.com/
344 El método de series de potencios [Copo 20
CAPÍTULO 21
Soluciones cerca de puntos
regulares singulares
P(x) = r(x)
x
donde r(x) es una función racional de x y r(x) existe en x = O. Sabemos
que tal función racional, esta r (x), tiene un desarrollo en serie de poten-
cias alrededor de x = O. Existe entonces el desarrollo
(3)
alrededor de x = O.
http://carlos2524.jimdo.com/
a cero los coeficientes de las potencias más grandes de x (x C-\ xc, xl+C, etc. )
en la identidad que está justamente antes de la ecuación (5).
Esta conclusión superficial es correcta si la diferencia de las raíces
Cl y C2 no es entera. Si la diferencia es entera, sin embargo, puede entrar
un término logarítmico en la solución. Las razones para este extraño
comportamiento se aclararán cuando desarrollemos un método para
obtener las soluciones.
(2)
~ ~
~ 2 (n + e) (n + e - 1) anXn+c- l +¿ n= (}
(n + e) anX,,+C-l
~ '"
+~ (n + e)anXMC - 2 ~ anX n+c = 0,
110 = 0 110=0
o
~ ~
\
http://carlos2524.jimdo.com/
352 Soluciones cerca de puntos regulares singulares [Cap. 21
(2n - 5)an-l
a" = -
2n(2n +
1) ,
así que el producto nos da, para n :2: 1,
( -1 )"[( -3 )( -1)(1) ... (2n - 5)]ao
( 6) n
a = [2 . 4 . 6 .. . (2n)] [3 . 5 . 7 ... (2n 1)]' +
La fórmula (6) puede simplificarse a la forma,
( 7)
_____-...:(~-_1.!.-
) n_,.:,..3a--=o---,--,--_ _
a" -
2"n! (2 n - 3) (2n - 1) (2n + 1) ,
Empleando ao = 1, la an de (7) ,. y el valor pertinente de e, el = i,
podemos ahora escribir una solución particular. Esta es
00 ( _ 1 )" 3x"+!
(8)
y1 = xi + "~1 2"n ! ( 2n - 3) (2n - 1)( 2n + 1) .
La notación y1 es para enfatizar que esta solución particular corres-
ponde a la raíz el de la ecuación indicial. Nuestro siguiente objetivo
será obtener una solución particular y2 correspondiente a la raíz más
pequeña e2. Entonces, la solución general, si se desea, puede escribirse
inmediatamente como
y = AY1 + By2
donde A y B son constantes arbitrarias,
Al volver a la relación de recurrencia que está justamente antes de
la ecuación indicial (5) con la intención de usar e = e2 = 0, es evi-
dente que las a serán diferentes de aquellas con e = el. Por tanto, es
deseable cambiar de notación. Emplearemos las b en lugar de las a. Con
http://carlos2524.jimdo.com/
b o arb.
( -2) bo
1.1
( -1) b1
2·3
_ ( O) b2
3·5
(n - 3)bn-l
n ( 2n - 1) •
EJERCICIOS
2. 4xy" + 3y' - 3y = O.
( -3 )"xn+t
X· + "=ln!5.9.13 .
ex)
SOL. Y1 = ~
... (4n + 1) '
ex) ( -3 )"x"
Y2 = 1 + ~ n! 3 .
1 7 . 11 . . . (4n - 1)'
3. 2X2(X + l ) y" + x (7x - +Y=
l)y' O.
ex)
Y2 = x; +! ~ (- 1) " (n + 1) (n + 2) xn+~ .
4. 2xy" + 5(1 + 2x )y' + 5y = O.
ex) 3( -5 )"x"
SOL. Y1 = 1 + ~1 n! (2n + 1) (2n + 3) ;
Y2 = x-~ - 10x-;.
5. 8x 2y" + 10xy' - (1 + x)y = O.
x"+.
ex)
y - x-; +~ .
2- "=12"n!I.5.9 ... (4n -3 )
6. 3xy" + (2 - x )y' - 2y = O.
(3n + 4 )Xn+i ex)
SOL. Y1 = ' ~
,,=0 4 . 3"n.I ;
(n + l)x" ex)
y -1+~ .
2 - ,,=12 . 5 . 8 .. . (3n - 1)
7. 2x(x + 3)y"- 3 (x + l)y' + 2y = O.
• ex)-1) "+1xn+~ . (
y2 = 1 +~x +iX2.
8. 2x2y" - x(2x + l )y' + (1 - 5x)y = O.
1 (2n + 3) (2n + 5 )Xn+1 . ex)
SOL. Y1 = X + 15 ~ I '
"=1 n.
2"-1 (n + 1) (n + 2) xn+!
Y2 = xi + ~ .
ex)
"=1 1· 3 . 5 . . . (2n - 1)
9. 9x 2y" + 3x(x + 3)y' - (1 + 4x)y = O.
SOL. h = xs + ~ X: ;
_ _11 10 ( - 1 ) "x"-~ . ex)
Y2 - X + ~13nn!(3n _ 5) (3n - 2)
10. 2xy" + (1 - 2X2) y' - 4xy = O.
http://carlos2524.jimdo.com/
Y2 = 1 + ~
<Xl
.
( -1) n2 n (n 1) X" +
n=1 1. 3·5··· (2n - 1)
14. 2xy" + (1 + 2x) y' - 5y = O.
SOL. Y1 = Xl; + ~X! + l\X!;
<Xl 15 ( -1) n+1 x"
Y2 = n~on!(2n - 5) (2n - 3) (2n - 1)'
15. 2x 2y" - 3x ( 1 - x) y' + 2y = O.
SOL. Y1 = X
2
+ ~
<Xl ( - 1) n3 n (n 1) X"+2 •+ ,
n =15 • 7 ·9 .. . (2n + 3)
<Xl ( _1 )n+l3 n (2n _ 1) xn+.
Y2 = x' + ~1 2nn! .
16. 2x 2y" + x(4x - 1) y' + 2 (3x - 1)y = O.
<lO (_ 1) n2nxn+2
SOL. Y1 = ~
n=O n.
I = x e-
2 2
",;
_ '" ( -1 ) n4nxn-!!
Y2 = X ,+ ~ d _ 3) ( - 1) . 1 ... (2n - 5)"
17. 2xy" - (1 + 2X2 )y' - xy = O.
<Xl 2kx 2 k",.~
SOu.
h = xl + k;l 7 . 11 . 15 . . . (4k + 3) ;
Y2 =1+ ::1
<lO x2k
2kk! = exp (iX2) .
18. La ecuación del ejercicio 17 tiene una solución enteramente particular
Y2 = exp (iX2) obtenida por el método de series. Hacer un cambio
de la variable dependiente en la ecuación diferencial, usando y = v
exp (iX2) (el artificio de la sección 84), y obtener entonces la so-
lución general en "forma cerrada".
SOL. Y = Cl exp (ix ) + C2 exp (ix )
2 2 f3!! exp ( - if32) df3. J:
En los ejercicios 19-22, usar el método de serie de potencias. ¿ Cuál
es la causa de que las relaciones de recurrencia degeneren en rela-
ciones de un término?
http://carlos2524.jimdo.com/
356 Soluciones cerca de puntos regulares singulares [Cap. 21
EJEMPLO b) : Si
e+n
y = e(c + l)(e + 2) ... (e +n - 1)
http://carlos2524.jimdo.com/
358 Soluciones cerca de puntos regulares singulares [Cap. 21
entonces
~~ = y t! n - ~- e ! 1- e ! 2- .
_1_¡}.
e+n-
EJEMPLO e): Si
(2)
L (y) = n=O
~ (n + e) (n + e -1 )anXn+c + ~ 3(n
n=O
+ c)anX 1HC
~ ~
+ ~ anXn+c - ~ 2anXn+c+1,
,,=0 .. =0
de lo que
~ ~
a - ------::-.
2an-l
n - (e + n + 1) 2'
2"ao
n >
-
1: a" = - - - -------..:--- - -- -
[( C + 2 )( c + 3 ) ... ( c + n + 1) ]2'
(6 ) y ( x, c) = XC + ¿"" a" ( c ) x M C
,
n= l
en la que
(7 ) 2"
n ~ 1: a,, ( c) = [( c + 2 )( c + 3 ) . .. (c + n + 1 ) ]2'
La y de la ecuación (6) ha sido determinada en tal forma que la y del
miembro derecho de la ecuación (4) deberá reducirse a un solo término,
el término n = O. Esto es, para la y (x, c) de la ecuación (6), tenemos
(8)
Una solución de la ecuación diferencial original es una función y
para la cual L (y) = O. Ahora vemos porqué la elección c -1 nos =
conduce a una solución i esta elección hace que el miembro derecho de la
ecuación (8) se anule.
Una consecuencia automática de la igualdad de las raíces de la ecua-
ción indicial es que el factor (c + 1) se encuentre al cuadrado en la
ecuación (8). Conocemos del cálculo elemental que si una función
contiene una potencia de un cierto factor, dependiente de c, entonces la
derivada con respecto a c de esa función contiene el mismo factor a una
potencia disminuida en una unidad con respecto a la original.
Para la ecuación (8) en particular, la diferenciación de cada miembro
con respecto a c nos da
Por tanto, calculamos a,, (c) y a',, (c) en c = -1. Sabemos que
2"
a,, (c ) = [(c +2) (c+3) ... (c+n+ 1))2
de donde por el método de la sección 108 obtenemos inmediatamente
H" ~ !..
= 1 + !2 + !3 + ... + !n = k=lk
Podemos ahora escribir a',,( -1) en forma más simple como
, 2"+1H,,
a,,( - I) = - (n!)2'
Y
(14)
SOL. Y1 = ~
(k')2 ;
"=0 .
00
2
( -1 ) "H"X k--1
~
Y2 = Y1 Inx - k=1 (k'• )2 .
5. x(l + x)y" + (1 +5x )y' + 3y = O.
00
2 _ 00 ( _1 )"H"x 2k
Y - y¡ln x - k-;l 22"( k!) 2 •
7. x 2y" - x(l + 3x )y' + (1 - 6x) y = O.
00 3"(n + 1) (n + 2 )x"+l
SOL. Y1 = ~ .
n=O 2 . n! '
00
3n (n + 1) (n + 2) (Hn+2 - 2Hn - 1 )X7H1
Y2 = Y1 lnx + ~ 2· n! .
n=l
SOL. Y1 =1- x; X2 =h In x +f •
t" -
00
~ 2n
(n
(
+ 1 )x" •
11=2 n n - 1)
10. Resuélvase la ecuación del ejercicio 9, alrededor del punto x = 2.
SOL . Y1 = 1 + (x - 2) ;
Y2 = Y1 In (x _ 2) _ ~ (x _ 2) _ i; (- 1 ) n(n + 1) (x - 2) n
11=2 2n n(n - 1) •
http://carlos2524.jimdo.com/
364 Soluciones cerca ae puntos regulares singulares [Cap. 21
Y2 = Yl In x + ~
'"
,,=1
( - 1 ) "[ ( - 1) . 1 . 3 . . . (2n - 3)]( 2H 2,,-2 - H "-1 - 2H" - 2) x"-!
(n !) 2 •
81 2( -3 ) "X"-l
'"
¿ '( -1 )( n -2)'
Y2 = Y1Inx -1 5--¡-x+ ,,=3n.nn
17. xy"+ (l - x )y' -y=O.
SOL. Y1 1+ ¿ =
'" x"
11=1 n.
eX; ,=
Y2 = Y1 lnx _ ¿'" H x"
n~l n.
_n_, .
18. Con respecto al ejercicio 17 anterior. Se encontró que una solución
es Y1 = eX. Úsese el cambio de variable dependiente, y = ve", para
obtener la solución general de la ecuación diferencial en la forma
http://carlos2524.jimdo.com/
365
(1 ) xy" (4 + x)y' + 2y = O.
(2) Y = ,¿
110=0
n=O
00
- ,¿ (n +c- 2) anXn+c
n=O
o
00 00
/
http://carlos2524.jimdo.com/
366 Soluciones cerca de puntos regulares singulares [Cap. 21
\
http://carlos2524.jimdo.com/
§110] Casa no logarítmico 367
(n - 3 ) a"-l
n (n-5) ,
de lo que
_ 3· 4 . 5 ... ( n - 3) as _ 3 . 4 . 5as
a" - - - --::--,---:------,-:--:------=-...,....-:
[6·7·8 . . . n]( n-5)! (n-2)(n-1)n (n-5)!'
Por tanto, con ao y as arbitrarias, la solución general puede escribirse
como
-
y - ao
( 1 + 1X
2
+ i21 X 2) + as [ x 5 ,,"~":s(n 60x" ]
_ 5) !n (n - 1)(n - 2) .
El coeficiente de as puede escribirse también con un corrimiento de Ín-
dice en la forma mostrada en seguida
ex) 60x"+s
n~o n! ( n + 5)(n + 4) (n + 3)'
Antes de proceder a hacer ejercicios, examinemos una ecuación en la
que no se encuentra la situación afortunada de que ao y as son arbitrarias.
Para la ecuación
(4) x2y" + X (1 - x) y' - ( 1 + 3x) y = O
el tanteo
nos conduce a
ex) ex)
n = 1: -a1 - 2ao = O,
n = 2: O . a2 - 3a1 = O,
_ ( n + 1 )an-1
n ~ 3: an - .
n(n - 2)
Se sigue que
a1 = -2ao,
O . a2 = 3a1 = - 6ao.
Estas relaciones no pueden satisfacerse excepto eligiendo ao = O. Pero,
si esto se hace, a2 será la única constante arbitraria y la única solución
que se obtendrá del trabajo será la correspondiente al valor mayor de
e, e = 1. Este es un caso donde se indica una solución logarítmica.
La ecuación será resuelta en la sección siguiente.
Una buena manera de desperdiciar tiempo es usar ao = O, a1 = O
Y a2 arbitraria para determinar an, n ~ 3, de la relación de recurrencia
anterior. En esta forma se puede hacer trabajo extra para obtener una
solución, la cual se reobtendrá automáticamente al resolver la ecuación
por el método de la sección siguiente.
EJERCICIOS
4. xy" - (3 + x) y' + 2y = O.
~ 24 (n - 3) x1I
y = ao( 1 + + ¡¡X +
2 1 2)
SOL. :¡X a4 ¿. ,
11=4 n.
5. x(1 + x)y" + (x + 5)y' -
4y = O.
SOL. y = ao(x- 4 + 4x- 3 + 5X-2) + a4 ( 1 + ix + ~X2).
6. Resuélvase la ecuación del ejercicio 5 alrededor del punto x = -1.
SOL. Y = ao[1 + (x + 1) + i(x + 1) 2] +
a 00
7. x2y" + x 2y' - 2y = O.
SOL. y = ao( x-1 - t )
8. x( 1 - x )y" - 3y' + 2y = O.
00
SOL. Y = ao(1 +!x + -kx2) + a. ~
n=4
(n - 3)xn.
9. Resuélvase la ecuación del ejercicio 8 alrededor del punto x = 1.
SOL. Y = ao[ (x - 1)-2 + 4(x - 1 )-1] + ~[1 + Hx - 1 ) + Hx - 1)2].
10. xy" + (4 + 3x )y' + 3y = O.
00 (_ 3 ) n-3 x n-3
SOL. y = ao(x-3 - 3x- 2 + ~X-1) + 6a 3 ~ I
n=3 n.
11. xy" - 2(x + 2 )y' + 4y = O.
00 60 . 2n - sx"
SOL. y = ao(1 + x + -kx2 + a s ,,~s (n _ 5) !n (n _ 1) (n _ 2) •
12. xy" + (3 + 2x )y' + 4y = O.
SOL.
De la forma supuesta
n=O
(e + 2) [e ( e + 1) . . . (e + n - 1) ]
y deberá seguirse que
( 3) L (y ) = ( e + 1) (e - 1) a ox c •
,
De la raíz grande, e = 1, sólo se puede obtener una solución. De la
raíz más pequeña, e = - 1, se pueden obtener dos soluciones siguiendo
la técnica de usar y (x) e) y ay ( x) e) / oe como en el caso de raíces iguales
si el miembro derecho de (3) contuviera el factor ( e + 1) 2 en lugar de
(e + 1) a la primera potencia. Sin embargo, a o todavía es arbitrario
de tal forma que podemos tomar
ao = (e + 1)
para obtener el factor cuadrado deseado en el miembro derecho de la
ecuación (3) .
http://carlos2524.jimdo.com/
(7) y2 -_ (oy(x,e)) .
oc C=-l
oy ( x, e)
(9) -,-,:,--"":'" = y ( x, e ) In x
oc
+ XC
+ (e + 1) (e + 3) Xl+ {_1 _+ _1_ __1__ ~}
C
(e + 2 )e e + 1 e +3 e +2 e
( e + 4) X 2+C { 1 1 1}
+ (e + 2)e c + 4 - c + 2 - ~
00
+~
n=3
EJERCICIOS
Encuéntrense las dos soluciones linealmente independientes válidas alre-
dedor de x = 0, a menos de que se indiqúe otra cosa.
1. xy" + y = o.
co ( -1 ) 1IX"
Y1=};
SOL.
11=1 n. ( n - 1)1;
' .
11
Y2 = Y1 in x + 1+ x _ ~ ( -1 )",( H n + H"(1 )X •
11=2 n.(n - 1) .
2. x 2y" - 3xy' + (3 + 4x)y = O.
co(-1) 1I+l411 xn+1
SOL. y¡ =};
n=2 n.' ( n - 2) ., ;
Y2 = Y1ln x + Y1 + X + 4X2 + ~ (-4) n(Hn + H,,-2) Xn+1.
11 =2 n!(n - 2)!
3. 2xy" + 6y' + Y = O.
co ( -1 ) 11+lxn-2
Y1 =
SOL.
"';2
211 n! (n - 2) ! ;
co (-1)n(H + H -'2)X 2 1l
1I n
Y2 = y¡ In x + Y1 + x- 2 + ix-1 + };
-
Y2 = y¡ In x
co
+ 2x-?¡ - 11;2 211-1t~1.
H
_
X 1l- '
~ ) !.
5. x 2y" - x( 6 + x) y' + lOy = O.
(n + 1) xn+2
SOL. Y1 = };
11=a
co
2 (n - 3)! '
.
Y2 = Y In x + '!'y + X2 _ x3 + .:!x4 + ~ [1 - (n + 1 ) Hn_3]Xn+~
1 < 1 2 11=3 2 (n - 3) ! . '
http://carlos2524.jimdo.com/
374 Soluciones cerca de puntos regulares singulares [Cap. 21
(-1)k(H k +Hk_1)X"2k-1
Y2_y
- 1
l nx +x -1
-
k=l
¿OC>
22kk!(k-1)! .
7. xy" + (3 + 2x)y' + 8y = O.
00 (_ 1) "+12" (n + 1) xn-2
SOL.
Y1 = ~2 (n - 2) ! ;
~ ( -2 )"[(n + 1)H"-2 - 1]X"-2
Y2 = Y1 In x + x- 2 + 4x- + n=2
..:..
1
(n - 2) ! .
8. x(1 - x)y" + 2(1 - x)y' + 2y = O.
YI = -2 + 2x;
SOL.
2xn-l
00
Y2 = YIln x + x- + 1
1 - 5x + ~3 (n _ 1) (n _ 2) .
112. SUMARIO
Confrontados por una ecuación lineal
(1 ) L (y) = O,
determinamos primero la localización y naturaleza de los puntos sin-
?;ulares de la ecuación. En la práctica, el uso para el que se van a
destinar los resultados, dictará cerca de cuál punto o puntos deberá
resolverse la ecuación. Al buscar soluciones válidas alrededor del punto
x = Xo, se trasladará el problema siempre al origen, haciendo la trans-
formación x - Xo = v.
Las soluciones válidas cerca de un punto ordinario x = O de la ecua-
ción ( 1) tomarán la forma
(2)
(3)
n=O
(4)
para la cual
(5 )
será determinada, donde k es cero o uno para las ecuaciones que están
siendo tratadas aquí y el y e2 son las raíces de la ecuación indicial.
Estas raíces pueden determinarse entonces de los coeficientes reales en
( 4) usando ya sea e = C1 o e = C2, en ambos casos se obtendrán dos
soluciones de la ecuación diferencial. Si C1 = C2 los resultados deberán
ser idénticos y es indicado el uso de 3y (x, c) loc. El otro caso logarí!mi-
ca será identificado por el hecho de que alguno o algunos de los coefi-
cientes fn( c) no existirán cuando C = C2 sea la raíz pequeña. Se necesita
entonces nuevamente el proceso de diferenciación después de la intro-
ducción de ao = C - C2.
El método aquÍ esbozado tiene la desventaja de que parece tender
al uso automático de reglas, lo cual es siempre un procedimiento peli-
groso en matemáticas. Sin embargo, cuando el método ha sido comple-
tamente entendido en cada uno de los cuatro casos posibles, puede
aplicarse con seguridad y facilita mucho los cálculos.
La extensión de los métodos de este y el capítulo precedente, es
directo a ecuaciones lineales de orden más alto. Como un ejemplo, será
tratada inmediatamente una ecuación de cuarto orden cuya ecuación
indicial tiene raíces c = 2, 2, 2, l Deberá determinarse una sene
para y(x, c),
y(x, c) = ao [xc + ~1 fn(c)xn+c]
para la cual el miembro izquierdo de la ecuación se reduzca a un tér-
mino, tal como
L[y(x, c)] = (c - 2) 3(2c - l)aox c.
Así, podrían obtenerse cuatro soluciones linealmente independientes
y
-
Y1-y(x,2),ya- a C)] ,y3 = [02 ()( x,2 C)] ,Y4=y(x,i).
- [Oy(X,
e C=2 c C= 2
http://carlos2524.jimdo.com/
377
EJERCICIOS VARIOS
CAPÍTULO 22
Ecuaciones de tipo hipergeornétrico
/
http://carlos2524.jimdo.com/
380
(3) r (x + 1) = xr ( x ) .
Empleando repetidamente la relación ( 3), encontramos que si n es
un entero,
r (a + n) ( a+n-1)r ( a+n - 1)
(a +n - 1) ( a +n - 2 )r ( a +n - 2)
(a +n - 1) ( a +n - 2) ... ( a )r(a)
(a)"r(a).
Por tanto la función factorial y la gamma están relacionadas por
(4) a - r (a
n -
+ n) ,
r(a)
n entero, n > o.
Realmente (4) ,aunque esto no se prueba aquí, es válida para cual-
quier número a complejo excepto cero o un entero negativo.
http://carlos2524.jimdo.com/
381
n=O
(a ) ( b ) x n
F (a, b; e; x ) = 1 + 00
~ n ~,
n=l (c )nn.
y = ~ fnxn+1-c
n=O
y = ~ ad'
k =o
k
(5)
- L ( )1 + ~ ( -n) k(H.....k - Hn - 2Hk)X
y2 - n X n X;:l (k!) 2
00 (_ 1) nn ! (k - 1) !x"+n
+ k~l [ ( k + n) 1]2 •
La solución (3 ) es válida para toda x finita; la solución (5) es válida
para 0< Ixl < oo.
1
J--n( X) = 2-nr ( 1 _ n) y2.
y = ~ aix;+c,
;=0
y procediendo con la técnica de la sección 111 , determinemos y (x, e) y
/
http://carlos2524.jimdo.com/
386 Ecuaciones de tipo hipergeométrico [Cap. 22
Para nJ un entero mayor que uno, las ecuaciones (5) y (6) pueden
usarse como el par fundamental de soluciones linealmente independien-
tes de la ecuación de Bessel (1).
y= ¿ ajx
j
;=0
y empleando los métodos del capítulo 20, obtenemos la solución general
_ [ 00 2" ( - n ) ( - n + 2) . . . ( - n + 2k - 2) X 2k]
(2) Y - ao 1 + k~l (2k) ! .
00 2k (1 - n) (1 - n + 2) . . . (1 - n + 2k - 2) X 2"+1]
+ aJ. [ x + ~l (2k + 1)! '
válida para toda x finita y con ao y al arbitrarios.
El interés de la ecuación (1) es mucho mayor cuando n es un entero
positivo o cero. Si n es un entero ,par, el coeficiente de ao en (2) tiene·
un término final, cada término para k ~ H n + 2) será cero. Si n es un
entero impar, el coeficientede al en (2) tiene un término final, cada
término para k ~ H n + 1) será cero. Por tanto, la solución de la
ecuación de Hermite siempre es polinomial, de grado nJ para n cero o
entero positivo. Es elemental pero tedioso obtener de (2) una expresión
única para esta solución polinomial. El resultado es
[!nJ ( _1 ) k n ! ( 2x ) n-21<
(3) Hn(x) = ~o k!(n - 2k)! '
en la que [in] simboliza la parte entera de in.
http://carlos2524.jimdo.com/
§120] Pol inomios de legendre 387
(5) y1 = 1 + ~ (- n ) k (n + 1) k (1 - X)
' k
k=1 ( 1) kk! 2
El miembro derecho de la ecuación (5) es un ejemplo de la función
hipergeométrica que tratamos en la sección 115. De hecho
/ http://carlos2524.jimdo.com/
388 Ecuaciones de tipo hipergeométrico [Cap. 22
(6) 1-
Y1=F ( -nJn+ 1;1;-2- . x)
Si n es cero o entero positivo, las series en (4), (5) o (6) tienen
un término final. En este caso estas series se llaman polinomios de
Legendre y se simbolizan por P,, (x) . Escribimos nuestra solución no
logarítmica de la ecuación de Legendre como
00 (a) n( b ) nX n
(3) F(a,b;c;x)=+2:
n=l (C)
nn.'
de la sección 115. En (2) la serie tiene solamente un parámetro en el
numerador, a, y uno en el denominador, c. En (3) hay dos parámetros
en el numerador, a y b, yun parámetro en el denominador, c. Se acos-
tumbra entonces usar para el miembro derecho de (2) una notación
como la de (3).
Escribimos
•• _ ro (a)nXn
(4) lFl(a, c, x) - 1 + n~l (C)nn!'
CAPÍTULO 23
Tópicos adicionales sobre
soluciones en
series de potencias
(2)
en la que ao = 1, al = 0, a2 = 0, y para
(n - 1) a"-l - an -3
n ~ 3: a" = ----~~--------
n (n - l) ,
en donde )\2 está dada por
00
(4) y2 = ~ b"x",
"=0
http://carlos2524.jimdo.com/
n Z 4:
00
Y2 = ~ bnx") en la que bo = O, b, = 1, b2 = O, ba =- ~,
1
n Z 4: b" = - n(n _ 1) (bn_2 + bn-a + b"-4)'
2. Demuéstrese que las series en las respuestas al ejercicio 1 principian
como sigue:
y' = 1 - tx2 - txa - + """ X 5 + /27 x 6 + 2!OX7 + . .
2\X4 Q
y
Y2 = X
'l
- '6 Xa - I 4- 2'4
TIx I X 5+ l20 1 x6+ 'i'44
1 x 7+ . . '.
3. (1 - x ) y" + (2 + x) y' - 2y = O.
ro
SOL. Y1 = ~ anx") en la que ao
710 = 0
= 1, al = O,
(n - 4)
n ~ 2: a" = n (n _ 1) [( n - 1) an-1 - an-2];
00
Y2 = ~ b"x") en la que bo
710=0
= O, b = 1,
1
n >_ 2: (n -_ 4)
b" -- n(n 1) [( n - 1). bn-1 - bn-.2
]·
http://carlos2524.jimdo.com/
~ ~
n=2
Por tanto la ecuación indicial es (e + 1) 2 = O.
Ya que las raíces de la ecuación indicial son iguales, e = -1, -1,
se sigue que existen las soluciones
(4)
~
(5 ) y2 = y, In x + ~ bnxn-\
n=l
(6)
'" anXn-\
y, = ~
n:::'Ü
en la que ao = 1, a, = O
(n - 1) an-l + an-2
n2
En seguida deseamos pedir que L(y2) = O. De
y
ca
y'; = y~ In x + 2x-Iy: - x-2y I + ~ (n - 1 ) (n - 2 ) b'11x'11-3.
El cálculo directo de L (y2 ) nos da
ca
L (Y2) = L (YI) In x + 2xy; - YI + XYI + 3YI + ~ n 2bnx'n-l
ca
'11=1 ca
+ ~ (n - 1) bn-1X'n-1 + ~ bn _2X'11 -1•
ca ca
- ~ 2nanx'11 -1 - ~ an _1X'n-1,
n=l n=l
( 1)
n='Ü
en la que ao = 1, al = O,
( 2)
( n - 1) an-l + an-2
n2
Sea En (x) el error hecho al calcular yl (x) usando los términos del
miembro derecho de ( 1) solamente e incluyendo el término anX,>-l.
Entonces En(x) es igual a la suma de los términos que quedan
00
(4 ) IEn (X ) 1 S ~ lasx8-ll·
8=n+l
(5)
an 1 [2"-1 (n - 1) 2n -2 ]
l l<n 2 (n - l )! + (n -2)!
< -2n-2
n 2
[2
(n -2)!
1 ]
+ (n-2)! ...,--------,::-c-:-
3 . 2n-2
< -----::--c:----=:--:--:-
2
n (n-2)!
< 3 (n - 1). 2"
4n n!
Por tanto
00 28+7H1 Ixl8+" _ 2n+1lxl" 00 281xl8
IE,, ( x) 1 < 8~O (n + 1) !s! - (n + 1) ! 8~ ~'
y, finalmente
( 10)
con un error menor que 0.0017. Concluimos que, para dos decimales,
Y1 (i) = 3.94.
Volvamos ahora a Y2(X). Se encontró en la sección 123 que esta
serie está dada por
00
en la que bl = -1 , b 2 = i,
(13) n > 3: b n = _ (n - 1) bn-l + bn--2 _ 2a n _ a"-l.
- n2 n n2
Supongamos
(14)
(15)
oc
Ahora, sea E" * (x) el error cometido al calcular ~ b"x"-l usando los
"=1
términos superiores a, y que incluyen a brlx n - 1 • Entonces
E,,*(x)
y por ( 15)", para n ~ 1,
IE,,* (x) I <
así
2
(16) IEn*(x)1 < 2n+2Ixlne
(n+1)!
/"/
.
La solución y2 (x) principiará corno sigue
)1,2 ( x) = y1 ( x) In x - 1 +~- 1 ~ SX2 - 1 ~ ! 2X3 +
Calculemos y2 ( i) con dos cifras decimales. Vernos primero que
http://carlos2524.jimdo.com/
402 Tópicos adicionales en series de potencias [Cap . 23
Entonces
EJERCICIOS
2a"
n ~ 4: b" =
n
2. 2x(1 - x ) y" + (1 - 2x)y' + (2 + x)y = O.
co
SOL. YI = n=O
~ anx1l+i, en la que' ao = 1, al = -i,
n ~ 2:
(2n + 1) (2n - 3) an - l - 2a n-2
an = 2n (2n + 1) ;
co
Y2 = ~ bnx n
, en la que bo = 1, bl = -2,
b _ 2n (n - 2) bn - l - bn _2.
n ~ 2:
"- n(2n - 1)
3. xy" + y' + x (1 + x) y = O.
co
SOL. YI =~ n ==O
anx", en la que ao = 1, a¡ = 0, a2 = -:1,
n ~ 3:
co
Y2 = yllnx + ~ bnx", en la que bl = 0, b2 =f, b3 =/7'
"=1
n ~ 4: bn =
4. xy' + x(1 + x)y' - (1 - 3x + 6X2)y = O.
co
SOL. Y = ~ anx"-\ en la que ao es arbitraria, al =
2ao, a2 es
n=O (n+l)a n _1-6an_2
arbitraria, n ~ 3: aro = - -'----=--'----'.:....::...--..:.:....:.
n(n - 2)
5. Demuéstrese que la serie, en la respuesta al ejercicio 4, principia
como sigue
http://carlos2524.jimdo.com/
n==O
Y2 '" b"x".
= ylln X + ~
El coeficiente bI (s = 1) vuelve a ser arbitrario y lo elegimos igual
a cero. Demuéstrese que las YI y Y2 indicadas son soluciones si
ao = 1,
n > 5: a
n
= _ (n + 1) a n - + an -
I 5
n(n+1 ) ,
y si las b están dadas por
bo = -1 , bI = O (así elegidas ) , b 2 = 1, b3 =- ¡, b. =;~,
b nb n- I + bn - 5 (2n - 1) a,,-t +
n > 5: n
__
- n (n - 1)
_
n(n - 1)
a " -2
.
7. Para la YI del ejercicio 6, probar que la a" alterna su signo. Calcular
también los términos de YI hasta el término x 7 •
SOL.= X - X2 + ~X3 - kx· + 2'.X5
YI - /.X6 + 1 ~ ~ sx 7 + .. .. .
8. x(x - 2 )2y" - 2(x - 2)y' + 2y = O.
SOL . JI = 1 - i x;\
x
Y2 = YI In x + ~ b"x", en la que bl =~, b z
71=1
= - ~, b3 =- .11!>
n -2
n>4 : b" = ~[2(2n - 1) b"-1 - (n - 3 ) b"-1!J.
n> 2:
(2n + 1) an - l + 2a n -2
a" = 2n (2n + 1) ;
Yz = '"
~ b"x" en la que bo = 1, bI = 1,
n=O
.... 2
n,;:;:-: b,, __ _n_bn-:-=-l_+'----b-"''t'-z
--=
n(2n - 1)
http://carlos2524.jimdo.com/
404 Tópicos adicionales en series d e potencias [Cap. 23
11. Demuéstrese que las respuestas al ejercicio 10 están dadas también por
Y2 = eX, Yl = eX r
x
exp ( - YJ2) d(3.
(5)
Y
2
d y = w 4 d y + 2w 3 dy
2
(6)
dX2 dw 2 dw
Con la ayuda de (5) y (6) la ecuación (4) se transforma en
y= ~ anW'HC,
n::::O
n+c-1'
http://carlos2524.jimdo.com/
se sigue que
( -1 )nao
.
an = ------~~--~---------
c( c + 1) (c +n - 1) .
Por tanto elegimos
1 ) nw n +c .
y( w, c) = Wc + 00 ( _
~ ---,..:----'.--~-
n"l c( C + 1) . : . (c +n - 1)'
y encontramos entonces que
ay(~: c) = y (w, c) In w
-f
( -1 ) nw n +c {.!c + _1
c+1
_ + ... + 1 }
c+n-1.
n::l c( C + 1) ... (c +n - 1)
Empleando la raíz c = 1, llegamos a las soluciones
00 (_ 1 ) nw n +1
yl = W + ~
n::l n.
I
(8)
Y
(9)
EJERCICtos
En los ejercicios 1-6, los puntos singulares en el plano finito han sido
clasificados y localizados. Para cada ecuación determinar si el punto en el
infinito es un punto ordinario, un punto regular singular o un punto
irregular singular. No resolver las ecuaciones.
1. x 3 (x - l )y"+ (x - l )y'+4xy=O. (Ejercicio 1 de la sección 104.)
SOL . Punto regular singular.
2. 2 2
X (X - 4 )y" + 2x 3y' + 3y = O. (Ejercicio 2 de la sección 104. )
SOL. Punto ordinario.
http://carlos2524.jimdo.com/
§125] Solución poro x grandes 407
00 00
= ~ (n
11=0
+ e + 1) 2anX1HC - ~ (n + e) anX1HC-1
n=O
00 00
= ~ (n + e) 2an_1X'HC-1 - ~ (n + e) an x 1HC-1.
11=1 n=O
http://carlos2524.jimdo.com/
§126] Ecuación con un punto irregular singular 409
(3)
(4) yl = 1 +¿ n !xn.
n=l
(5) yl = f
n =O
Joo e-ttn dt xll.
o
o
yl = r e- t n~o (xt)n dt,
00 e-t
(6) yl = Jo -1-- dxt t .
Probaremos ahora que yl de (6) es una solución de la ecuación dife-
rencial (1) sobre el rango semi-infinito x ::;; O. Si sucede que x es com-
pleja, se puede usar la misma clase de discusión para 'demostrar que
yl es una solución en la región en la cual la parte real de x es no
positiva.
En (6) la variable de integración nunca es negativa, t ¿ O. Para
x ::;; 0, tenemos por tanto
-xt ¿ °
tal que (1 - xt) ¿ 1, de lo cual concluimos que
(7) 0<-1_1_<l.
- xt-
Debido a (7) la integral de (6) converge:
9 , Joo te- t dt
( ) Y = o (1 - xt )2'
http://carlos2524.jimdo.com/
§126] Ecuación con un punto irregular singular 411
XYI
I
+ JII -
-
Jooo e-t [ 1 _1 xt + xt J
d .
(1 _ xt) 2 t,
esto es,
(12)
De (12)
(13)
x dt /
( 1 - xt) 2
1 - xt
Obtenemos por tanto
X ?-YII + XYI = [
- e-
--
t
J"" + J"" -e- t--
dt ,
1 - xt o o 1 - xt
o
x2y~ + Xyl = -1 + yI
Entonces la yI de (6) satisface la ecuación
(14) x2y; + (x - l)Yl = -1.
http://carlos2524.jimdo.com/
412 Tópicos adicionales en series de potencias [Cap. 23
(17) 1 + ~ n!x
n::::l
n
•
( 1)
(7) x ::; o.
Hasta ahora hemos demostrado que O < (1 - xt) -1 ::; I para x en el
rango x ::; O. Por tanto, de (7), podemos concluir que
(1 )
\,
t
en las cuales n, gn, y ao se conocen, debiéndose encontrar an cuando
n > O.
Un caso especial de (1) es fácil de resolver. Supóngase que conocemos
b o y tn, y que
b1 = bo +t 1,
b 2 = b1 +t 2,
bn = bn-l + tn,
de tal forma que sumando los miembros respectivos obtenemos
..
bro = bo + ~ t k.
k=1
donde la 'f" se elegirá de manera tal que se elimine In. Usando (3) en
la relación (1) llegamos a
Ahora deseamos elegir 'fro tal que 'f.. = lro 'fn~1. Siempre podemos hacer
esto utilizando métodos que ya hemos discutido y usado y que serán
exhibidos en los ejemplos que aparecerán a continuación. Una vez que
'f.. se ha elegido, la relación (4) es del tipo (2) y se puede resolver
fácilmente para bn. Entonces aro está dada por la ecuación (3).
Hagamos a.. = 'fnb .. con 'fo = 1 de tal forma que ao = bo. Entonces
de (5) se obtiene que
n ;;::: 1: 'f"b" = n 2 'fn-lbn-1 + (n + 1),
debemos ahora escoger 'f.. de tal forma que satisfaga
(6)
http://carlos2524.jimdo.com/
416 Tópicos ad ici onales en series de potencias [Cap. 23
Ahora
'P1 = 1 2'P(J,
'P2 = 22'P1,
y, finalmente
12 ( -1 )"ao ~ ( -1 ) -k ( k+2 ) ! ( k+l ) !
n >
- 1: an = ( n + 3) ! (n + 2 ) ! + k~l (n + 3) ! (n + 2 ) ! .
Con un poco de práctica la 'P" puede obtenerse por inspección.
EJERCICIOS
1- x
Ahora nuestra ecuación diferencial podría representarse como
<Xl
(5 ) xy" +y = ~ X1H'2,
n=O
(6)
n=O n=l
Dos series de potencias que son idénticas en una región tienen sus
coeficientes correspondientes iguales, por tanto concluimos que las a"
están dadas por las relaciones
n = O: 6ao = 1,
n ~ 1: (n + 3 ) ( n + 2) a" + an-1 = 1.
Es bastante razonable parar en esta etapa de nuestro desarrollo y es-
cribir la solución general de la ecuación diferencial original como
y = AY1 + By2 + y.p,
en la cual la y1 y y2 están dadas por las ecuaciones (3) Y (4) Y
00
(7 ) Yv = ¿ a nXn+a,
n:=Q
donde ao -t y
1 - an - l
an = .
(n+3 )( n+2 )
Se pueden calcular tantas a como se deseen de la relación dé recu-
rrencia anterior.
Puede encontrarse por el método de la sección precedente una forma
explícita para la ano En efecto, de acuerdo con el ejemplo b) de esa
sección
n~1: 12 ( -1 )nao ~ ( - 1 )n+k( k+2 ) ! ( k+1)!
a" = ( n + 3) !( n + 2) ! + k=l (n + 3) !( n + 2) ! .
Por tanto, ya que ao =t , podemos escribir la solución particular in-
dicada en ( 7) como
( 8 ) Yv = ~X3 + ~
00
1 + ::1 (-
[
2
" ] ( - 1) "Xn+3
1) k ( k + 2) ! ( k + 1)! ( n + 3) ! ( n + 2) !.
/
http://carlos2524.jimdo.com/
420 Tópicos adicionales e n series de potencias [Cap. 23
ex)
xy; + yP = ~ X'H2.
n=O
EJEMPLO b): Resolver la ecuación
(10) 2x2y" + x (4x - l)y' + 2(3 x - l )y =:r + x 4
cerca de x = O.
La ecuación homogénea correspondiente tiene las soluciones lineal-
mente independientes
( -2) "X"+2 00
(12) y.2 - X - ¿ ,
(- 3)( - 1) ·1·· . (2n - 5)
como se indicó en el ejercicio 16 de la sección 107.
Si usamos
ex)
y = ¿ a"xn +y,
n=O
o
_ 1 3 1 1
00
~
( _ 1 ) n-12 n xn+3
- ' - _ - - ' -_ __
y - 'i
X
- 12"67>=1 (n + 1) ! .
o
(14)
lo cual se puede verificar fácilmente por sustitución directa.
Si en este ejemplo elegimos y = 2, una elección menos aparente
quizá, entonces las a deben determinarse de las relaciones
n = O: O· ao = O, (ao arbitraria),
n = 1: 7al + 14ao = 1,
n=2: 18a2 + 18a1= 1,
n ~ 3: na,. + 2an-1 = O.
Así forzamos a a2 a ser cero, an = O para n > 2. Si a2 = O,
al = 1\'
ao = 1\ - ~a1 = 2\\ .
De aquí que las elecciones de y = 2 Y ao nos conducen a an = 2\\ =O
para n > 2 dándonos otra vez la solución (14).
EJERCICIOS
y= ¿ a"x"
n=O
8. Verificar, por medio del cálculo de L(y) que la ecuación del ejercicio
7 tiene la solución particular
x
y = -ln- x
?
+ In x »=1
¿
00 (-1 ) "xn+l
I + »=1
¿
00 (-1) ,,+1 G+ H,,)
l .
x ..+1
2 n · n. n · n.
http://carlos2524.jimdo.com/
CAPÍTULO 24
Métodos numéricos
(2)
_ _- l_ _ _L -_ _ _ _ _ _....I..-_ _ 'x
o
FIGURA 49
D.X = 0.1
x y y2 X2 (y2 _ X2) dy
0.0 1.00 1.00 0.00 1.00 0.10
0.1 1.10 1.21 0.01 1.20 0.12
0.2 1.22 1.49 0.04 1.45 0.14
0.3 1.36 1.85 0.09 1.76 0.18
0.4 1.54 2.37 0.16 2.21 0.22
0.5 1.76
CUADRO 2
D.x = 0.05
x y y2 X2 (y2 _ x2 ) dy
0.00 1.000 1.000 0.000 1.000 0.050
0.05 1.050 1.102 0.002 1.100 0.055
0.10 1.105 1.221 0.ü10 1.211 0.061
0. 15 1.166 1.360 0.022 1.338 0.067
0.20 1.233 1.520 0.040 1.480 0.074
0.25 1.307 1.708 0.062 1.646 0.082
0.30 1.389 1.929 0.090 1.839 0.092
0.35 1.481 2.193 0.122 2.071 0.104
0.40 1.585 2.512 0.160 2.352 0.118
0.45 1.703 2.900 0.202 2.698 0.135 /
0.50 1.838
http://carlos2524.jimdo.com/
426 Métodos numéricos [Cap. 24
EJERCICIOS
1. Úsese el método de incrementos con f:>.x = 0.1 para calcular, en el
rango O :s; x :s; 1, la solución de la ecuación
y'=x+y
para que la y =
1 cuando x =
O. También, resolver exactamente el
problema de valores a la frontera por métodos elementales y com-
parar los valores de yen intervalos de 0.1 en x.
2. Empléese f:>.x = 00.5 en el problema del ejercicio 1.
3. Empléese el método de incrementos en el problema de valores a la
frontera
y' = e- ZY , cuando x O, Y = O, =
pnmero con f:>.x = 0.2 y después con f:>.x = 0.1, calculando y para
el intervalo O :s; x :s; 2.
http://carlos2524.jimdo.com/
§132] Método de aproximación sucesiva 427
Y1(X) = 1 + r: ( 1 - t 2 ) dt,
y, ( x) = 1 +x - h 3
•
y2()
X -1+
- X + X 2- 1
¡¡x4 -
2
15x 5+ 1
6"3x1 .
CUADRO 4
EJERCICIOS
yl (x ) = 1 r[(
+ o 1 + tr - t 2 ] dt,
yl (x) = 1 + x + X2;
y
Y.2(x) 1 + J: [( 1 + t + t r - t dt,
2 2
]
CUADRO 5
EJERCICIOS
430
Yl(X) = b + r
a
f[t, yo (t )] dt,
(2) y= ¿ anX''',
»=0
http://carlos2524.jimdo.com/
o
~
'" nanX1l-l = - X2 + '"
~
fr.=l n =O
y(x) = ~ G.nx'"
n=O
converge al menos para [xl < 1.
Para el rango en que estamos interesados, °: ;
x ::; i, denotaremos el
error debido a que sólo consideramos hasta el término G.nx" mediante
E .. (x); esto es, sea
Ya que 0< a n ::; 1, es cierto para toda x del rango O ::; x ::; t que
CAPÍTULO 25
Ecuaciones diferenciales
parciales
o sus equivalentes
Xl = HX + Y - z) , yl = HX - Y - z) , Zl = Z,
se puede llegar a menudo a conclusiones incorrectas como resultado del
uso de un subíndice como en Zl, aun cuando Z = Z1. Por ejemplo,
del cambio de variables bajo discusión se sigue que
av
{lz
De aquí que
all'1 cuando Z = Zl .
Retornemos a la ecuación de Laplace. En coordenadas esféricas p,
O, 'P, se relacionan con X, y, Z mediante las ecuaciones
(5)
(6)
REFERENCIAS
(1)
!!:...f' (t ) =o
dt f (t ) ,
df = kf
dt '
de la que es solución general
SOL. U = [Al cos a[3t + A,2 sen a,Bt][Bl cos,Bx + B 2 sen ,Bx],
u = (A3 + A4t) (B 3 + B4X) , .
y
u = [A5 cosh a,Bt + A6 senh a,Bt][B 5 cosh,Bx + B6 senh ,Bx],
donde ,B y las A y las B son constantes arbitrarias.
02V 02V .
2. ~
ex
+ -::;--z = O.
ay
SOL. V = [A 1 ef3 " + A'2e-¡3"][Bl cos,By + B2 sen ,By], y
v = (A3 + A4X) (Ba + B4Y)'
donde [3, las A y las B son constantes arbitrarias; los papeles desempe-
ñados por x y y en la primera solución pueden intercambiarse.
(Fu OU 02U '
3. :::;-;¡ + 2b -;- = a2 ::.;-z.
ot ut uX
Demuéstrese que estas ecuaciones tienen las soluciones
u = g(t) [Bl cos kx + B 2 sen kx],
donde g(t) puede tomar una de las siguientes formas
e-bt[Aleyt + A 2e-yt ], 1'2 = b2 - a2 k 2 , si b2 _ a2k 2 > 0,
o
e-bt[Aa cos 8t + A4 sen 8t], 82 = a 2k 2 - b 2, si a2k 2 - b 2 > 0,
o
e- bt [A 5 + A t], si a k b = 0,
6
2 2
-
2
y las soluciones
u = (A7 + A se- 2bt ) (B3 + B4X).
Hállense todas las soluciones conteniendo el'" y e-kx •
4. OW = yCW. SOL. W = A exp {k (2x + y2)], k Y A arbitrarias.
oy ox
5.
ow ow
x - = W + y- SOL. W = AXkJ'-t, k y A arbitrarias.
OX oy
6. Sométase la ecuación' diferencial parcial del ejercicio 5 al cambio de
variable dependiente W = v/y y demostrar que la ecuación resul-
tante para v es
http://carlos2524.jimdo.com/
§138] Métodos de separación de variables 441
ov
x-= y-o
ov
OX oy
7. Demuéstrese que el método de separación de variables no altera la
ecuación
22U (Fu 02U
ox2 + 4 OX ot + 5 at2 = o.
8. Para la ecuación del ejercicio 7, búsquese una solución de la forma
u = ektf(x)
y por tanto obtener las soluciones
u = exp [k (t - 2x)][Al. cos kx +A 2 sen kx],
donde k, Al, A 2 son arbitrarias. Demuéstrese que también la ecuación
del ejercicio 7, tiene las soluciones u =
Ax + B Y u Ct + D, con =
A, B, C, D arbitrarias.
9. Para la ecuación
22U eu e2u
ex2 + 4x ex + ey2 = 0,
11. Demuéstrese por cálculo directo que si fl (Y) y f2(Y ) son funciones
cualesquiera con segundas derivadas continuas, f~ (y) y 1';
(y), en-
tonces
u = fl(X - at ) + f2(X + at )
satisface la ecuación simple de onda (ejercicio 1),
¡Fu ¡Fu
-2- a2 -
ot - ox 2 •
------~-----------+------------- x
o e
FIGURA 50
(7)
con A B a constantes arbitrarias.
J J
U = ~ u ..
n-1
o
( 13) u (x, t ) ~
00
Bnexp [(n7rh)2
- - - t] sen-
n7rX
n-1 C e
es también una solución de la ecuación diferencial (2).
La u ( x, t ) de la ecuación ( 13 ) satisface la ecuación ( 2) y las condi-
ciones a la frontera (4 ) Y (5 ) . Si u ( x, t) satisface la condición (3),
entonces para cada x en el intervalo O < x < c, el miembro derecho
de la ecuación ( 13 ) deberá tender a f (x) cuando t ~ 0+. Suponemos
que podemos intercambiar el orden del límite (cuando t ~ 0+) y la suma
y concluir que la condición ( 3 ) formalmente requiere que
00 n7rX
(4) f (x ) = ~ Bn sen-, para O < x < c.
11 -1 C
CAPÍTULO 26
Conjuntos ortogonales
140. ORTOGONALIDAD
Lw ( x)fn (x)fm ( x) dx =
9
O para m =F n,
=F O para m = n.
La ortogonalidad es una propiedad que se encuentra con
mucha frecuencia en ciertas ramas de la matemática. También
se emplea a menudo la representación de funciones en series
de la forma
n- O
(1) L b
w (x)x"fn (x) dx = 0, k = 0, 1,2,' .. ', (n - 1) ,
=1= 0, k = n.
PRUEBA. Supóngase que se satisface (1). Ya que x" forma un con-
junto simple, podemos escribir
(2) fm (X)
Si m < n) se sigue que
J
b
a w (x ) fm( x)f... (x) dx
m
= k~O ak Lw(x )x"fn( x) dx = °
b
= an Lw (x)xnfn(x) dx =1= o.
b
http://carlos2524.jimdo.com/
452 Conjuntos ortogonales [Copo 26
( w(x)xnfn(x) dx = mt bm J: w(x)fm(X)fn(X) dx
= bn J: w(x)fn (x) dx =1= O.
2
J: w(x)fn(x) dx = O
y w(x) no puede cambiar de signo en a < x < b.
Supóngase que fn(x) cambia de signo en a < x < b precisamente en
los puntos distintos x = al, <1<2, • • • , as. Las a son precisamente los ceros
de multiplicidad impar de fn(x) en el intervalo. Ya que fn(x) es de
grado n, tiene n ceros, tomada en cuenta la multiplicidad. Sabemos en-
tonces que s ~ n.
http://carlos2524.jimdo.com/
§144] Ortogonol idad de los polinomios de Legendre 453
De la función
( 1)
Entonces, en a < x < b, if¡ (x ) cambia de signo en x = al, « 2, . . . , as
y en ningún otro lado. .
Si s < n, if¡ (x) es de grado < n, entonces
( 1) 1-
P,,( x)=F ( -n, n+1 ; 1;-2- x)
de la sección 120 fueron obtenidos resolviendo la ecuación diferenéial
(2) ( 1 - x2 )y" - 2xy' + n (n + l)y = O.
Los Pn(X) forman un conjuntto simple de polinomios para los cuales
obtendremos ahora una propiedad de ortogonalidad.
De ( 2) tenemos
( 1 - x2 )P;: (x) - 2xP~(x) + n(n + l)P,, (x ) = O,
(8) (n - m) (n + m + 1) r-1
P",(x)P,,(x) dx = O.
(9) f1 Pm(X)P.,(x) dx = O.
( 11 ) r
t
1
p.,2(X) dx = 2n
2
+ 1'
l
(12) P (x)
.,
= -2 1n! D10 (X2 -
n
1) 1o. D
'dx'
= -d
(13) xP~ ( x) = nP.,(x) + P'n-l(X) ,
(14 ) (X2 - 1) P~ ( x) = nxP.. ( x) - nPn-1 ( X) ,
(15) nP.. (x) = (2n - l)xPn--l(x) - (n - I)P"_2 (x).
CAPÍTULO 27
Series de Fourier
c 2n7Tx]e
= [ x - 2m¡sen-C- o = c.
Entonces
re n7TX
J sen2 --dx = c para n = 1,2,3, .
-e C
[5
re n7TX
= J cos 2 - - dx =
-e C
c 2n x]e 7T
= [x + -2
n7T
sen - - =
c o
c.
[6 = r
e
1 . dx = 2c.
http://carlos2524.jimdo.com/
460 Series de Fourier [Cap. 27
Así que
c n7rX
J-c
2
COS - -
c
dx = c para n = 1,2, 3, .
= 2c para n = O.
Hemos demostrado que el conjunto
sen (n7rx/c), n = 1,2,3,·
(A) { cos (m7rx/c), m = 0, 1,2, .
es ortogonal con respecto a la función peso w (x) = 1 sobre el intervalo
- c ::; x ::; c. Hemos también calculado las integrales de los cuadrados
de las funciones del conjunto (A).
110-1
[ (
a" J cos-sen- dx
-e C C
+
+ r n7rX k7rx ]
C
(3)
re
J cos -sen -
-e
n7rX
C
k7rx
C
dx = °
para toda k y n,
y
(4) J( sen -n7l"X sen -
k7l"X
- dx = O para k =1= n; k ,n = 1, 2, 3,.. ..
-e C C
tenemos
bk = ! Je f(x) sen k7l"X dx, k = 1,2,3, .
C -e C
( 6) bn = -1 Je f (x) sen-
n7l"X dx, n = 1,2,3, ....
e -c C
(7) e k7l"x
cos- dx = iao Je k7l"x
+
J-c f(x) e -c
cos- dx
C
+ 00
~
[Je
an
n7l"X k7l"x
cos-cos- dx +
11= 1 -c C C
n7l"X k7l"x ]
+ bn J-ee sen-cos-dx·
c e
El coeficiente de b ll en (7) es cero para toda n y k. Si k =1= O, sa-
bemos que
e n7l"X k7l"x
J
cos -
-e
cos - dx = O para n =1= k,
C e
= c para n = k,
Y también el coeficiente de iao es cero. De aquí que para k =1= O la
ecuación (7) se reduce a
re
J f (x) cos-dx =
k7l"X
ak
Je k7l"X
cos2 -dx,
-c C -e C
http://carlos2524.jimdo.com/
L f (x) dx = iao L dx +
+ n~1 [ an ( cos n;x dx + b,. [ sen n;x dx ]-
Los términos que contengan n ~ 1 son cero cada uno. De aquí que
(12) b,. =-
1 Je I (x) sen--
n7rX dx, n = 1,2,3, ....
e -c e
Antes de proceder con los ejemplos específicos y las aplicaciones nos
comprometeremos a establecer condiciones bajo las cuales la igualdad
( 10) tiene sentido.
Cuando a,. y b,. están dadas por (11) y (12) el miembro derecho
de la ecuación (10) se llama serie de FourierJ sobre el intervalo
-c <;x <; CJ para la función f (x). Se establecerán de inmediato las
condiciones suficientes para asegurar que la serie de Fourier en (10)
representa la función f( x) de una manera razonable.
Sea f(x) continua y diferenciable en cada punto del intervalo
- c <; x <; c excepto para a lo más un número finito de puntos, y en
http://carlos2524.jimdo.com/
§147] Series de Fourier: un teorema de desarrollo 463
f(x)
--~~--------~~--------------x
-e O e e
T
FIGURA 51
f(x)
•
----~~-+----~--+---~---+----~--~--~---x
o
---
-2c 3c
-2 •
e 3. 2c
"2 "2
>----<
FIGURA 52
repite tal estructura una y otra vez fuera del intervalo. Para la función
mostrada en la figura 51, la serie de Fourier correspondiente deberá
converger a la función periódica mostrada en la figura 52. Nótese que
la convergencia tiende al valor promedio de las discontinuidades, la
periodicidad y la forma en que las dos juntas determinan el valor al que
la serie converge en x = c y x = -c.
Esta exposición será ampliamente ilustrada con los ejemplos numé-
ricos y los ejercicios de la siguiente sección.
I(x)
-------~2-------0~------~------x
FIGURA 53
y
= i Lr 2 f (x) sen-
2
mrx
(3 ) bn -dx;
2
n = 1,2,3, ....
Como la descripción de f (x); está dada por dos fórmulas una para el
intervalo - 2 < x < O y la otra para O < x < 2, es conveniente separar
las integrales ( 2) y (3) en las partes correspondientes. Así que susti-
tuyendo f (x) de ( 1) en la integral (2) llegamos a la forma
rO n7rX r2 n7rX
(4) an = f J 2 cos - dx + i J x cos - dx.
-2 2 o 2
Para estas integrales, el método de integración diferirá de acuerdo con
que n = O o n =F O.
Si n =F O, entonces
an = 2
- [n7rxJo
sen-- 2 n7rX ( -2
+ i [ -xsen-+ )2cos -2-
n 7rxJ2
o'
n7T 2 -2 n7r 2 n7r
o
an = n
2
7r[0 - OJ + f [O + (n:Y cos n7r - O - (n:Yl
De donde, para n =F 0, la a.. está dada por la fórmula
_ -2 (1 - cosn7r)
(5) an - , n = 1, 2,3, ....
n 27T 2
Para n = O las integraciones anteriores no son válidas (división entre
n), pero regresando a (4), hacemos n = O, y obtenemos
2
ao = i rO 2 dx
L2
+i r x dx
Jll
de donde
b .. =~
mI'
[- 7r
cos n xJo
2 -2
+ i [- (~) x cos n7rX + (~)2 sen n 7rxJ2,
n7T 2 n7r 2 o
pero como cos ( -n7r) = cos n7r,
o
(6) b" = - ~;
n7r
n = 1, 2,3, ....
( 10)
I(x)
-6 -4 -2 O 2 4 6
FIGURA 54
Por tanto
1=i!---~
4 00 1 __2"'0~_
2 7T2k =O(2k + 1)2 7T n =l n '
o
'" 1 7T 2
( 11 )
?~o (2k + 1) 2 = 8
Para x = 1, la serie en (10) vale uno otra vez. Usando x =1 en
( 10 ), llegamos a
Para n impar, n = 2k + 1,
sen (2k + 1)71' = sen (b + t7l') = cosk7l' = ( - 1)7<.
2
Así llegamos a la ecuación
1 = ~ _ ~~ ( ._ 1)7<,
2 71' 7<=0 2k + 1
o
""( _ 1)7<_71'
~
+ 1 - -4'
(12)
7<=02k
que puede verificarse siempre por el hecho de que el miembro izquierdo
representa arctan 1.
f(x)
---+----~~----+---~~~---+----~~----4---X
FIGURA 55
EJERCICIOS
" =1
+ n7!'
sen-,
c
o
f (x) ,. . ., ! +
~ i¡ sen [(2k +
l)7rx/c].
7r k =o 2k + 1
6. Intervalo, -c < x < c; función f (x) = x 3 •
SOL. f(x) ,.....,2~3 f ( _1 ) " +1 (n 7!'2 - 6) ~en (n7rx/ c).
2
71' "=1 n
7. Intervalo, - 7r <X< 71'; función
f (x) = 37r + 2x, -71' < X < 0,
= 7r + 2x, O< x < 7r.
00 sen 2kx
SOL. f (x) ,....., 27r - 2 "~1 -k-'
8. Intervalo, -c < x < c;
función f (x) = x(c x), - c < x < O, +
= (c - x) 2, O < X < e,
2 2
c c 1 [ n7!'X n7rx]
f (x),....., 12 +"2 ~"2 {3 + (-1)"} cos- + n7rsen- .
00
SOL.
. 7r "=1 n c c
9. Intervalo, -2 < x < 2; función f (x) = x 1, -2 < x < O, +
= 1, O ~ x < 2.
SOL. j(x),.....,!
2 1
+ "2 + n7!'( -1) 1>+1 sen !n7rx].
00
~ "2 [{l - (-1) ,,} cos !n7rx
7r 11=1 n
http://carlos2524.jimdo.com/
472 Series de Fourier [Cap. 27
X
< 0,
< 7r.
~ 7r ~ (- 1) " cos nx +
2
SOL. f (X ) ~ 6 + 2 "'=1
~
n
2
1 100 .
+- .~ _g[(_1 )n+1n 27r 2 - 2 + 2(-1)"]sennx.
7r "=1 n
12. Intervalo, -7r < X < 7r; función f (x) = cos 2x.
SOL. f(x),-..' cos 2x.
13. Intervalo, -7r < X < 7r; función f (x) = cos (x/2).
2 4 00 ( - 1) n+) cos nx
f(x ) ,-..'; + ;"~1 (2n _ 1) (2n + 1)
SOL.
14. Intervalo, :"-'7r < x < 7r; función f(x ) = sen 2 x .
SOL. sen 2 x ,-..' i - i- cos 2x.
15. Intervalo, -e < x < e; función f (x) = eX.
senh e
SOL. f(x ),-..' --+
e
~ 2 ( -1 )"senhe[ecos (n7rx/ e) - n7rsen (n7rx/e)]
+ e2 + n 27r 2 '
16. Intervalo, - e < x < e; función f (x) = 0, - e < x < 0,
= e-x, < x < e.
1 - e-e
°
SOL. f(x ),-..' 2e +
~ 1- ( -1 ) "e-e (n7rx n7rx)
+ "'=1
..:;.
e
2 + 2
n 7r
2 eCos -
e
+ n7r sen -
e
.
17. Intervalo, -e <x< e; función f(x) = 0, -e < x < ie,
= 1, ie < x < e.
SOL. f( x),-..' i -
- -1 ~ 00
-1 [ n
sen T
i-n7r coso
- x + (cos n7r - cos in7r) sen-
n7rx] .
-:r "=1 n e c
18. Intervalo, - e < x < e; función f (x) = 0, -e < x < 0,
= x,
f(x ) ,-..' ie -
°< SOL.
x < e.
e 00 1 [ n7rX + n7T ( -1 ) " sen -nTox]
- 2 ~ -; {1- ( - 1) ,,} cos - .
Ti n .::;l n c c
= x (c - x),
19. Intervalo, -4 < x < 4; función f(x) = 1, - 4 < x < 2,
< x < c. °
= 0, 2 < x < 4.
20. Intervalo, - e < x < e; función f(x ) = 0, -c < x < 0,
= x(c - x), < x < e. °
http://carlos2524.jimdo.com/
§149] Series senoidoles de Fourier 473
21. Intervalo, -c
22. Intervalo,
< x < c;
-c < x < c;
función f(x) = c + x,
= 0,
función f (x ) = x 4 •
-c
°<< xx << c. 0,
c4 00 n 27T 2 - 6 n7TX
SOL. f(x) ,..,-+
5
8c 4
~
7>=1
(-1 ) "
n
4
7T
4 cos-.
c
23. Empléese la solucióh del ejercicio 3 para demostrar que
00 ( _ 1 ) n+l 7T 2
~ =-.
7>=1 n2 12
~ ~=
2
24. Empléese la solución del ejercicio 8 para demostrar que 71"6 •
n=l n
25. Empléese la solución del ejercicio 22 para demostrar que
00 1 71"4
~-=-.
4
7>=1 n 90
26. Empléese x = O en la solución del ejercicio 15 para sumar la serie
00 (_1)" c - senh c
~1 c + n 7T2 2 2 SOL.
2c 2 senh c •
27. Supóngase que c ~ O en el resultado del ejercicio 26, y compárese
con el ejercicio 23.
e jo
C
(2) b" =- f (x) sen-dx, n = 1, 2, 3, . . . .
e
La representación (1) se llama serie senoidal de Fourier para f(x) sobre
el intervalo O < x < c.
Debe hacerse notar que el artificio de introducir la función g (x )
fue una herramienta para llegar a (1) Y ( 2 ) ; no hay necesidad de
repetirlo en problemas específicos, o sea, aquellos que resolveremos por
uso directo de p) y (2) definidos anteriormente.
(3)
n:::l
en la qlle
(4) b" = 2 J: X2 sen n~x dx
(5)
/
______ --~~~------~~~----~--~~~--~x
-3 3
/
FIGURA 56
EJERCICIOS
13. Intervalo O < x < 7T; función f(x) = cos 2x. Nótese el procedi-
miento especial necesario para el cálculo de b2 •
4 ~ (2k + 1) sen [(2k + 1)x]
SOL.
f()x~:;;: k":O (2k _ 1) (2k + 3)
14. Intervalo, O < x < 7r; función f (x) = cos X.
8 00 k sen 2kx
SOL. f(x) ~::;;: k~l 4k 2 _ 1 .
15. Intervalo, O < x < e; función f(x) = e-z.
~ 2n7r[1 -
( _1 )ne-c] sen (n7rx jc)
f () ~ ~
e + n-7r
SOL. X 2' 2 •
n=l
"' = 1
(2"2sen---cos-
n
n7T
2 n
n7r) sen-o
2
n7rx
c
7T
co n7TX
(1 ) f ( x) "'""' iao + ~ a n cos - - , O < X < e,
n=l e
en la que
2 JC n7TX
(2) an =- f(x) cos-dx.
e o e
I(x)
FIGURA 57
Inmediatamente tenemos
n7TX
f(x) "'""' itl{¡ + n~co= l an cos-,
e
en la que
an =-2 JC xcos-dx.
n7TX
e o e
Para n =1= O, la an puede calcularse como sigue
2 [ -xsen--+
an =-
e n7T
e n7TX (e
e
-
n7T
)2 cos--
n7T xJC
e o
http://carlos2524.jimdo.com/
§150] Series cosenoidales de Fourier 479
EJERCICIOS
SOL. f(x)
4
3
e 6c ~ - -1)+
,.....-+-
-rr 2
3
n
"
-2 . 1 -
00
-rr
n=1
[( -
2 2
( - 1)
n4
"J n-rrx
cos-·
e
16. Intervalo, O < x < e ; función f( x) = x 4
•
( 13) u (x , t) = ~
00
Bnexp [(n7rh)2
n=1
- - n7rX
t sen-
C
J C
152. ELECTROCARDIOGRAFíA
Puesto que las series de Fourier proporcionan una representación
periódica para una amplia gama de fun ciones, es razonable pensar que
las series de Fourier, tengan conexión con funciones que sean verdade-
ramente periódicas en su estado natural. Un cierto tipo de tales funcio-
nes que son naturalmente periódicas es ampliamente empleado en
medicina, particularmente como una ayuda en el diagnóstico de enferme-
dades del corazón. Se descubrió hace mucho tiempo que existe una
corriente eléctrica causada por la actividad propia del corazón y que
era posible registrar gráficamente el voltaje debido a dicha corriente.
El producto final se llamó un electrocardiograma. Las figuras 58-61
muestran partes del electrocardiograma.*
,¡
¡-.¡
:1 ;
; ;':
!;
FIGURA 58
g uía "
FIGURA 59
guía 11
FIGURA 60
FIGURA 61
( 1)
i: (- 1 )"n cos nx
n=l n3 +7 .
La serie de la ecuación (1) converge absolutamente puesto que su
00
término general es menor en valor absoluto que 1/n2 y ~ l/n?' con-
verge. Llamemos I" (x) la suma de la serie (1). k=l
(2) X
2 -
-
7T
-
2
+4 ~ ( -1 )" 2cos nx , para
~
-7T _< X ::::;
_ 7T.
3 "=1 n
Por tanto,
2
(3) i (X2 _ 7T ) = ~ (-1)" cos nx -7T ::::; X ::::; 7T.
3 n =l n2
Puesto que los coeficientes de la serie (3) y los de la serie para 1" ( x )
_ ~ (- 1 ) "n cos nx
(4 ) 1"
(X ) - ,,'7:1 n3 + 7 ' O< x < 7T,
son casi iguales para n grande, es obvio que la diferencia de dichos
coeficientes debe ser pequeña. De tal manera que si restamos los miem-
bros de la ecuación (3) de sus correspondientes en la ecuación (4)
obtenemos
1" (x) - i (X 2 - 7T )
3
2
= i: (- 1 ) ,,[n
n=l
3 n
+7
- l.J
n2
cos nx
00 7 ( - 1 ).,...1 cos nx
= 1>~= 1 n 2( n 3 + 7 ) .
http://carlos2524.jimdo.com/
§ 153] Artificios para obtener mayor rapidez en la convergencia 487
CAPÍTULO 28
Problemas de valores a la frontera
ou _ 2 ( 02 U 02U 02U)
( 1) a¡ - h ox2+ oy2 + {lz2 '
No es difícil, pero requiere conceptos que están más allá del nivel
de este curso, verificar que (6) con los coeficientes Bn dados por (8), es
realmente una solución, esto es, que (6 ) posee, para una elección ade-
cuada de f (x), las propiedades de convergencia requeridas, además de
que satisface formalmente la ecuación diferencial (2) Y las condiciones
a la frontera (3), (4) Y (5).
La cantidad de calor que fluye a través de un elemento de superficie
en un tiempo especificado, es proporcional al grado de cambio de la
temperatura en la dirección normal (perpendicular) a esa superficie.
Por tanto, se tiene que el flujo de calor en la dirección x (a través de
una superficie normal a la dirección x) es
-K ox
ou'
http://carlos2524.jimdo.com/
§ 154] Ecuación térmica unidimensional 491
C() n7l"X
(18) f ( X) ro-' n~l b n sen -2-' para O < x <2
donde
f (x) = Uox, para O < x < 1,
Y f (x ) se ha elegido de tal forma que para el intervalo 1 < x < 2 en
( 18) se hayan eliminado los términos pares en n.
. Físicamente, no es difícil ver que deseamos extender de alguna forma
la anchura de la losa más allá de x = 1, para evitar que el calor fluya
a través de la cara aislada x = 1. Una vez que este hecho se ha
llevado al cabo, se sigue inmediatamente que todas las condiciones de
temperatura sean simétricas con respecto a la cara aislada x = 1.
((x) f(x)
Uo
o o 2
FIGURA 62 FIGURA 63
EJERCICIOS
SOL.
1 8uo 00
u = 60
280 ~ (-1 ) "
+ -;- ,,~ 2k + 1 exp
[
-
(2k +1
1000
Ft].
7. Dos losas hechas del concreto mencionado en el ejercicio 6 (con
h 2 7r 2 = 0.4 pie 2 /hr), una de 15 pies de grueso y la otra de 5, están
colocadas una al lado de la otra. La más gruesa está inicialmente
a la temperatura de 120°F mientras que la más delgada está
a 30°F. Las caras exteriores se mantienen a 30°F para t > O. En-
cuéntrese la temperatura en toda la losa para t > O. Medir x desde
la cara exterior de la losa más gruesa.
SOL.
180 ~ 1 - cos (3n7r/4)
u = 30 +- ¿;.,
( n2
exp - - - sen--·
n7rx t)
7r n=l n 1000 20
8. Dos losas del mismo material, una de dos pies de grueso y la otra
de uno, están puestas una al lado de la otra. La más gruesa está
inicialmente a la temperatura A y la más delgada a cero grados °F
Las caras exteriores se mantienen a una temperatura cero para
t > O. Encontrar la temperatura en el centro de la losa de 2 pies
de grueso. 1 4A
3 n7r [ (hn7r)2]
00
SOL. - ~
-sen -exp - - - t·
n 3 7r 3n=l
00 (-I)"(2k + 1) 27r
[:::~:::::::::::I::::::::::::::]
O L 2L
u =B
FIGURA 64
FIGURA 65
A2 en el segundo mes
etcétera.
La función A (t) puede expresarse, A(I¡
A3
empleando to = O, por :-~
I A2 : I :
(5 ) A (t) = An,
:: ~'
para tn-1 < t < tn ; I
Al::
I
I
l
::
I I
n = 1,2,3, . I
I
I
I
"
I
I
I
--~
4 ~ 1 [(2k
ap -
+ 1)2Prr2t] Kn-. .(2_k_+_1-,---)7T_X
:. .
7T k=o 2k +1 c2 c'
para O < t, O < x < c. Directamente, o por comparación con el ejer-
cicio 1, sección 149, la función F (x, t) de (6) puede verse que es una
solución de la ecuación térmica (1) Y tiene las siguientes propiedades:
(7) Cuando t ~ 0+, F(x, t) ~ O, para O < x < c;
( 8) Cuando x ~ O+, F(x, t) ~ 1, para 0< t;
(9) Cuando x ~ c-, F(x, t) ~ 1, para 0< t.
http://carlos2524.jimdo.com/
500 Problemas de valores a la frontera [Cap. 28
debido a (8), (9) yel hecho de que F (x, t - tn-l) ==0 para t::; tn-l.
En la práctica real, los cálculos con la solución (11) se simplifican
muchísimo porque la F (x) t) de (6) es esencialmente la misma para
todas las difusividades y todas las anchuras de losa. En la ecuación (6),
hagamos h27r 2 t/ c2 = T Y 7rx/ C = {;. La forma como se eligen los nuevos
argumento T y {; será discutida en la sección 167.
Ahora podemos escribir
t
Sea R el radio de la esfera y ( p) la temperatura inicial. Entonces el
problema que se nos presenta es
h2 ( 02 + ~ au)
U
(4) au = para O < t, O < p < R;
ot ap2 p op ,
(5) Cuando p -,) R-, u -,) O, para O < t;
(6) Cuando t-,)O+, u-,)f(p), para O:S; p < R.
El estudiante puede comprobar fácilmente que el cambio de variable
dependiente,
(7) u =-,v
p
transforma el problema (4) a ( 6) en el nuevo problema
2
EJERCICIOS
1. Resolver el problema (4) a (6) por el método arriba señalado.
SOL. u = ; n~l bn exp [ - (h~7rYtJ sen n~p, en donde
2 rE n7rp
bn = RJ pf ( p) sen¡¡ dp•
o
2. Una esfera de radio R tiene inicialmente una temperatura constante
Uo en toda ella, para t > O mantiene su superficie a la temperatura
Ul. Encontrar la temperatura en toda la esfera para t > O Y en
particular la temperatura U c en el centro de la misma.
SOL. u = Ul +
2R (uo - Ul ) ~ (_l)n+l
¿;. exp --' - - t sen-;
[(hn7r)2 ] n7T.p
p7r n =l n R R
Uc = Ul + 2(uo - Ul) ~l ( - l ) n+l exP [ _(h~7ry tJ
para t > O.
158. ECUACIóN SIMPLE DE ONDA
Si una cuerda elástica que se mantiene fija en dos puntos está tensa,
se desplaza entonces de su posición de equilibrio soltándose después, los
desplazamientos subsecuentes a partir de su posición de equilibrio pueden
determinarse resolviendo un problema de valores a la frontera. La figura
68 muestra un desplazamiento representativo de la cuerda que se man-
tiene fija en x = O Y x = c. El desplazamiento y para O < x < c y O < t
se encuentra a partir del desplazamiento inicial conocido, (x ), la veloci- t
dad iniciall" (x) y el hecho de que y debe satisfacer la ecuación de onda
unidimensional
(1 )
~o~----------------~c--x
FIGURA 68
EJERCICIOS
En los ejercicios 1 a 5 encuéntrese el desplazamiento para t > O del pro-
blema de la curva vibrante de esta sección, bajo la condición de que la
velocidad inicial sea cero y que el desplazamiento inicial esté dado por
la f (x) descrita.
1. f (x) = x, para O ~ x ~ ej2,
= e - x,paraej2 ~ x ~ e.
4e ~ ( -1 )k cos[(2k + l )a7Ttje] sen [(2k + l)7Txje]
SOL. y = 7T2k~O (2k + 1 )2 •
((x)
----~----~------------~--~~-x
o e 3c
"4
e
4
FIGURA 69
f(x'
FIGURA 70
(1)
u:f(x)
b~------~----~
--~-------_L-_--.x
O u=O a
FIGURA 71
02U 02U
(2) OX 2 + oy2 = O, para O < x < a, O< Y < bj
EJERCICIOS
SOL.
2 ce
u = i + ;; k~O R
(r
)2k+1 sen (2k -f-' 1)8
2k+ 1 .
11. Una placa con fronteras circulares concéntricas r = a, y r = b,
O < a < b, tiene su frontera interna a la temperatura A y la externa
a la temperatura B. Encuéntrese la temperatura en toda la placa.
B In (r/a) - A In (r/b )
SOL . u= In (bja) •
12. La placa del ejerclClO 11 tiene su orilla interior r = a aislada, su
orilla externa está a la temperatura unitaria para O < O < '7T Y a cero
para '7T < 8 < 2'7T. Encuéntrese la temperatura en toda la placa.
2 00 (r jb ) 2k+1 + (ajr) 2k+1( ajb) 2k+1 sen (2k + 1)0
SOL. u=i+;eo 1+ (ajb) 41<+2 • 2k+l .
13. Un prisma plano está definido en coordenadas polares por la región
O < r < R, O < O < ¡J,. Encuéntrese la temperatura a todo lo largo
del prisma si las aristas O = O Y O = f3 se mantienen a la temperatura
cero y la arista curvada r = R está a la temperatura unitaria.
SOL. u = ±f (:"-)(2k+1
l
1r//3sen [(2k + 1)7rO/f3J.
'7Tk=o R 2k + 1
14. Para el prisma plano del ejercicio 13, encontrar la temperatura si la
arista O = O está a cero, la arista O = f3 está a la temperatura unitaria,
y la arista curvada r = R se mantiene a la temperatura f (O) para
O < O < f3.
SOL. U = Ojf3 ~+
00 r)"
7r
CAPÍTULO 29
Propiedades adicionales de la
transforIllada de Laplace
00 x"
( 2) e"'= ~- para toda Xj
n!'
n=O
511
http://carlos2524.jimdo.com/
512 Propiedades adicionales de la transformada de Laplace [Cap. 29
00 X2n+l
(6) ~ (2
senh X = n=O • n +
1) .l' para toda X;
00 (-1 ) "X21H1
( 7) ~ 2n + 1 '
Arctan x = n=O JxJ < 1;
1 )"x1H1
+ x) = ~
00 (
(9) In (1 , JxJ < 1;
"=0 n +1
1+x 00 x 2n +l
( 10) ln--=2~
1- x n=O 2n + 1, JxJ < 1;
Al buscar la transformada de Laplace o la transformada inversa
de una función da4a, pueden .e!1contrarse serios inconvenientes y difi-
cultades, o inclusive ser imposible obtener el resultado deseado haciendo
uso directo de los teoremas de los capítulos 11 y 12 en un número
finito de pasos. En este caso, volveremos a usar frecuentemente las
series infinitas. Si desarrollamos nuestra función en una serie tal que
podamos conocer la forma de obtener la transformada deseada o la
transformada inversa de cada término, podremos resolver nuestro pro-
blema original. .
L-1 { f(s) }
senh (es) .
Sabemos que senh z = i (e Z - e-Z ) . Entonces
f(s) 2f (s)
(11 )
senh (es) - eC8 _ e-C8 •
Para h > 0, s > 0, sabemos cómo calcular L-l{e-h8 f(s)} por el teorema
19, sección 65. En efecto,
(12) L-l{e-h8 f (s)} = FU - h)a (t - h), h > 0, s> O.
Por tanto, representamos ( 11) como
f (s) 2f (s ) e-C8
(13)
senh (es) 1 - e- 2CS
debido a que podemos emplear la serie de potencias ( 1) para desarro-
llar (1 - e- 2C8 ) -1 en una serie de exponenciales con argumentos nega-
tivos. De (1) obtenemos
http://carlos2524.jimdo.com/
§ 160] Series de potencia s y transformadas inversas 513
1 00
(14)
sen
f~?es ) = 2 }; f (s) exp ( -2nes - es ).
n=O
e- t = 00
~
( _ 1 ) nf'
=1+ 00
~
( - 1)
-'---:'----
nt n
n=O n! n=1 n!
Por tanto, podemos escribir
1 - e- t 1 ) n+1['>--l
t
~
n=1
00 ( _
n.
,.
Un corrimiento en el índice, de n a (n + 1) nos da
1 - e- t 00 ( - 1 ) nt n
t =~ (
n=O n +
1) .1"
tn} -1. De donde
Sabemos que L {-n! = sn+1
L{1 - e- t
} = (- 1)"
~o (n + 1) s"+1'
00
t
comparando con la ecuación (9) anterior, se tiene
(16) s> O.
http://carlos2524.jimdo.com/
514 Propiedades adicionales de la transformada de Laplace [Cap. 29
{ 1} = (2nt ) !
2n
1
Ahora L- - - - o De donde
S2n+1
L - 1 {In s +
s- 1
1} = 2~ n=O (2n
t
+
2
"
1) !
,
(17)
S +
L-1 { l n
s- 1 t
1} 2
- - = - senh t.
EJERCICIOS
,
1. Cálculese L - t - . {sen kt} SOL.
,
Arctan - ,
k
s
s> O.
2. Cálcúlese L {
kt} •
1 - cos
t SOL.
SOL. F(t ) = ~
'" (-1) "(t - 4n -2)2a(t - 4n - 2); F( 12 ) = 68.
n ::::: O
7. Dado rp(t) = L-1 {s4sel (3S) } Calcúlese rp( lO). SOL. 344.
8. Dado e> O, s> O, Y L-l{f(S ) } = F(t) . Pruébese que
L-l {_f~} = 2 i (-l)nF(t - 2na - c)a(t - 2nc - e).
cosh (es ) "=0
http://carlos2524.jimdo.com/
L-1 { 1 }_~
Vs +1 - V7rt'
Entonces el teorema de convolución nos da
L
-1 { 1 } 2
s v' s + 1 = v' 7r
Jvo t exp ( -"1 2
) dy.
Esto es
( 3) L -1 { 1 } = fer ( Yt) .
sv' s + 1
En nuestro trabajo son útiles algunas propiedades básicas de fer x
que obtendremos inmediatamente. Directamente de la definición (1)
se sigue que la derivada de fer x está dada por
(5) fer x = 2
-= ~
00 ( -1 )"x 2n
-'--_-'---__ +1
v' 7r " =0 ( 2n + 1) n!
En cálculo elemental encontramos que
oo v'-;
( 6)
Jo exp ( - W ) df3 = 2'
De (6 ) se obtiene
(7) Lim fer x = 1.
o
Hagamos (3 = k en ( 12) de tal forma que los límites de integra-
vv
ción se transformen en v = t Y v = O. Ya que d/3 = -1kv-~ dv obte-
nemos, empleando el signo menos para invertir el orden de integración,_
Sea ahora
(15)
(16)
Y
(17)
o
(18)
L {ir! exp ( - ~2) + k2r~ exp ( _ ~2)} = sL {tt eKp ( _ ~2)} - O.
Debido a (15), (16) Y (17), la ecuación (18) puede escribirse como
2
_idA + k2A = sd A.
ds ds 2
Por tanto la función deseada A (s) es una solución de la ecuación dife-
rencial.
(19)
r
8 ... 0+ 8 ... 0+
(20) Jor
Y
(k2)
v-~ exp - -;- dv = k)
Ty; eríc ( yy'
De donde
. ,r;. (k)
Llm A (s) = -k Llm ferc _r = -k ferc O = -k.
y; y-:;;'
8 .... 0+ y .... '" YY
De donde
d 2A 1 d 2A 1 dA
J - = -2- - - -
ds 2 4 dz 4z dz
y la ecuación (19) se transforma en
2
(21) d A _ 4k 2 A = O.
dz 2
La solución general de (21) es
A = blexp (-2kz) + b2exp (2kz),
entonces la solución general de (19) es
(22) A = blexp (-2k Ys) + b2exp (2k ys).
Debemos ahora determinar las constantes b1 y b2 de las condiciones
que indican que A ~ O cuando s ~ 00 y A ~ y; / k cuando s ~ 0+.
Cuando s ~ 00, A no tiende a un límite a menos que b2 = O. Entonces,
dejando que s ~ O+, obtenemos
y-:;;'
T=b 1 •
Por tanto
{ , (k
2
A(s) =L r 2 exp -t )} = y;
T exp (-2kYs).
* Tal cambio de variable está dictado por la prueba de la página 16 del libro de
E. D. Rai!lville, Intermediate Differential Equations. 2" Edici6n. Nueva York: The
Macmillan Co., 196.4.
http://carlos2524.jimdo.com/
520 Propiedades adicionales de la transformada de laplace [Cap. 29
Emplearemos ( 23 ) en la forma
EJERCICIOS
1. Demuéstrese que para toda x real, Ifer xl l. <
2. Demuéstrese que fer x es una función impar de x.
'
3. D emuestrese que L'1m - fer- 2
x = -=.
,, ~ O x V 71'
4. Aplíquese el método de integración por partes para demostrar que
Jo ferydy = xferx _
I/) 1 [1- exp ( _ X2 ) ].
V7r
5. Obténgase la ecuación (11), sección 161.
6. Principiar con la serie de potencias para fer x, ecuación (5), sec-
cción 161, y demuéstrese que
1
9. Calcúlese L-l {.¡:;: 1+ 1}' SOL. _ r-; - e t ferc (yt).
y 7rt
dA =
ds -2 Joo
o exp ( - sf32 - P[3-2) df3
= -2 exp (-2k Ys) r exp [- ([3 ys - kf3-1 ) 2] d[3.
dA -
Ys-
ds
-kA = - 2 Y7r exp (-2k Ys)
Obtener de aquí la ecuación diferencial A (s) .
http://carlos2524.jimdo.com/
523
(2) L
-1 {.!.
s exp
(_
s
!!.)} -_~o ( -1k !k)"x"t"
00
!
Cuando n = O en (1) obtenemos, ya que r (k + 1) = k!,
00 ( -1 )" ( tz)2"
(3) Jo( z) = k~O k!k! •
Comparando (2) y (3) se tiene
De ( 1) obtenemos
( _1 )k(fx)2kt2k
OC>
Jo( xt ) = :¿ -'-----'-:~:-'--
70=0 k!k! .
Entonces
L{f ( t ) }=
o x 70=0
i:
( -1)k (fx)2k(2 k ) !
k!k !S27Hl
Pero ( 2k ) ! = 2kk![1 ·3·5· .. ( 2k - 1)J. Por tanto
o
L{fc(xt) } =-1 ( 1
S
+ ~S22)-11 .
Por tanto
(6) L{fo(xt ) } = 1
YS2 + X2
De (1) es fácil concluir que
d
dz Jo (z) = - Jl (Z) .
Entonces
y obtenemos
(7)
EJERCICIOS
(x)} = (t)!
1 exp S
L-l { sn+l -; 1I In ( 2 V-xt ) .
http://carlos2524.jimdo.com/
525
CAPÍTULO 30
Ecuaciones diferenciales parciales;
Métodos de transformación
a
renciación con respecto x son conmutativas, así como las transformadas
de Laplace lo son con respecto a t.
Debido a que la ecuación (1) tiene coeficientes constantes, no apare-
cen las derivadas con respecto a la variable transformada s. La ecuación
diferencial parcial (1) será transformada en una ecuación diferencial
ordinaria con variable independiente x y donde s aparece como un
parámetro. En vista de lo estipulado por la relación (6), la aplicación
del operador L transforma (1), (2) Y (3) en
(7) x> O.
( 10)
EJERCICIOS
1. -oy oy =
+ 4-;:;- - 8t> parat> 0, x> O,.
ox ut
t -7 O+, Y -7 0, para x > O;
X -7 ()+, Y -7 2t 3, para t > O.
SOL. y(x, t) = _t 2 + 3(t ~ 4x )2a(t - 4x).
oy
2. ox + 2 oy
ot = 4t, para t > 0, x > O;
t -7 O+, y -7 0, para x> O;
X -7 O+, Y -7 2t 3 , para t > O.
SOL. y(x>t ) = t 2 + [2(t - 2X )3 - (t - 2x) 21a( t - 2x).
3. Resolver el ejercicio 1 con la condición cuando t -7 0+ reemplazada
por t -7 O+, Y -7 )O.
SOL . y(x,t) = x - ¡t - t 2 + [3 (t - 4X)2 + Ht - 4x)]a(t - 4x).
4. Resolver el ejercicio 2 con la condición cuando t -7 0+ reemplazada
por t -7 O+, Y -7 2x .
SOL. y(x, t ) =2x-t+t 2+ [2 (t-2x ) 3 _ (t- 2x )2 + (t -2x) la(t-2x).
o2y o2y
5. ox 2 = 16 ot 2 ' para t > O, x > O;
t -7 0+, Y -7 0, para x > O;
t -7 0+ o~ -7 - 2 para x > O;
, ot '
X -7 O+, Y -7 t, para t > O;
Lim y(x> t ) existe para t > O.
SOL. Y = 3(t - 4x)a(t - 4x) - 2t. ~
02y ¡)2y
6. ot 2 = 4 ox 2 ' para t > O, x > O;
t -70+, y -7 O, para x> O;
t -7 0+ Oy -7 2 para x > O,·
, ot '
X-7 O+, y -7 sen t, para t > O;
Lim y(x, t ) existe para t > O.
2>->00
SOL. Y 2t + [sen (t - ix) -
- 2 (t - ix)]a (t - ix).
Y =O y ay =O
at '
para x ~ O.
FIGURA 72
forma que F' (O) Y F" (O) se anulen cuando lo hace F (O), de tal forma
que la continuidad de y y sus derivadas no se interrumpa a lo largo de
la recta x = ato Se deja al lector como ejercicio el completar la veri-
ficación de la solución.
En la sección 158 estudiamos el desplazamiento tranversal de una
cuerda de longitud finita fija en ambos extremos; para tales problemas
parece más apropiado el método de series de Fourier que el de la trans-
formada de Laplace. Trátese, por ejemplo, el método de la transfor-
mada en el ejercicio 1 de la sección 158.
EJERCICIOS
1. Interpretar y resolver el problema
(J2y a2y
ift2 = ax 2 ' para t > O, O < x < 1;
t ~ O+, Y ~ x - x 2
, para O < x < 1;
t ~ 0+ ay ~ O para O < x < 1,·
, ot '
x~ O+, y~ O, para t > O;
x ~ l-,y ~O, parat > O.
Verificar la solución directamente
00
SOL. Y = x - X2 - ( l. + n=O
¿ (- 1) n[ (t - n - x) 2a (t - n-x ) +
+ (t - n - 1 + x ) 2a (t - n - 1 - x)].
2 d w
2
sw - A = h dr' x> 0,
o
(7) x> O.
Y
B
(9) x~ O+, w ~ -s [1 - exp ( - tos)] .
y la w de (10) se transfonna en
(12) w = clexp ( -
X
-h-
-..Is) +"5'
A
(16) L-
l
H exp ( - x :s) exp (-toS)} = fercCh lt ~ tol~) a (t - t o),
donde los signos de valor absoluto han sido insertados para permitir
que t sea empleada en el rango O a t o en donde la función a forzará al
miembro derecho de la ecuación (16) a anularse.
Estamos ahora en posición de escribir la transformada inversa de
la w de la ecuación (14). Para x > O Y t > O,
( 1) au = h2 au
2
para t > O, O< x < c;
at ax 2
'
u = temperatura (OF)
t = tiempo (horas )
x = . coordenada espacial (pies)
h 2 = difusividad térmica (pie 2 por horas)
e = longitud (pies )
A = temperatura inicial (OF )
Buscaremos ahora las nuevas variables sin dimensiones ?:, T, "', pro-
porcionales a las variables físicas x, t y u. Por el momento supongamos
( 5) x = /3t,
en donde /3, y, S son constantes positivas que deben ser determinadas de
tal forma que las nuevas variables sean sin dimensiones. Los cambios
de variables (5) transforman a las expresiones de (1) a (4) en
S o'" _ h S 0
2 2
",
para T > O, 0< /3t < e;
:y OT - 7F ot 2
'
1 h2
---
Y ¡32'
de donde
(17)
(2 ) t O+, u ~ A,
~ para O < x < c;
(3) x~ O+, u ~O, para t > O;
(4) x~ c-, u~O, para t > O.
Resolveremos el problema correspondiente en variables canónicas.
Esto es, en las ecuaciones ( 1) a (4) hacemos
x
, =-,
(5)
C
tf ~ 1 - ~
ca [ ferc + 1) -
(2n2 y; ferc (2n2 Y+ l)J -
n=O T T
_ ~
n=o'
[ferc (_ ~) - ferc
V T
(n V+_l)J.
T
tf ~1- Lim
11o---i-OO
~
k =o
[ferc (_ "r-) - ferc (k -~)J,
y T V T
(2) t ~
O+, u ~ 1, para x> 0, y> O;
(3) x ~ O+, u ~ 0, para t > 0, y > O;
(4) y ~ O+, u ~ 0, para t > 0, x > O;
(5) Lim u(x, y, t) existe para t y Y fijas y positiva;;
,,-> 00
y
(16) Lim w(y, t) existe para t fija y positiva.
1I~ c<J
d 2g
(18) sg - 1 = dx2l
(24) v ( x, t) = fer(~).
2 vt
Por t anto la solución de nuestro problema original ( 1) a (6) es
( 25)
1. Demuéstrese que para la u de (25 ), O < u < 1, para toda x, y, t > ().
2. ,Supóngase que el punto con coordenadas (x, y, t ) está en el primer
..octante del espacio rectangular x, y, t; Y que este punto se acerca al
origen a lo largo de la curva
X2 = 4a 2 t,
y2 = 4a2t,
en la que a es positiva pero arbitraria. Demuéstrese que cuando
x, y, t ~ 0+ en la forma aescrita anteriormente, puede hacerse que u
se aproxime a cualquier número deseado que esté entre cero y la
unidad.
http://carlos2524.jimdo.com/
http://carlos2524.jimdo.com/
,
Indice analítico
Ecuación: Ecuación:
hipergeométrica, 381-382 con coeficientes variables, 15-122,417-
indicial, 349 422
con diferencia entera de raíces, con coeficientes constantes, ' 147-
365-368, 369-373 181, 222-227, 279-291
con diferencia no entera de raíces, métodos de la transformada de
350-353 Laplace, 226
con raíces iguales, 358-362, 396- métodos del operador diferencial,
397 165.-179, 287-289
integrales especiales, 299-234 reducción de orden, 287-291
lineal, 16 solución en series, 417-422
cambio de variables en, 289-290 variación de parámetros, 279-286
coeficientes indeterminados, 153- no lineal, 17
156 de orden dos O más, 321·-328
de orden uno y de mayor grado,
de orden n, 115-123
295-297, 299-318
de orden uno, 49-52
de orden uno y grado uno, 29-47,
de orden uno, forma estándar, 51
69-91
homogénea con coeficientes cons- indicaciones sobre la solución ge-
tantes, 133-146 neral de, 331 .
homogénea con coeficientes varia- resuelta mediante series de poten-
bles, 115-121 , 329-412 cia, 430-432
métodos de la transformada de punto regular singular de, 345
Laplace, 222-227 soluciones cerca de, 350-379/ 395-
métodos de series de potencia, 329-
~0 1
422 sistema de, 261-267
métodos operacionales, 165-181, térmica, 442, 489
222-227 bidimensional, 544-547
no homogénea con coeficientes en coordenadas esféricas, 500
constantes, 147-180 unidimensional, 442-448, 489-502
no homogénea con coeficientes va- validez de, 496-498
riables, 15-122, 417-422 variables canónicas, 538-540
punto regular singular de, 345 Electrocardiografía, 482-486
punto ordinario de, 333 E} ctrocardiograma, 482
punto ordinario de, soluciones cer- Eliminación de constantes arbitrarias,
ca de, 335-341, 391-393 17-25
punto singular de, 345 Eliminando la variable:
punto singular de, soluciones cerca dependiente, 309-312
de, 350-379, 395-402 independiente, 316-318
sistemas de, 261-267 Envolvente, 304
solución cerca de, 408-412 Esfera, conducción de calor en una,
variación de parámetros, 279-285 500-502
lineal homogénea, 115 Estrafoides:
con coeficientes constantes, 13~ - ecuación diferencial de, 27
146 Existencia de soluciones, 293-297
con coeficientes variables, 115, Extensión periódica, 464, 475, 479
329-412
lineal no homogénea, 115 Factor de amortiguamiento, 240, 249
cambio de variables en, 289-291 Factor integrante, 49, 69-76
coeficientes indeterminados, 153- para ecuaciones con coeficientes ho-
156 mogéneos, 77
http://carlos2524.jimdo.com/
552 fndice analítico
OTROS TÍTULOS
Introducción a la matemática
superior ,~ ,
Bush - Obreanu
Es un texto dirigido a los estudiantes p rsdniver-
sitarios y a los primeros grados de profesional.
Examina el problema del análisis de los argumen-
tos matemáticos, problema que pertenece en su
totalidad al campo de la lógica matemática o sim-
'bólica, esto es, el enfoque tradicional y el moder-
no, aunque se encamina más hacia 1& matemática
moderna.
Matemática finita
Louis O. Katlsoff - Albert J., Simone
Esta es una obra muy adaptable que se presta a
ser utilizada en cursos de varios niveles y de tipo
,diverso,.y8¡ q,;ie 'niá~~a debidamente las lecturas
.~" 9-ue los estuqtantesi' de primeros lírados pueden
, omitir sin perdér la unidad de su curs'o , No obs-
tante, su adaptabilidad no le resta en A.bsoluto
coherencia en sus exposiciones ni en'torpece
la secuencia de sus capítulos. Cuenta fcon un
... ........ extenso material práctico que ofrece numerosos
ejemplos, problemas y secciones enteras sobre
sus aplicaciones.