Fausto Ecuaciones PDF

Integracion de
ecuaciones diferenciales
ordinarias
Fausto Cervantes Ortiz

.
Universidad Autnoma de la Ciudad de Mxico
Integracin de ecuaciones
diferenciales ordinarias
Academia de Matemticas
Colegio de Ciencia y Tecnologa
ECUACIONES DIFERENCIALES ORDINARIAS

Ciclo Bsico Ingenieras
Integracin de ecuaciones diferenciales ordinarias
Primera edicin, 2008
D.R. Universidad Autnoma de la Ciudad de Mxico

Av. Divisin del Norte 906, Col. Narvarte Poniente,
Delegacin Benito Jurez, C.P. 03020, Mxico, D.F.
Academia de Matemticas del Colegio de Ciencia y Tecnologa de la UACM

Ciclo Bsico para Ingenieras
ISBN: 978-968-9259-45-9
Fotografa de la portada:
Vista satelital del huracn Katrina
sobre el Golfo de Mxico en agosto del 2005.
Tomada de en.wikipedia.org
Diseo de Portada: Aarn Aguilar

Diagramas del texto elaborados por el autor
Material educativo universitario de distribucin gratuita. Prohibida su venta.

Hecho e impreso en Mxico / Printed in Mexico
Correo electrnico:
bibliotecadelestudiante@hotmail.com
Introduccion
El libro que ahora llega a manos de los estudiantes, tiene su origen en las notas que el autor
preparo para impartir distintos cursos en la Universidad Autonoma de la Ciudad de Mexico.
Esta herramienta pretende ser un material intermedio entre los extensos tratados de Ecua-
ciones Diferenciales que encontramos en las bibliotecas y los breves formularios sin un sustento
conceptual. No es un Curso de Ecuaciones Diferenciales en el sentido tradicional. Busca ofrecer
al lector un enfoque diferente, mas practico. Por lo tanto, se plasman los metodos de integracion
sin demostraciones: en su lugar, se espera que cuando al estudiante de Ciencias e Ingenieras se le
presente una ecuacion, tenga todos los recursos necesarios para integrarla, y no necesariamente
deba demostrar que la solucion existe y esta acotada.
Sin embargo, no es suficiente con que el estudiante tenga las formulas; es necesario que
aprenda a usarlas correctamente, por lo cual se hace hincapie en los ejemplos. Ademas, es
fundamental que practique la resolucion de ejercicios, para lo que se incluye una buena cantidad
de ellos al final de cada seccion de estudio, con sus respectivas soluciones, permitiendo as que
el universitario mida su avance por s mismo. Finalmente, puesto que en las evaluaciones el
enfasis se hace en la habilidad para resolver ejercicios, es primordial trabajar la mayor cantidad
posible de ellos, durante el curso y antes de cada evaluacion.
A fin de que el aprendizaje sea significativo para los estudiantes, se incluyen ejercicios de apli-
cacion que permiten un acercamiento interdisciplinario a otras asignaturas de su carrera, aunque
ha de aceptarse que, cuando se avanza a mayores niveles de abstraccion en la Matematica, re-
sulta cada vez mas difcil hallar conexion con situaciones de la vida cotidiana: el futuro ingeniero
debera comprender que su vida diaria, si bien sera apasionante, tambien sera muy diferente a
la de otros profesionistas.
El autor desea agradecer a aquellos estudiantes y profesores que leyeron versiones anterio-
res, con lo que ayudaron a localizar inconsistencias en el texto y los ejercicios, ademas de
haber hecho sugerencias interesantes y muy utiles. En particular se agradece a Oscar Garca
y Humberto Mondragon (UIA). Sin embargo, las erratas que aun persistan, seran la absoluta
responsabilidad del autor y se agradecera a los lectores que se sirvan senalarlas, as como sus su-
gerencias y aportaciones. Para ello pueden contactar al autor a traves de la direccion electronica
fausto.cervantes@uacm.edu.mx , y/o en el cubculo E-256 del Plantel San Lorenzo Tezonco.
v
vi Introduccion
Muchas personas ayudaron a que este libro llegara a ser publicado despues de un largo
camino. En particular se agradece a Ana Beatriz Alonso por todo el apoyo editorial brindado.
Se espera que este material sea util a todos los lectores; pero sobre todo a los estudiantes de
Ingeniera de la UACM, ya que fue en ellos en quienes, en primer lugar, se penso al escribirlo.
Nada humano me es ajeno

San Lorenzo Tezonco, D. F.
Febrero de 2008
Indice
1. Ecuaciones de primer orden 1

1.1. Ecuaciones separables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2. Ecuaciones homogeneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3. Ecuaciones exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4. Factores integrantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.5. Ecuaciones lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.6. Ecuacion de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.7. Ecuacion de Riccati . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.8. Aplicaciones de modelado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.8.1. Trayectorias isogonales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.8.2. Movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.8.3. Circuitos electricos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
1.8.4. Ley de enfriamiento de Newton . . . . . . . . . . . . . . . . . . . . . . . 27
2. Ecuaciones de segundo orden 33

2.1. Ecuaciones lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.1.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.1.2. Reduccion de orden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.2. Coeficientes constantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.3. Ecuacion de Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.4. Ecuaciones no homogeneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.4.1. Metodo de los coeficientes indeterminados . . . . . . . . . . . . . . . . . 42
2.4.2. Metodo de variacion de parametros . . . . . . . . . . . . . . . . . . . . . 45
2.5. Aplicaciones de modelado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.6. Ecuaciones lineales de orden superior . . . . . . . . . . . . . . . . . . . . . . . . 56
3. Sistemas de ecuaciones lineales 57

3.1. Metodo de eliminacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.2. Metodos matriciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.3. Sistemas homogeneos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4. Sistemas no homogeneos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.4.1. Coeficientes indeterminados . . . . . . . . . . . . . . . . . . . . . . . . . 64
3.4.2. Variacion de parametros . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
viii Indice
4. Transformadas de Laplace 69
4.1. Calculo de la transformada de Laplace . . . . . . . . . . . . . . . . . . . . . . . 69
4.2. Algunas propiedades de la transformada de Laplace . . . . . . . . . . . . . . . . . 71
4.2.1. Traslacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.2. Funcion escalon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.3. Derivadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.3. Metodos para invertir la TL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
4.3.1. Fracciones parciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
4.3.2. Convolucion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.4. Solucion de ecuaciones con la TL . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.5. Funcion de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.6. Solucion de sistemas usando TL . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
5. Soluciones usando series de potencias 85

5.1. Solucion de ecuaciones alrededor de puntos ordinarios . . . . . . . . . . . . . . . 85
5.2. Solucion de ecuaciones alrededor de puntos singulares regulares . . . . . . . . . . 90
5.2.1. Races diferentes por un numero no entero . . . . . . . . . . . . . . . . . 91
5.2.2. Races diferentes por un numero entero . . . . . . . . . . . . . . . . . . . 96
5.2.3. Races iguales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
Tablas 107
Bibliografa 109
Captulo 1
Ecuaciones de primer orden
Una ecuacion diferencial es una ecuacion en la cual intervienen derivadas de funciones

desconocidas (incognitas). Es importante distinguir una ecuacion diferencial de una identidad
diferencial. Por ejemplo, sabemos que la igualdad:
d dg(x) df (x)
[f (x)g(x)] = f (x) + g(x)
dx dx dx
siempre se cumple, por lo tanto no es una ecuacion diferencial, sino una identidad diferencial;
una ecuacion solo sera valida para algunas funciones que deseamos encontrar. Resolver (o
integrar) una ecuacion diferencial es encontrar la funcion incognita a partir de conocer la
relacion dada entre sus derivadas y otras funciones. A menudo a tal funcion se le denomina
curva integral.
Las ecuaciones diferenciales se clasifican como ordinarias o parciales, dependiendo de si
intervienen en ellas derivadas de funciones de una sola o de varias variables independientes.
Tambien se distinguen por el orden de la derivada de mayor rango. En este captulo se estu-
dian las ecuaciones diferenciales ordinarias (solo hay una variable independiente) de primer
orden (solo hay derivadas de orden uno). Estas se clasifican de diferentes maneras, segun las
caractersticas que se indican en cada caso.
1.1. Ecuaciones separables

Una ecuacion se llama separable si posee la forma
M(x)dx + N(y)dy = 0. (1.1)

Algunas ecuaciones separables estan escritas en la forma
dy
= f (x, y), (1.2)
dx
pero es posible pasarlas a la forma (1.1) para resolverlas, lo cual se logra tratando a la derivada
como si fuera 1 un cociente de la diferencial de y entre la diferencial de x; y entonces, al pasar al
1
En realidad no lo es. Aunque la definicion de derivada involucra un cociente, al tomar el lmite obtenemos
una nueva funcion. La definicion de diferencial es lo que nos permite manipular la derivada como si fuera un
cociente de diferenciales.
2 1 Ecuaciones de primer orden
otro miembro de la ecuacion las funciones dependientes de x, tambien pasamos multiplicando

la diferencial de x.
Las ecuaciones separables se resuelven integrando en ambos terminos diferenciales con
respecto a la variable indicada por la diferencial. El resultado de ambas integraciones se iguala a
una constante arbitraria, que podemos considerar como la suma de las constantes de integracion
que resultan de las integraciones realizadas.
Ejemplo
Resolver la ecuacion
dy
= sen x
dx
Solucion
La ecuacion dada es equivalente a
dy = sen x dx,
que podemos integrar con respecto a cada variable independiente dada
Z Z
dy = sen x dx,
esto nos da
y = cos x + C.
Ejemplo
dy 1+y
=
dx 1x
Solucion
Esta ecuacion es equivalente a
dy dx
= ,
1+y 1x
que podemos integrar con respecto a cada variable independiente dada
dy dx
Z Z
= ,
1+y 1x
esto nos da
ln |1 + y| = ln |1 x| + ln C,
1.1 Ecuaciones separables 3
esto es,
(1 x)(1 + y) = C,
o bien
C
y= 1.
1x
Para eliminar la constante arbitraria en la solucion de una ecuacion diferencial de primer

orden, se necesita alguna condicion extra, llamada frecuentemente condicion inicial, que nos
permita determinar el valor de dicha constante.
Ejemplo
dy 3x2 + 4x + 2
= , y(0) = 1.
dx 2(y 1)
Solucion
Primero separamos variables
2(y 1)dy = (3x2 + 4x + 2)dx,

lo que al integrar se transforma en
y 2 2y = x3 + 2x2 + 2x + C;
sustituyendo la condicion inicial, que y = 1 cuando x = 0
3 = C,
con lo que finalmente
y 2 2y = x3 + 2x2 + 2x + 3.
Esta solucion se puede dejar en forma implcita, o despejar con la formula cuadratica.
Ejemplo
Resolver la ecuacion dada, sujeta a la condicion inicial indicada.
dy y cos x
= , y(0) = 1.
dx 1 + 2y 2
Solucion
Separando variables obtenemos
1 + 2y 2
dy = cos x dx,
y
lo que integrando da
ln y + y 2 = sen x + C;
sustituyendo la condicion inicial y = 1 cuando x = 0 encontramos
1 = C.
Con este valor de C obtenemos finalmente
ln y + y 2 = sen x + 1.
Esta solucion se debe dejar en forma implcita, ya que no se puede despejar con la formula
cuadratica por la presencia del termino ln y.
Ejercicios
Integrar las ecuaciones diferenciales siguientes
1. (1 + y 2 ) dx + (1 + x2 ) dy = 0 R: x + y = C(1 xy)
2. (1 + y 2 ) dx + xy dy = 0 R: x2 (1 + y 2 ) = C

1+x
3. (y 2 + xy 2 )y + x2 x2 y = 0 R: (x + y)(x y 2) + 2 ln 1y =C
4. (1 + y 2 ) dx = x dy R: y = tg ln Cx
p p
5. x 1 + y 2 + yy 1 + x2 = 0 R: 1 + x2 + 1 + y 2 = C
p p
6. x 1 y 2 dx + y 1 x2 dy = 0, y(0) = 1 R: 1 x2 + 1 y 2 = 1, y = 1
7. ey (1 + y ) = 1 R: ex = C(1 ey )
8. y ln y dx + x dy = 0, y(1) = 1 R: y = 1
9. y = ax+y , (a > 0, a 6= 1) R: ax + ay = C
10. 2x(1 + ey ) dx + ey (1 + x2 ) dy = 0 R: 1 + ey = C(1 + x2 )

x
11. (1 + ex )yy = ey , y(0) = 0 R: (1 + y)ey = ln 1+e
2 +1x
p
12. (1 + y 2 )(e2x dx ey dy) (1 + y) dy = 0 R: 12 e2x ey ln 1 + y 2 arc tg y = C
13. (xy 2 y 2 + x 1) dx + (x2 y 2xy + x2 + 2y 2x + 2) dy = 0 R: (x2 2x + 2)(y 2 + 1)e2 arc tg y = C
y = sen(x y) R: x = ctg yx

14. 2 + 4 +C
15. y = ax + by + c (a, b, c constantes) R: a + b(ax + by + c) = Kebx
(x + y)2 y = a2 R: x + y = a tg ya + C

16.
1.2 Ecuaciones homogeneas 5
y2 ey
17. (1 y)ey y + x ln x =0 R: y = ln | ln x| + C
2
p 1+y
18. (1 + y 2 ) dx = (y 1 + y 2 )(1 + x2 )3/2 dy R: ln = x +C
|y+ 1+y 2 | 1+x2
19. xy 2 (xy + y) = a2 R: 2x3 y 3 = 3a2 x2

1 1
20. (x2 y 2 + 1) dx + 2x2 dy = 0 (Sugerencia: z = xy) R: y = 1xy + 2 ln x = C
1
21. (x2 y 2 + 1)y + (xy 1)2 xy = 0 R: exy xy = Cy 2
22. (x2 y 3 + y + x 2) dx + (x3 y 2 + x) dy = 0 R: 3x2 12x + 2x3 y 3 + 6xy = C

x3 y3
23. (x6 2x5 + 2x4 y 3 + 4x2 y) dx + (xy 2 4x3 ) dy = 0 R: 3 x2 + 2x + 3x3 4 xy = C
(x+y)m (x+y)nm+1 (x+y)pm+1

24. y + 1 = (x+y)n +(x+y)p R: x = nm+1 + pm+1 +C
25. (ln x + y 3 ) dx 3xy 2 dy = 0 R: y 3 = Cx 1 ln x
26. (xy + 2xy ln2 y + y ln y) dx + (2x2 ln y + x) dy = 0 R: 2x2 + (2x ln y + 1)2 = C, x = 0

Cx
27. y xy = a(1 + x2 y ) R: y = a + ax+1

(a2 + y 2 ) dx + 2x ax x2 dy = 0, R: y = a tg xa 1
p
28. y(a) = 0
x+y xy
R: 2 sen x2 + ln tg y4 = C

y + sen

29. 2 = sen 2
1.2. Ecuaciones homogeneas

En esta seccion, llamar a una ecuacion homogenea quiere decir que al escribirla en la forma
dy
= f (x, y),
dx
la funcion f (x, y) es homogenea2 (de cualquier grado). En este caso, la ecuacion se resuelve
facilmente si se hace el cambio de variable v = xy . Con este cambio de variable, se tendra que
dy dv
y = vx, dx = v + x dx .
Ejemplo
dy y 2 + 2xy
= .
dx x2
Solucion
La ecuacion se puede escribir en la forma

2
Una funcion f (x, y) es homogenea de grado n si se cumple que f (x, y) = n f (x, y).
dy y 2 y
= +2 ,
dx x x
y al hacer el cambio de variable indicado se transforma en
dv
v+x = v 2 + 2v,
dx
dv
x = v 2 + v,
dx
dv dx
= ,
v2 +v x
dx dv dv dv
= = ,
x v(v + 1) v v+1
ln x = ln v ln(v + 1) ln C,
v
Cx = .
v+1
Regresando a las variables originales se tiene la solucion
y
v x
Cx = = y ,
v+1 x +1
y
Cx + Cy = ,
x

1
Cx = C y,
x
Cx Cx2
y= 1 = .
x C 1 Cx
Ejercicios
Resolver las siguientes ecuaciones homogeneas
1. 4x 3y + y (2y 3x) = 0 R: y 2 3xy + 2x2 = C

p
2. xy = y + y 2 x2 R: 2Cy = C 2 x2 + 1, y = x
3. 4x2 xy + y 2 + y (x2 xy + 4y 2 ) = 0 R: (x + y)(x3 + y 3 ) = C
4. 4x2 + xy 3y 2 + y (5x2 + 2xy + y 2 ) = 0 R: (y x)8 (y 2x)9 = C(y + 2x)5

2xy
5. y = 3x2 y 2 R: C(y 2 x2 ) = y 3
2
6. 2xy (x2 + y 2 ) = y(y 2 + 2x2 ) R: y 2 = Cxe xy2

p p y(y+ y2 x2 )
7. xy = y 2 x2 R: y + y 2 x2 = Cx3 e x2
1.3 Ecuaciones exactas 7
8. ax2 + 2bxy + cy 2 + y (bx2 + 2cxy + f y 2 ) = 0 R: ax3 + 3bx2 y + 3cxy 2 + f y 3 = C

9. (y 4 3x2 ) dy = xy dx R: x2 = y 4 + Cy 6
10. y 3 dx + 2(x2 xy 2 ) dy = 0 R: y 2 = x ln Cy 2
11. (y xy )2 = x2 + y 2 R: C 2 x2 = 1 + 2Cy, C 2 x2 = 2Cy
12. 3x + y 2 + y (x 1) = 0 R: (x 1)(3x + 2y 1) = C
2x+y
13. 2x + 2y 1 + y (x + y 2) = 0 R: x + y + 1 = Ce 3
14. (3y 7x + 7) dx (3x 7y 3) dy = 0 R: (x + y 1)5 (x y 1)2 = C

p p
15. (y + y x2 y 4 + 1) dx + 2x dy = 0 R: x2 y 4 + 1 = Cx2 y 2 1
16. 4xy 2 dx + (3x2 y 1) dy = 0 R: y(x2 y 1)2 = C
3
17. (x + y 3 ) dx + (3y 5 3xy 2 ) dy = 0 R: arc tg yx = 1
ln(x2 + y 6 ) + ln C
2
p p
18. 2(x2 y + x2 y 4 + 1) dx + x3 dy = 0 R: x2 (x2 y + 1 + x4 y 2 ) = C
19. (2x 4y) dx + (x + y 3) dy = 0 R: (y 2x + 3)3 = C(y x + 1)2
20. (x y + 3) dx + (3x 6y + 2) dy = 0 R: x + 3y ln |x 2y| = C
2x+2
21. (x y + 3) dx + (3x + y + 1) dy = 0 R: y = 1 x + Ce x+y1
22. (x + y) dx + (x + y 1) dy = 0 R: (x + 3y 1)2 + 2x = C
23. y cos x dx + (2y sen x) dy = 0 R: sen x + 2y ln |y| Cy = 0
y y y
24. (x y cos x) dx + x cos x dy = 0 R: x = Ce sen x
p
25. y 3 dy + 3y 2 x dx + 2x3 dx = 0 R: C x2 + y 2 = 2x2 + y 2
x
26. y dx + (2 xy x) dy = 0 R: ye y = C
1.3. Ecuaciones exactas

Una ecuacion de la forma
M(x, y)dx + N(x, y)dy = 0

es exacta si se cumple que:
M N
= .
y x
Con esto, al integrar M con respecto a x y N con respecto a y se obtiene una solucion
implcita de la forma
f (x, y) = C.
Es importante hacer notar que la solucion f (x, y) no es la suma de las integrales halladas,
sino que es la expresion mas simple que quede al eliminar los terminos que se repiten.
Ejemplo
2xy 3 dx + 3x2 y 2 dy = 0.
Solucion
La ecuacion es exacta, ya que

(2xy 3 ) = 6xy 2 = (3x2 y 2 ),
y x
entonces, al integrar:
Z Z
3 3 2
2xy dx = y x , 3x2 y 2 dy = x2 y 3 ,
por lo tanto, la solucion es
x2 y 3 = C.
Ejemplo
(y cos x + 2xey )dx + (sen x + x2 ey + 2)dy = 0.
Solucion
La ecuacion es exacta, puesto que
M N
= cos x + 2xey = .
y x
Al integrar obtenemos
Z Z
M dx = (y cos x + 2xey )dx = y sen x + x2 ey
Z Z
N dy = (sen x + x2 ey + 2)dy = y sen x + x2 ey + 2y.
De aqu, la solucion es
y sen x + x2 ey + 2y = C.
1.4 Factores integrantes 9
1.4. Factores integrantes

A veces se tiene alguna ecuacion de la forma
M(x, y)dx + N(x, y)dy = 0 (1.3)

en la que no se cumple que
M N
= .
y x
A pesar de esto, es posible hallar una solucion para esta ecuacion si se puede encontrar
alguna funcion (llamada factor integrante) que al multiplicarla por la ecuacion 1.3 nos de una
ecuacion exacta. Esto es, si en la ecuacion
(x, y)(M(x, y)dx + N(x, y)dy = 0)

se cumple que
(M) (N)
= ,
y x
entonces la funcion (x, y) es un factor integrante de la ecuacion dada.
Ejemplo
(y 2 + xy)dx x2 dy = 0.
Solucion
Esta ecuacion no es exacta, ya que
M
= 2y + x,
y
mientras que
N
= 2x.
x
Pero supongamos que existe un factor integrante de la forma
(x, y) = xm y n ,
entonces, al sustituir en la ecuacion original, tendremos que
(xm y n+2 + xm+1 y n+1 )dx xm+2 y n dy = 0

y las derivadas parciales son
M
= (n + 2)xm y n+1 + (n + 1)xm+1 y n ,
y
N
= (m + 2)xm+1 y n .
x
O sea que para que la ecuacion sea exacta es necesario que
n + 1 = m 2, n+2=0
lo cual nos da, al resolver el sistema de ecuaciones
1
= .
xy 2
Entonces, al sustituir el factor integrante en la ecuacion, esta se transforma en

1 1 x
+ dx dy = 0,
x y y2
cuya solucion es
x
ln x + = C.
y
A veces el factor integrante solo depende de x, en cuyo caso se puede hallar resolviendo la
ecuacion
M
d y
N
x
= .
dx N
Si en cambio el factor integrante depende solo de y, se puede hallar resolviendo
Z N M !
x y
= exp dy .
M
Si el factor integrante es funcion de x2 + y 2 , se encuentra con

Z N M !
1 x y
= C exp d(x2 + y 2 ) .
2 My Nx
En general un factor integrante satisface la ecuacion:
M N
M + =N +
y y x x
sin embargo, si suponemos que depende de una forma especfica (por ejemplo, x y, x2 y 2 ,
etc), se puede simplificar la ecuacion y su solucion sera mas facil de encontrar.
1.5 Ecuaciones lineales 11
Ejercicios
1. (2x + 3) dx + (2y 2) dy = 0 R: x2 + 3x + y 2 2y = c
2. (3x2 2xy + 2) dx + (6y 2 x2 + 3) dy = 0 R: x3 x2 y + 2x + 2y 3 + 3y = c
3. (2xy 2 + 2y) dx + (2x2 y + 2x) dy = 0 R: x2 y 2 + 2xy = c
4. (ex sen y 2y sen x) dx + (ex cos y + 2 cos x) dy = 0 R: ex sen y + 2y cos x = c
5. (yexy cos 2x 2exy sen 2x + 2x) dx + (xexy cos 2x 3) dy = 0 R: exy cos 2x + x2 3y = c
6. ( xy + 6x) dx + (ln x 2) dy = 0 R: y ln x + 3x2 2y = c

x dx y dy
7. (x2 +y 2 )3/2
+ (x2 +y 2 )3/2
=0 R: x2 + y 2 = c
c3x
8. (2xy + 3) dx + (x2 1) dy = 0 R: y = x2 1
9. (cos x cos y + 2x) dx (sen x sen y + 2y) dy = 0 R: y = sen x cos y + x2 y 2 = c
10. (1 + ln y) dx + ( xy ) dy = 0 R: x ln y + x = c
11. cos x dy (y sen x ex ) dx = 0 R: y = (c ex ) sen x
12. ( 1+xy2 y2 + 2x) dx + ( 1+xx2 y2 2y) dy = 0 R: x2 y 2 + arc tg(xy) = c
13. (xy 2 + 3x2 y) dx + (x + y)x2 dy = 0 R: x2 y 2 + 2x3 y = c
14. (ye2xy + x) dx + xe2xy dy = 0 R: e2xy + x2 = c
15. (3x2 y + 2xy + y 3 ) dx + (x2 + y 2 ) dy = 0, = e3x R: e3x (3x2 + y 3 ) = c
16. y = e2x + y 1, = ex R: y = cex + 1 + e2x

e2y
17. y dx + (2xy e2y ) dy = 0, = y R: xe2y ln y = c
18. ex dx + (ex ctg y + 2y csc y) dy = 0, = sen y R: ex sen y + y 2 = c

3
19. 4 xy2 + y3 dx + 3 yx2 dy = 0 R: x4 + 3xy + y 4 = c
2
20. 3x + y6 + xy + 3 xy y = 0 R: x3 + 3x2 + y 3 = c
1.5. Ecuaciones lineales

La ecuacion lineal de primer orden mas general es
y + p(x)y = g(x).
R
p(x)dx
Esta ecuacion tiene como factor integrante a = e . Con esto, la solucion general es:
"Z #
R R
p(x)dx p(x)dx
y = e e g(x)dx + C
Ejemplo
Resolver la ecuacion dada, con la condicion inicial especificada.
y 2xy = x, y(0) = 1.
Solucion
Primero calculamos la integral de p(x)

Z Z
p(x)dx = 2xdx = x2
con lo que la solucion sera de la forma

Z
x2 x2
y=e e xdx + C .
Para resolver la integral procedemos como sigue. Sea u = x2 , entonces, du = 2xdx con lo cual
1 1
Z Z
2 2
ex xdx = eu du = ex .
2 2
O sea que
1 2
y = + Cex .
2
Evaluando en cero para cumplir con la condicion inicial impuesta
1
y(0) = + C = 1,
2
de donde obtenemos que
3
C= ,
2
y finalmente, la solucion es
1 3 2
y = + ex .
2 2
Ejemplo
Resolver el problema
y 2xy = 1, y(0) = 1
Solucion
Aqu tenemos que:

Z
p(x)dx = x2
Z
2
y = ex ex dx + C
Z
x2 2 2
y=e ex dx + Cex .
Si en la formula que nos da la solucion de la ecuacion lineal general las funciones tienen
discontinuidades esenciales, la solucion tendra discontinuidades en los mismos puntos y solo
en ellos. En todos los demas es continua, por lo que se obtiene informacion sobre la solucion
con solo identificar los puntos de discontinuidades de p(x) y g(x). En algunos casos se pueden
elegir las constantes en la solucion de tal modo que se eliminen las discontinuidades.
Ejemplo
xy + 2y = 4x2 , y(1) = 2.
Solucion
El ejercicio planteado es equivalente a
y + 2y/x = 4x, p(x) = 2/x, g(x) = 4x,

donde vemos que p(x) es discontinua en x = 0, por lo cual la solucion sera discontinua all en
algunas condiciones. Primero resolvemos la integral de p(x)
2
Z
dx = 2 ln x = ln x2
x
Z
ln x2 ln x2
y=e e (4x)dx + C
Z
1 C
= 2 4x dx + C = x2 +
3
x x2
y(1) = 2 = 12 + C
C=1
1
y(x) = x2 +
x2
que es discontinua en x = 0, por lo cual la solucion sera discontinua all. Pero si se tuviera que
y(1) = 1, se tendra que C = 0, y la solucion sera
y(x) = x2
que es continua en todos los reales.
Si p(x) o g(x) tienen discontinuidades de salto, por ejemplo en x0 , se resuelve por separado
cuando x > x0 y cuando x < x0 . Posteriormente se eligen las constantes de forma que y sea
continua en x0 .
Ejemplo
y + 2y = g(x), y(0) = 0,
con:

1 0x1
g(x) =
0 x>1
Solucion
Aqu x0 = 1, entonces, para x < x0 :
y + 2y = 1
Z
2dx = 2x
Z
2x 2x
y=e e dx + C1 =

1 2x
= e2x e dx + C1 =
2
1
+ C1 e2x
2
Ahora, para x > x0 :
y + 2y = 0
Z
y = e2x e2x (0)dx + C2 = C2 e2x .
Para que haya continuidad en x0 = 1, procedemos del siguiente modo.

Como y(0) = 0:
1
y(0) = + C1 e2(0) = 0
2
1
C1 =
2
Tambien
1
+ C1 e2(1) = C2 e2(1) .
2
Entonces
1
2 21 e2 e2 1
C2 = =
e2 2 2
y con esto
e2 1

1 1
y = e2x + 2 e2x
2 2 2
Ejercicios
x 1
1. y + 3y = x + e2x R: y = ce3x + 3 9 + e2x
2x
2. y 2y = x2 e2x R: y = ce2x + x3 e3
y + y = xex + 1 R: y = cex + 1 + x2 e 2
x
3.
c 3 cos 2x 3 sen 2x
4. y + y/x = 3 cos 2x R: y = x + 4x + 2
5. y y = 2ex R: y = cex + 2xex

cx cos x+sen x
6. xy + 2y = sen x R: y = x2
2 2 2
7. y + 2xy = 2xex R: y = x2 ex + cex
1 arc tg x+c
8. (1 + x2 )y + 4xy = (1+x2 )2 R: y = (1+x2 )2
9. y y = 2xe2x , y(0) = 1 R: y = 3ex + 2(x 1)e2x
y + 2y = xe2x , R: y = (x2 1) e 2
2x
10. y(1) = 0
1 3x4 4x3 +6x2 +1
11. xy + 2y = x2 x + 1, y(1) = 2 R: y = 12x2
2y cos x sen x
12. y + x = x2 , y() = 0 R: y = x2
13. y 2y = e2x , y(0) = 2 R: y = (x + 2)e2x

c+sen x
14. y + y/x = sen x R: y = x cos x
15. y + y tg x = x sen 2x R: y = (c 2x cos x + 2 sen x) cos x
c+(x1)ex
16. xy + 2y = ex R: y = x2
1 3x2 4x+6+x2
17. xy + 2y = x2 x + 1, y(1) = 2 R: y = 12
ex +1e
18. xy + y = ex , y(1) = 1 R: y = x
19. y + y cot x = 2 csc x, y( 2 ) = 1 R: y = 2x+1

sen x
1 sen xx cos x
20. xy + 2y = sen x, y() = R: y = x2
21. y + x2 y = 1
x2 R: y = c
x2
e2x

2 0x1 0x1
22. y + p(x)y = 0, y(0) = 1, con: p(x) = R: y = (x+1)
1 x>1 e x>1
1.6. Ecuacion de Bernoulli

La ecuacion:
y + p(x)y = q(x)y n ,
llamada ecuacion de Bernoulli, se convierte en una ecuacion lineal al hacer la sustitucion
v = y 1n .
Ejemplo
x2 y + 2xy y 3 = 0
Solucion
Dividiendo entre x2 :
2 1
y + y = 2 y3,
x x
ahora hacemos
v = y 2
v = 2y 3 y
y v 1
3
= , v=
y 2 y2
y 2 1 1
+ = 2
y3 x y2 x
v 2 1
+ v= 2
2 x x
4 2
v v= 2
x x
4 1
Z
dx = 4 ln x = ln 4
x x
Z
1 1 2
v = e ln x4 eln x4 dx + C =
x
Z
2 2
v = x4 6 dx + C = + Cx4
x 5x
1 2 2 + 5Cx5
2
= + Cx4 =
y 5x 5x
r
5x
y=
2 + 5Cx5
1.7 Ecuacion de Riccati 17
1.7. Ecuacion de Riccati

La ecuacion de Riccati es
dy
= q1 (x) + q2 (x)y + q3 (x)y 2 .
dx
Si se conoce una solucion particular y1 (x), se puede obtener una mas general con la susti-
tucion
1
y = y1 (x) + .
v(x)
La funcion v(x) satisface la ecuacion
dv
= (q2 + 2q3 y1 )v q3 .
dx
Ejemplo
y = 1 + x2 2xy + y 2 , y1 (x) = x
Solucion
Dada la solucion particular y1 , la ecuacion a resolver para v es
dv
= (2x + 2x)v 1 = 1,
dx
que tiene la solucion
v = x + C,
y la solucion general es
1
y =x .
xC
Ejercicios
Resolver las siguientes ecuaciones
ex
1. y = y y 2 , > 0, > 0. R: y = ex +C
x
2. y = y y 3 , > 0, > 0 R: y = ee
x +C
iex
3. y = y f (x)y 3 , >0 R: y = R
C2 f (x)e2x dx
1
4. y + y 2 sen x = 0 R: y + cos x = C
x2
5. y = 1 + x + y 2 + xy 2 R: arc tg y x 2 =C
6. y = (cos2 x)(cos2 2y) R: 2 tg 2y 2x sen 2x = C

7. (2x + 3)dx + (2y 2)dy = 0 R: y = x2 + 3x + y 2 2y = C
8. ( xy + 6x)dx + (ln x 2)dy = 0, x>0 R: y ln x + 3x2 2y = C
9. (ye2xy + x)dx + xe2xy dy R: e2xy + x2 = C
1
10. x2 y 3 + x(1 + y 2 )y = 0 R: x2 + 2 ln |y| y2 =C
11. ydx + (2x yey )dy = 0 R: xy 2 (y 2 2y + 2)ey = C

12. ydx + xy sen y dy = 0 R: xy + y cos y sen y = C
p
13. xy = x2 y 2 + y R: y = x sen(ln Cx)

x2 +y 2
14. xdx + ydy + xdy ydx = 0 R: arc tg xy + ln x2 =C
2
dy y 2 xy

15. dx = R: =C
x+ x +y 2 x4 + x +y 2
2
dy y4x
16. dx = xy R: xy + 2x2 = C
x+y
17. y = x R: y = x ln x + Cx
x2 +3y 2

18. y = 2xy R: y = Cx3 x2
y
19. y = x12 x + y2, y1 (x) = 1
x R: y = 1
x + 2x
Cx2
dy 2 cos2 xsen2 x+y 2 1

20. dx = 2 cos x , y1 (x) = sen x R: y = sen x + C cos x 12 sen x
2
21. y + xy = x3 y 3 R: y 2 (x2 + 1 + Cex ) = 1

22. (1 x2 )y xy 9xy 2 = 0 R: (C 1 x2 9)y = 1
23. 3y 2 y 2y 3 x 1 = 0 R: 4y 3 = Ce2x 2(x + 1) 1
y2 y2
h i
24. y (x2 y 3 + xy) = 1 R: x (2 y 2 )e 2 + C = e 2
25. (y ln x 2)ydx = xdy R: y(Cx + ln x + 1) = 1

2 tg x+sec x
26. y y cos x = y cos x(1 sen x) R: y = sen x+C
1.8. Aplicaciones de modelado

Para resolver problemas usando las ecuaciones diferenciales es necesario traducir la situacion
fsica en terminos matematicos. Suele llevarse a cabo esto al establecer hipotesis, acerca de lo
que esta sucediendo, que sean coherentes con los fenomenos observados. Por ejemplo, se ha
observado que los materiales radiactivos decaen con una rapidez proporcional a la cantidad de
material presente, que el calor pasa de un cuerpo caliente hacia uno mas fro con una rapidez
proporcional a la diferencia de temperaturas, que los objetos se mueven segun las leyes de
Newton del movimiento, etc. Cada una de estas proposiciones comprende una razon de cambio
(derivada) y, en consecuencia, al expresarse matematicamente, toma la forma de una ecuacion
diferencial.
Una vez que el fenomeno se ha modelado matematicamente, el problema se reduce a resolver
una o mas ecuaciones diferenciales.
1.8 Aplicaciones de modelado 19
1.8.1. Trayectorias isogonales

En una familia de curvas dependiente de un parametro C, digamos y = (x, C), las lneas
de la forma y = g(x, K) que cortan a todas las curvas de tal familia en un mismo angulo se
llaman trayectorias isogonales. En el caso particular de que el angulo de corte sea recto, se
denominan trayectorias ortogonales. Para que se cumpla lo dicho anteriormente respecto al
angulo de corte, las pendientes de las rectas tangentes a las familias de curvas de la familia y
las trayectorias isogonales en un angulo deben cumplir la relacion
dyt
dy tg
= dx , (1.4)
dx 1 + tg dy
dx
t
donde el subndice t indica que se toma la derivada de las curvas tangentes. La ecuacion anterior
proviene de la identidad trigonometrica para la tangente de la resta de dos angulos. En el caso
particular de las trayectorias ortogonales, la relacion 1.4 se reduce a
dy 1
= dyt . (1.5)
dx dx
Ejemplo
Encontrar las trayectorias ortogonales a la familia de curvas
y = cx2 .
Solucion
La derivada de la familia de curvas (parabolas) anterior es
dy
= 2cx,
dx
que al combinar con la ecuacion de la familia para eliminar el parametro c nos da
dy y 2y
= 2x 2 = .
dx x x
Aplicando la condicion 1.5, obtenemos la ecuacion
dy x
= ,
dx 2y
que se resuelve muy facilmente por separacion de variables, dando como solucion
x2
+ y 2 = k.
2
La solucion encontrada define una familia de elipses. En la figura 1.1 se ven ambas familias de
curvas.
Figura 1.1: Trayectorias ortogonales
Ejemplo
Hallar las trayectorias isogonales de la familia de rectas
y = Cx
que cortan las lneas en un angulo de 45o .
Solucion
Para la familia dada se tiene que
dy
= C,
dx
que al eliminar el parametro C (despejando de la ecuacion original) se convierte en
dy y
= .
dx x
De la ecuacion 1.4 encontramos que la ecuacion diferencial a resolver es
dy
dx 1 y
dy
= .
1+ dx
x
Despejando la derivada de la ecuacion anterior obtenemos la ecuacion
y
dy 1+ x
= y .
dx 1 x
Esta es una ecuacion homogenea. Al resolverla por la tecnica estudiada anteriormente se obtiene
la solucion
p y
ln x2 + y 2 = arc tg + k.
x
La ecuacion anterior define una familia de espirales logartmicas, cuya grafica se da en la figura
1.2.
Figura 1.2: Trayectorias isogonales
Ejercicios
Encontrar las trayectorias ortogonales a las familias de curvas dadas a continuacion
1. y = cx3 R: x2 + 3y 2 = k
y 2 = 4(x c)
x
2. R: y = ke 2
k
3. x2 y 2 = c R: y = x
4. x2 + y 2 = 2cx R: y = k(x2 + y 2 )
x3
5. y2 = 2cx R: (x2 + y 2 )2 = k(y 2 + 2x2 )
6. (x2 + y 2 )2 = c2 (x2 y 2 ) R: (x2 + y 2 )2 = kxy
7. y 2 + 2cx = c2 R: y 2 2kx = k 2
8. y = cx3 R: x2 + 3y 2 = k
9. y = 5ecx R: 2x + 5ay 2 = k
10. cos y = cex R: sen y = kex
11. x2 + 12 y 2 = c2 R: y 2 = 2kx
12. y = axn R: x2 + ny 2 = c
13. y = aex R: 2x + ay 2 = c
14. x2 13 y 2 = a2 R: xy 3 = c
Hallar las trayectorias isogonales en el angulo dado para las familias de curvas siguientes
1. x 2y = c, = 45o R: y 3x = k

2. x2 = 2c(y 3x), = 60o R: y 2 = k(x 3y)

6 2yx
3. y 2 = 4cx, = 45o R: y 2 xy + 2x2 = k exp 7 arc tg
7x

= 30o R: x2 + y 2 = k exp 2 3 arc tg xy

4. y = cx,
= 45o R: x2 + y 2 = k exp 2 arc tg xy

5. y = cx,
x2 + y 2 + arc tg xy = k
p
6. x2 + y 2 = c2 , = 45 R: ln
1.8.2. Movimiento
La segunda ley de Newton establece que cuando una fuerza actua sobre una partcula, sufre
una aceleracion proporcional a la fuerza que se aplica, e inversamente proporcional a la masa
de tal partcula. Esto se escribe como
F = ma, (1.6)
donde F es la fuerza total aplicada (esto es, la suma de todas las fuerzas que actuan sobre la
partcula), m es la masa de la partcula y a es la aceleracion que experimenta. Si recordamos
que a = dv(t)
dt
y que v(t) = dx(t)
dt
, con v(t) la velocidad de la partcula en el instante t y x(t) la
posicion de la partcula en el mismo instante, tendremos las ecuaciones
dv(t) F
= (1.7)
dt m
y
dx(t)
= v(t), (1.8)
dt
que nos permitiran determinar el estado de movimiento de la partcula en un instante dado.
A continuacion se tratan un par de casos de movimiento, en donde se aplican las leyes de
Newton para determinar la cinematica de las partculas involucradas.
Movimiento con resistencia

Generalmente, en los problemas de cada libre, se supone que la unica fuerza que actua
sobre el objeto que cae es la gravedad. Pero cuando vemos caer a un paracaidista, despues de
un tiempo ya no lo hace con aceleracion, sino que cae con velocidad constante. Esto es porque
hay una fuerza causada por la friccion viscosa del aire.
En este caso, al aplicar la segunda ley de Newton, al peso del sistema hay que agregar una
fuerza proporcional al cuadrado de la velocidad de cada, que es la ley que sigue dicha fuerza
viscosa. Con esto, la fuerza total durante la cada libre esta dada por
F = mg + cv 2 , (1.9)
siendo c una constante que depende de diversos factores como tamano y forma del objeto,
medio en que se desarrolla la cada (aire o agua, por ejemplo), etc.
Al alcanzarse el equilibrio entre la fuerza gravitacional y la fuerza viscosa se llega a la
llamada velocidad terminal. Como en este caso F = 0, se tiene que
r
mg
vt = . (1.10)
c
Despues de esto, la velocidad de cada es constante.
Ejemplo
Un paracaidista cuya masa es de 70 kg cae desde una altitud de 1400 m. a) Encontrar la ley de
movimiento. b) Hallar su velocidad terminal. Tomar c =0.25 kg/m.
Solucion
La ecuacion de movimiento es
ma = mg + cv 2 ,
o, equivalentemente
dv
m = mg + cv 2 .
dt
Para integrar esta ecuacion, primero dividimos entre c, y usando la ecuacion 1.10 tenemos la
ecuacion equivalente
vt2 dv
= vt2 + v 2 .
g dt
Separando variables en esta ecuacion obtenemos (al integrar entre valores indeterminados de v y
t) las integrales
v t
dv g dt
Z Z
= ,
0 v vt2
2
0 vt2
que al resolverlas nos dan el resultado

1 vt + v g
ln = t.
2vt vt v vt2
Despejando v de la ecuacion anterior obtenemos
2gt
1 e vt
v = vt 2gt .
1 + e vt
Sustituyendo en esta ultima ecuacion los valores de m y c dados, hallamos
1 e0.03742t
v = 52.38 ,
1 + e0.3742t
en unidades de m/s.
Para la velocidad terminal calculamos el lmite
lm v = vt = 52.38 m/s
t
El signo menos en la anterior igualdad indica que la velocidad esta dirigida hacia abajo.
Movimiento con masa variable

Cuando un cuerpo cuya masa vara conforme pasa el tiempo (como un cohete), la segunda
ley de Newton tiene el aspecto
dv dM
M(t) =F +u , (1.11)
dt dt
siendo dM/dt la razon a la que sale (o entra) la masa del sistema bajo estudio, u la velocidad
a la que sale (o entra) masa (medida desde el mismo sistema), y F las fuerzas externas que
actuan sobre el. As, el movimiento de dicho cuerpo se describe al resolver la ecuacion (1.11).
Ejemplo
Un cohete tiene una masa inicial m0 = 25000 kg, de los cuales, 20000 kg son de combustible. La
velocidad de expulsion de los gases es de u = 400 m/s, durante lo cual hay perdida de a = 1000
kg de combustible. El cohete parte desde el reposo en la Tierra y viaja verticalmente hacia arriba,
siendo la unica fuerza externa su peso. Encontrar la velocidad del cohete despues de 15, 20 y 30
segundos. Hallar la altura que alcanza cuando se ha quemado la mitad del combustible.
Solucion
Aplicamos la ecuacion (1.11) a este problema, notando que M (t) = M0 at, tenemos que la
ecuacion para el cohete es
dv
(M0 at) = (M0 at)g ua,
dt
que para resolver dividimos entre (M0 at) y multiplicamos por dt, quedando
dt
dv = g dt ua ,
M0 at
as que integrando tenemos
v t t
dt M0
Z Z Z
dv = g dt ua = gt + u ln .
0 0 0 M0 at M0 at
Sustituyendo los valores dados, tenemos que v(t = 15) = 220 m/s, v(t = 25) = 448 m/s y
v(t = 30) = 350 m/s.
Para encontrar la altura, recordamos que v(t) = dy/dt, as que la altura se encontrara integrando
la velocidad
Z t Z t Z t
y(t) = gt dt + u ln(M0 ) dt ln(M0 at) dt,
0 0 0
lo cual nos da la ley de movimiento
1 2 u
y(t) = gt + u ln M0 t + {(M0 at)[ln(M0 at) 1] M0 [ln M0 1]}.
2 a
Sabiendo que M = 15000 kg se alcanza cuando t = 10 s, encontramos que h = 445 m.
Ejercicios
1. Una masa de 25 g cae desde el reposo por accion de la gravedad. Encontrar la distancia viajada y la
velocidad alcanzada 3 segundos despues de iniciar su movimiento.
R: 4.41 cm, 2.94 cm/s
2. Se lanza hacia arriba una masa de 200 g, con una velocidad de 2450 cm/s.
a) Encontrar las distancias desde el punto de partida y las velocidades alcanzadas 2 s y 4 s despues
de iniciar el movimiento
b) Encontrar el punto mas alto que alcanza y el tiempo requerido
c) Cuales son las distancias totales recorridas despues de 2 s y 4 s?
R: a) 2.94 cm, 490 cm/s ; 1.96 cm, -1.47 cm/s b) 3062.5 cm, 2.5 s c) 2.94 cm, 4.165 cm
3. Una pequena gota de aceite, de 0.2 g de masa, cae en el aire desde el reposo. Para una velocidad de 40
cm/s, la fuerza debida a la resistencia del aire es de 160 dinas. Suponiendo que la fuerza de resistencia
es proporcional a la velocidad,
a) encontrar la velocidad y la distancia como funciones del tiempo

b) encontrar la velocidad lmite
R: a) v = 49(1 e20t ), x = 49t + 2.45e20t 2.45 b) 49 cm/s

4. Un paracaidista que pesa 192 libras tiene una velocidad lmite de 16 ft/s cuando cae en el aire, el cual
genera una resistencia proporcional al cuadrado de la velocidad instantanea. Si cae desde el reposo,
a) encontrar la velocidad del paracaidista despues de 1 s

b) en que momento alcanza una velocidad de 15 ft/s?
R: a) 15.4 ft/s b) 0.86 s

5. Una partcula se mueve a lo largo del eje x estimulada solo por una fuerza opuesta proporcional a su
velocidad instantanea. La partcula empieza en el origen con una velocidad de 10 ft/s, la cual se reduce
a 5 ft/s despues de moverse 2.5 ft. Encontrar la velocidad cuando esta a 4 ft del origen.
R: 2 ft/s
1.8.3. Circuitos electricos

Cuando se tienen componentes electricas en un circuito al cual se le suministra una fuerza
electromotriz (fem), el voltaje suministrado disminuira al pasar por cada componente. Las
relaciones entre las cadas de voltajes, las cargas y las corrientes involucradas en sus elementos
son:
1. En una malla cerrada, la suma de cadas de voltaje es igual a la fem (E) suministrada.
2. En un nodo, la suma algebraica de las corrientes es cero.
3. La cada de voltaje de una resistencia es VR = iR, con i la corriente y R la resistencia

de esa componente.
Q
4. La cada de voltaje en un capacitor es VC = C
, con Q la carga acumulada y C la
capacitancia de la componente.
di
5. La cada de voltaje de una bobina (inductor) es VL = L dt , con L la inductancia de la
componente.
Ejemplo
Un generador con una fem de 100 V se conecta en serie con una resistencia de 10 y una bobina
de 2 H. Si en t = 0 se cierra el circuito, encontrar la funcion que da la corriente como funcion del
tiempo.
Solucion
Como el circuito esta en serie, solo hay una malla, donde la suma de voltajes es igual a la fem, lo
cual nos da la ecuacion
VR + VL = E,
o bien
di
iR + L = E.
dt
Sustituyendo los valores dados obtenemos la ecuacion
di
2 + 10i = 100,
dt
que es equivalente a la ecuacion

di
+ 5i = 50.
dt
Como en t = 0 se cierra el circuito, la condicion inicial es i(0) = 0.
Resolviendo esta ecuacion (que es lineal) obtenemos la solucion general
i(t) = 10 + Ce5t ,
que al sustituir la condicion inicial nos da la solucion particular
i(t) = 10(1 e5t ).
Ejercicios
1. En un capacitor, la carga electrica sale a una tasa que es proporcional a la carga instantanea del capacitor.
Inicialmente la carga es de 5 Coulomb, y en 20 minutos sale un tercio de la carga inicial. En cuanto
tiempo quedara solo un Coulomb en el capacitor?
R: 79 min
2. Un circuito consta de un inductor de 2 H en serie con una resistencia de 40 . En t = 0 (cuando la
corriente vale cero) se aplica al circuito una fem dada por 100 sen 10t. Encontrar el valor de la corriente
para cualquier instante.
R: i(t) = 2 sen 10t cos 10t + e20t
3. En un circuito hay un capacitor de 5 mF en serie con una resistencia de 25 y una fem 50 cos 6t. Si la
carga en el capacitor es cero al inicio, determinar la carga y la corriente en cualquier instante.
R: Q(t) = 0.16 cos 6t + 0.12 sen 6t 0.16e8t, i(t) = 0.72 cos 6t 0.69 sen 6t + 1.28e8t
4. Una resistencia de 4 y un inductor de 1 H se conectan en serie, suministrandoles un voltaje 100e4t cos 50t.
Encontrar i(t) si i(0) = 0.
R: i(t) = 2e4t sen 50t
5. Se conecta una resistencia de 20 con un capacitor de 0.01 F y una fem 40e3t + 20e6t . Si Q(0) = 0,
hallar el voltaje maximo en el capacitor.
R: 0.25 C
1.8.4. Ley de enfriamiento de Newton

Cuando un cuerpo (o sustancia) esta a una temperatura mayor que la del ambiente, el calor
se disipa al medio hasta que el cuerpo quede a la misma temperatura. La rapidez con que esto
ocurre es mayor entre mas grande sea la diferencia de temperaturas, esto es, la velocidad a que
se enfra un cuerpo es proporcional a la diferencia de temperaturas entre el cuerpo y el medio
que lo rodea. Esto se escribe como
dT (t)
= k [T (t) Ta ] , (1.12)
dt
siendo Ta la temperatura ambiente. El enunciado anterior constituye la ley de enfriamiento
de Newton. Cabe senalar que aunque se le llama ley de enfriamiento, tambien funciona para
describir el calentamiento de un cuerpo que esta mas fro que el medio. Tambien hay que
aclarar que dicha ley supone que la temperatura ambiente es constante, lo cual solo sucede para
intervalos de tiempo de pocas horas, por lo cual se debe tener cuidado durante su aplicacion.
Ejemplo
Supongase que en el instante t = 0 se descubre un cadaver y que su temperatura es 85 F. Se
supone que en el instante del deceso td la temperatura del cuerpo Td tena el valor normal de
98 F. Dos horas mas tarde su temperatura es de 74 F. La temperatura ambiente es de 68 F,
encontrar el tiempo de deceso.
Solucion
La ecuacion a resolver es
dT (t)
= k [T (t) 68] ,
dt
sujeta a las condiciones
T (0) = T0 = 98, T (td ) = Td = 85, T (2 + td ) = 74.
Por separacion de variables se encuentra que
T (t) = 68 + 30ekt .
De las condiciones dadas encontramos el valor de k como

1 74 68
k = ln 0.5207h
2 85 68
y finalmente,

1 98.6 68
td = ln 1.29h.
0.5207 85 68
De donde se concluye que el cuerpo se descubrio aproximadamente 1 hr. 8 min. despues del
fallecimiento.
Ejemplo
Un caluroso da de verano (la temperatura ambiente es de 30 C), Gonzalo llega al bar y pide
una cerveza. Mientras tanto, llega su amigo Gabriel, as que cuando le traen su cerveza (a 4 C)
a Gonzalo, este decide esperar a que le traigan su cerveza a su amigo. A los 5 minutos, tiempo
que tarda el mesero en traer la cerveza de Gabriel, la cerveza de Gonzalo ya esta a 7 C. En
ese momento se les une Ramiro y pide otra cerveza, la cual tardan en traer otros 5 minutos. A
que temperatura toma su cerveza cada quien?
Solucion
Sustituyendo en la ecuacion 1.12, obtenemos la ecuacion
dT
= k [T 30] ,
dt
cuya solucion es
T (t) = 30 + C2 ekt .
Sustituyendo la condicion inicial T (0) = 4, encontramos que

T (t) = 30 26ekt ,
y aqu usamos la condicion de que T (5) = 7, lo cual nos da k = 0.02452, o sea que la ley de
enfriamiento es
T (t) = 30 26e0.02452t .
Sustituyendo en esta funcion t = 10, que es el tiempo en que empiezan a beber, obtenemos para
la cerveza de Gonzalo
T (10) = 30 26e10(0.02452) 9.6 C.
Entonces las temperaturas de las cervezas fueron 4 C para Ramiro, 7 C para Gabriel y 9.6 C
para Gonzalo.
Ejercicios
1. En un cuarto a temperatura de 25 C se deja enfriar agua originalmente a 100 C. Despues de 10
minutos, su temperatura es de 80 C.
a) Encontrar la temperatura del agua despues de 20 minutos.

b) En que momento sera la temperatura de 40 C?, 20 C?
R: a) 65.3C b) 52 min, 139 min
2. Un refresco que estaba a 10 C, se calienta hasta 20 C en un lugar con temperatura ambiente de 40 C.
a) Encontrar su temperatura despues de 20 minutos y despues de media hora.

b) En que momento sera la temperatura de 25 C?
R: a) 34.1C, 37.41C b) 8.5
3. Si una taza de cafe tiene una temperatura de 200 F cuando acaba de servirse y un minuto despues se
ha enfriado hasta 190 F en un cuarto cuya temperatura es de 70 F, determinar el momento en que el
cafe alcanza una temperatura de 150 F.
R: 6.07 minutos
4. A medianoche, se descubre un cuerpo con una temperatura de 85 F, donde la temperatura ambiente es

constante y es de 70 F. El cuerpo se enva rapidamente (suponga que instantaneamente) a la morgue,
en donde la temperatura ambiente se mantiene a 40 F. Al cabo de una hora se encuentra que la
temperatura del cuerpo es de 60 F. Estimar el momento de la muerte.
R: 11:12 pm
5. La temperatura maxima que puede medir cierto termometro es 110 F. Cuando el termometro marca
36 F, se coloca en un horno. Despues de 1 y 2 minutos, las temperaturas respectivas son de 60 F y 80

F. Cual es la temperatura del horno?
R: 324 F
Decaimiento radiactivo
Cuando se tiene un material radiactivo, sus atomos se desintegran con una rapidez propor-
cional a la cantidad presente de los mismos en ese instante. El tiempo que transcurre para que
la masa de material radiactivo decaiga a la mitad de su valor original se llama vida media del
material. La ecuacion diferencial que gobierna el decaimiento radiactivo es
dQ
= rQ, (1.13)
dt
en donde la constante r > 0 se conoce como razon de decaimiento.
Ejemplo
Se tienen 100 mg de isotopo radiactivo de torio 234. Si la masa de este material se reduce a 82.04
mg en una semana, encontrar una expresion para la cantidad presente en cualquier instante.
Tambien, hallar el intervalo de tiempo que debe transcurrir para que la masa decaiga hasta la
mitad de su valor original.
Solucion
Sea Q(t) la cantidad de torio 234 presente en cualquier instante t, en donde Q se mide en
miligramos y t, en das. Se busca la solucion de la ecuacion 1.13 que tambien satisfaga la condicion
inicial
Q(0) = 100,
as como la condicion:
Q(7) = 82.04.
La ecuacion es lineal y tambien separable; su solucion general es:
Q(t) = Cert
en donde C es una constante arbitraria. La condicion inicial requiere que C = 100 y, por lo tanto,
Q(t) = 100ert.
A fin de satisfacer la ecuacion original, se hace t = 7 y Q = 82.04; esto da:
82.04 = 100e7r ,
entonces,
r = 0.02828/da
Por tanto, se ha determinado la razon de decaimiento r. Si se usa este valor de r en la solucion

se obtiene:
Q(t) = 100e0.02828t mg,

con lo que se obtiene el valor de Q(t) en cualquier instante.

Encontremos la vida media del material. Sea el tiempo en que Q(t) es igual a 50 mg. Entonces,
por la ecuacion 1.13,
50 = 100er ,
o bien,
r = ln 2.
La relacion entre la razon de decaimiento y la vida media es valida no solo para el torio 234, sino
para cualquier material que obedezca la ecuacion diferencial 1.13. Al usar el valor de r dado por
la ecuacion hallada antes, se encuentra que, para el torio 234:
ln 2
= 24.5 das.
0.02828
Ejercicios
1. El isotopo radiactivo plutonio 241 decae de forma que se satisface la ecuacion diferencial
dQ
= 0.0525Q
dt
en donde Q se mide en miligramos y t en anos. a) Determinar la vida media r del plutonio 241. b) Si en
este momento se cuenta con 50 mg de plutonio, cuanto quedara en 10 anos?
R: a) 13.20 anos b) 29.6 miligramos
2. El einstenio 253 decae con una rapidez proporcional a la cantidad que se tenga. Determinar la vida
media si este material pierde un tercio de su masa en 11.7 das. R: 20 das
3. El radio 226 tiene una vida media de 1 620 anos. Encontrar el periodo en el que un cuerpo de este
material se reduce a tres cuartas partes de su tamano original. R: 672.4 anos
4. La vida media del carbono 14 es aproximadamente de 5 568 anos. a) Hallar una expresion para Q(t)
en cualquier instante t, si Q(0) = Q0 . b) Suponer que se descubren ciertos restos en los que la cantidad
residual presente de carbono 14 es el 20 % de la cantidad original. Determinar la antiguedad de estos
restos. R: a) Q(t) = Q0 e1.244810 t b) 12 929 anos
4
5. Se encuentran 100 mg de torio 234 en un recipiente cerrado y a este se le agrega torio 234 con una rapidez
constante de 1 mg/da. a) Hallar la cantidad Q(t) de torio 234 que hay en el recipiente en cualquier
instante. b) Hallar la cantidad lmite Q(t) de torio 234 que existira en el recipiente cuando t . R:
a) Q(t) = 35.36 + 64.64e0.02828t, b) 35.36 mg
Captulo 2
Ecuaciones de segundo orden
2.1. Ecuaciones lineales

2.1.1. Introduccion
Las ecuaciones diferenciales lineales (es decir, que no hay potencias ni productos entre la
funcion o sus derivadas) de segundo orden tienen dos soluciones linealmente independientes
(esto es, que no son multiplos una de la otra). En las ecuaciones de segundo orden hay dos
constantes arbitrarias y la solucion general se expresa como:
y = c1 y1 (x) + c2 y2 (x). (2.1)

Las soluciones y1 (x) y y2 (x) son linealmente independientes si el siguiente determinante
(llamado Wronskiano) vale cero

y (x) y2 (x)
W (y1 (x), y2 (x)) = 1

= y1 y y y2 = 0 (2.2)
2 1
y1 (x) y2 (x)
Primero veremos el caso en que la ecuacion es homogenea. Aqu homogenea significa que
en la forma general
y + p(x)y + q(x)y = g(x), (2.3)

g(x) vale cero.
Las dos constantes arbitrarias de una ecuacion de segundo orden se eliminan si conoce-
mos dos condiciones iniciales, una para la funcion y otra para la derivada; pero ambas en el
mismo punto. Cuando se tienen condiciones en dos puntos (los extremos del intervalo donde
esta definida la funcion) se les llama condiciones de frontera. Este tipo de problemas no se
estudiara en este texto.
Una ecuacion de segundo orden de la forma
P (x)y + Q(x)y + R(x)y = 0 (2.4)

es exacta si se puede reescribir en la forma
[P (x)y ] + f [P (x), Q(x), R(x)]y = 0. (2.5)

Al integrar esta ecuacion se obtiene otra de primer orden. Para esto es necesario que
34 2 Ecuaciones de segundo orden
P (x) + Q(x) + R(x) = 0. (2.6)

Si una ecuacion homogenea no es exacta, se puede convertir en exacta al multiplicarla por
un factor integrante (x) que cumpla que:
P + (2P Q) + (P Q + R) = 0. (2.7)
Esta es la ecuacion adjunta de la original y su adjunta es la original. Si la ecuacion original
y su adjunta son la misma, se llama ecuacion autoadjunta, para lo cual se necesita que P (x) =
Q(x).
2.1.2. Reduccion de orden

Si tenemos la ecuacion
y + p(x)y + q(x)y = 0 (2.8)

y conocemos una solucion y1 (x), podemos hallar la segunda solucion por medio de
Z
y2 (x) dx
Z
= exp p(x)dx . (2.9)
y1 (x) [y1 (x)]2
Esto reduce el problema a solo encontrar una solucion y determinar la otra por medio de
la formula anterior.
Ejemplo
(1 x2 )y 2xy + 2y = 0,
sabiendo que una solucion es y1 = x.
Solucion
De la formula 2.9 tenemos que
dx R
Z
2x
dx
y2 = x e 1x2
x2
as que tenemos que resolver dos integrales. Primero hacemos u = 1 x2 , con lo cual du = 2xdx
y tenemos que
2x du
Z Z
= = ln u = ln(1 x2 ).
1 x2 u
As pues tendremos que
dx
Z
y2 = x .
x2 (1 x2 )
Resolviendo esta otra integral por descomposicion en fracciones parciales

2.2 Coeficientes constantes 35
dx dx 1 dx 1 dx 1 1
Z Z Z Z
= + + = + ln(1 x2 ).
x2 (1 x2 ) x2 2 1+x 2 1x x 2
Con esto tenemos que la segunda solucion es

1 1 x
y2 = x + ln(1 x2 ) = 1 + ln(1 x2 ).
x 2 2
Ejercicios
Dada una solucion de las siguientes ecuaciones, encuentre una segunda solucion para cada una de ellas
1. y 4y 12y = 0, y1 = e6x R: y2 = e2x
2. y + 2y + y = 0, y1 = ex R: y2 = xex
1
3. x2 y + 2xy = 0, y1 = 1 R: y2 = x
1
4. x2 y + 2xy 2y = 0, y1 = x R: y2 = x2
ln x
5. x2 y + 3xy + y = 0, x > 0, y1 = 1/x R: y2 = x
6. x2 y x(x 2)y + (x + 2)y = 0, x > 0, y1 = x R: y2 = xex
7. (1 x cot x)y xy + y = 0, 0 < x < , y1 = x R: y2 = sen x

cos
x
8. x2 y + xy + (x2 1/4)y = 0, y1 = x1/2 R: y2 = x

3x 3x2 1 1x
9. (1 x2 )y 2xy + 6y = 0, y1 = 3x2 1 R: y2 = 4 + 8 ln 1+x
1 C2
10. 2x2 y + 3xy y = 0, y1 = x R: y2 = 32 C1 x + x
2.2. Coeficientes constantes

Para ecuaciones lineales es frecuente proponer soluciones de la forma y(x) = eax . Esto hace
que al derivar dichas soluciones y sustituir en la ecuacion original nos quede solo una relacion
entre las constantes puesto que, como sabemos, la derivada de una funcion exponencial es un
multiplo de ella misma.
La ecuacion
L[y] = ay + by + cy = (aD 2 + bD + c)y = 0, (2.10)

tiene asociada la siguiente ecuacion algebraica, llamada ecuacion caracterstica
ar 2 + br + c = 0. (2.11)
Dependiendo de sus races sera la solucion de esta ecuacion, siendo posibles tres casos para
races reales diferentes, reales iguales o complejas conjugadas.
Races reales diferentes

y = C1 er1 x + C2 er2 x (2.12)
Races reales iguales

y = C1 erx + C2 xerx (2.13)
Races complejas conjugadas
y = C1 ex cos x + C2 ex sen x, (2.14)
donde r = i.
Este ultimo caso resulta de aplicar la formula de Euler para la exponencial imaginaria que
multiplica a la solucion con valores reales.
Ejemplo
y + 5y + 6y = 0, y(0) = 0, y (0) = 1.
Solucion
La ecuacion caracterstica de esta ecuacion es
r2 + 5r + 6 = 0,
r1 = 3, r2 = 2,
y = C1 e3x + C2 e2x .
Al sustituir esta ecuacion en las condiciones iniciales obtenemos
y(0) = C1 + C2 = 0,
o sea, C1 = C2 . Por otra parte, la derivada es
y (x) = 3C1 e3x + 2C2 e2x ,

que nos da lo siguiente
y (0) = 3C1 + 2C1 = 1,

o sea que C1 = 1, y finalmente
y(x) = e2x e3x .
Ejemplo
2.2 Coeficientes constantes 37
y + 4y + 4y = 0.
Solucion
La ecuacion caracterstica de esta ecuacion es
r2 + 4r + 4 = 0,
r1 = r2 = 2,
y = C1 e2x + C2 xe2x .
Ejemplo
Resolver la ecuacion:
y + y + y = 0.
Solucion
La ecuacion caracterstica es
r2 + r + 1 = 0,
cuya solucion es

1 3
r i ,
2 2
con lo que la solucion general es
!
21 x 3 3
y=e C1 cos x + C2 sen x .
2 2
Ejercicios
1. y + 2y 3y = 0 R: y = C1 ex + C2 e3x
1
2. 4y + 4y + y = 0 R: y = (C1 x + C2 )e 2 x
x x
3. 6y y y = 0 R: y = C1 e 2 + C2 e 3
x
4. 2y 3y + y = 0 R: y = C1 e 2 + C2 ex
5. y y = 0 R: y = C1 ex + C2 ex
6. y 2y + y = 0 R: y = (C1 x + C2 )ex

7. y + 5y = 0 R: y = C1 cos 5x + C2 sen 5x

93 5 9+3 5
8. y 9y + 9y = 0 R: y = e 2 x
+e 2 x

9. y 2y 2y = 0 R: y = C1 e(1+ 2)x
+ C2 e(1 2)x
10. y + 2y + y = 0 R: y = (C1 x + C2 )ex
11. y + y 2y = 0, y(0) = 1, y (0) = 1 R: y = ex + 23 e2x
12. y 6y + 9y = 0, y(0) = 0, y (0) = 2 R: y = 12 xe3x

9 x1 1 9(x1)
13. y + 8y 9y = 0, y(1) = 1, y (1) = 0 R: y = 10 e + 10 e
14. y 2y + 2y = 0 R: y = ex (C1 cos x + C2 sen x)

15. y 2y + 6y = 0 R: y = C1 e(1+ 5)x
+ C2 e(1 5)x
16. y + 2y 8y = 0 R: y = C1 e2x + C2 e4x

x
17. 9y 6y + y = 0 R: y = (C1 x + C2 )e 3
18. y + 6y + 13y = 0 R: y = e3x (C1 cos 2x + C2 sen 2x)

1
19. y + 4y = 0, y(0) = 0, y (0) = 1 R: y = 2 sen 2x
20. y + 4y + 5y = 0, y(0) = 1, y (0) = 0 R: y = e2x sen x
2.3. Ecuacion de Euler

La ecuacion
x2 y + xy + y = 0, x > 0, (2.15)
se transforma en una con coeficientes constantes si se hace z = ln x. Al resolver para z, se tiene
que regresar a la variable original x.
Ejemplo
x2 y + xy y = 0
Solucion
Esta ecuacion es de la forma 2.15 por lo que, haciendo el cambio de variable
z = ln x,
obtenemos la ecuacion:
d2 y
= y(z),
dz 2
cuya solucion es
y(z) = A senh z + B cosh z,

2.3 Ecuacion de Euler 39
o bien
y(x) = A senh ln x + B cosh ln x,

que, despues de aplicar la definicon de las funciones hiperbolicas, se simplifica como
C2
y(x) = C1 x + .
x
La solucion de esta ecuacion se encuentra mas rapidamente si se propone una solucion de

la forma
y(x) = xr (2.16)
lo cual nos dara para las derivadas
y (x) = rxr1 , (2.17)
y = r(r 1)xr2 , (2.18)

las cuales al sustituirse en la ecuacion 2.15 y eliminar xr , nos daran una ecuacion algebraica
en r llamada ecuacion caracterstica. Al resolver esta ecuacion caracterstica, tendremos los
tres casos para la solucion, dependiendo si las races de la ecuacion caracterstica son reales
diferentes, reales iguales o complejas conjugadas.
Races reales diferentes La solucion es de la forma:
y(x) = C1 xr1 + C2 xr2 . (2.19)
Races reales iguales La solucion es de la forma:
y(x) = xr (C1 + C2 ln(x)). (2.20)
Races complejas conjugadas En este caso, las races tienen la forma:
r1,2 = a ib,
con esto, la solucion es de la forma:
y(x) = xa [C1 cos(b ln x) + C2 sen(b ln x)]. (2.21)
Ejemplo
x2 y + xy y = 0
Solucion
Esta ecuacion es de la forma 2.15, por lo que proponiendo una solucion de la forma 2.16 y
sustituyendo las ecuaciones 2.17 y 2.18 obtenemos
x2 r(r 1)xr2 + xrxr1 xr = 0.
Eliminando xr obtenemos la ecuacion caracterstica
r2 1 = 0,
cuya solucion es
r1 = 1, r2 = 1.
Con esto, la solucion es
y(x) = C1 x + C2 /x.
Ejercicios
1. x2 y + xy + y = 0 R: y = C1 cos(ln x) + C2 sen(ln x)
2. x2 y + 4xy + 2y = 0 R: y = C1 x2 + C2 x1

3. x2 y + 3xy + 5y/4 = 0 R: y = x1 [C1 cos(ln x) + C2 sen(ln x)]
4. x2 y 4xy 6y = 0 R: y = C1 x1 + C2 x6
5. x2 y 3xy + 4y = 0 R: y = x2 (C1 ln x + C2 )
1
6. x2 y + 2xy + y/4 = 0 R: y = x 2 x (C1 + C2 ln x)
7. x2 y + 3xy + 5y = 0 R: y = x1 [C1 cos(2 ln x) + C2 sen(2 ln x)]
8. x2 y xy + y = 0 R: y = x(C1 ln x + C2 )
9. x2 y 5xy + 9y = 0 R: y = x3 (C1 + C2 ln x)
1

15 15
10. x2 y + 2xy + 4y = 0 R: y = x 2 [C1 cos( 2 ln x) + C2 sen( 2 ln x)]
11. x2 y 4xy + 4y = 0 R: y = C1 x + C2 x4

12. x2 y + 6xy y = 0 R: y = C1 x(5+ 29)/2
+ C2 x(5 29)/2
3

3 3
13. 2x2 y 4xy + 6y = 0 R: y = x 2 [C1 cos( 2 ln x) + C2 sen( 2 ln x)]
14. (x 1)2 y + 8(x 1)y + 12y = 0 R: y = C1 (x 1)3 + C2 (x 1)4
15. (x + 1)2 y + 3(x + 1)y + 34 y = 0 R: y = C1 (x + 1)1/2 + C2 (x + 1)3/2
16. (x 2)2 y + 5(x 2)y + 8y = 0 R: y = (x 1)2 [C1 cos(ln(x 2)) + C2 sen(ln(x 2))]
2.4 Ecuaciones no homogeneas 41
2.4. Ecuaciones no homogeneas

En esta seccion, no homogenea significa que es de la forma
y + p(x)y + q(x)y = g(x), (2.22)
con g(x) 6= 0.
Para resolver ecuaciones no homogeneas, primero se debe hallar una solucion de la ecuacion
homogenea asociada (o sea, la misma ecuacion, pero sustituyendo g(x) por 0), para despues
hallar una solucion particular que satisfaga a la ecuacion no homogenea. La solucion general
constara de la solucion homogenea mas la solucion particular.
Ejemplo
y + 4y = 1 + x + sen x.
Solucion
La ecuacion homogenea asociada es
y + 4y = 0,
cuya solucion es
y = C1 cos 2x + C2 sen 2x.
Ahora necesitamos soluciones particulares de las ecuaciones
y + 4y = 1, y + 4y = x, y + 4y = sen x.
Las soluciones particulares para cada una de esas ecuaciones son
1 1 1
y= , y= x, y= sen x
4 4 3
respectivamente, segun se puede ver al derivar.
Entonces, la solucion general es
1 1 1
y= + x + sen x + C1 cos 2x + C2 sen 2x.
4 4 3
Enseguida se dan dos metodos para hallar la solucion particular de la no homogenea.

2.4.1. Metodo de los coeficientes indeterminados

Para hallar la solucion particular de la ecuacion:
y + p(x)y + q(x)y = g(x), (2.23)

se busca una funcion que tenga una forma tal que al derivar nos de la funcion g(x), por
ejemplo, g(x) multiplicada por un multiplo de x; o en el caso de funciones trigonometricas,
una combinacion lineal de funciones trigonometricas apropiada y se dejan las constantes como
parametros libres a determinar una vez que se sustituyan en la ecuacion diferencial.
Es importante tener en cuenta que si g(x) no es solucion de la homogenea asociada, fre-
cuentemente es suficiente con un multiplo de ella para hallar la solucion particular; pero si
g(x) es solucion de la homogenea asociada, sera necesario multiplicarla por un polinomio. Esto
se ilustra en los ejemplos.
Ejemplo
Hallar una solucion particular para la ecuacion
y 3y 4y = 2 sen x.
Solucion
Para que al derivar y hacer la operaciones dadas en la ecuacion podamos obtener 2 sen x, la
solucion particular debe contener sen x o cos x. Entonces proponemos la siguiente solucion par-
ticular
yp = A cos x + B sen x,
que al derivarla nos da
yp = A sen x + B cos x, yp = A cos x B sen x.
Sustituyendo estas derivadas en la ecuacion obtenemos
(A cos x + B sen x) 3(A sen x + B cos x) 4(A cos x + B sen x) = 2 sen x.
Para que esta ecuacion se cumpla es necesario que la suma de coeficientes de cosenos sea 0, y la
de senos sea 2, o sea
5A 3B = 0, 3A 5B = 2.
3 5
Este sistema tiene como solucion A = 17 , B = 17 , por lo que la solucion particular es
3 5
yp = cos x sen x.
17 17
Ejemplo
y 3y 4y = 4x2 .
Solucion
Aqu proponemos como solucion particular un polinomio del mismo grado que g(x)
yp = Ax2 + Bx + C,
que al derivarla nos da
yp = 2Ax + B, yp = 2A.
Sustituyendo estas derivadas en la ecuacion obtenemos
2A 3(2Ax + B) 4(Ax2 + Bx + C) = 4x2 .
De la ecuacion anterior obtenemos las condiciones
2A 3B 4C = 0, 6A 4B = 0, 4A = 4.
3 13
La solucion de este sistema es A = 1, B = 2 y C = 18 , por lo que la solucion particular es
13 3
yp = + x x2 .
18 2
Ejemplo
y 3y 4y = ex .
Solucion
Podramos pensar que una solucion particular yp = Aex nos servira; pero resulta que esta es
una solucion general de la homogenea asociada. Entonces propondremos la solucion
yp = Axex ,
y derivamos para obtener
yp = Aex Axex , yp = 2Aex + Axex .
Sustituyendo en la ecuacion diferencial
(Ax 2A)ex 3(A Ax)ex 4Axex = ex ,

o sea:
1
5A = 1, A= .
5
La solucion particular es
1
yp = xex .
5
En general entonces, si g(x) tiene la forma indicada, debemos probar yp (x) de la forma
sugerida en la siguiente tabla
g(x) yp (x)
Pn = an xn + ... + a0 s
x Qn (x)
Pn ex xs Qn ex
sen x
Pn ex xs Qn ex sen x + xs Rn ex cos x
cos x
siendo s = 0, 1, dependiendo de si g(x) es o no solucion de la ecuacion homogenea asociada.
Ejercicios
3
1. y 2y 3y = 3e3x R: y = C1 ex + C2 e3x + 16 (4x 1)e3x
1
R: y = e 2 x C1 cos 19 19
12 6
y + y + 5y = 3 sen 2x

2. 4 x + C2 sen 4 x + 5 sen 2x 5 cos 2x
3. y 2y 3y = 3xex R: y = C1 e3x + C2 ex + 12 x2 ex + 14 xex

1 1
4. y + 2y = 3 + 4 sen 2x R: y = C1 + C2 e2x + 3x 2 sen 2x 2 cos 2x
2 1 2 3x
5. y + 9y = x2 e3x + 6 R: y = C1 cos 3x + C2 sen 3x + 3 + 18 x e
6. y + 2y + y = 2ex R: y = (C1 x + C2 )ex + 32 x2 ex

x
7. 2y + 3y + y = x2 + 3 sen x R: y = C1 ex + C2 e 2 + x2 3 sen x 6x + 14
8. y + y = 3 sen 2x + x cos 2x R: y = C1 cos x + C2 sen x 59 sen 2x 91 cos 2x 12

x
sen 2x x3 cos 2x
1
y + y + 4y = 2 senh x R: y = e 2 x C1 cos 15 15 5 1

9. 4 x + C2 sen 4 x + 12 senh x 12 cosh x
3
10. y y 2y = cosh 2x R: y = C1 ex + C2 e2x + 2 cosh 2x + senh 2x
11. y + y 2y = 2x, y(0) = 0, y (0) = 1 R: y = 43 ex 13 e2x x 1

19 3
12. y + 4y = x2 + 3ex , y(0) = 0, y (0) = 2 R: y = 40 cos 2x 10 sen 2x + 14 x2 1
+ 35 ex
8
y 2y + y = xex + 4, y(0) = 1, y (0) = 1 R: y = ex 21 x2 + x 3 + 4

13.
14. y 2y 3y = 3xe2x , y(0) = 1, y (0) = 0 R: y = e3x + 23 ex xe2x 23 e2x
15. y + 4y = 3 sen 2x, y(0) = 2, y (0) = 1 R: y = 2 cos 2x 18 sen 2x 34 x cos 2x
y + 2y + 5y = 4ex cos 2x, y(0) = 1, y (0) = 0 R: y = ex cos 2x + 21 sen 2x + x sen 2x

16.
17. y 5y + 6y = 2ex R: y = 32 xex
18. y y 2y = 2ex R: y = 23 x2 ex
19. y + 2y + y = 3ex R: y = C1 (1 + x) + C2 ex + 12 e2x (x 1)

2.4.2. Metodo de variacion de parametros

La solucion particular de la ecuacion:
y + p(x)y + q(x)y = g(x), (2.24)

se puede encontrar si conocemos las soluciones y1 y y2 de la homogenea asociada, por medio
de la formula:
y2 ()g() y1 ()g()
Z Z
yp (x) = y1 (x) d + y2 (x) d. (2.25)
W (y1 , y2 ) W (y1, y2 )
Esta formula es general, pero puede llevar a integrales difciles de resolver, por lo cual se
recomienda utilizarla solo cuando no se puede aplicar el metodo de coeficientes indeterminados.
Ejemplo
y + y = sec x.
Solucion
La homogenea asociada es:
y + y = 0,
cuyas soluciones son:
y1 = A cos x, y2 = B sen x,
con lo cual calculamos el wronskiano:

A cos x B sen x
W (y1 , y2 ) = = AB cos2 x + AB sen2 x = AB.
A sen x B cos x
Esto nos lleva a tener que resolver las integrales:
A cos x B sen x
Z Z
yp = B sen x sec xdx + A cos x sec xdx
AB AB
Z Z
= cos x tan xdx + sen x dx
= cos x ln(cos x) + x sen x.
Entonces la solucion general es:
y = A cos x + B sen x + cos x ln(cos x) + x sen x.

Ejercicios
1. y + y = tan x R: y = C1 cos x + C2 sen x cos x [sen x + ln | sec x + tg x]
2.
y + 2y + y = e x
ln x R: y = (C1 x + C2 )ex + x2 ex ln x 43 x2 ex
1
3. y + 4y = sec x R: y = C1 cos 2x + C2 sen 2x + cos x cos 2x + sen x sen 2x 2 sen 2x ln | sec x + tg x|
1 ex e2x
4. y + 3y + 2y = 1+e2x R: y = C1 ex + C2 e2x + 2 ln |1 + ex | + 2 ln |1 + e2x |
5. y 2y + 2y = e tan x
x
R: y = ex (C1 cos x + C2 sen x) ln | sec x + tg x|ex cos x
1 1
6. y y = 1+ex R: y = 2 [ex ln(1 + ex ) ex ln(1 + ex ) 1]
1
7. y + y = 1+sen x R: y = sen x ln(1 + sen x) + 1 x cos x + cos x tg x
x
e 1
8. y y = 1+e2x R: y = C1 ex + C2 ex 2 ln(1 + e2x ) cosh x
9. y + y = cot x R: y = C1 sen x + C2 cos x + sen x ln(csc x ctg x)
10. y + y = sec x R: y = C1 cos x + C2 sen x + cos x ln(cos x) + x sen x
1 x
11. y + 4y = csc 2x R: y = C1 sen 2x + C2 cos 2x + 4 sen 2x ln sen 2x 2 cos 2x
4x 2x
e xe e e4x
y y = ex R: y = C1 e2x + C2 e (3+ 3x (3 3)x 2x
12. + C3 e + 1742 88 32 + 192
13. y + 3y + 2y = 3e2x + x R: y = C1 ex + C2 e2x + x2 34 3xe 2x
x
1 e e e e2x
y + y 2y = ln x R: y = C1 e2x + C2 ex
R x 2x R
14. 2 ln x + 3 x dx + 6 x dx
2 2
y y = ex R: y = C1 ex + C2 ex + 2e x
exx dx
R
15.
1 3x x
16. 2y + 3y + y = e3x R: y = C1 ex + C2 e 2 + 10 e

y 4y + 4y = x R: y = C1 e2x + C2 xe2x 12e2x x3/2 e2x dx + 12x xe2x dx
R R
17.

18. y + 9y = 9 sec2 3x, 0<x< 2 R: y = sen 3x ln | tg x + sec x| 1
e
y + 4y + 4y = R: y = e2x ln x
2x
19. x2 , x>0
3
20. y + 4y = 3 csc 2x, 0<x< 2 R: y = 4 sen 2x ln(sen 2x) 32 x cos 2x
2.5. Aplicaciones de modelado
Oscilaciones mecanicas
La segunda ley del movimiento de Newton establece que la aceleracion que experimenta
un cuerpo es proporcional a la suma de las fuerzas que actuan sobre ella, e inversamente
proporcional a su propia masa. Esto se escribe generalmente como
F~ = m~a,
X
(2.26)
o en forma escalar
X X
Fx = max , Fy = may . (2.27)
Para un sistema mecanico vibratorio podemos tener, ademas de las fuerzas provocadas por
el estiramiento del resorte, fuerzas de friccion con el medio o alguna fuerza externa constante
o variable.
Oscilaciones libres
Cuando no hay fuerzas externas (incluida la friccion), las oscilaciones solo estan asociadas
con la rigidez del resorte y la inercia de la partcula. En este caso la ecuacion de movimiento
es
mx + kx = 0, (2.28)
donde x es la posicion de la partcula, y las primas significan derivadas con respecto al tiempo.
La solucion general de esta ecuacion es:
x(t) = A cos 0 t + B sen 0 t, (2.29)

donde 02 = k/m es la llamada frecuencia natural de la oscilacion. Los valores de las constantes
A y B se determinan de las condiciones iniciales del problema dado. La oscilacion tiene un
periodo
r
m
T = 2 .
k
A veces se utiliza mas la notacion
x(t) = R cos cos 0 t + R sen sen 0 t, (2.30)

donde R = A2 + B 2 y = tan1 (B/A). R representa el desplazamiento maximo de la
partcula y la fase, o sea, el corrimiento de la onda con respecto al tiempo.
Ejemplo
Una masa de 0.1 kg se cuelga de un resorte y lo alarga 20 cm. Si la masa se desplaza 2 cm de la
posicion de equilibrio y se le da una velocidad inicial de 0.1 m/s, encontrar la posicion x de la
partcula como funcion de t.
Solucion
De la ley de Hooke F = kx, encontramos que
2
0.1kg(9.8m/s )
k = F/m = = 4.9N/m.
0.2m
La ecuacion de movimiento es
0.1x + 4.9x = 0,
que tiene como solucion
x(t) = A cos(7t) + B sen(7t).
Las condiciones iniciales son x(0) = 0,2, x (0) = 0,1, lo cual nos da los valores de las constantes
A = 1/50, B = 1/70. Con esto, la solucion final es
1 1
x(t) = cos(7t) + sen(7t),
50 70
x(t)
Figura 2.1: Oscilaciones libres
o bien
x(t) = 0.0246 cos(7t + 0.6202).
La grafica de esta funcion se ve en la figura 2.1.
Oscilaciones amortiguadas
Si el medio en el que se mueve el oscilador es viscoso, la ecuacion de movimiento es
mx + x + kx = 0, (2.31)
con la constante de amortiguamiento. Las races de la ecuacion caracterstica son
p
2 4km
r1,2 = . (2.32)
2m
Dependiendo del valor de = 2 4km se tienen tres posibles formas de la solucion, que
presentan tres formas de movimiento diferentes
a) > 0 nos da la solucion
x(t) = Aer1 t + Ber2 t . (2.33)

b) = 0 nos da la solucion
x(t) = (A + Bt)ert . (2.34)

c) < 0 nos da la solucion
x(t) = et (A cos t + B sen t). (2.35)

El caso a) nos da una solucion que decae exponencialmente segun se ve en la figura 2.2.
x(t)
Figura 2.2: Oscilaciones sobreamortiguadas
El caso b) nos da una solucion que comienza a oscilar, pero conforme transcurre el tiempo,
decae exponencialmente, como se ve en la figura 2.3. A estas soluciones se le llama crticamente
amortiguadas.
El caso c) da una solucion que oscila con su amplitud decreciendo constantemente. Este es
el caso mas general de oscilaciones. Esto se ilustra en la figura 2.4.
Oscilaciones forzadas
Cuando en un oscilador actua una fuerza externa periodica, se tiene un movimiento difer-
ente. Tomemos el caso en que F = F0 cos t, en este caso la ecuacion de movimiento es
mx + x + kx = F0 cos t. (2.36)
El caso mas simple de movimiento forzado es cuando = 0, o sea, cuando no hay amor-
tiguamiento. Con esto, la ecuacion se reduce a
mx + kx = F0 cos t. (2.37)
Cuando 6= 0 , la solucion es
F0
x(t) = C1 cos 0 t + C2 sen 0 t + cos t. (2.38)
m(02 2)
x(t)
Figura 2.3: Oscilaciones con amortiguamiento crtico
Dependiendo de las condiciones iniciales, las constantes C1 y C2 representan diferentes

clases de movimiento. Por ejemplo, si x(0) = 0 y x (0) = 0, los valores de C1 y C2 son
F0
C1 = , C2 = 0,
m(02 2)
y la solucion es
F0
x(t) = (cos t cos 0 t), (2.39)
m(02
2)
que nos da un comportamiento como el que se muestra en la figura 2.5. A esta clase de
movimiento se le llama oscilaciones moduladas.
Por otro lado, si = 0 , la solucion es de la forma
F0
x(t) = C 1 cos 0 t + C2 sen 0 t + t sen 0 t. (2.40)
2m0
El ultimo termino domina la solucion conforme crece el tiempo y, como se ve, crece lineal-
mente. A esta clase de movimiento se le llama oscilaciones con resonancia y el comportamiento
se ilustra en la figura 2.6.
Finalmente, cuando 6= 0, la solucion es de la forma
F0
x(t) = C1 er1 t + C2 er2 t + p 2
cos(t ). (2.41)
m (0 2)2 + 2 2
2
Aqu tambien el ultimo termino domina cuando t , lo cual significa que no importa
cuales sean las condiciones iniciales, el movimiento tiende a la frecuencia que habra unicamente
bajo la accion de la fuerza externa, como se ve en la figura 2.7. Antes de que su frecuencia sea
igual a la de la fuerza externa, se dice que esta en un estado transitorio.
x(t)
Figura 2.4: Oscilaciones amortiguadas
Circuitos electricos
En un circuito electrico donde haya elementos tpicos como lo son resistencias, capacitores
y bobinas, as como fuentes de voltaje, tambien puede haber oscilaciones. De acuerdo con
las leyes de Kirchhoff, la suma de los voltajes en los elementos del circuito debe ser igual al
voltaje aplicado por la fuente, llamado fuerza electromotriz (FEM). Por esto, la ecuacion de
un circuito con una resistencia de valor R, un capacitor de capacitancia C y un inductor de
inductancia L al que se le aplica una FEM de valor E(t) es
dI 1
+ RI + Q = E(t).
L (2.42)
dt C
En esta ecuacion hay dos incognitas, I y Q. Sin embargo, como conocemos la relacion entre
estas cantidades, o sea, I = dQ/dt, la ecuacion se transforma en
d2 Q dQ 1
L 2
+R + Q = E(t). (2.43)
dt dt C
Al resolver esta ecuacion, se tienen situaciones analogas a las que se presentan en las
oscilaciones mecanicas, tan solo cambiando x por Q, m por L, por R, k por 1/C y F (t) por
E(t). En particular, las oscilaciones moduladas son la base para la transmision radiofonica.
Ejemplo
Un circuito en serie tiene un capacitor de 0.25 106F y un inductor de 1 H. Si la carga inicial en
el capacitor es de 106 C y no hay corriente inicial, encuentre la carga en el capacitor en cualquier
instante t.
Solucion
La ecuacion diferencial es (sabiendo que no hay resistencia ni voltaje aplicado)
x(t)
Figura 2.5: Oscilaciones moduladas
d2 Q 1
L + Q = 0,
dt2 C
cuya solucion es
1 1
Q(t) = A cos t + B sen t.
LC LC
Usando las condiciones de frontera obtenemos Q(0) = 106 C = A, I(0) = 0, o sea B = 0.

As pues, la solucion es
Q(t) = 103 cos 2 103 t.
Ejercicios
1. Una masa que pesa 2 lb alarga 6 pulg un resorte. Si se tira hacia abajo la masa otras 3 pulg mas y
luego se suelta, y si no hay resistencia del aire, determinar la posicion de la masa en cualquier instante
t. Encontrar la frecuencia, el periodo y la amplitud del movimiento.
R: u(t) = cos 8t, = 8 rad/s, T = /4 s, R = 41 pie.
2. Una masa de 100 g alarga 5 cm un resorte. Si la masa se pone en movimiento desde su posicion de
equilibrio con una velocidad hacia abajo de 10 cm/s y no hay resistencia del aire, determinar la posicion
de la masa en cualquier instante t. Cuando regresa por vez primera la masa a su posicion de equilibrio?

R: u(t) = 57 sen 14t, t = 14 s.
3. Una masa que pesa 3 lb estira 3 pulg un resorte. Si la masa se empuja hacia arriba, contrayendo el
resorte una distancia de 1 pulg y luego se pone en movimiento con una velocidad hacia abajo de 2
pies/s, y no hay resistencia del aire, encontrar la posicion de la masa en cualquier instante t. Determinar
la frecuencia, el periodo, la amplitud y la fase del movimiento.
x(t)
Figura 2.6: Oscilaciones con resonancia
1
1

R: u(t) = 4 sen 8 2t 12 cos 8 2t
2
= 8 2rad/s,
T = /4
q 2 s,
11
R = 288 pie,
= arc tg 32 2.0113.
4. Un circuito en serie tiene un capacitor de 0.25 106 F y un inductor de 1 H. Si la carga inicial en el

capacitor es de 106 C y no hay corriente inicial, encontrar la carga en el capacitor en cualquier instante
t.
R: Q(t) = 106 cos 2000t Coulomb.
5. Una masa de 20 g estira 5 cm un resorte. Suponga que la masa tambien esta sujeta a un amortiguador
viscoso cuya constante de amortiguamiento es de 400 dinas-s/cm. Si se tira hacia abajo de la masa 2
cm o mas y luego se suelta, encontrar su posicion en cualquier instante t. Determinar la cuasifrecuen-
cia y el cuasiperiodo, as como la razon del cuasiperiodo al periodo de movimiento no amortiguado
correspondiente.

R: u(t) = e10t ( 56 sen 4 6t + 2 cos 4 6t),

=4 6 rad/s,
Td = /2 6 s,
Td 7
T = 2 6 1.4289.

6. Una masa que pesa 10 lb estira un resorte 3 pulg. La masa esta sujeta a un amortiguador viscoso con
constante de amortiguamiento de 2 lb-s/pie. Si la masa se pone en movimiento desde su posicion de
equilibrio con una velocidad hacia abajo de 3 pulg/s, encontrar su posicion en cualquier instante t.
Determinar cuando la masa vuelve por primera vez a su posicion de equilibrio.
1

R: u(t) = 8 3
e2t sen 2 31t, t = /2 3 s.
7. Un resorte se alarga 10 cm por la accion de una fuerza de 3 N. Del resorte se cuelga una masa de 2 kg
y se sujeta tambien a un amortiguador viscoso que aplica una fuerza de 3 N cuando la velocidad de la
masa es de 5 m/s. Si se tira hacia abajo de la masa 5 cm por debajo de su posicion de equilibrio y se
x(t)
Figura 2.7: Oscilaciones forzadas amortiguadas
le imprime una velocidad inicial hacia abajo de 10 cm/s, determine su posicion en cualquier instante t.
Encontrar la cuasifrecuencia y la razon de a la frecuencia natural del movimiento no amortiguado
correspondiente.
R: u(t) = 0.057198e0.15t cos(3.87008t 0.50709),
= 3.87008rad/s,
3.87008
0 =
15
1.4289.
8. Un circuito en serie tiene un capacitor de 105 F, un resistor de 3 102 y un inductor de 0.2 H. La
carga inicial en el capacitor es de 106 C y no hay corriente inicial. Encontrar la carga en el capacitor
en cualquier instante t.
R: Q(t) = 106 (2e500t e1000t ).
9. Cierto sistema vibrante satisface la ecuacion u + u + u = 0. Encontrar el valor del coeficiente de
amortiguamiento para el que el cuasiperiodo del movimiento amortiguado es 50 % mayor que el periodo
del movimiento no amortiguado correspondiente.
p
R: = 20/9 1.4907.
10. Una masa que pesa 8 lb alarga 1.5 pulg un resorte. La masa tambien esta sujeta a un amortiguador con
coeficiente . Determinar el valor de para el que el sistema esta crticamente amortiguado; proporcionar
las unidades .
R: = 8 lb s/pie.
11. Si un circuito en serie tiene un capacitor de C = 0.8 106 F y un inductor de L = 0.2 H, encontrar la
resistencia R de modo que el circuito esta crticamente amortiguado.
R: R = 103 .
12. Suponer que el sistema descrito por la ecuacion mu + u + ku = 0 esta crticamente amortiguado y
que las condiciones iniciales son u(0) = u0 , u (0) = u0 . Si u0 = 0, demuestre que
lm u = 0
t0
pero que u nunca es cero. Si u0 es positiva determine una condicion para u0 que asegure que la masa
pasa por su posicion de equilibrio despues de que se suelta.
R: u0 < u0 /2m.
13. Una masa que pesa 4 lb alarga 1.5 pulg un resorte. La masa se desplaza 2 pulg en la direccion positiva
y se suelta sin velocidad inicial. Si se supone que no hay amortiguamiento y que sobre la masa actua
una fuerza externa de 2 cos 3t lb, formular el problema con valor inicial que describe el movimiento de
la masa.
R: u + 256u = 16 cos 3t, u(0) = 61 , u (0) = 0.
14. Una masa de 5 kg estira 10 cm un resorte. Sobre la masa actua una fuerza externa de 10 sen(t/2) N. La
masa se mueve en un medio que le imparte una fuerza viscosa de 2 N cuando su velocidad es de 4 cm/s.
Si la masa se pone en movimiento a partir de su posicion de equilibrio con una velocidad inicial de 10
cm/s, plantear las ecuaciones que describen el movimiento de la misma.
R: u + 10u + 98u = 2 sen 2t , u(0) = 0, u (0) =0.08.
15. Si un sistema resorte-masa no amortiguado, con una masa que pesa 6 lb y una constante de resorte
de 1 lb/pulg, se pone repentinamente en movimiento en t = 0 por una fuerza externa de 4 cos 7t lb,
determinar la posicion de la masa en cualquier instante y trazar una grafica del desplazamiento contra
t.
64
R: u(t) = 45 (cos 7t cos 8t).
16. Una masa que pesa 8 lb alarga 6 pulg un resorte. Sobre el sistema actua una fuerza externa de 8 sen 8t
lb. Si se tira hacia abajo de la masa 3 pulg y luego se suelta, determinar la posicion de la masa en
cualquier instante t. Determinar los instantes de las primera cuatro veces en las que la velocidad de la
masa es cero.
R: u(t) = 41 (cos 8t + sen 8t 8t cos 8t).
17. Un resorte se alarga 6 pulg por una masa que pesa 8 lb. La masa esta sujeta a un mecanismo de
amortiguacion viscosa que tiene una constante de amortiguamiento de 0.25 lb/pie, y sobre ella actua
una fuerza externa de 4 cos 2t lb. a) Determinar la respuesta de estado estable de este sistema. b) Si
la masa dada se sustituye por una masa m, determinar el valor de m para el que la amplitud de la
respuesta de estado estable es maxima.
8
R: a) u(t) = 901 (30 cos 2t + sen 2t) pies, b) m = 4 slug.
18. Un circuito en serie tiene un capacitor de 0.25106 F, una resistencia de 5 103 y un inductor de 1
H. La carga inicial del capacitor es cero. Si se conecta una batera de 12 V al circuito, y este se cierra
cuando t = 0, determinar la carga en el capacitor cuando t =0.001 s, cuando t =0.01 s y en cualquier
instante t. Determinar tambien la carga lmite cuando t .
R: Q(t) = 106 (e4000t 4e1000t + 3) Coulomb,
Q(0.001) 1.52 106 Coulomb,
Q(0.01) 3 106 Coulomb,
lm Q(t) = 3 106 Coulomb.

t
19. Demostrar que el periodopdel movimiento de una vibracion no amortiguada de una masa suspendida de
un resorte vertical es 2 l/g, en donde l es el alargamiento del resorte debido a la masa m y g es la
aceleracion debida a la gravedad.
20. Demostrar que la solucion del problema con valor inicial
mu + u + ku = 0, u(t0 ) = u0 , u (t0 ) = u0
puede expresarse como la suma u = v + w, en donde v satisface las condiciones iniciales v(t0 ) = u,
v (t0 ) = 0, w satisface las condiciones w(t0 ) = 0, w (t0 ) = u0 y tanto v como w satisfacen la misma
ecuacion diferencial que u. Este es otro caso de superposicion de soluciones de problemas mas simples
para obtener la solucion de un problema mas general.
2.6. Ecuaciones lineales de orden superior

Cuando se tienen ecuaciones lineales de orden superior al segundo, podemos generalizar los
conceptos dados en este captulo. Por ejemplo, en ecuaciones de tercer orden hay tres soluciones
linealmente independientes, el wronskiano de tercer orden es el determinante

y1 y2 y3

W = y1 y2 y3 ,
y y y
1 2 3
con lo que la ecuacion caracterstica es de tercer grado. Tambien debe haber tres condiciones
iniciales (una para la funcion, una para la primera derivada y una para la segunda derivada).
Ademas, mientras mayor sea el orden de la ecuacion, mas combinaciones se pueden tener
en las races de la ecuacion caracterstica. Por ejemplo, si es de tercer grado, puede haber tres
races reales diferentes, dos iguales y una diferente, una real y dos complejas conjugadas, etc.
Como en la mayora de las aplicaciones practicas de las ecuaciones diferenciales es suficiente
trabajar con ecuaciones de primer y segundo ordenes, no se tratara el tema aqu, pero es util
que el lector sepa que hacer en caso necesario.
Captulo 3
Sistemas de ecuaciones lineales
3.1. Metodo de eliminacion

Para resolver un sistema de ecuaciones diferenciales se puede usar el metodo de eliminacion,
que consiste en transformar el sistema en una sola ecuacion de orden superior. Esto se consigue
derivando una (o mas) de las ecuaciones y sustituyendo en la(s) otra(s). Una vez que se tiene
la solucion de la ecuacion de orden superior, recuperamos las otras soluciones de cada una de
las incognitas sustituyendo a la inversa.
Ejemplo
Reslover el sistema
dx
=y (a)
dt
dy
= x (b)
dt
Solucion
Derivando la ecuacion (a) obtenemos
d2 x dy
= .
dt2 dt
Sustituyendo en la ecuacion (b) llegamos a la ecuacion
dx
+ x = 0,
dt2
cuya solucion es
x(t) = A cos t + B sen t.
Para obtener y(t), segun la ecuacion (a), tan solo derivamos x(t), con lo cual
y(t) = A sen t + B cos t.

58 3 Sistemas de ecuaciones lineales
3.2. Metodos matriciales

Para resolver un sistema de ecuaciones diferenciales lineales es muy util hacer uso de matri-
ces por lo que, en lo que sigue, unicamente se mencionaran las operaciones necesarias, as que
se recomienda al lector repasar el algebra de matrices1 .
Si se tiene un sistema de n ecuaciones lineales, como
x1 = p11 (t)x1 + ... + p1n (t)xn + g1 (t)

.. .. ..
. . .

xn = pn1 (t)x1 + ... + pnn (t)xn + gn (t),
cuya solucion es
x1 = 1 (t), ..., xn = n (t);
podemos escribirlo en forma matricial como
x = Px + g(t), (3.1)
siendo

x1 p11 p12 . . p1n

.

p21 p22 p2n

x =
. ,
P=
. . . ,
. . . .
xn pn1 pn2 . . pnn

x1 g1

.

.

x=
. ,
g =
. .

. .
xn gn
Aunque en general la matriz de coeficientes puede ser funcion de t, nos ocuparemos unica-
mente del caso en que todos los coeficientes son constantes, esto es, P(t) = A.
3.3. Sistemas homogeneos

Para resolver este sistema, es necesario encontrar los eigenvalores de la matriz de coefi-
cientes, y dependiendo de los valores que tomen estos, sera la forma de la solucion.
1
Vease por ejemplo: Howard Anton Introduccion al algebra lineal. 2a edicion. Limusa. 2002.
3.3 Sistemas homogeneos 59
Eigenvalores reales diferentes

En este caso, la solucion es de la forma
x = C1 ~1 er1 t + C2 ~2 er2 t + ...Cn ~n ern t , (3.2)

donde r1 , r2 , ... rn son los eigenvalores de la matriz de coeficientes y ~1 , ~2 , ... ~n , son los
eigenvectores asociados a los eigenvalores encontrados.
Ejemplo
Resolver el sistema

1 1
x = x
4 1
Solucion
Para encontrar los eigenvalores resolvemos la ecuacion
det(A rI) = 0,
o sea

1r 1

4 = 0 = (1 r)2 4,
1r
cuyas soluciones son r1 = 1 y r2 = 3. Estos son los eigenvalores, cuyos eigenvectores asociados
son

1 1
~1 = , ~2 = .
2 2
Con esto, la solucion es

1 3t 1
x = C1 e + C2 et .
2 2
Ejercicios
Resolver los siguientes sistemas

3 2 1 t 2
1. x = x R: x = C1 e + C2 e2t
2 2 2 1

4 3 3 1
2. x = x R: x = C1 + C2 e2t
8 6 4 2

2 1 1 1
3. x = x R: x = C1 t
e + C2 et
3 2 1 3

1 1 1 1
4. x = x R: x = C1 e3t + C2 e2t
4 2 4 1

2 1 1 1
5. x = x R: x = C1 e3t + C2 et
1 2 1 1

1 i 1 1
6. x = x R: x = C1 + C2 e2t
i 1 i i

2 10 5 2
7. x = x R: x = C1 + C2 e3t
1 5 1 1

3 6 2 3
8. x = x R: x = C1 + C2 et
1 2 1 1

1 2 1 t 2
9. x = x R: x = C1 e + C2 e2t
3 4 1 3

5 1 2 1 1
10. x = x, x(0) = R: x = 23 e2t + 72 e4t
3 1 1 3 1

2 1 1 1 1 t 1 1
11. x = x, x(0) = R: x = 2 e +2 e3t
5 4 3 1 5
Eigenvalores reales repetidos

Cuando hay eigenvalores repetidos, se buscan soluciones adicionales que incluyan productos
de polinomios y funciones exponenciales, en la forma
~ rt
x1 = e
~ rt + ~ert
x2 = te
~ 2 ert + ~ tert + e
x3 = t ~ rt (3.3)
. .
. .
. .
~ siendo r el eigenvalor repetido y

donde ~ es un vector que satisface la ecuacion (A rI)~ = xi,
~ el eigenvector asociado a dicho eigenvalor. Debemos senalar que la solucion para ~ tendra un
~ En la solucion general podemos prescindir de dicho multiplo.
multiplo de .
Ejemplo
Resolver el sistema

1 1
x = x
1 3
Solucion
Para encontrar los eigenvalores tenemos que

1r 1

1 = r2 4r + 4 = 0,
3r
cuyas soluciones son r1 = 2 y r2 = 2, esto es, hay eigenvalores repetidos. Con esto, una solucion
de la ecuacion dada es

1
x1 (t) = e2t .
1
Para hallar una segunda solucion, hacemos

~ 2t + ~ e2t ,
x2 (t) = te
cuya derivada es
~ 2t + e
x2 (t) = 2te ~ 2t + 2~ e2t .
Sustituyendo esto en la ecuacion original tenemos que

~ 2t + e
2te ~ 2t + 2~ e2t = A te
~ 2t + ~ e2t ,
o bien, agrupando terminos tenemos que
2~ = A,
~
~ + 2~ = A~ .
Si en la segunda de estas ecuaciones despejamos ~ en terminos de ~ tenemos que

~ 0
(A 2I)~ = = ,
1
lo cual nos da como solucion el siguiente vector

0 1
~ = +k ,
1 1
as que la solucion general es

1 1 0
x(t) = C1 e2t + C2 te2t + e2t .
1 1 1
Ejercicios

3 4 2 2 1
1. x = x R: x = C1 et + C2 tet + et
1 1 1 1 0

4 2 1 t 1 t 0
2. x = x R: x = C1 e + C2 te + et
8 4 2 2 1/2

3/2 1 1 1 0
3. x = x R: x = C1 + C2 t
1/4 1/2 2 2 1/2

3 5/2 2 2 0
4. x = x R: x = C1 et + C2 tet + et
5/2 2 1 1 2

1 1 1
5. x = 2 1 1 x
0 1 1

3 0 0 1
R: x = C1 4 et + C2 1 e2t + C3 1 te2t + 0 e2t
2 1 1 1

0 1 1 1 1 0
6. x = 1 0 1 x R: x = C1 1 e2t + C2 0 et + C3 1 et
1 1 0 1 1 1

1 4 3 3 + 4t
7. x = x, x(0) = R: x = e3t
4 7 2 2 + 4t

1 0 0 1 1 0 0
8. x = 4 1 0 x, x(0) = 2 R: x = 2 et + 4 1 tet + 3 0 e2t
3 6 2 30 33 6 1

1 1 1
9. x = 2 1 1 x Ayuda: r1 = r2 = r3 = 2
3 2 4

0 0 1 0 1 1
R: x = C1 1 e2t +C2 1 te2t + 0 e2t +C3 1 t2 e2t + 2 0 te2t + 2 0 e2t
1 1 1 1 1 2
Eigenvalores complejos
Si los eigenvalores son complejos, se tienen eigenvectores de la forma
~i = a + ib,
con lo que los vectores xi son de la forma
xi = et (a cos t b sen t) + iet (a cos t + b sen t) = u + iv,
y los vectores
u = et (a cos t b sen t) y v = et (a cos t + b sen t)
son soluciones del sistema (observese que y son las partes real e imaginaria respectivamente
de los eigenvalores encontrados).
Ejemplo

1/2 1
x = x
1 1/2
Solucion
Los eigenvalores de la matriz de coeficientes son r1 = 21 + i y r2 = 21 i, por lo que los
eigenvectores son
~ 1 1 0
1 = = +i
i 0 1
y
1 1 0
~2 = = i
i 0 1
lo cual nos da como vector x el vector

1 0 1 0
x = et/2 cos t sen t + iet/2 sen t cos t
0 1 0 1
con lo que las soluciones son

t/2 1 0
u=e cos t sen t
0 1
y
1 0
v = et/2 sen t + cos t ,
0 1
que se pueden simplificar como

t/2 cos t t/2 sen t
u=e , v=e .
sen t cos t
Ejercicios

3 2 t cos 2t t sen 2t
1. x = x R: x = C1 e + C2 e
4 1 cos 2t + sen 2t cos 2t + sen 2t

1 1 2 cos 2t 2 sen 2t
2. x = x R: x = C1 et + C2 et
2 1 sen 2t cos 2t

2 5 5 cos t 5 sen t
3. x = x R: x = C1 + C2
1 2 2 cos t + sen t cos t + 2 sen t
5 cos 23 t 5 sen 23 t

2 5/2 t/2 t/2
4. x = x R: x = C1 e 3 3 + C2 e
9/5 1 3(cos 2 t + sen 2 t) 3( cos 32 t + sen 32 t)

1 1 2 cos 3t 2 sen 3t
5. x = x R: x = C1 + C2
5 3 cos 3t + 3 sen 3t 3 cos 3t + sen 3t

5/4 3/4 cos t sen t
6. x = x R: x = C1 et + C2 et
3/4 5/4 2 cos t + sen t cos t + 2 sen t

1 0 0 2 0 0
7. x = 2 1 2 x R: x = C1 et 3 + C2 et cos 2t + C3 et sen 2t
3 2 1 2 sen 2t cos 2t

3 0 2 2 2 sen 2t
8. x = 1 1 0 x R: x = C1 e2t 2 + C2 et cos 2t
+
2 1 0 1 cos 2t 2 sen 2t

2 cos 2t
+C3 et sen
2t

2 cos 2t sen 2t

1 5 1 cos t 3 sen t
9. x = x, x(0) = R: x = et
1 3 1 cos t sen t

3 2 1 cos t 5 sen t
10. x = x, x(0) = R: x = e2t
1 1 2 2 cos t 3 sen t
3.4. Sistemas no homogeneos

Hay varios metodos de solucion de sistemas no homogeneos. Describiremos los casos mas
simples, que son el de coeficientes indeterminados y el de variacion de parametros.
3.4.1. Coeficientes indeterminados

Analogamente a lo que sucede en ecuaciones de segundo orden (captulo 2), suponemos
una solucion de forma conocida para el sistema y la sustituimos, dejando como incognitas a
despejar las constantes, que en este caso son vectores.
Ejemplo
2et

2 1
x = x+ = Ax + g(t)
1 2 3t
Solucion
Con los metodos descritos anteriormente se puede comprobar que la solucion de la homogenea es

1 t 1
x(t) = C1 e + C2 e3t .
1 1
Para hallar una solucion particular del sistema dado, primero lo reescribimos como

2 1 2 0
x = x+ et + t,
1 2 0 3
as que proponemos una solucion particular de la forma
xp = atet + bet + ct + d,
que al derivar nos da

xp = aet atet bet + c,
y sustituyendo en la ecuacion original obtenemos
= aet atet bet + c =

t t 2 t 0
= Aate + Abe + Act + Ad + e + t.
0 3
De la ecuacion anterior agrupamos terminos semejantes como sigue
a = Aa

2
a b = Ab +
0

0
0 = Ac +
3
c = Ad.
Si en las ecuaciones vectoriales-matriciales anteriores se despejan los vectores constantes obten-

emos
3.4 Sistemas no homogeneos 65

1 0 1 3/4
a= , b= , c= , d= ,
1 1 2 5/3
as que la solucion particular es

1 t 0 t 1 3/4
xp = te + e + t+ .
1 1 2 5/3
Finalmente, la solucion general es

1 t 1 3t 1
x(t) = C1 e + C2 e + tet +
1 1 1

0 t 1 3/4
+ e + + t+ .
1 2 5/3
3.4.2. Variacion de parametros

Este metodo consiste en construir una solucion para el sistema no homogeneo a partir de
la solucion homogenea. Para esto se forma una matriz que tiene como columnas a los vectores
solucion del sistema homogeneo. Esta nueva matriz, llamada , se invierte y se multiplica por
el vector que nos da el termino no homogeneo. El producto resultante se integra (ignorando
las constantes de integracion), y el resultado final se multiplica nuevamente por la matriz. En
resumen, la solucion no homogenea es
Z
xp = 1 g(t) dt. (3.4)
Ejemplo
Reslover el sistema

3 1 3t
x = x+
2 4 et
Solucion
La solucion homogenea de este sistema es
e2t e5t

x1 = , x2 =
e2t 2e5t
lo que nos da la matriz
e2t e5t

= .
e2t 2e5t
Al invertir la matriz obtenemos
2 2t 1 2t

1
= 3e 3e ,
1 5t
3e 31 e5t
con lo que la solucion particular buscada sera
2 2t 1 2t
e2t e5t
Z
xp = 3e 3e 3t
dt =
1 5t
e2t 2e5t 3e 13 e5t et
e2t e5t 2te2t + 13 et

Z
= dt =
e2t 2e5t te5t 13 e4t
e2t e5t te2t 12 e2t + 31 et

= 1 5t 1 5t 1 4t =
e2t 2e5t 5 te 25 e 12 e
6 27
+ 14 et 6 27 1

= 5t 50 = 5 t 50 + 4 et .
3 21
5t 25 + 12 et 3
5
21
50
1
2
Con lo anterior, la solucion final es

6 27 1

1 1
x = C1 e2t + C2 e5t + 5
3 t 50
21 + 4
1 et .
1 2 5 50 2
Ejercicios
t
2 1 e
1. x = x+
3 2 t

1 1 1 1 1 0
R: x = C1 et + C2 et + 32 tet 14 et + t
1 3 1 3 2 1

et

2. x = 1 3
x + t
3 1 3e

3 1
2/3
1

R: x = C1 e2t + C2 e2t + et + et
1 3 1/ 3 2/ 3

2 5 cos t
3. x = x+
1 2 sen t

5 cos t 5 sen t
R: x = 15 (2t 23 sen 2t 21 cos 2t+C1 ) + 51 (t 12 sen 2t+ 23 cos 2t+C2 )
2 cos t + sen t cos t + 2 sen t
2t
1 1 e
4. x = x+
4 2 2et

1 1 0 1
R: x = C1 e3t + C2 e2t e2t + 12 et
4 1 1 0
3
4 2 t
5. x = x+
8 4 t2

1 t 1 1 0 1 2 1 1/2
R: x = C1 e + C2 t+ 2 2 ln t + t t2
2 2 1 2 5 0
t1

4 2
6. x = x+
2 1 2t1 + 4

1 2 5t 1 8 1 4 2
R: x = C1 + C2 e + ln t + 5 t + 25
2 1 2 2 1
3.4 Sistemas no homogeneos 67

1 1 2
7. x = x+ et
4 1 1

1 1 1
R: x = C1 e3t + C2 et + 41 et
2 2 8

2 1 1
8. x = x+ et
3 2 1

1 t 1 t 1 t 1
R: x = C1 e + C2 e + e +2 tet
1 3 0 1

5/4 3/4 2t
9. x = x+
3/4 5/4 et

1 t/2 1 2t 1 5 1 17 t 1 1
R: x = C1 e + C2 e +2 t 4 e +6 et
1 1 3 15 3

3
2 1
10. x = x+ et
2 2 1

1 t 2 4t 1 2 1 t 1 2 + 2
R: x = C1 e + C2 e 3 te + 9 et
2 1 2 2 1 2
Captulo 4
Transformadas de Laplace
Un metodo muy sencillo y efectivo para resolver ecuaciones diferenciales es el que utiliza la
transformada de Laplace. En este captulo se explica el calculo de la transformada de Laplace
y su inversion, para su posterior aplicacion a la solucion de ecuaciones y sistemas.
4.1. Calculo de la transformada de Laplace

Se define la transformada de Laplace (Abreviado TL) de una funcion f (t) con la ecuacion:
Z
L{f (t)} = F (s) = est f (t)dt. (4.1)
0
Recordemos de los cursos de calculo que una integral que contenga algun lmite de inte-
gracion infinito se llama integral impropia. Para evaluar la integral impropia dada se usa:
Z a
lm est f (t)dt.
a 0
De aqu observamos que la TL solo existe en el caso de que la integral dada sea conver-
gente. Frecuentemente se escribe f () para indicar el lmite de f (t) cuando t , y as se
usara aqu.
Ejemplo
Calcular la transformada de Laplace de la funcion f (t) = 1.
Solucion
De la ecuacion (4.1), tenemos que:
Z
L{1} = est dt.
0
Para resolver esta integral hacemos u = st, y du = sdt, o bien, dt = (1/s)du, con lo cual:

1 1 u 1
Z
L{1} = eu du = e = .
s 0 s 0 s
70 4 Transformadas de Laplace
Ejemplo
Calcular la transformada de Laplace de la funcion f (t) = t.
Solucion
De la ecuacion (4.1), tenemos que:
Z
L{t} = test dt.
0
Para resolver esta integral hacemos u = t, y dv = est dt (integracion por partes), lo cual nos da
Z
est

1
L{t} = t + est dt = 2 .
s 0 0 s
Ejercicios
Encontrar la transformada de Laplace de cada una de las siguientes funciones. En todas ellas a es una
constante arbitraria.
2
1. f (t) = t2 R: F (s) = s3
6
2. f (t) = t3 R: F (s) = s4
n!
3. f (t) = tn R: F (s) = sn+1
s
4. f (t) = cos at R: F (s) = s2 +a2
a
5. f (t) = sen at R: F (s) = s2 +a2
s
6. f (t) = cosh at R: F (s) = s2 a2
a
7. f (t) = senh at R: F (s) = s2 a2
2
8. f (t) = sen2 t R: F (s) = s(s2 +4
s(s2 +m2 +n2 )

9. f (t) = cos mt cos nt R: F (s) = (s2 +m2 +n2 )4m2 n2
s2 2a2
10. f (t) = cosh2 at R: F (s) = s(s2 4a)2
1
11. f (t) = eat R: F (s) = sa
1
12. f (t) = teat R: F (s) = (sa)2
n!
13. f (t) = tn eat R: F (s) = (sa)n+1
s2 a2
14. f (t) = t cos at R: F (s) = (s2 +a2 )2
2as
15. f (t) = t sen at R: F (s) = (s2 +a2 )2
1
16. f (t) = et sen t R: F (s) = (s1)2 +1
4.2 Algunas propiedades de la transformada de Laplace 71
4.2. Algunas propiedades de la transformada de Laplace

4.2.1. Traslacion
Al calcular la TL de la funcion eat f (t) notamos que
Z Z
at at st
L{e f (t)} = e f (t)e dt = e(sa)t f (t)dt, (4.2)
0 0
para resolver la integral, hacemos u = s a, con lo cual tenemos
Z
eut f (t)dt = F (u) = F (s a). (4.3)
0
Esto quiere decir que, si tenemos alguna funcion cuya TL ya conocemos, multiplicada por
una exponencial, es suficiente con tomar la TL ya conocida y sustituir s por sa. Esto equivale
a trasladar a F (s) en a unidades en el espacio s.
Ejemplo
Hallar la TL de
f (t) = e2t cos 3t.
Solucion
Sabemos que
s
L{cos 3t} = ,
s2 + 9
por lo tanto, usando la ecuacion 4.3 tenemos que
(s 2) s2
L{e2t cos 3t} = 2
= 2 .
(s 2) + 9 s 4s + 13
4.2.2. Funcion escalon

En muchos problemas de circuitos electricos se usa la funcion u(t a), llamada escalon
unitario o funcion de Heaviside, definida como

0 t<a
u(t a) = . (4.4)
1 ta
Su grafica se muestra en la figura 4.1.
Cuando una funcion f (t) definida en [0, ), como la que se muestra en la figura 4.2, se
traslada a unidades; cambia la forma de la funcion debido al cambio en el dominio de definicion,
como se puede ver en la figura 4.3. Para evitar esto, se multiplica la funcion f (t a) por el
escalon unitario, con lo cual se obtiene algo como lo que se muestra en la figura 4.4.
La TL de la funcion u(t a)f (t a) es
Z
L{u(t a)f (t a)} = u(t a)f (t a)est dt, (4.5)
a
f(t)
0
t
Figura 4.1: Funcion de Heaviside
haciendo v = t a se tiene que

Z
u(v)f (v)esteas dt = eas F (s). (4.6)
0
Entonces concluimos que
L{u(t a)f (t a)} = eas F (s), (4.7)

la cual establece que si tenemos una funcion y la trasladamos a unidades, su TL es la que
tena antes de la traslacion multiplicada por eas .
Ejemplo
Hallar la TL de
f (t) = u(t 5) cos(2(t 5)).
Solucion
De la formula 4.7
s
L{cos 2t} = ,
s2 + 4
por lo tanto, usando la ecuacion 4.3 tenemos que
s
Lu(t 5) cos 2(t 5)} = e5s L{cos 2t} = e5s .
s2 + 4
1
0.8
0.6
0.4
0.2
1 2 3 4 5
Figura 4.2: Funcion inicial
0.8
0.6
0.4
0.2
1 2 3 4 5
Figura 4.3: Funcion recorrida
4.2.3. Derivadas
TL de las derivadas
Calculemos la transformada de Laplace de la derivada de una funcion, por ejemplo, f (t)
Z

L{f (t)} = f (t)est dt. (4.8)
0
Integrando por partes, con u = est y dv = f (t), tenemos que

Z Z
st st
f (t)e dt = e f (t) f (t)(sest )dt

0 0
Z 0
= f (0) + s f (t)est dt,
0
1
0.8
0.6
0.4
0.2
1 2 3 4 5
Figura 4.4: Funcion recorrida y multiplicada por el escalon unitario
y en esta ultima ecuacion reconocemos la ultima integral como la TL de f (t), lo cual nos lleva
a la igualdad
L{f (t)} = sF (s) f (0). (4.9)

Para la TL de la segunda derivada de una funcion tenemos, de manera analoga a la TL
anterior
L{f (t)} = s2 F (s) sf (0) f (0). (4.10)

Se puede demostrar que, en el caso general
L{f (n) (t)} = sn F (s) sn1 f (0) sn2 f (0) ... f (n1) (0). (4.11)
La formula para la TL de una derivada sera muy util en la solucion de ecuaciones diferen-
ciales con condiciones iniciales.
Derivada de las TL
Tenemos que para alguna funcion de t
Z
F (s) = est f (t)dt. (4.12)
0
Como la integral anterior depende de t, pero no de s, se puede derivar con respecto a s de
la siguiente manera
dF d d st
Z Z
st

= e f (t)dt = e f (t) dt =
ds ds ds
Z
= (t)est f (t)dt = L {tf (t)} .
0
De la ecuacion anterior, podemos establecer la igualdad
dn
F (s) = L{tf (t)}. (4.13)
dsn
Al aplicar reiterativamente la derivada anterior obtenemos, para el caso general
dn
L{tn f (t)} = (1)n
F (s), (4.14)
dsn
formula que establece que, si tenemos que calcular la TL de una funcion cuya TL conocemos,
multiplicada por alguna potencia de t, solo necesitamos derivar dicha TL el mismo numero de
veces que la potencia de t y colocar el signo apropiado.
Ejemplo
Hallar la TL de
f (t) = t2 cos 2t
Solucion
De la formula 4.14 obtenemos que
d2 d2

s
L t2 cos 2t = (1)2 2 L {cos 2t} = 2 2

=
ds ds s + 4
4 s2 2s(s2 12)

d
= 2 2
= .
ds (s + 4) (s2 + 4)3
Ejercicios
Encontrar la TL de las funciones siguientes
5
1. f (t) = e3t sen 5t R: F (s) = s2 6s+34
120
2. f (t) = t5 e4t R: F (s) = (s4)6
2
3. f (t) = e5t senh 2t R: F (s) = (s+3)(s+7)
s2
4. f (t) = t e2t sen 3t R: F (s) = 2 (s2 2s+5) 2
es
5. f (t) = u(t 3) sen(t 3) R: F (s) = 1+s
e3s
6. f (t) = u(t 1)e(t1) R: F (s) = 1+s2
e 1
7. f (t) = u(t + 1) senh(t + 1) R: F (s) = 2(s2) 2e(s+1)
1
8. f (t) = t e2t R: F (s) = (s2)2
4(3s2 4)
9. f (t) = t2 sen 2t R: F (s) = (s2 +4)3

1 1
10. f (t) = t3 cosh 4t R: F (s) = 3 (s4)4 + (s+4)4
4.3. Metodos para invertir la TL

Ya vimos como, dada una funcion f (t), obtener otra funcion F (s) a la que llamamos
transformada de Laplace. Ahora nos plantearemos el problema inverso: dada una funcion F (s),
de que funcion f (t) es TL? A la funcion f (t) se le llama transformada inversa de F (s) y se
simboliza como L1 {F (s)} = f (t). Aunque existe una formula para invertir la TL de cualquier
funcion, no se estudiara en este texto porque hace uso de la teora de variable compleja. Si
elaboramos una tabla de transformadas (en el apendice se da una breve tabla de TL de las
funciones mas frecuentes), tendremos de inmediato las funciones originales. Para otros casos,
habra que simplificar las TL obtenidas, para reducirlas a TL de funciones conocidas. Por lo
tanto veremos un par de metodos de simplificacion para poder realizar la inversion.
4.3.1. Fracciones parciales

Un metodo para invertir una TL consiste en descomponerla en fracciones parciales, e iden-
tificar cada una de estas con la TL de alguna funcion conocida.
Ejemplo
Encontrar la transformada inversa de
s1
Y (s) = .
s2 s2
Solucion
Para hallar la inversa, notamos primero que el denominador se puede factorizar como
s2 s 2 = (s 2)(s + 1),
con lo cual podemos escribir
s1
Y (s) =
(s 2)(s + 1).
Esta ultima fraccion se puede descomponer en fracciones parciales en la forma
s1 A B
= + ,
(s 2)(s + 1) s2 s+1
lo cual nos da el siguiente sistema de ecuaciones para A y B
A + B = 1, A 2B = 1,
cuya solucion es A = 13 , B = 23 . Esto implica que
1 1 2 1
Y (s) = + ,
3s2 3s+1
y las ultimas son las transformadas de las funciones
1 2t 2 t
e , e ,
3 3
de donde concluimos que
4.3 Metodos para invertir la TL 77
1 2t 2 t
y(t) = L1 {Y (s)} = e + e .
3 3
4.3.2. Convolucion
Otro metodo para invertir una TL consiste en usar la convolucion. Para dos funciones,
digamos f (t) y g(t), se define su convolucion como
Z t Z t
(f g)(t) = f (t x)g(x) dx = f (x)g(t x) dx. (4.15)
0 0
La igualdad anterior establece ademas que la convolucion es una operacion conmutativa. El

uso de esta operacion es util dado que la TL de la convolucion de dos funciones es el producto
de las TL de cada una de ellas. Esto generalmente se usa en sentido inverso: la transformada
inversa de un producto de las transformadas de dos funciones conocidas es la convolucion de
esas dos funciones, esto es
L1 {F (s)G(s)} = (f g)(t). (4.16)
Ejemplo
Hallar la transformada inversa de
a
H(s) = .
s2 (s2 + a2 )
Solucion
Esta funcion se puede considerar como el siguiente producto

1 a
H(s) = F (s)G(s) = ,
s2 s2 + a 2
en donde
1 a
F (s) = = L {t} , G(s) = = L{sen at}.
s2 s2 + a 2
Por lo tanto, usando la formula 4.16 tenemos que
h(t) = L1 {H(s)} = (t sen at) =

Z t Z t
= xu sen a(t x) dx = (t x) sen ax dx.
0 0
Como la convolucion es conmutativa, las dos integrales anteriores deben dar el mismo resultado,
por lo que existe la libertad de elegir la integral a resolver. Si optamos por la segunda, obtendremos
Z t Z t Z t
(t x) sen ax dx = t sen ax dx x sen ax dx
0 0 0
t h x cos ax it h sen ax it
= ( cos ax)t0 +
a a 0 a 0
t t t cos at sen at
= cos at + + 0 +0
a a a a
at sen at
= .
a2
Ejercicios
Encontrar las funciones de t de las cuales son TL las siguientes funciones.
3
1. s2 +a2 R: f (t) = 3 sena at
4
2. (s1)3 R: f (t) = 2t2 et
2
3. s2 +3s4 R: f (t) = 54 e4t + 45 et
3s 1
q
3 (3 15)t
h i
4. s2 s6 R: f (t) = 2 5e 3 + 15 + (3 + 15)e2 15t
2s+2
5. s2 +2s+5 R: f (t) = et [cos 2t + sen 2t]
2s3
6. s2 4 R: f (t) = 47 e2t + 41 e2t
2s+1
7. s2 2s+2 R: f (t) = 2et [2 cos t + 3 sen t]
8s2 4s+12
8. s(s2 +4) R: f (t) = 3 + 5 cos 2t 2 sen 2t
12s
9. s2 +4s+5 R: f (t) = e2t (3 sen t 2 cos t)
2s3 5
10. s2 +2s+10 R: f (t) = et (2 cos 3t 3 sen 3t)
4.4. Solucion de ecuaciones con la TL

Para resolver una ecuacion diferencial con condiciones iniciales por medio de la TL, primero
obtenemos la TL de la ecuacion, usando la formula 4.9. Observese que es necesario tener
condiciones iniciales, ya que intervienen en la TL de la derivada. Una vez que se tiene una
ecuacion en donde interviene la TL de la funcion buscada, se despeja y se busca su inversa.
Esta funcion es la solucion buscada. Notese que con este metodo no se obtienen constantes
arbitrarias que despues haya que determinar, ni es necesario resolver primero la homogenea.
4.4 Solucion de ecuaciones con la TL 79
Ejemplo
Resolver la ecuacion con condiciones iniciales
y + y = sen 2t, y(0) = 0, y (0) = 1.
Solucion
Sea Y (s) la TL de y(t). Entonces la TL de la ecuacion es
2
s2 Y (s) sy(0) y (0) + Y (s) = .
s2 + 4
Despejando Y (s) obtenemos
s2 + 6
Y (s) = ,
(s2 + 1)(s2 + 4)
y ahora hay que invertir Y (s) para hallar la solucion y(t). Descomponiendo en fracciones parciales
tenemos
as + b cs + d
Y (s) = + ,
s2 + 1 s2 + 4
lo cual nos da los valores a = 0, b = 5/3, c = 0 y d = 2/3, que nos dan para Y (s) la ecuacion
5 1 1 2
Y (s) = .
3 s + 1 3 s2 + 4
2
La funcion anterior la reconocemos como la TL de la funcion
5 1
y(t) = sen t sen 2t
3 3
que es la solucion buscada.
Ejercicios
Resolver las siguientes ecuaciones diferenciales con las condiciones iniciales dadas usando la TL.
1. y y 6y = 0, y(0) = 1, y (0) = 1 R: y = 51 (e3t + 4e2t )
2. y + 3y + 2y = 0, y(0) = 1, y (0) = 0 R: y = 2et e2t
3. y 2y + 2y = 0, y(0) = 0, y (0) = 1 R: y = et sen t
4. y 4y + 4y = 0, y(0) = 1, y (0) = 1 R: y = e2t te2t

5. y 2y 2y = 0, y(0) = 2, y (0) = 0 R: y = 2et cosh 3t 23 et senh 3t
6. y + 2y + 5y = 0, y(0) = 2, y (0) = 1 R: y = 51 e3t + 45 e2t
y + 2 y = cos 2t, 2 6= 4, y(0) = 1, y (0) = 0 R: y = 214 ( 2 5 cos t + cos 2t

7.
8. y 2y + 2y = cos t, y(0) = 1, y (0) = 0 R: y = 15 (cos t 2 sen t + 4et cos t 2et sen t)
9. y 2y + 2y = et , y(0) = 0, y (0) = 1 R: y = 15 (et et cos t + 7et sen t)
10. y + 2y + y = 4et , y(0) = 2, y (0) = 1 R: y = et (2 + t + 2t2 )

4.5. Funcion de Dirac

A veces, en la ingeniera es necesario tratar con fenomenos de naturaleza impulsiva, por
ejemplo fuerzas o voltajes de gran magnitud que actuan en intervalos de tiempo muy cortos.
Tales fenomenos generan ecuaciones de la forma
y + py + qy = f (t), (4.17)
siendo f (t) la fuerza (o voltaje, etc.) que actua en el intervalo [t0 , t0 + ], y es cero en
cualquier otro instante t. Una medida de la intensidad de tal fuerza es
Z t0 + Z
I( ) = f (t)dt = f (t)dt, (4.18)
t0
cantidad llamada impulso total. Supongamos que t0 = 0 y que f (t) es

1
2
< t <
f (t) = (4.19)
0 t ot .
Con esta definicion, I( ) = 1, independientemente del valor de (a excepcion de = 0,
por supuesto). Para intervalos de tiempo cada vez menores, tendremos
lm f (t) = 0, t 6= 0 (4.20)
0
pero tambien
lm I( ) = 1. (4.21)
0
La funcion especial que cumple los lmites anteriores se llama delta de Dirac y esta definida
como
(t) = 0, t 6= 0 (4.22)
Z
(t)dt = 1. (4.23)

Otra propiedad muy usual es que
Z
(t t0 )f (t)dt = f (t0 ). (4.24)

La TL de la funcion se calcula de la siguiente manera
L{(t t0 )} = lm L{f (t t0 )}. (4.25)

0
Z
L{(t t0 )} = lm est f (t t0 )dt =
0 0
t0 +
1
Z
lm est dt =
0 2 t0
4.5 Funcion de Dirac 81
1 st t=t0 + 1 st0 s
= lm e = lm e e es =
0 2 s t=t0 0 2 s
senh s st0
lm e = est0 (4.26)
0 s
(4.27)
Por lo tanto
L{(t t0 )} = est0 . (4.28)

Ahora calculemos la TL del producto de una funcion continua g(t) por la funcion delta
Z
L{(t t0 )g(t)} = (t t0 )g(t)estdt.
0
Para encontrar esta integral, sea f (t) = g(t)est. De la ecuacion 4.28 obtenemos
Z
(t t0 )f (t)dt = f (t0 ).
0
As que, finalmente
L{(t t0 )g(t)} = g(t0 )est0 . (4.29)

Esta ultima igualdad es muy usual al resolver ecuaciones con condiciones iniciales.
Ejemplo
Resolver la ecuacion dada, sujeta a las condiciones iniciales indicadas
y + 2y + 2y = (t ), y(0) = 0, y (0) = 0.
Solucion
Calculando la TL de la ecuacion obtenemos
s2 Y (s) sy(0) y (0) + 2sY (s) 2y(0) + 2Y (s) = es .

Usando las condiciones iniciales dadas, esto se reduce a
es es
Y (s) = = .
s2 + 2s + 2 (s + 1)2 + 1
La transformada obtenida nos recuerda a s21+1 , cuya inversa es sen t. El termino (s + 1)2 nos
indica el producto por la exponencial et , lo cual da

1 1
L = et sen t.
(s + 1)2 + 1
Finalmente, el factor es indica una traslacion en unidades con un producto por u(t ), lo
cual nos da la solucion
es

y(t) = L1 = = u(t )e(t) sen(t ).
(s + 1)2 + 1
Ejercicios
Resolver las siguientes ecuaciones diferenciales con las condiciones iniciales dadas usando la TL.
1. y + 2y + 2y = (t ), y(0) = 1, y (0) = 0 R: y = et cos t + et sen t u(t )e(t) sen t
2. y + 4y = (t ) (t 2), y(0) = 0, y (0) = 0 R: y = 21 u(t ) sen 2t 12 u(t 2) sen 2t
3. y + 2y + y = (t) + u(t 2), y(0) = 0, y (0) = 1

R: y = 2tet + u(t 2) 1 e(t2) (t 2)e(t2)

4. y y = 2(t 1), y(0) = 1, y (0) = 0 R: y = cosh t + 2u(t 1) senh(t 1)
5. y + 2y + 3y = sen t + (t ), y(0) = 0, y (0) = 1

R: y = y 12 2et sen 2t + 14 et cos 2t + 14 (sen t cos t) + 21 2 u(t )e(t) sen 2(t )
6. y + y = 2(t ), y(0) = 1, y (0) = 0 R: y = cos t u(t ) sen t
7. y + y = 2(t ) cos t, y(0) = 0, y (0) = 1 R: y = [1 + u(t )] sen t
8. y + 4y = 2(t /2), y(0) = 1, y (0) = 0 R: y = cos 2t 12 u(t 2 ) sen 2t
9. y + 9y = 2(t /3), y(0) = 1, y (0) = 0 R: y = cos 3t 13 u(t 3 ) sen 3t
10. y + 2 y = 2(t /), y(0) = 1, y (0) = 0 R: y = cos t 1 u(t ) sen t
11. y 3y = (t 2), y(0) = 0 R: y = e3(t2) u(t 2)
12. y + y = (t 1), y(0) = 2 R: y = 2et + e(t1) u(t 1)
13. y + y = (t 2), y(0) = 0, y (0) = 1 R: y = sen t + sen t u(t 2)

1
14. y + 16y = (t 2), y(0) = 0, y (0) = 0 R: y = 4 sen 4(t 2) u(t 2)
15. y + y = (t /2) + (t 3/2), y(0) = 0, y (0) = 0 R: y = cos t u(t 2 ) + cos t u(t 3

2 )
16. y + y = 1 + (t 2) + (t 4), y(0) = 1, y (0) = 0

R: y = cos t + cos(t 2) u(t 2) + cos(t 4) u(t 4)
y + 2y + y = (t 1), y(0) = 0, y (0) = 0 R: y = 12 21 e2t + 12 12 e2(t1) u(t 1)

17.
1 1
18. y 2y = 1 + (t 2), y(0) = 0, y (0) = 1 R: y = 2 cosh 2t + senh 2t 2 cosh 2(t 2) u(t 2)
19. y + 4y + 5y = (t 2), y(0) = 0, y (0) = 0 R: y = e2(t2) sen t u(t 2)
20. y + 2y + y = (t 1), y(0) = 0, y (0) = 0 R: y = te(t1) u(t 1)
4.6. Solucion de sistemas usando TL

La Transformada de Laplace se puede usar para resolver un sistema lineal con condiciones
iniciales. Para ello se toma la TL de cada ecuacion, reservando un smbolo convencional para la
transformada de cada una de las funciones involucradas, y se despeja el sistema de ecuaciones
algebraicas por cualquier metodo conocido. Despues se invierte cada una de estas transformadas
para obtener la solucion final. Notese que tambien en este caso es irrelevante si el sistema es
homogeneo o no.
4.6 Solucion de sistemas usando TL 83
Ejemplo
Resolver el sistema
y + 6y = x , 3x x = 2y , x(0) = 2, y(0) = 3.
Solucion
Llamemos X(s) a la TL de x(t) y sea Y (s) la TL de y(t). Tomando las TL de cada ecuacion
obtenemos
[sY (s) y(0)] + 6Y (s) = [sX(s) x(0)],
3X(s) [sX(s) x(0)] = 2[sY (s) y(0)],

lo cual, despues de sustituir las condiciones iniciales nos da el sistema
sX (s + 6)Y = 1,
(3 s)X 2sY = 8,
cuya solucion es
2s + 16 4 2
X(s) = = ,
(s 2)(s + 3) s2 s+3
3s 1 1 2
Y (s) = = + .
(s 2)(s + 3) s2 s+3
Invirtiendo las transformadas encontradas obtenemos

4 2
x(t) = L1 = 4e2t 2e3t ,
s2 s+3

1 1 2
y(t) = L + = e2t + 2e3t ,
s2 s+3
que constituyen la solucion buscada.
Ejercicios
Resolver los siguientes sistemas dados usando la TL.
1. y = y, y = x, x(0) = 2, y(0) = 1 R: x = 2 cos t sen t, y = 2 sen t cos t
2. x + x 5y = 0, y + 4x + 5y = 0, x(0) = 1, y(0) = 2 R: x = e3t (2 sen 4t cos 4t), y = 2e3t cos 4t
3. x + 3y + y = et , y x = y, x(0) = 0, y(0) = 1 R: x = 3tet , y = 41 et + 34 et + 32 tet
4. x 3x 6y = 27t2 , x + y 3y = 5et x(0) = 5, y(0) = 1 R: x = 3et 9t2 + 6t + 2, y = et 6t
5. x = 2y, y = y x , x(0) = 0, x (0) = 10, y(0) = 5 R: x = 10(1 et ), y = 5et
6. y = x 2, x = y + 2, x(0) = 0, x (0) = 1, y(0) = 2, y (0) = 3 R: x = 2 sen t 3 cos t + et + 2,

y = 2 sen t + 3 cos t + et 2
1 1 1
7. x + y = cos t, x + y = 2, x() = 2, y(0) = 0, y (0) = 2 R: x = 2 + 2 sen t, y = 2 sen t
8. x = y + 4e2t , y = x e2t , x(0) = y(0) = 0, x (0) = y (0) = 0

1
R: x = sen t 2 cos t + e2t + 41 et 34 et , y = 1
2 cos t
9. x = x + y + z, y = 2x + 5y + 3z, z = 3x + 9y + 5z, x(0) = 2, y(0) = 1, z(0) = 3

R: x = 2, y = 1, z = 3
10. x = x + y, y = x y, z = 2y, x(0) = 2, y(0) = 0, z(0) = 0

R: x = 1 + 22 e 2t + 1 22 e 2t , y = 2
2
e 2t
e 2t
,z=e 2t
+ e 2t
2
Captulo 5
Soluciones usando series de potencias
En este captulo se estudia un metodo para resolver algunas ecuaciones cuya solucion es
imposible de hallar con los metodos estudiados anteriormente. Se supone que el lector conoce los
conceptos basicos de las series de potencias as que, en lo que sigue, unicamente se mencionaran
las operaciones necesarias. Se recomienda al lector repasar series de potencias1 .
5.1. Solucion de ecuaciones alrededor de puntos ordinarios

Para resolver la siguiente ecuacion con coeficientes variables
P (x)y + Q(x)y + R(x)y = 0, (5.1)

alrededor del punto x0 usando series de potencias, primero debemos verificar que P (x0 ) 6= 0. Si
esto no sucede (esto es, si P (x0 ) = 0), el metodo a utilizar es diferente. Esto se estudiara pos-
teriormente.
Suponiendo que ya se verifico lo anterior, ahora vamos a proponer una solucion de la forma

X
y(x) = an xn , (5.2)
n=0
con an valores a determinar para cada ecuacion particular. Para la serie propuesta las derivadas
son de la forma

X
X
n1

y (x) = an nx ,
y (x) = an n(n 1)xn2 . (5.3)
n=0 n=0
Al sustituir estas series en la ecuacion original, se tendran que hacer los ajustes necesarios
a las series involucradas y encontrar una formula de recurrencia que nos permita obtener los
coeficientes para cualquier valor de n. Esto nos permitira fijar los valores de las constantes an
y dar la solucion definitiva en forma de la serie propuesta con los valores an encontrados, como
se ilustra en el ejemplo siguiente.
Ejemplo
1
Vease por ejemplo el libro de G. B. Thomas, et al. Calculo de una variable, citado en la bibliografa.
86 5 Soluciones usando series de potencias
y + y = 0
empleando series.
Solucion
Aunque en esta ecuacion no hay coeficientes variables, la resolveremos para ilustrar el metodo.
Proponemos una solucion de la forma

X
y= an xn ,
n=0
cuya segunda derivada es

X
y = n(n 1)an xn2 ,
n=0
la cual se puede reescribir como

X
X
y = n(n 1)an xn2 = (n + 2)(n + 1)an+2 xn ,
n=2 n=0
puesto que los dos primeros coeficientes son cero.

Ahora, sustituyendo las series anteriores en la ecuacion dada obtenemos

X
X
(n + 2)(n + 1)an+2 xn + an xn = 0.
n=0 n=0
Para que estas series sumen cero, se necesita que la suma de los coeficientes valga cero para cada
valor de n, esto es
an+2 (n + 2)(n + 1) + an = 0,
que al despejar an+2 nos da la siguiente relacion
an
an+2 = ,
(n + 2)(n + 1)
llamada formula de recurrencia. Para hallar la solucion a la ecuacion vamos a dar valores a n.
Para n = 0,
a0
a2 =
12
n = 1,
a1
a3 =
23
n = 2,
a2 a0
a4 = =
34 1234
n = 3,
5.1 Solucion de ecuaciones alrededor de puntos ordinarios 87
a3 a1
a5 = = .
45 2345
Continuando con este proceso, encontramos que para n par (siendo n = 2k) se cumple
a0
a2k = (1)k ,
(2k)!
y para n impar (siendo n = 2k + 1) se cumple
a1
a2k+1 = (1)k .
(2k + 1)!
Lo anterior nos dice que la solucion general de la ecuacion que deseamos resolver es

X (1)k 2k
X (1)k 2k+1
y = a0 x + a1 x .
(2k)! (2k + 1)!
k=0 k=0
Si hubieramos resuelto la ecuacion por los metodos del captulo 2, habramos hallado como
solucion
y = A cos x + B sen x.
Podemos ver la equivalencia de ambas soluciones al recordar el desarrollo de las funciones seno y
coseno por medio de series de Taylor.
Ejemplo
Resolver la ecuacion con condiciones iniciales
y + xy + y = 0, y(0) = 0, y (0) = 1,
empleando series.
Solucion
Proponemos una solucion de la forma

X
y= an xn ,
n=0
cuya derivada es

X
y = nan xn1 ,
n=0
y la segunda derivada

X
y = n(n 1)an xn2 ,
n=0
las cuales se pueden reescribir como


X
y = nan xn1
n=1
y

X
y = n(n 1)an xn2 .
n=2
Sustituyendo en la ecuacion

X
X
X
n(n 1)an xn2 + x nan1 xn + an xn = 0,
n=2 n=1 n=0
que es equivalente a la ecuacion

X
X
X
n n
(n + 2)(n + 1)an+2 x + nan x + an xn = 0.
n=0 n=0 n=0
Para que estas series sumen cero necesitamos que
an+2 (n + 2)(n + 1) + an (n + 1) = 0,
que al despejar an+2 nos da la siguiente formula de recurrencia
an
an+2 = ,
n+2
y ahora le damos valores a n.
Para n = 0,
a0
a2 =
2
n = 1,
a1
a3 =
3
n = 2,
a2 a0
a4 = =
4 24
n = 3,
a3 a1
a5 = =
5 35
n = 4,
a4 a0
a6 = =
6 246
n = 5,
5.1 Solucion de ecuaciones alrededor de puntos ordinarios 89
a5 a1
a7 = =
7 357
n = 6,
a6 a0
a8 = =
8 2468
n = 7,
a7 a1
a9 = = .
9 3579
Continuando con este proceso, encontramos
a a a a
0 1 0 1
y = a0 + a1 x + x2 + x3 + x4 + +
2 3 24 35
a0 a1 a0 a1
x6 + x7 + x8 + x9 + ...
246 357 2468 3579
que se puede reescribir como

1 2 1 4 1 6 1 8
y = a0 1 x + x x + x ...
2 24 246 2468

1 3 1 5 1 7 1 9
+a1 x x + x x + x + ... .
3 35 3 5 67 3579
La solucion general de la ecuacion que deseamos resolver es

X x2k 2k X 2k k!x2k+1
y = a0 (1)k k
x + a1 (1)k+1 .
2 k! (2k + 1)!
k=0 k=0
Para encontrar la solucion particular que satisface las condiciones iniciales dadas, encontramos
y(0), lo cual nos da

X x2k X k
k+1 2 k!(0)
2k+1
y(0) = a0 (1)k (0)2k
+ a 1 (1) = a0 ,
2k k! (2k + 1)!
k=0 k=0
que nos indica que a0 = 0, puesto que y(0) = 0.
Para la otra condicion, primero hallamos la primera derivada

X 2k k!x2k
y = a1 (1)k+1 ,
(2k)!
k=0
que al sustituir nos da

X 2k k!(0)2k
y (0) = a1 (1)k+1 = a1 ;
(2k)!
k=0

como la condicion es y (0) = 1, concluimos que a1 = 1. La solucion buscada es

X 2k k!x2k
y = (1)k+1 .
(2k)!
k=0
Ejercicios
Resolver las siguientes ecuaciones usando series.
1. y y = 0, x0 = 0.
x2n x2n+1
P P
R: y = n=0 (2n)! + n=0 (2n+1)!
2. y xy y = 0, x0 = 0.
x2n 2n n!x2n+1
P P
R: y = n=0 2n n! + n=0 (2n+1)!
3. y + xy + 2y = 0, x0 = 0.
P (1)k x2k P (1)k 2k k! 2k+1
R: y = k=0 2k k! + k=0 (2k)! x
4. xy + y + xy = 0, x0 = 0, y(0) = 1, y (0) = 0.
P (1)k x2k
R: y = k=0 (k!)2 22k
5. (1 + x2 )y + 2xy = 0, x0 = 0, y(0) = 0, y (0) = 1.

P (1)k x2k+1
R: y = k=0 2k+1
6. (1 x2 )y xy = 0, x0 = 0 y(0) = 0, y (0) = 1.
P (2k)!x2k+1
R: y = k=0 22k k!(2k+1)
7. (1 x)y + y = 1 + x, x0 = 0, y(0) = 0.
xk+1
P
R: y = 1 + k=1 k(k+1)
8. y + k 2 x2 y = 0, x0 = 0, k una constante.
P (1)m+1 (k2 x4 )m+1
R: y = 1 + m=0 3478...(4m+3)(4m+4) +
h
(1)m+1 (k2 x4 )m+1
i
+x 1 +
P
m=0 4589...(4m+4)(4m+5)
9. (1 x)y + y = 0, x0 = 0.
x2 x3 x4 x3 x4
R: y = 1 2 6 24 + ... + x 6 24 + ...
10. (2 + x2 )y xy + 4y = 0, x0 = 0.
x4 x6 x3 x 5
19x7
R: y = 1 x2 + 6 30 + ... + x 4 + 7 160 1920 + ...
5.2. Solucion de ecuaciones alrededor de puntos singulares

regulares
Cuando en la ecuacion
d2 y dy
2
P (x)
+ Q(x) + R(x)y = 0,
dx dx
se requiere una solucion alrededor de un punto x0 en que P (x0 ) = 0, decimos que x0 es un
punto singular. Si en dicho punto singular se cumple que
Q(x)
lm (x x0 ) es finito
xx0 P (x)
R(x)
lm (x x0 )2 es finito
xx0 P (x)
5.2 Solucion de ecuaciones alrededor de puntos singulares regulares 91
se dice que tal punto singular es regular. El metodo que se describe a continuacion es obra
del matematico aleman F. G. Frobenius, por lo que tambien es conocido como metodo de
Frobenius. Existen tres casos del metodo de Frobenius, que se estudian por separado.
5.2.1. Races diferentes por un numero no entero

Consideraremos ecuaciones de la forma
d2 y dy
P (x) + Q(x) + R(x)y = 0, (5.4)
dx2 dx
restringiendonos al caso en que P (x) = x2 y Q(x) = xq(x); esto es, ecuaciones de la forma
2
2d y dy
x 2
+ xq(x) + r(x)y = 0, (5.5)
dx dx
siendo q(x) y r(x) funciones para las cuales existen expansiones en series de potencias en x,
digamos

X
q(x) = qm xm , (5.6)
m=0

X
r(x) = rm xm , (5.7)
m=0
convergentes en algun intervalo que contenga al punto x = 0.

La idea fundamental en el metodo de Frobenius es probar una solucion para la ecuacion
5.5 de la forma

X
s 2 n
z(x, s) = x (a0 + a1 x + a2 x + ... + an x + ...) = an xn+s , (5.8)
n=0
con a0 6= 0. Esto siempre es factible, puesto que si a0 = 0, simplemente tendramos que ahora
a1 juega el papel de a0 , con otro valor de s.
De la solucion propuesta (ecuacion 5.8) tendremos que

dz X
= an (n + s)xn+s1 ,
dx n=0
y

d2 z X
= an (n + s)(n + s 1)xn+s2 ,
dx2 n=0
o sea que si queremos que se satisfaga la ecuacion 5.5, se debe cumplir que
d2 z dz
x2 2
+ xq(x) + r(x)z = 0,
dx dx
lo cual, al sustituir z y sus derivadas implica que se debe cumplir la ecuacion

X
X
X
n+s2 n+s1
an (n + s)(n + s 1)x + an (n + s)x + an xn+s = 0,
n=0 n=0 n=0
la cual, luego de utilizar los desarrollos 5.6 y 5.7 para q(x) y r(x) y dividir entre xs se
transforma en

X
X
X
n+s2 n+s1
an (n + s)(n + s 1)x + an (n + s)x + an xn+s = 0. (5.9)
n=0 n=0 n=0
Para que la suma de las series infinitas nos de cero, es necesario que todos los coeficientes
involucrados en cada potencia de x den cero por separado. Esto da origen a un conjunto de
ecuaciones como se vera enseguida.
Si queremos que el coeficiente de x0 sea cero, debemos elegir m = 0 y n = 0, lo cual da la
ecuacion
a0 s(s 1) + a0 q0 s + a0 r0 = 0. (5.10)
Ahora, para que el coeficiente de x1 sea cero, necesitamos que en la primera serie n = 1,
mientras que en la segunda y tercera series se requiere que n + m = 1, lo cual puede suceder
si m = 0 y n = 1, o bien si m = 1 y n = 0. Esto nos da la ecuacion
a1 (s + 1)s + a1 q0 (s + 1) + a0 q1 s + a1 r0 + a0 r1 = 0. (5.11)
En general, para que el coeficiente de cualquier potencia de x, digamos xi , sea cero, es
necesario que n + m = i, lo cual puede suceder para diferentes valores de n y m (por ejemplo,
n = i y m = 0, o bien n = i 1 y m = 1, etc.). Con esta condicion se obtiene que
ai (s + i)(s + i 1) + {ai1 q1 (s + i 1) + ai2 q2 (s + i 2) + ...+
+a0 qi s} + ai r0 + ai1 r1 + ai2 r2 + ... + a0 ri = 0, (i 1). (5.12)

Aqu juntaremos todos los coeficientes de ai e identificaremos a su coeficiente con
f (s + i) = (s + i) (s + i)(s + i 1) + q0 (s + i) + r0 , (5.13)
con lo cual, la ecuacion 5.12 queda como
ai f (s + i) + {{ai1 q1 (s + i 1) + ai2 q2 (s + i 2) + ...+
+a0 qi s} + ai1 r1 + ai2 r2 + ... + a0 ri , (i 1). (5.14)

Regresando a la ecuacion 5.10, tras una simplificacion queda como
a0 s2 + (q0 1)(s + r0 ) = 0,
y recordando que a0 6= 0, tendremos que
s2 + (q0 1)(s + r0 ) = 0. (5.15)

Esta ecuacion, llamada ecuacion indicial, es cuadratica en s, as que tiene dos races, las
cuales simbolizaremos como s1 y s2 , ademas de convenir en ordenarlas de tal forma que s2 s1 .
Aqu hemos supuesto que las races de la ecuacion indicial son reales.
La ecuacion 5.15 es simplemente la ecuacion f (s) = 0, as que tenemos que
f (s) = (s s1 )(s s2 ). (5.16)

El conjunto de ecuaciones 5.14 se puede utilizar ahora para determinar los coeficientes a1 ,
a2 , ... en terminos de a0 . Por ejemplo, la ecuacion 5.14, con i = 1 nos da
a1 f (s + 1) + a0 q1 s + a0 r1 = 0,
as que tenemos que
a0 (q1 s + r1 )
a1 =
f (s + 1)
a0 h1 (s)
= , (5.17)
f (s + 1)
donde h1 (s) = (q1 s + r1 ) es un polinomio de primer grado en s.
La ecuacion 5.14 con i = 2 da
a2 F (s + 2) + a1 q1 (s + 1) + a0 q2 s + a1 r1 + a0 r2 = 0
que, sustituyendo en la ecuacion 5.17 para determinar a2 nos da

a0 h1 (s) a0 h1 (s)
a2 f (s + 2) + q1 (s + 1) + a0 q2 s + r1 + a0 r2 = 0,
f (s + 1) f (s + 1)
dandonos
a0 q1 (s + 1)h1 (s) + q2 sf (s + 1) + r1 h1 (s) + r2 f (s + 1)

a2 =
f (s + 1)f (s + 2)
h2 (s)
= a0 , (5.18)
f (s + 1)f (s + 2)
donde h2 (s) = q1 (s + 1)h1 (s) + q2 sf (s + 1) + r1 h1 (s) + r2 f (s + 1) es un polinomio que de-
pende de s.
En general, para el coeficiente ai deberemos tener que
hi (s)
ai = a0 , (5.19)
f (s + 1)f (s + 2)...f (s + i)
siendo hi (s) un polinomio de s.
Si en la ecuacion 5.19 sustituimos s = s1 obtenemos expresiones para ai en terminos de a0 ,
lo cual nos da una solucion que involucra a a0 como una constante arbitraria de acuerdo a lo
que sabemos que funciona para ecuaciones de segundo orden. Al hacer esto mismo para s = s2
deberamos tener otra solucion, que nos dara la segunda solucion linealmente independiente.
Sin embargo, esto solo es cierto si el denominador en la ecuacion 5.19 es diferente de cero, esto
es, si sucede que
f (s1 + i) 6= 0
y
f (s2 + i) 6= 0, (5.20)
siendo i cualquier numero entero positivo.
Ahora bien, la ecuacion 5.16 nos dice que
f (s) = (s s1 )(s s2 ),
de manera que
f (s + i) = (s + i s1 )(s + i s2 ),
y entonces
f (s1 + i) = i(s1 s2 + i),

y
f (s2 + i) = (s2 s1 + i)i, (5.21)

o sea que la ecuacion 5.20 se satisface si, y solo si s2 s1 no es un entero positivo, esto es, si las
races de la ecuacion indicial no difieren por un numero entero (tambien, como caso particular,
si no son iguales).
As pues, cuando s2 s1 no es un numero entero (o cero) se tienen dos soluciones lin-
ealmente independientes de la ecuacion 5.5. En caso de que s2 s1 sea un numero entero o
cero (races iguales), necesitamos otro metodo para construir la segunda solucion linealmente
independiente.
Ejemplo
2xy + y + y = 0.
Solucion
Suponiendo que

X
X
X
z= an xn+s , z = an (n + s)xn+s1 , z = (n + s)(n + s 1)an xn+s2
n=0 n=0 n=0
tenemos la ecuacion

X
X
X
2an (n + s)(n + s 1)xn+s1 + an (n + s)xn+s1 + an xn+s = 0,
n=0 n=0 n=0
para la cual, si queremos que los coeficientes se hagan cero, primero debemos tener en fase las series, esto es,
debemos hacer que los exponentes sean iguales, prefiriendo el de orden menor (en este caso, n + s 1), con lo
cual tendremos

X
X
X
2an (n + s)(n + s 1)xn+s1 + an (n + s)xn+s1 + an1 xn+s1 = 0.
n=0 n=0 n=1
0
Aqu, si queremos que el coeficiente de x sea cero, se debe cumplir la ecuacion anterior para n = 0, o sea:
2a0 (0 + s)(0 + s 1) + a0 (0 + s) = 0,
la cual, despues de dividir entre a0 (que nunca es cero) y reducir se convierte en
2s2 s = 0.
La anterior es la ecuacion indicial, que tiene como races s1 = 0 y s2 = 12 , las cuales difieren por un numero no
entero. Tenemos que para n 1 se cumple que
2an (n + s)(n + s 1) + an (n + s) + an1 = 0,

con lo cual obtenemos la formula de recurrencia
an1
an = .
(n + s)[2(n + s) 1]
Sustituyendo algunos valores de n tenemos que
a0
a1 = ,
(1 + s)(1 + 2s)
a1 a0
a2 = = ,
(2 + s)(3 + 2s) (1 + s)(2 + s)(1 + 2s)(3 + 2s)
a2 a0
a3 = = ,
(3 + s)(5 + 2s) (1 + s)(2 + s)(3 + s)(1 + 2s)(3 + 2s)(5 + 2s)
con los cuales podemos generalizar como sigue
a0
an = (1)n .
(1 + s)(2 + s)...(n + s)(1 + 2s)(3 + 2s)...(2n + 2s 1)
Si en la formula anterior tomamos s = s1 = 0 obtenemos la solucion
x2 x3

x
z(x, 0) = a0 1 + + ...+
1 1 1 2 1 3 1 2 3 1 3 5
xn

(1)n + ... .
1 2 3...n 1 3 5...(2n 1)
Si ahora tomamos s = s2 = 12 , tenemos la segunda solucion
x2 xn

x n
z(x, 0) = a0 1 3 + 3 5 + ... + (1) 3 5 2n+1 + ... .
2 2 2 2 24 2 2 ... 2 2 4 6..,2n
Finalmente, la solucion general es la suma de las dos soluciones anteriores, que se puede escribir en forma
mas compacta como
n
X
n (2x) 1/2
X (2x)n
y(x) = A (1) + Bx (1)n ,
n=0
(2n)! n=0
(2n + 1)!
siendo A y B constantes arbitrarias.
5.2.2. Races diferentes por un numero entero

Si las races de la ecuacion indicial difieren por un numero entero, el procedimiento anterior
solo funciona para la raz mayor, esto es, para s = s2 , segun veremos. Si s2 s1 = r, de las
ecuaciones 5.21 tenemos que
f (s1 + i) = i(r + i) (5.22)
f (s2 + i) = i(r + i), (5.23)

y vemos que f (s1 + i) = 0 para i = r, mientras que f (s2 + i) 6= 0 para cualquier entero positivo
i, de manera que el procedimiento solo sirve para encontrar una solucion.
Para encontrar una segunda solucion de la ecuacion considerada, usamos la relacion 5.19
antes de fijar el valor de s, con lo cual tenemos que

s h1 (s) h2 (s)
z(x, s) = a0 x 1 + x+ + ...+
f (s + 1) f (s + 1)f (s + 2)x2

hi (s) i
+ x + ... . (5.24)
f (s + 1)f (s + 2)...f (s + i)
Sabemos por el metodo de construccion de a1 , a2 ,... que, si todos los terminos de el lado
derecho estan bien definidos, deben cumplir la ecuacion
d2 z dz
x2 2
+ xq(x) + r(x)z = a0 f (s)xs .
dx dx
Esto es porque, como las ecuaciones 5.14 se satisfacen, los coeficientes de cada potencia de xs+1
con i 1 se deben anular; y el coeficiente de xs es, por las ecuaciones 5.10 y 5.13, simplemente
a0 f (s). Esta ecuacion se puede reescribir en la forma
2

2 d d
x + xq(x) + r(x) z(x, s) = a0 f (s)xs . (5.25)
dx2 dx
Por desgracia, los terminos en el lado derecho de la ecuacion 5.25 no estan definidos para
s = s1 , dado que los denominadores se anulan para i r. Sin embargo, si multiplicamos z(x, s)
por f (s + r) cancelamos el factor en los denominadores de los ultimos terminos que se anulan
cuando s = s1 . Y como f (s + r) = (s + r s1 )(s + r s2 ) = (s + s2 2s1 )(s s1 ) simplemente
debemos tambien multiplicar por (s s1 ). Tambien, como este factor es independiente de x,
la ecuacion 5.25 toma la forma
2

2 d d
x + xq(x) + r(x) (s s1 )z(x, s) = a0 (s s1 )f (s)xs
dx2 dx
= a0 (s s1 )2 (s s2 )xs , (5.26)
usando la ecuacion 5.16.
Usando s = s2 , tenemos
2

2 d d
x + xq(x) + r(x) (s2 s1 )z(x, s2 ) = 0,
dx2 dx
lo cual muestra que z(x, s2 ) es una solucion.
Similarmente, con s = s1 tenemos [(s s1 )z(x, s)]s=s1 como solucion. Sin embargo, resulta
que esta serie no es independiente de z(x, s2 ), sino que es un multiplo de ella. Esto se puede
ver como sigue. El factor s s1 cancela los ceros en el denominador para los terminos con
i r, pero dara ceros en el numerador para todos los terminos con i < r; entonces, la primera
potencia en la serie es simplemente xs1 +r = xs2 , que es precisamente el primer termino en
z(x, s2 ). Y como usamos la misma regla para calcular cualquier coeficiente en terminos de los
precedentes en ambos casos, obtenemos la misma serie en ambos casos, salvo por una constante
multiplicativa.
Si derivamos ambos lados de la ecuacion 5.26 con respecto a s, obtenemos
2

2 d d d
x + xq(x) + r(x) (s s1 )z(x, s) =
dx2 dx ds

2 d s 2
= a0 (s s1 ) (s s2 )x + 2(s s1 )(s s2 )x (5.27)
ds
y vemos que el lado derecho de esta ecuacion es cero cuando s = s1 , as que debemos tener
[(d/ds)(s s1 )z(x, s)]s=s1 como solucion de la ecuacion diferencial. Para verificar que esta se-
gunda solucion es linealmente independiente de la primera, es suficiente con calcular el wron-
skiano y observar que no se anula.
Esto implica que las soluciones linealmente independientes de la ecuacion dada son
[(s s1 )z(x, s)]s=s1 (5.28)

y

d
(s s2 )z(x, s) . (5.29)
ds s=s1
Cuando las races difieren por un numero entero puede surgir otra situacion. Como f (s+i) =
0 para s = s1 e i = r, puede suceder que hr (s1 ) = 0. En este caso ar esta indeterminada (puesto
que tiene cero tanto en el numerador como en el denominador), as que podemos usarla como
otra constante arbitraria, obteniendo as las dos soluciones independientes de la primer serie.
Esto se vera en un ejemplo particular.
Ejemplo
x2 y 3xy + (3 x)y = 0.
Solucion
Suponiendo que

X
X
X
z= an xn+s , z = an (n + s)xn+s1 , z = (n + s)(n + s 1)an xn+s2
n=0 n=0 n=0
tenemos la ecuacion

X
X
X
X
an (n + s)(n + s 1)xn+s 3 an (n + s)xn+s + 3 an xn+s an xn+s+1 = 0,
n=0 n=0 n=0 n=0
que al enfasar las series nos da

X
X
X
X
an (n + s)(n + s 1)xn+s 3 an (n + s)xn+s + 3 an xn+s an1 xn+s = 0.
n=0 n=0 n=0 n=1
0
Si queremos que el coeficiente de x sea cero, necesitamos que
a0 (0 + s)(0 + s 1) 3a0 (0 + s) + 3a0 = 0,

la cual, al dividir entre a0 y simplificar se convierte en
s2 4s + 3 = 0.
Esta es la ecuacion indicial, que tiene como races a s1 = 1 y s2 = 3. Estas races difieren por un entero,
as que necesitamos encontrar la solucion para s = s2 , y luego construir la segunda solucion en alguna de las
formas indicadas anteriormente.
Para n 1 tenemos que se debe cumplir que
an (n + s 1)(n + s 3) an1 = 0,
lo cual nos da la formula de recurrencia
an1
an = .
(n + s 1)(n + s 3)
En la ecuacion anterior se ve que todas las an para n 1 han quedado bien determinadas, por lo cual la
segunda solucion se debera hallar derivando con respecto a s la solucion z(x, s). Mientras tanto, demos algunos
valores a n, a fin de encontrar la primera solucion
a0
a1 = ,
s(s 2)
a1 a0
a2 = = ,
(s + 1)(s 1) s(s + 1)(s 2)(s 1)
a2 a0
a3 = = ,
(s + 2)s s(s + 1)(s + 2)(s 2)(s 1)s
y en general, para cualquier n tendremos
a0
an = ,
s(s + 1)(s + 2)...(s + n 1)(s 2)(s 1)s...(s + n 3)
que nos da para z(x, s) la expresion

( )
s
X xn
z(x, s) = a0 x 1 + .
n=1
s(s + 1)(s + 2)...(s + n 1)(s 2)(s 1)(s...(s + n 3)
La primera solucion linealmente independiente es [(s 1)z(x, s)]s=1 , esto es

X xn X xn+1
y1 (x) = a0 x = a0 .
n=2
1 2...n(1) 1 2...(n 2) n=2
n!(n 2)!
Si ahora tomamos la derivada con respecto a s de (s 1)z(x, s), y evaluamos para s = s1 = 1, tendremos
la otra solucion buscada

d d s (s 1)
[(s 1)z(x, s)] = a0 x (s 1) + x+
ds ds s(s 2)

)
xn

X d (s 1)
+ a0 xs 1 + x+
n=2
s(s + 1)...(s + n 1)(s 2)s(s + 1)...(s + n 3) ds s(s 2)

)
X d 1
+ xn .
n=2
ds s(s + 1)...(s + n 1)(s 2)s(s + 1)...(s + n 3)
Para calcular la derivada que aparece dentro de la suma usaremos la formula de derivacion logartmica en
la forma siguiente. Llamemos al termino a derivar g(s), calculemos su logaritmo y derivemos
d 1 d
ln g(s) = g(s),
ds g(s) ds
as que, despejando la derivada de g(s) en la ecuacion anterior, tenemos que
d
g(s) = g(s) ln g(s).
ds
Por lo tanto, para calcular la derivada es necesario calcular el logaritmo de g(s), que nos da
ln g(s) = ln s ln(s + 1) ln(s + 2)... ln(s + n 1)
ln(s 2) ln s ln(s + 1)... ln(s + n 3),

ahora derivemos esta expresion

d 2 2 2 1 1 1
ln g(s) = + + ... + + + + ,
ds s s+2 s+n3 s2 s+n2 s+n1
y evaluando en s = s1 = 1 obtenemos

d 1 1 1 1 1
ln g(s) = 2 1 + + + ... + 1+ +
ds s=1 2 3 n2 n1 n
n2
( )
1 1 X
= + 1+2 .
n n1 m=1
Tambien evaluamos g(s) en s = 1, obteniendo
1
g(s = 1) =
1 2...n (1) 1 2 ...(n 2)
1
= ,
n!(n 2)!
que, con el resultado anterior, nos da

n2
!
d 1 1 1 X 1
g(s) = + 1+2 .
ds s=1 n!(n 2)! n n1 m=1
m
Entonces, con el resultado anterior y las derivadas de los otros terminos, tendremos que la expresion
d

ds {(s
1)z(x, s)} s=1 es igual a
n2
( ! )
1 X 1 1 X 1
n
a0 (ln x)y1 (x) + a0 x 1 x + + 1+2 x .
n=2
n!(n 2)! n n1 m=1
m
Finalmente, la solucion general es la suma de ambas soluciones encontradas.
Ejemplo
(x2 1)y + 3xy + xy = 0.

Solucion
Suponemos

X
X
z= an xn+s , z = an (n + s)xn+s1 ,
n=0 n=0

X
z = (n + s)(n + s 1)an xn+s2 ,
n=0
con lo cual tenemos la ecuacion

X
X
an (n + s)(n + s 1)xn+s an (n + s)(n + s 1)xn+s2 +
n=0 n=0

X
X
3 an (n + s)xn+s + an xn+s+1 = 0,
n=0 n=0
que al enfasar las series nos da

X
X
an2 (n + s 2)(n + s 3)xn+s2 an (n + s)(n + s 1)xn+s2 +
n=2 n=0

X
X
+3 an2 (n + s 2)xn+s2 + an3 xn+s2 = 0.
n=2 n=3
0
Para que el coeficiente de x sea cero, necesitamos que
a0 (0 + s)(0 + s 1) = 0,
que al dividir entre a0 y simplificar nos da
s2 s = 0.
Esta es la ecuacion indicial, que tiene como races s1 = 0 y s2 = 1, y vemos que difieren por un numero
no entero. Tambien notamos que para n = 1 solo hay un termino, lo cual no permite determinar el valor de la
constante a1 en terminos de a0
a1 (s + 1)s = 0,
lo que hace que podamos utilizar a a1 como constante arbitraria para la segunda solucion linealmente
independiente. Para n = 2 tenemos que
a0 s(s + 2) a2 (s + 2)(s + 1) = 0,
mientras que para valores de n 3 se cumple
an3 + an2 (s + n 2) an (s + n)(s + n 1) = 0,

lo cual nos lleva a la formula de recurrencia
an3 + n(n 2)an2

an = , (n 3),
n(n + 1)
donde ya se ha sustituido s = s1 = 0, puesto que solo necesitamos usar una relacion para fijar ambas soluciones.
Dando valores a n en las formulas encontradas tenemos que
a2 = 0,
1 1
a3 = a0 + a1 ,
6 2
1 2 1
a4 = a1 + a2 = a1 ,
12 3 12
1 3
a5 = a0 + a1 .
8 8

1 1 1 1 3
y(x) = a0 1 + x3 + x5 + ... + a1 x + x3 + x4 + x5 + ... .
6 8 2 12 8
Notese que en este caso no es sencillo encontrar una forma mas compacta para escribir la solucion, por lo
cual solo escribimos los primeros terminos de la serie correspondiente.
5.2.3. Races iguales

El ultimo caso a considerar es cuando ambas races de la ecuacion indicial son iguales, esto
es, s1 = s2 , y si utilizamos la ecuacion 5.25 (puesto que este caso es un caso particular del
anterior) obtenemos
2

2 d d
x + x(q(x) + r(x) z(x, s) = a0 (s s1 )2 xs . (5.30)
dx2 dx
Si usamos s = s1 veremos que el lado derecho vale cero, as que z(x, s) es una solucion.
Tambien, si derivamos ambos miembros de la ecuacion con respecto a s, obtenemos
2

2 d d d
x + xq(x) + r(x) z(x, s) =
dx2 dx ds

s 2 d s
= a0 2(s s1 )x + (s s1 ) x . (5.31)
ds
Aqu nuevamente el lado derecho es cero cuando s = s1 , as que tenemos como solucion
[(d/ds)z(x, s)]s=s1 . Esta solucion es independiente de z(x, s) as que, para este caso, las dos
soluciones son
z(x, s1 ) (5.32)
y

d
z(x, s) . (5.33)
ds s=s1
Ejemplo
x(1 x)y + (1 x)y y = 0.

Solucion

X
X
n+s
z= an x , z = an (n + s)xn+s1 ,
n=0 n=0

X
z = (n + s)(n + s 1)an xn+s2
n=0
con lo cual tenemos la ecuacion

X
X
an (n + s)(n + s 1)xn+s1 an (n + s)(n + s 1)xn+s +
n=0 n=0

X
X
+ an (n + s)xn+s1 an (n + s)xn+s
n=0 n=0

X
an xn+s = 0,
n=0
enfasando las series obtenemos

X
X
an (n + s)(n + s 1)xn+s1 an1 (n + s 1)(n + s 2)xn+s1 +
n=0 n=1

X
X
+ an (n + s)xn+s1 an1 (n + s 1)xn+s1
n=0 n=1

X
an1 xn+s1 = 0.
n=1
0
Para que el coeficiente de x sea cero, necesitamos que
a0 s(s 1) + a0 s = 0,
que al dividir entre a0 y simplificar nos da
s2 = 0.
Esta es la ecuacion indicial, que tiene como races s1 = s2 = 0, as que procederemos como se indico antes.
Primero vamos a encontrar z(x, s). Cuando n 1 tenemos
an (s + n)2 an1 [(s + n 1)2 + 1],

lo que lleva a la formula de recurrencia
[(s + n 1)2 + 1]
an = an1 .
(s + n)2
Si damos valores a n obtenemos
(s2 1)
a1 = a0 ,
(s + 1)2
[(s + 1)2 + 1] [s2 + 1][(s + 1)2 + 1]

a2 = a1 = a 0 ,
(s + 2)2 (s + 1)2 (s + 2)2
[(s + 2)2 + 1] [s2 + 1][(s + 1)2 + 1][(s + 2)2 + 1]

a3 = a2 2
= a0 ,
(s + 3) (s + 1)2 (s + 2)2 (s + 3)2
as que para el caso general tenemos
[s2 + 1][(s + 1)2 + 1][(s + 2)2 + 1]...[(s + n 1)2 + 1]

an = a0 .
(s + 1)2 (s + 2)2 ...(s + n)2
Con esta ecuacion, obtenemos para z(x, s)

" #
s
X [s2 + 1][(s + 1)2 + 1]...[(s + n 1)2 + 1] n
z(x, s) = a0 x 1 + x ,
n=1
(s + 1)2 (s + 2)2 ...(s + n)2
de donde obtenemos de inmediato la primer solucion

" #
s
X 1[12 + 1]...[(n 1)2 + 1] n
y1 (x) = z(x, 0) = a0 x 1 + x =
n=1
12 22 ...n2

" #
2
X 1 2 5 10 17...[(n 1) + 1]
a0 xs 1 + 2
xn .
n=1
(n!)
d
Para la segunda solucion necesitamos hallar ds z(x, s). Esto nos lleva a lo siguiente

" #
[s2 + 1][(s + 1)2 + 1]...[(s + n 1)2 + 1] n

d d s X
z(x, s) = a0 x 1+ x +
ds ds n=1
(s + 1)2 (s + 2)2 ...(s + n)2

d [s2 + 1][(s + 1)2 + 1]...[(s + n 1)2 + 1] n
X
+a0 xs x ,
n=1
ds (s + 1)2 (s + 2)2 ...(s + n)2
y nuevamente tenemos que utilizar la tecnica de derivacion logartmica para el termino dentro de la suma g(s),
que segun vimos antes, nos da la derivada en la forma
d 1 d
g(s) = ln g(s).
ds g(s) ds
Calculando el logaritmo natural de g(s) obtenemos
n
X n
X
ln g(s) = ln[(s + m 1)2 + 1] ln(s + m)2 ,
m=1 m=1
y su derivada es
n n
d X 2(s + m 1) X 2
ln g(s) = 2
=
ds m=1
(s + m 1) + 1 m=1 s + m
n
X s+m1 1
=2 ,
m=1
(s + m 1)2 + 1 s + m
as que
n
2[s2 + 1]...[(s + n 1)2 + 1] X

d s+m1 1
g(s) = .
ds (s + 1)2 ...(s + n)2 m=1
(s + m 1)2 + 1 s + m
Al evaluar en s = 0 la derivada anterior obtenemos
n
2(1)(12 + 1)(22 + 1)...[(n 1)2 + 1] X

d m1 1
g(s) = ,
ds s=0 12 22 ...n2 m=1
(m 1)2 + 1 m
con lo cual tenemos que

" #
[1][12 + 1][22 + 1]...[(n 1)2 + 1]

d X
y2 (x) = z(x, s) a0 ln x 1 + +
ds s=0 n=1
12 22 32 ...n2
n
X 2 1 2 5 10 17...[(n 1)2 + 1] X m2
+a0 .
n=1
(n!)2 m=1
m[(m 1)2 + 1]
y(x) = Ay1 (x) + By2 (x).
Ejercicios
Resolver las ecuaciones siguientes usando el metodo de Frobenius
1. 4xy + y y = 0
h P i h P i
1 n 1
R: y(x) = A 1 + n=1 n!159...(4n3) x + Bx3/4 1 + n=1 n!71115...(4n+3) x
n
2. xy + (x + 1)y + y = 0
h
(1)n
i
n xn n
P P P
ln x n=0 (1)n xn! n=0 1 1 1 1
xn

R: y(x) = A n=0 (1) n! + B n! 1+ 2 + 3 + 4 + ... + n
3. xy 2y + y = 0
hP i h
x2 x3
n 1 n
P 1 x
R: y(x) = A n=3 (1) 2n!(n3)! x + B ln x n=3 (1)n 2n!(n3)! +1+ 2 + 4 36 +
n
P o i
(1) 1 1 1 1 1 1 3
xn

+ n=4 2n!(n3)! (2) 1 + 2 + 3 + ...+ + n3 n2 n1 n + 2
4. 9x(1 x)y 12y + 4y = 0

h P i h n P oi
147...(3n2) n 81114...(3n+5) n
R: y(x) = A 1 + n=1 3n n! x + B x7/3 1 + n=1 101316...(3n+7) x
5. (1 x2 )y + 2xy + y = 0
h
(1)(3)...(4n2 14n+9) 2n
i h
(3)(1)...(4n2 10n+3) 2n+1
i
R: y(x) = A 1 + +B x+
P P
n=1 (2n)! x n=1 (2n+1)! x
6. x2 y + (x2 3x)y + (4 2x)y = 0

(1)n n
h P i
R: y(x) = Ax2 + B x2 ln x + x2 n=1 nn! x
7. (1 x2 )y + xy y = 0
h
32 52 ...(2n3)2 n
i
x2 x4
P
R: y(x) = Ax + B 1 + 2! + 4! + n=3 (2n)! x
8. x2 y 3xy + (x3 5)y = 0

hP i h
n+1 1 3n1
P 1
R: y(x) = A n=2 (1) 32n1 (n2)!n! x + B ln x n=2 (1)n+1 32n1 (n2)!n! x3n1 + x1 +
o i
1
+ 91 x2 + 324 x5 + n+1 1 1 1 1 1 3n1
P
n=2 (1) 32n1 (n2)!n! {1 2 1 + 2 + ... + n2 n1 n x
9. y x2 y y = 0
R: y(x) = A 1 + 12 x2 + 1 4 1 5 1 6 13 7
+ ... +B x + 61 x3 + 1 4 1 5 7 6 41 7

24 x + 20 x + 720 x + 2520 x 12 x + 120 x + 360 x + 5040 x + ...
10. x2 y + 2xy + xy = 0
(1)n xn (1)n xn
hP i h P P i
1 Hn +Hn1 n n
R: y(x) = A n=0 n!(n+1)! + B ln x n=0 n!(n+1)! + x 1 n=1 n!(n1)! (1) x ,
1 1 1
con Hn = n + n1 + ... + 2 +1
Tablas
Integrales
R
R sen xdx = cos x + C
R cos xdx = sen x + C
sec2 xdx = tg x + C
csc2 xdx = ctg x + C
R
R
R sec x tg xdx = sech x + C
R csch xctg xdx = csch x + C
R tg xdx = ln(cos x) + C
R ctg xdx = ln(sen x) + C
R sec xdx = ln(sec x + tg x) + C
csc xdx = ln(csc x ctg x) + C
1
R
1x2
dx = arc sen x + C
1
R
1x2 dx = arc cos x + C
1
R
1+x2
dx = arc tg x + C
1
R
1+x2 dx = arcctg x + C
1
R
x x2 1
dx = arcsec x + C
xx12 1 dx = arccsc x + C
R
Transformadas de Laplace
L{1} = 1s
n!
L{tn } = sn+1
1
L{eat } = sa
n!
L{tn eat } = (sa) n+1
a
L{sen at} = s2 +a2
s
L{cos at} = s2 +a 2
b
L{senh bt} = s2 b2
s
L{cosh bt} = s2 b 2
dn
L{tn f (t)} = (1)n ds n F (s)
L{u(t a)f (t a)} = eas F (s)

L{eat f (t)} = F (s a)
L{(t a)} = eas
L{f (t)} = sF (s) f (0)
L{f (t)} = s2 F (s) sf (0) f (0)
108 5 Tablas
Desarrollos en serie de algunas funciones

2 3 n xn
ex = 1 + x + x2! + x3! + ... + xn! + ... =
P
n=0 n!
3 5 x2n1 n1 x2n1
sen x = x x3! + x5! ... + (1)n1 (2n1)! + ... =
P
n=0 (1) (2n1)!
x2 x4 n x2n n x2n
P
cos x = 1 2! + 4! ... + (1) (2n)! + ... = n=0 (1) (2n)!
2 3 n1
n xn1
ln(1 + x) = x x2 + x3 ... + (1)n xn1 + ... =
P
n=0 (1) n1
(1 + x) = 1 + x + (1) 2!
x2 + ... + (1)...(n+2)
(n1)!
xn1 + ... =
P (1)...(n+2) n1
n=0 (n1)!
x
Bibliografa
1. Amelkin, V. V. Ecuaciones diferenciales en la practica. URSS. Moscu, 2003.
2. Boyce, W. E. y Richard C. DiPrima Ecuaciones diferenciales y problemas con

valores en la frontera. 4a ed. Limusa. Mexico, 1998.
3. Coddington, E. A. Ecuaciones diferenciales. CECSA. Mexico, 1974.
4. Elsgoltz, L. Ecuaciones diferenciales y calculo variacional. Mir. Moscu, 1969.
5. Ince, E. L. Ordinary differential equations. Dover. Nueva York, 1956.
6. Ince, E. L. Integration of ordinary differential equations. Oliver & Boyd. Edimburgo,

1956.
7. Kaplan, W. Ordinary differential equations. Addison Wesley. Nueva York, 1958.
8. Krasnov, et al. Curso de Matematicas Superiores. 2a ed. URSS. Moscu, 2003.
9. Kiseliov, et al. Problemas de ecuaciones diferenciales ordinarias. Mir. Moscu, 1979.
10. Ledder, G. Ecuaciones diferenciales. McGraw Hill. Mexico, 2006.
11. Marcus, D. A. Ecuaciones diferenciales. CECSA. Mexico, 2004.
12. Martin, R. Ordinary differential equations. McGraw Hill. Nueva York, 1985.
13. Nagle, R. K. et al. Ecuaciones diferenciales. 4a ed. Pearson - Addison Wesley.

Mexico, 2005.
14. Piskunov, N. Calculo diferencial e integral. Mir. Moscu, 1977.
15. Spiegel, M. Ecuaciones diferenciales aplicadas. Prentice Hall. Mexico, 1981.
16. Zill, D. G. Ecuaciones diferenciales con problemas de valores en la frontera. 6a ed.

Thompson. Mexico, 2006.
.
Esta obra se termin de imprimir
en el mes de mayo de 2008
en el taller de impresin de la
Universidad Autnoma de la Ciudad de Mxico
con un tiraje de 2000 ejemplares.

Fausto Ecuaciones PDF

Cargado por

Copyright:

Formatos disponibles

Fausto Ecuaciones PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Fausto Ecuaciones PDF

Cargado por

Copyright:

Formatos disponibles

Integracion de

Fausto Cervantes Ortiz

Fausto Cervantes Ortiz

ECUACIONES DIFERENCIALES ORDINARIAS

Fausto Cervantes Ortiz

D.R. Universidad Autnoma de la Ciudad de Mxico

Academia de Matemticas del Colegio de Ciencia y Tecnologa de la UACM

Diseo de Portada: Aarn Aguilar

Material educativo universitario de distribucin gratuita. Prohibida su venta.

Nada humano me es ajeno

Fausto Cervantes Ortiz

1. Ecuaciones de primer orden 1

2. Ecuaciones de segundo orden 33

3. Sistemas de ecuaciones lineales 57

5. Soluciones usando series de potencias 85

Ecuaciones de primer orden

Una ecuacion diferencial es una ecuacion en la cual intervienen derivadas de funciones

1.1. Ecuaciones separables

M(x)dx + N(y)dy = 0. (1.1)

otro miembro de la ecuacion las funciones dependientes de x, tambien pasamos multiplicando

La ecuacion dada es equivalente a

Esta ecuacion es equivalente a

Para eliminar la constante arbitraria en la solucion de una ecuacion diferencial de primer

Primero separamos variables

2(y 1)dy = (3x2 + 4x + 2)dx,

Resolver la ecuacion dada, sujeta a la condicion inicial indicada.

Separando variables obtenemos

Con este valor de C obtenemos finalmente

10. 2x(1 + ey ) dx + ey (1 + x2 ) dy = 0 R: 1 + ey = C(1 + x2 )

13. (xy 2 y 2 + x 1) dx + (x2 y 2xy + x2 + 2y 2x + 2) dy = 0 R: (x2 2x + 2)(y 2 + 1)e2 arc tg y = C

15. y = ax + by + c (a, b, c constantes) R: a + b(ax + by + c) = Kebx

19. xy 2 (xy + y) = a2 R: 2x3 y 3 = 3a2 x2

22. (x2 y 3 + y + x 2) dx + (x3 y 2 + x) dy = 0 R: 3x2 12x + 2x3 y 3 + 6xy = C

(x+y)m (x+y)nm+1 (x+y)pm+1

25. (ln x + y 3 ) dx 3xy 2 dy = 0 R: y 3 = Cx 1 ln x

26. (xy + 2xy ln2 y + y ln y) dx + (2x2 ln y + x) dy = 0 R: 2x2 + (2x ln y + 1)2 = C, x = 0

1.2. Ecuaciones homogeneas

La ecuacion se puede escribir en la forma

1. 4x 3y + y (2y 3x) = 0 R: y 2 3xy + 2x2 = C

3. 4x2 xy + y 2 + y (x2 xy + 4y 2 ) = 0 R: (x + y)(x3 + y 3 ) = C

4. 4x2 + xy 3y 2 + y (5x2 + 2xy + y 2 ) = 0 R: (y x)8 (y 2x)9 = C(y + 2x)5

8. ax2 + 2bxy + cy 2 + y (bx2 + 2cxy + f y 2 ) = 0 R: ax3 + 3bx2 y + 3cxy 2 + f y 3 = C

14. (3y 7x + 7) dx (3x 7y 3) dy = 0 R: (x + y 1)5 (x y 1)2 = C

1.3. Ecuaciones exactas

M(x, y)dx + N(x, y)dy = 0

La ecuacion es exacta, ya que

por lo tanto, la solucion es

(y cos x + 2xey )dx + (sen x + x2 ey + 2)dy = 0.

La ecuacion es exacta, puesto que

1.4. Factores integrantes

M(x, y)dx + N(x, y)dy = 0 (1.3)

(x, y)(M(x, y)dx + N(x, y)dy = 0)

Esta ecuacion no es exacta, ya que

(xm y n+2 + xm+1 y n+1 )dx xm+2 y n dy = 0

y las derivadas parciales son

O sea que para que la ecuacion sea exacta es necesario que

Entonces, al sustituir el factor integrante en la ecuacion, esta se transforma en

Si el factor integrante es funcion de x2 + y 2 , se encuentra con

En general un factor integrante satisface la ecuacion: