Ecuaciones Diferenciales Tecnicas de Solucion y Ecuaciones

ECUACIONES DIFERENCIALES
TCNICAS DE SOLUCIN Y APLICACIONES

JOS VENTURA BECERRIL ESPINOSA
DAVID ELIZARRARAZ MARTNEZ
UNIVERSIDAD
AUTONOMA
METROPOLITANA
Casa abierta al tiempo
Azcapotzalco
DERECHOS RESERVADOS 2004, Universidad Autnoma Metropolitana (Mxico). Prohibida la reproduccin de esta obra as como la distribucin y venta fuera del mbito de la UAM. E-libro Bibliomedia Bibliomedia@mail.com
JOS VENTURA
BECERRIL ESPINOSA
Curs la Licenciatura en Fsica y Ma-
temticas en la ESFM del IPN, titu-
lndose en 1985. Obtuvo el grado de
Maestro en Ciencias en la Seccin
de Matemtica Educativa del CIN-
VESTAV en 1987 con la tesis "Algunos
resultados clsicos de la integracin
de funciones elementales" Ingres de
forma definitiva a la Universidad Au-
tnoma Metropolitana, Azcapotzalco
en 1985, en el Departamento de Cien-
cias Bsicas. Desde entonces ha pu-
blicado materiales diversos para la
docencia, as como reportes de inves-
tigacin que se han presentado en con-
gresos nacionales tanto de Matemti-
cas, como de investigacin educativa.
Actualmente es miembro del Grupo
de Investigacin de Matemtica Edu-
cativa, donde trabaja en la aplicacin
de la tecnologa para la enseanza de
las matemticas.
COLECCIN / LIBROS DE TEXTO Y MANUALES DE PRCTICA
SERIE / MATERIAL DE APOYO A LA DOCENCIA
transformando el dihgo por la razn
UNIVERSIDAD AUTNOMA METROPOLITANA
3O AOS
Jos Ventura Becerril Espinosa
David Elizarraraz Martnez
UNIVERSIDAD
AUTONOMA
METROPOLITANA
Casa abierta al tiempo
A z c a p o t z a l c o
Dr. Luis Mier y Tern Casanueva
RECTOR GENERAL
Dr. Ricardo Sols Rosales
SECRETARIO GENERAL
UNIDAD AZCAPOTZALCO
Mtro. Vctor Manuel Sosa Godnez
RECTOR
Mtro. Cristian Eduardo Leriche Guzmn
SECRETARIO
Mtra. Mara Aguirre Tamez
COORDINADORA GENERAL DE DESARROLLO ACADMICO
DCG. Ma. Teresa Olalde Ramos
COORDINADORA DE EXTENSIN UNIVERSITARIA
DCG. Silvia Guzmn Bofill
JEFA DE LA SECCIN DE PRODUCCIN Y DISTRIBUCIN EDITORIALES
ECUACIONES DIFERENCIALES. TCNICAS DE SOLUCIN Y APLICACIONES
Primera edicin, 2004
D.R. 2004 Universidad Autnoma Metropolitana
Unidad Azcapotzalco
Av. San Pablo 180, Col. Reynosa Tamaulipas
C. P. 02200, Mxico, D. F.
e.mail: secedi@correo.azc.uam.mx
Diseo y produccin editorialnopase. Eugenia Herrera/Israel Ayala
Ilustracin de portada. Israel Ayala.
Fotografa de autores. Roberto Cano
ISBN 970-31-0230-1
Impreso en Mxico/Printed in Mxico
Prlogo
Este libro est diseado para un curso trimestral de ecuaciones diferenciales or-
dinarias. Presentamos los teoremas y tcnicas de solucin que consideramos bsicos en
un estudio introductorio de sta importante disciplina de las Matemticas. Aunque no
hemos puesto nfasis en las demostraciones, proporcionamos una buena cantidad de ejer-
cicios resueltos, de modo que un estudiante de Ingeniera podra obtener, mediante su
anlisis, un nivel satisfactorio en los diferentes mtodos de solucin de ecuaciones dife-
renciales y sus aplicaciones elementales ms comunes. Para el curso sugerimos seguir
el orden previsto, no obstante, si algn lector considera que es demasiado material se
pueden omitir las secciones 2.8, 2.9 y 2.10.
Los diferentes temas se exponen en forma clara y sencilla para su inmediata com-
prensin. En las primeras secciones los desarrollos se hacen de manera exhaustiva. Ms
adelante, lo que ya es conocido no se desarrolla completamente, sino que se dejan al lector
los detalles que en ese momento ya est en capacidad de realizar. En consecuencia, el
texto se puede utilizar para un curso tradicional o bien, para un curso en el Sistema de
Aprendizaje Individualizado, en el que el alumno estudia por su propia cuenta.
En el captulo uno ilustramos la aplicacin de las ecuaciones diferenciales. Ms que la
solucin y el entendimiento de los problemas planteados, buscamos motivar el inters en
estudiar ecuaciones diferenciales, a travs de problemas actuales como son el crecimiento
de poblaciones, el impacto de la publicidad y las curvas de persecucin, entre otros.
El captulo 2 est dedicado a la solucin de ecuaciones diferenciales ordinarias de
primer orden. Se analizan, uno por uno, varios de los mtodos ms usuales e incluimos
la seccin 2.7 en la que se aborda el problema de resolver una ecuacin dada empleando
el mtodo ms conveniente. En esta seccin se proponen 80 ejercicios, debido a que
usualmente los libros de texto estudian la solucin de ecuaciones diferenciales por tema
y el alumno sabe que las ecuaciones que se le plantean son del tema estudiado; pero en
la prctica y en sus cursos posteriores las ecuaciones con que se encuentra las tiene que
resolver sin conocer de antemano de que tipo son.
8
El captulo 3 muestra la aplicacin de las ecuaciones diferenciales ordinarias de primer
orden, entre las cuales estn las ms accesibles para los estudiantes, ya que se espera que
ellos sean capaces de entenderlas y resolver problemas relacionados con stas.
En el captulo 4 nos concentramos en resolver ecuaciones diferenciales ordinarias linea-
les de orden dos. Pensamos que no nes difcil entender los resultados para orden mayor que
dos y aplicarlos a casos sencillos de ecuaciones homogneas con coeficientes constantes,
razn por la cual aparece la seccin 4.5.
Finalmente, el captulo cinco contiene diversas aplicaciones de las ecuaciones diferen-
ciales ordinarias lineales de segundo orden con coeficientes constantes. Proponemos slo
como ejercicios los correspondientes al movimiento vibratorio de una masa sujeta a un
resorte y de circuitos elctricos, para no caer en un exceso de material, aunque en el
captulo aparecen otras aplicaciones.
Queremos destacar que en las pginas 220 a la 243 se encuentran las respuestas
de todos los ejercicios propuestos, lo cual ser de gran ayuda para que el estudiante
compruebe sus conocimientos.
El libro incorpora el producto del trabajo realizado por los autores al impartir en
varias ocasiones el curso de ecuaciones diferenciales ordinarias en la UAM-A. Tambin
agradecemos a la M. en M. Marina Salazar Antnez y al M. en C. Jos Luis Huerta
Flores por sus valiosos comentarios y problemas aportados.
Al final, presentamos algunas referencias bibliogrficas que esperamos sean tiles para
los lectores interesados en profundizar en el estudio de algn tema o de conocer otras
aplicaciones.
Captulo 1
Introduccin
"El lenguaje para entender a la naturaleza es la matemtica."
Galileo Galilei.
1.1 Ecuaciones Diferenciales y Modelos Matemticos
Una gran cantidad de leyes en la Fsica, Qumica y Biologa tienen su expresin natural
en ecuaciones diferenciales ordinarias o parciales. Tambin, es enorme el mundo de las
aplicaciones de las ecuaciones diferenciales en Ingeniera, Economa, Ciencias Sociales,
Astronoma y en las mismas Matemticas. La causa es simple, si un fenmeno se puede
expresar mediante una o varias razones de cambio entre las variables implicadas entonces
correspondientemente tenemos una o varias ecuaciones ecuaciones diferenciales.
El ejemplo ms simple de una ecuacin diferencial proviene de la segunda ley de
Newton F = raa, ya que si un cuerpo cae bajo la influencia de la fuerza de gravedad
entonces
y como a
donde y(t) denota la posicin del cuerpo al tiempo , tenemos
que es una ecuacin diferencial ordinaria, cuya solucin es la funcin de posicin y (i).
Si adems suponemos que sobre el cuerpo acta una fuerza de friccin con el medio
que lo rodea, cuya magnitud es proporcional a la velocidad instantnea dy/dt, se sigue
que
de donde
10 Captulo 1. Introduccin
Otros ejemplos son las famosas ecuaciones en derivadas parciales del calor, de onda
y de Laplace, que tienen la forma
Ao
Poblacin (millones de hab.)
1900
13.61
1910
15.16
1920
14.33
1930
16.53
1940
19.65
1950
25.78
1960
34.92
Con base en los datos de la tabla y ubicndonos en el ao de 1960, se podra haber
hecho una estimacin para la poblacin de los aos 1970 y 1980 ?
Solucin. Una suposicin razonable es que la rapidez de variacin de la poblacin con
respecto al tiempo es proporcional a la poblacin, es decir si P(t) denota la poblacin al
tiempo t entonces
donde a es una constante positiva.
As, para conocer la poblacin en cualquier tiempo hay que resolver la ecuacin
anterior. La solucin es P(t) = ce
at
, con c una constante arbitraria. Para determinar c
tenemos la condicin inicial que en t = 0 (correspondiendo al ao de 1950) la poblacin
es 25.78, de donde P(t) = 25.78e
at
.
Para encontrar la constante de proporcionalidad podemos usar que P{IO) = 34.92.
En consecuencia
respectivamente, que han sido fuente inagotable de diversos trabajos de investigacin.
En nuestro caso nos restringiremos al estudio de las ecuaciones diferenciales ordinarias.
A continuacin presentaremos algunos problemas que motivan el inters en el estudio de
estas ecuaciones.
EJEMPLO 1. Se puede predecir la poblacin de un pas?
La siguiente tabla muestra el nmero de millones de habitantes que haba en toda la
Repblica Mexicana, de acuerdo al censo del ao que se indica.
Ahora para 1970 la poblacin aproximada sera P(20), que da por resultado
La poblacin para 1980 se estimar en F(30) = 64.07.
1.2. Ecuaciones Diferenciales y Modelos Matemticos 11
Es interesante comparar los valores calculados con los que se reportaron en los censos
respectivos. Los censos realizados mostraron que la poblacin en 1970 y 1980 fue de 48.22
y 67.41 millones de habitantes, respectivamente.
Con base en los dos ltimos datos qu poblacin se esperar para los aos 2000 y
2010?
EJEMPLO 2. Es posible medir el impacto de la publicidad?
Cierta compaa produce un artculo destinado a una poblacin en la que hay un
nmero M de potenciales compradores. La compaa decide establecer una campaa de
publicidad para promocionar su producto. Los propietarios de la compaa han solicitado
a su departamento de publicidad una medida del impacto de la publicidad. Se puede
ayudar a los publicistas?
Solucin. Hay varias maneras de medir el impacto de la publicidad, una es la siguiente.
Sea y(t) el nmero de personas que conocen el producto al tiempo t. Supongamos que la
velocidad con que vara el nmero de personas que conocen el producto es proporcional
tanto al nmero de personas que conocen el producto, como al de las que todava no lo
conocen. Entonces
(1.1)
donde k es una constante positiva. En la seccin 3.3.3, ejemplo 2, se muestra como
resolver (1.1). Su solucin es la funcin
(1.2)
con c una constante.
En la literatura econmica a la ecuacin (1.2) se le conoce como ecuacin de la curva
logstica, la cual nos da una medida del nmero de personas que conocen el producto al
tiempo t. La forma general de su grfica se muestra en la figura 1.1.
Figura 1.1: Curva Logstica
EJEMPLO 3. Espejos Parablicos.
Supngase que seales luminosas (o de algn otro tipo) viajan en el plano xy y
paralelamente al eje y, chocando con la curva cuya ecuacin es y f(x) y reflejndose
de tal manera que todas ellas concurren en el punto F(0,p) con p una constante positiva.
Vase la figura 1.2. Comprobar que la curva y = f(x) es una parbola y que adems en
caso de que pase por (0, 0) su ecuacin es Apy = x
2
.
Figura 1.2: Seales luminosas
Solucin. Escribiremos primeramente el problema en lenguaje matemtico. Para un
punto cualquiera P(x,y) en la curva y = /(x), la derivada es igual a la pendiente de la
recta tangente a la grfica de / en dicho punto, es decir
(1.3)
En la figura 1.3 se puede ver la configuracin de la curva y = f(x), la recta tangente l
t
en el punto P(x, y), el punto F(0,p) y los ngulos ay 0. Esta configuracin se obtuvo de
los principios bsicos de la Geometra Euclideana y del hecho fsico de que la magnitud del
ngulo de incidencia es igual a la magnitud del ngulo reflejado. Del tringulo rectngulo
APQS se obtiene la igualdad
(1.4)
y del tringulo rectngulo APFV la igualdad
Recordando la identidad
(1.5)
(1.6)
Figura 1.3: Espejo parablico
y sustituyendo en sta (1.3),(1.4) y (1.5) obtenemos finalmente la ecuacin
La ecuacin (1.7) se resuelve en el ejemplo 7 de la seccin 2.2. Su solucin es
A
2
x
2
+ 2Ap-l
(1.7)
(1.8)
donde A es una constante. Ntese que (1.8) es en efecto la ecuacin de una parbola. Si
sabemos que la curva y = f(x) pasa por el punto (0, 0) es decir y = 0 cuando x = 0, de
acuerdo con (1.8) se sigue que,
de donde
Al sustituir el valor de A en (1.8) nos lleva a
o bien
4py = x
2
,
como se afirma en el enunciado.
EJEMPLO 4. Curvas de persecucin 1.
En medio de una gran bruma, un destructor y un submarino de dos naciones distintas,
se descubren estando a 4 km de distancia. Inmediatamente el submarino se sumerge y
avanza en una direccin fija a toda velocidad. Qu trayectoria debe seguir el destructor
para asegurarse de que estar exactamente encima del submarino si tiene una velocidad
igual a tres veces la del submarino?
Solucin. Utilizaremos coordenadas polares para la trayectoria que debe seguir el sub-
marino. Para estas coordenadas se cumple que el diferencial de arco es
El destructor avanzar a toda velocidad 3 kilmetros en direccin al punto donde
localiz al submarino, a partir de ah denotemos por r = f(9) la trayectoria que debe
seguir el destructor. Ver figura 1.4
Figura 1.4: Trayectorias del submarino y del destructor
La distancia ds que recorre el submarino hasta el punto donde las trayectorias se
cortan es
y la distancia d^ que recorre el destructor hasta la interseccin es
por ser su velocidad el triple de la del destructor.
Submarino
Destructor
Luego, tenemos que
o sea que
Derivando con respecto a 6 ambos lados de la igualdad anterior resulta
de donde
As, para conocer qu trayectoria debe seguir el destructor debemos resolver la ecuacin
diferencial
es decir, hay que encontrar una funcin r(6) cuya derivada sea
Una funcin que satisface la ecuacin (1.9) es
como puede comprobarse inmediatamente, por sustitucin directa de dicha funcin y su
derivada con respecto a 6 en (1.9).
De todo lo anterior podemos concluir que si el destructor sigue la trayectoria r e
e
^
s
despus de haber avanzado 3 kilmetros en lnea recta en direccin a donde localiz al
submarino, puede estar seguro de que pasar por encima del submarino sin importar que
direccin elija este ltimo.
Finalmente, es claro que no es la nica trayectoria que puede seguir el destructor y
a los que viajan en el submarino habra que recomendarles que no avancen en una sola
direccin.
EJEMPLO 5. Curvas de persecucin 2.
Dos equipos de ftbol americano se encuentran en juego. En un momento determinado
un jugador q recibe el baln y corre en lnea recta a la portera contraria con una velocidad
v
q
. En ese mismo instante otro jugador p (perseguidor) corre con velocidad v
p
en direccin
de q para tratar de interceptarlo, (ver figura 1.5). Nos interesa saber bajo que condiciones
p alcanza a q, es decir, si por ejemplo v
p
> v
q
, p alcanzar a q?
Figura 1.5: El jugador p persigue al q
Solucin. Como el corredor p va tras el corredor g, la direccin del vector velocidad de
p siempre apunta hacia q (figura 1.6).
Trayectoria del
corredor q
Figura 1.6: Trayectorias de los jugadores
En general, al tiempo , p se encuentra en (x(), j/()) y como q corre en lnea recta,
l se encontrar a un distancia s(t) = v
q
t del eje x, es decir, q se encuentra en el punto
(a,v
q
t), (ver figura 1.7).
Figura 1.7: Posicin de los jugadores al tiempo t
Usemos la figura 1.7 para obtener la ecuacin diferencial que decribe esta situacin.
Se tiene que
donde a, es una constante.
Esta ecuacin tiene tres variables x, t / y , pero puede reducirse a una que contenga
solamente dos de ellas. Si observamos un tramo infinitesimal de la trayectoria de p
tenemos
de donde
Al derivar ambos lados de (1.10) con respecto a x, aparecer En efecto resulta
Captulo 1. Introduccin
De acuerdo con (1.10), la ltima igualdad se reduce a
y usando (1.11)
As que
La ecuacin (1.12) se resolver en la seccin 2.1 de ecuaciones de variables separables
en el ejemplo 11. Su solucin es
con k
donde A y B son constantes arbitrarias. Ahora bien, como al tiempo t = 0 el corredor p
est en el origen, esto es y(0) = 0, la ecuacin (1.13) nos da
Hay otra condicin inicial al tiempo t 0, a saber y'(0) = 0 pues el vector velocidad
de p es horizontal. Derivando (1.13)
y sustituyendo x 0 y 0, se sigue que
Sustituimos este valor de A en (1.14), para obtener el valor de B
Finalmente obtenemos
1 y por tanto k + 1 > 0. Evaluando (1.15) en x = a resulta
es decir, el corredor p alcanza al corredor q en el punto
Por supuesto, dependiendo de los valores de a y fc se puede saber si p alcanza a q
dentro de la cancha.
EJEMPLO 6. Por qu un reloj de pndulo es impreciso?
Consideremos un modelo idealizado de reloj de pndulo, formado por una cuerda de
longitud / y un peso de masa m en su extremo, como se muestra en la figura 1.8.(a).
Inicialmente el peso se desva un ngulo a y luego se deja libre (ver figura 1.8.(b)).
Sea 9(t) el ngulo en radianes al tiempo t entre la cuerda y la posicin vertical, de la
figura 1.8.(a). Adoptamos la convencin de que 9 > 0 cuando la masa est a la derecha de
la posicin de equilibrio y 9 < 0 cuando est a la izquierda, lo cual esencialmente significa
que escogemos direcciones a lo largo del arco apuntando a la derecha como positivas y a
la izquierda como negativas.
La relacin entre la longitud del arco s y el ngulo 9 es s = 19, de donde
Por otra parte la fuerza causante del movimiento es el peso w mg. Esta fuerza se
descompone en sus componentes F\ y F2 en la direccin de la tangente a la trayectoria y
Figura 1.9: Movimiento del pndulo
La segunda ley de Newton dice que
ma = ^2( fuer zas),
y aplicada al pndulo, usando (1.16) y (1.17), conduce a
o equivalentemente
Para resolver (1.18) sea w De la regla de la cadena
con lo cual (1.18) toma la forma
cuya solucin es
con A una constante, as que
Usando que el valor de A resulta ser A
modo que
de donde
El signo menos en (1.19) toma en cuenta el hecho de que 9 decrece con el crecimiento de
t. Adems (1.19) implica que
y una segunda integracin nos lleva a
Captulo 1. Introduccin
donde B es una constante. Pero, inicialmente t 0, 9 = a, as que
En consecuencia
Calculemos el periodo a partir de (1.20). Ya que el periodo T es el tiempo que tarda el
pndulo en dar una oscilacin completa, al tiempo t T/4 la cuerda estar por primera
vez en posicin vertical, es decir haciendo t = T/4 y 9 = 0 en (1.20) resulta
o bien
Como se observa en (1.21) el periodo de las oscilaciones del pndulo depende del
ngulo inicial a. Precisamente este hecho es la causa principal por la que el reloj de
pndulo es impreciso, ya que en la prctica el peso cada vez se desva hacia su posicin
extrema en un ngulo distinto de a.
1.2 Conceptos Bsicos de Ecuaciones Diferenciales
Una igualdad que contenga una o ms derivadas de una funcin desconocida se llama
ecuacin diferencial.
Las siguientes igualdades son ejemplos de ecuaciones diferenciales.
con a, b constantes
1.2. Conceptos Bsicos de Ecuaciones Diferenciales 23
Las ecuaciones anteriores modelan cierto fenmeno. La ecuacin (1.22) es llamada
ecuacin logstica y describe el crecimiento de una poblacin. Las ecuaciones (1.23) y
(1.24) corresponden al movimiento armnico simple y forzado amortiguado, respectiva-
mente, que estudiaremos en el captulo 5. Las ecuaciones (1.25), (1.26) son las ecuaciones
del calor y de onda en una dimensin y finalmente, la ecuacin (1.27) es la ecuacin de
Laplace en dos dimensiones.
Se dice que una ecuacin diferencial es ordinaria, si la funcin incgnita depende de
una sola variable. Si la funcin incgnita depende de ms de una variable, entonces la
ecuacin se llama una ecuacin diferencial parcial.
Las ecuaciones (1.22), (1.23) y (1.24) son ordinarias, mientras que las ecuaciones
(1.25), (1.26) y (1.27) son parciales.
En todo lo que sigue consideraremos nicamente ecuaciones diferenciales ordinarias,
por lo que en ocasiones omitiremos mencionar explcitamente que se trata de ecuaciones
de esta clase.
El orden de una ecuacin diferencial es el de la derivada de mayor orden que aparece
en la ecuacin.
As, la ecuacin (1.22) es de primer orden, las ecuaciones (1.23) y (1.24) son de
segundo orden. Veamos otro ejemplo.
EJEMPLO 1. Determine el orden de la ecuacin diferencial dada.
Solucin. La ecuacin (a) es de primer orden, las ecuaciones (b) y (d) de segundo orden
y la ecuacin (c) es de cuarto orden.
Simblicamente, una ecuacin diferencial ordinaria de orden n, puede expresarse en
la forma
donde F es una funcin de n + 2 variables. Para nuestros propsitos supondremos que
(1.28) tambin admite la representacin
para alguna funcin / de n + 1 variables.
Una solucin de una ecuacin diferencial en un intervalo / es cualquier funcin
definida en / que satisface a la ecuacin, es decir, que al sustituirla la reduce a una
identidad.
EJEMPLO 2. Comprobar que la funcin f(x) = e~
3x
+ 5 es solucin de la ecuacin
diferencial dada en el intervalo (oc, oo).
Solucin. Es claro f{x) y f'(x) se definen para todo x en R. Sustituyendo sus expresiones
en la ecuacin diferencial resulta
La ltima identidad establece que efectivamente la funcin / es una solucin de la
ecuacin diferencial en R.
EJEMPLO 3. Verificar que la funcin g(x)
dada en el intervalo (0, oo).
es solucin de la ecuacin diferencial
Solucin. Sea x > 0. Derivando g(x) obtenemos
Sustituyendo g(x) y g"{x) en la ecuacin diferencial se sigue que
Lo cual comprueba que g es una solucin en (0, oo).
EJEMPLO 4. Sea c una constante arbitraria. Probar que la funcin definida por la
ecuacin xy
2
y
3
= c, es solucin de la ecuacin diferencial
Solucin. Derivando implcitamente la igualdad xy
2
y
3
= c, con respecto a x, se tiene
que
1.2. Conceptos Bsicos de Ecuaciones Diferenciales
En particular
de donde se observa que la ecuacin xy
2
y
3
= c define una solucin de la ecuacin
diferencial para todo valor de c. En este caso se dice que la solucin est en forma
implcita.
EJEMPLO 5. Determine si la funcin h(x) = e~
3x
es una solucin en R de la ecuacin
diferencial: y" + 6y' + 9y = 0.
Solucin. Derivando dos veces la funcin /i, tenemos que para x en R
Sustituyendo en la euacin diferencial encontramos que
As, la funcin dada es una solucin de la ecuacin diferencial.
Un problema de valores iniciales es aqul en el que se busca determinar una
solucin a una ecuacin diferencial, sujeta a condiciones sobre la funcin desconocida y
sus derivadas, dadas para un valor de la variable independiente. Tales condiciones se
llaman condiciones iniciales. En smbolos, un problema de valores iniciales de orden n,
puede representarse por la ecuacin (1.29) sujeta a las n condiciones
En los captulos 2 y 3 resolveremos algunos problemas de valor inicial, de primer
orden, que en forma general pueden expresarse como
Al encontrar un problema de este tipo, es natural preguntarse si tiene solucin y en
tal caso, si dicha solucin es nica. La respuesta a estas cuestiones viene dada por el
siguiente teorema.
Teorema 1.2.1 (Teorema de existencia y unicidad.) Sea R [a, b] x [c, d] un
rectngulo en el plano que contiene al punto (x
o
,y
o
) en su interior. Si las funciones f y
df/dy son continuas en R, entonces existe un intervalo I con centro en XQ y una funcin
nica y(x) = (p(x) definida en I que satisface el problema de valor inicial definido por
(1.31) y (1.32).
Captulo 1. Introduccin 26
Trataremos de explicar el significado del teorema anterior mediante los siguientes dos
ejemplos.
EJEMPLO 6. Demuestre que el problema de valor inicial
tiene solucin nica.
Solucin. Aplicamos el teorema 1.2.1. En primer lugar comprobaremos que se cumple
la hiptesis. En este caso Ambas funciones / y
df/dy son continuas en todo rectngulo R del plano xy. La condicin inicial y(2) = 1
significa que xo = 2 y yo 1, adems es claro que el punto (2, 1) est contenido en el
interior de algn rectngulo R. As que todas las hiptesis del teorema 1.2.1 se satisfacen
y por lo tanto la conclusin se cumple; es decir, existe una solucin nica ip de la ecuacin
7
diferencial definida en algn intervalo / con centro en XQ = 2, que satisface
la condicin inicial, esto es, que es tal que </?(2) = 1.
EJEMPLO 7. Considere los dos problemas de valor inicial
Estas funciones son continuas en el conjunto A de puntos del plano xy definido por
En el problema (a), x
0
= 2, yo = 0. El cuadrado R\ de lado 1 con centro en el punto
(2,0) est contenido en el conjunto A, de modo que las funciones / y df/dy satisfacen
las hiptesis del teorema 1.2.1 en R\. Dado que adems el punto (2,0) est en el interior
de R\
)
la conclusin del teorema 1.2.1 se aplica al problema (a) y sabemos que tiene una
solucin nica definida en algn intervalo alrededor de x
0
= 2.
Veamos el problema (b). Ahora xo =. TT/2, yo 1. En este punto df/dy no es continua.
Luego, el punto (TT/2, 1) no puede estar incluido en un rectngulo R donde las hiptesis
del teorema 1.2.1 se satisfagan. Entonces no podemos concluir, del teorema 1.2.1, que el
1.2. Conceptos Bsicos de Ecuaciones Diferenciales 27
problema (b) tenga una solucin nica. Con esto no afirmamos que no tenga solucin o
tenga varias; el teorema 1.2.1 simplemente no da informacin en un sentido u otro. De
hecho puede verificarse con un clculo sencillo que las funciones
con x en [7r/2,7r/2], son soluciones del problema (b), as que ste tiene al menos dos
soluciones definidas en el intervalo indicado.
Se llama solucin general de una ecuacin diferencial de primer orden y
1
= / (#, y) a
la funcin y = <p(x, c) que depende de una constante arbitraria c y satisface las siguientes
condiciones:
1. Es solucin de la ecuacin diferencial para cualquier valor de c.
2. Dada una condicin inicial arbitraria y(xo) yo, siempre es posible determinar un
valor de c = c
0
tal que la funcin y = (p(x, CQ) satisface la ecuacin diferencial y la
condicin inicial. La funcin y ip(x, CQ) se llama una solucin particular.
Geomtricamente la solucin general y = </?(#, c) representa una familia de curvas en
el plano xy. Estas curvas se llaman curvas integrales. Si las condiciones del teorema
de existencia y unicidad se satisfacen, estas curvas integrales no se intersectan.
EJEMPLO 8. Muestre que la funcin
donde c es una constante, es la solucin general de la ecuacin diferencial
Solucin. En primer lugar, podemos ver que (1.33) es solucin de (1.34) para todo c.
En efecto, derivando (1.33) con respecto a x se obtiene
y sustituyendo en (1.34) resulta
lo cual es una identidad.
Adems, haciendo x en (1.33), se sigue que
de donde
por lo cual, la funcin
es solucin de (1.33) y satisface la condicin inicial y(x
0
) = y
0
?
c o n x
o > 0 y yo arbitrario.
La figura 1.10 muestra algunas de las curvas integrales con c > 0.
y
Figura 1.10: Curvas integrales
1.2.1 Mtodo de Isoclinas
La ecuacin diferencial
(1.35)
donde la funcin f(x,y) est definida en algn conjunto D del plano xy, determina en
cada punto (x, y) de D, el valor de y', o sea, la pendiente de la recta tangente a la curva
integral en este punto. Luego, podemos interpretar la ecuacin diferencial (1.35) como
un conjunto de pendientes llamado campo de direcciones. La terna de nmeros (x,y,y
f
)
determina la direccin de una recta que pasa por el punto (x,y). El conjunto de los
segmentos de estas rectas es la representacin geomtrica del campo de direcciones.
El problema de resolver la ecuacin diferencial (1.35) puede entonces interpretarse
como sigue: encontrar una curva cuya tangente en cada punto tenga la misma direccin
que el campo en este punto. Para construir las curvas integrales introducimos las isocli-
nas. Se llama isoclina al lugar geomtrico de los puntos en los que las rectas tangentes a
las curvas integrales consideradas tienen una misma direccin. La familia de las isoclinas
de la ecuacin diferencial (1.35) viene dada por la condicin
f{x,y) = c, (1.36)
1.2.1. Mtodo de Isoclinas 29
donde c es una constante. Usando valores de c cercanos podemos dibujar una red bastante
compacta de isoclinas, a partir de las cuales es posible trazar aproximadamente las curvas
integrales de la ecuacin (1.35).
EJ EMPLO 1. Mediante las isoclinas, trace aproximadamente las curvas integrales de
la ecuacin diferencial
(1.37)
Solucin. En este caso f(x,y) x y las isoclinas estn dadas por la ecuacin x = c,
donde c es una constante. En consecuencia las isoclinas son rectas verticales. Para c 0
se obtiene la isoclina x = 0, el eje y. Este eje divide al plano en dos partes iguales, en
cada una de las cuales la derivada y' tiene un mismo signo. Las curvas integrales son
decrecientes si x < 0 y crecientes si x > 0, de modo que sobre esta recta se encuentran
sus puntos mnimos.
Las tangentes trazadas a las curvas integrales en los puntos de interseccin con las
isoclinas x 1 y x = 1, forman con el eje OX, ngulos de 45 y 135, respectivamente.
El campo de direcciones se muestra en la figura 1.11.
Figura 1.11: Campo de direcciones
Adems, derivando en (1.37) con respecto a x, se sigue que y" = 1. Por consiguiente
y" > 0 para todo x y las curvas integrales son cncavas hacia arriba. Tomando en cuenta
todo lo anterior un esbozo de la familia de curvas integrales de la ecuacin (1.37) se
muestra en la figura 1.12.
/ /
Figura 1.12: Curvas Integrales del Ejemplo 1
EJEMPLO 2. Represente el campo de direcciones e indique algunas curvas integrales
de la ecuacin diferencial
(1.38)
x x
Solucin. Ahora f(x,y) = y las isoclinas son rectas de la forma = c o bien
y y
Si c 1 tenemos la isoclina y = x a lo largo de la cual la inclinacin de las tangentes
a las curvas integrales es de 45. Con c = 1 resulta la isoclina y = x sobre la cual las
tangentes forman un ngulo de 135 con respecto al eje OX.
Adems, a partir de la ecuacin diferencial misma podemos concluir lo siguiente. Las
tangentes trazadas a las curvas integrales en los puntos de interseccin con el eje x (y = 0)
y con el eje y (x 0) son verticales y horizontales, respectivamente, con excepcin del
punto (0,0).
El campo de direcciones y algunas curvas integrales se muestran en la figura 1.13.
1.2.1. Mtodo de Isoclinas
31
Figura 1.13: Campo de direcciones y Curvas Integrales
EJEMPLO 3. Utilizando el campo de direcciones y las isoclinas determine aproximada-
mente la curva integral de la solucin del problema de valor inicial
(1.39)
Solucin. Tenemos que / (x, y) = x
2
+y
2
y las isoclinas / (x, y) c son las circunferencias
con centro en el origen definidas por
Con los valores de c = 1/4, c 1, c = 9/4 y c = 4 resultan las circunferencias de
radios 1/2,1,3/2 y 2, respectivamente, que se muestran en la figura 1.14. A lo largo de
cada una de ellas las pendientes de las tangentes a las curvas integrales es igual al valor
particular de c.
Examinando la figura 1.14 parece razonable que la grfica de la solucin del problema
(1.39) tenga la forma que se indica en la figura 1.15. Cabe destacar el valor de este anlisis,
en vista de que el problema (1.39) no puede resolverse empleando tcnicas elementales.
Figura 1.14: Campo de direcciones e Isoclinas
Figura 1.15: Solucin del problema de valor inicial
Ejercicios del Captulo 1. 33
EJERCI CI OS 1.2
En cada ejercicio indique el orden de la ecuacin diferencial y determine si la funcin
dada es una solucin de la ecuacin.
1. x- y = 2x
2
y, y = Ax^.
dx
y
y> y
2. y' + ycosx = eos x sen x, y ce~
senx
+ senx + 1.
3. y' + yt anx = 0, y = Acosx.
y + x o
9
4. y = , x^ - y
2
+ 2xy = c.
y x
5. y
1
y + 1, y = xl nx.
x
6. (x
2
+ y
2
)dx + (x
2
- xy)dy = 0, A(x + y)
2
= xe^.
Captulo 2
Ecuaciones Diferenciales de Primer
Orden
2.1 Ecuaciones Diferenciales de Variables Separables
Definicin 2.1.1 Se dice que una ecuacin diferencial ordinaria es de variables sepa-
rables si se puede escribir en la forma
La ecuacin (2.1) se expresa en forma diferencial como
y se resuelve integrando ambos miembros de (2.2).
EJEMPLO 1. Resolver
Solucin. Separando las variables resulta
e integrando
Resolviendo las integrales
35
(2.1)
(2.2)
(2.3)
36 Captulo 2. Ecuaciones Diferenciales de Primer Orden
Ya que la diferencia de constantes es una constante, podemos escribir c = C2 C\, obte-
niendo
As, al momento de integrar slo consideraremos una constante de integracin.
EJEMPLO 2. Resolver
Solucin. Separando variables la ecuacin se escribe como
integrando
y calculando las integrales, se sigue que
Como el producto de constantes es una constante tenemos
EJEMPLO 3. Resolver
Solucin. Tenemos que
2.1. Ecuaciones Diferenciales de Variables Separables
Ya que e
u
y lnt son funciones inversas,
37
Como c\ es una constante, e
ci
tambin es una constante, la cual podemos escribir sim-
plemente como c; de modo que
EJEMPLO 4. Resolver
Solucin. Procediendo como en los ejemplos anteriores, resulta
(2.6)
En este caso la solucin queda en forma implcita.
EJEMPLO 5. Resolver
Solucin. Tenemos
38
Captulo 2. Ecuaciones Diferenciales de Primer Orden
EJEMPLO 6. Resolver
Solucin. Separando variables se sigue que
Por lo tanto
es la solucin dada en forma implcita.
EJEMPLO 7. Resolver
Solucin. Para separar variables es de gran ayuda factorizar donde sea posible, en este
caso tenemos
Finalmente, al integrar encontramos que
2.1. Ecuaciones Diferenciales de Variables Separables 39
EJEMPLO 8. Resolver el problema de valor inicial
Solucin. Separando variables e integrando obtenemos
Haciendo x = 0 y ? / = l e n l a ltima igualdad, concluimos que
Por lo tanto la solucin de la ecuacin diferencial que cumple la condicin dada es
Solucin. Primero separamos variables
En segundo lugar integramos, usando fracciones parciales para la integral respecto de y.
Obtenemos
obtenemos as la solucin explcita
Ahora, despejamos y en la ltima igualdad
Si hacemos x = 2 y y = 4 en (2.12), tenemos
Finalmente, sustituyendo el valor de c en (2.12), llegamos a la solucin particular
La siguiente ecuacin diferencial es de segundo orden, sin embargo, mediante un
cambio de variable se reduce a una de primer orden. Corresponde a la ecuacin diferencial
(1.12) del ejemplo 5 del captulo 1.
Solucin. Separando variables e integrando se sigue que
Ahora despejamos y para expresar la solucin en forma explcita
es decir
con c = c\ 1.
Se quiere una solucin que cumpla con J/(TT) = 0, entonces
de donde
Sustituyendo en (2.14) obtenemos la solucin del problema de valor inicial
EJEMPLO 11. Resolver
donde a, k son constantes tales que
Solucin. Si hacemos
la ecuacin (2.15) se reduce a
Separando variables e integrando resulta
donde A es una constante. As que
Para obtener y necesitamos hacer una segunda integracin, pero antes escribimos (2.16)
en la forma
Ahora podemos integrar fcilmente, encontrando la solucin explcita siguiente
EJERCI CI OS 2.1
Mediante separacin de variables resuelva las siguientes ecuaciones diferenciales.
2.2 Ecuaciones Diferenciales Homogneas
Definicin 2.2.1 La funcin f(x,y) se llama homognea de grado n con respecto a las
variables x. y si para todo t se verifica que
lo cual muestra que la funcin f(x,y) 2x
3
5xy
2
+ 4j/
3
es una funcin homognea de
grado tres.
b) Se tiene que
Solucin.
a) En este caso
EJEMPLO 1. Diga si la funcin dada es homognea y determine su grado.
2.2. Ecuaciones Diferenciales Homogneas 45
t
2
y
2
d
t
2
(x
2
-y
2
)
t*xy
9 9
t
o% -y
xy
por lo cual
Lo cual muestra que / (#, y) si es una funcin homognea y de grado uno.
f) Ahora tenemos
no es homognea.
concluimos que f(x,y)
e) Se tiene que
Luego, f(x.y) si es una funcin homognea y de grado uno.
c) Tenemos ahora que
As, f(x,y) es homognea de grado cero.
d) Como
as que f(x,y) es homognea de grado cero.
Definicin 2.2.2 Se dice que la ecuacin diferencial
es homognea si las funciones M y N son homogneas y del mismo grado.
Note que la ecuacin diferencial
ser homognea si / es una funcin homognea de grado cero.
Mtodo de Solucin. Una ecuacin diferencial homognea Af (#, y)dx + N(x, y)dy = 0,
se resuelve reducindola a una ecuacin de variables separables, usando cualquiera de las
sustituciones v = y/x o bien v x/y, donde v es una nueva variable.
Nota: Aunque en teora cualquiera de las dos sustituciones anteriores reduce una ecuacin
homognea a una separable, en la prctica sugerimos utilizar
y = xv si N es de estructura "ms simple" que M. y
x = yv si M es de estructura "ms simple" que N.
El tomar en cuenta esta observacin, conduce a integrales ms fciles de calcular al
resolver la ecuacin diferencial separable que se obtiene.
EJEMPLO 1. Resolver
Solucin. Como
son funciones homogneas ambas de grado tres, la ecuacin dada es homognea. Adems
TV es de estructura algebraica ms simple que M, por lo cual, la sustitucin ms conve-
niente es y xv para reducir (2.17) a una ecuacin de variables separables.
Hacemos
Sustituyendo en (2.17) obtenemos
Integrando
de donde
Reemplazando v y simplificando encontramos que
EJEMPLO 2. Resolver
Solucin. Ya que
son funciones homogneas ambas de grado uno, nuestra ecuacin a resolver es homognea.
Como en el ejemplo anterior TV es de estructura algebraica mas simple que M. Por lo
cual hacemos
Sustituyendo en (2.18) obtenemos
Integrando
reemplazando v y simplificando, encontramos que
donde c
EJEMPLO 3. Resolver
Solucin. En este caso la estructura algebraica de M(x, y) es ms simple que la de
N(x,y), por lo cual proponemos
Entonces
o bien
Integrando
y usando v
obtenemos como solucin implcita
EJEMPLO 4. Resolver
Solucin. Escribimos primero la ecuacin diferencial como
(2.20)
2.2. Ecuaciones Diferenciales Homogneas
Ponemos x vy en la ecuacin diferencial, de lo cual se obtiene
o equivalentemente
Integramos
En trminos de la variable original, la solucin es
EJEMPLO 5. Resolver
(2.21)
Solucin. Escribimos la ecuacin diferencial como
Hacemos
Entonces
Sustituyendo
Separando variables
e integrando
La solucin, en forma implcita es
Proposicin. Las ecuaciones diferenciales de la forma
49
donde a^b^Ci G R (i = 1,2), se denominan cuasi-homogneas y se reducen a ho-
mogneas, haciendo el cambio de variables
siendo h y k las soluciones del sistema
Si el sistema no tiene solucin, es decir
la ecuacin diferencial puede escribirse como
la cual se reduce a separable con la sustitucin
EJEMPLO 6. Resuelve la ecuacin diferencial cuasi-homognea
(2.22)
Solucin. Hacemos
Entonces
y sustituyendo en la ecuacin, se tiene
Para que esta ecuacin diferencial sea homognea es necesario que
La solucin de este sistema de ecuaciones lineales es
Con estos valores de h y k la ecuacin diferencial se reduce a
que es homognea.
Hacemos ahora
Sustituyendo, se obtiene
o bien, separando variables
e integrando resulta
Regresando a las variables originales
obtenemos como solucin general de la ecuacin diferencial
El siguiente ejemplo es una ecuacin diferencial de segundo orden, pero una vez ms,
utilizando un cambio de variable se transforma en una ecuacin diferencial de primer
orden que adems es cuasi-homognea. Corresponde a la ecuacin (1.7) que apareci en
el ejemplo 3 del captulo 1.
EJEMPLO 7. Resolver
donde p es una constante.
Solucin. Poniendo w = (2.24) se escribe como
Esto nos hace recordar el trinomio cuadrtico, por lo que despejando w, obtenemos
o bien,
la cual como sabemos es una ecuacin cuasi-homognea. Si ponemos v
y p y naturalmente,
Ahora bien, si usamos (2.25) y seleccionamos el signo ms, la ecuacin (2.26) se reduce
a
(2.26)
Separamos variables e integramos.
La grfica de la solucin obtenida es una parbola. Para verlo, despejemos la raz
elevamos al cuadrado y simplificamos para obtener
de modo que la solucin de la ecuacin diferencial df 0 est dada implcitamente por
f{x,v) =c.
2.3. Ecuaciones Diferenciales Exactas 53
EJERCI CI OS 2.2
Resolver:
2.3 Ecuaciones Diferenciales Exactas
Definicin 2.3.1 Si z f(x,y) es una funcin con derivadas parciales de primer or-
den continuas en una regin rectangular R del plano xy, entonces su diferencial total,
denotada por dz o df, se define como
Ahora bien si f{x,y) = c, donde c es una constante, entonces
EJ EMPLO. Si f(x,y) entonces es decir
o bien
Note que la ecuacin diferencial (2.27) no es separable ni tampoco homognea, decimos
que es exacta y su solucin es x
4
+ 3x
2
y
2
+ y
3
c.
(2.27)
De manera general hacemos la siguiente definicin.
Definicin 2.3.2 Se dice que una ecuacin diferencial
(2.28)
es exacta si puede escribirse en la forma df = 0, es decir
Equivalentemente, la ecuacin diferencial (2.28) es exacta si existe una funcin f tal que
Note que, la solucin de una ecuacin diferencial exacta est dada implcitamente por
la ecuacin f(x,y) = c, donde c es una constante arbitraria.
A la funcin / que cumple las ecuaciones (2.29) se le denomina funcin potencial,
mientras que a la funcin vectorial
se le llama campo vectorial conservativo. En este contexto, resolver la ecuacin diferencial
exacta (2.28) es equivalente a encontrar la funcin potencial del campo (2.30). Puede
consultarse [8], para estudiar esta clase de campos vectoriales.
El siguiente teorema proporciona un criterio simple para determinar si una ecuacin
diferencial es exacta. Su aplicacin queda clara en los ejemplos posteriores.
Teorema 2.3.1 Sean las funciones M, N, M
y
y N
x
continuas en la regin rectangular
R. Entonces la ecuacin
es exacta en R si y slo si
para todo punto (x, y) en R.
(2.31)
EJEMPLO 1. Resolver
(2.32)
Solucin. Verifiquemos, primero, que (2.32) es una ecuacin diferencial exacta. Aqu
tenemos que
y como
afirmamos que (2.32) es una ecuacin diferencial exacta. Luego, existe una funcin / (#, y)
tal que la ecuacin (2.32) se puede escribir en la forma df(x,y) = 0. Es decir, que
Para determinar / integramos (2.33) con respecto de x, resulta
o bien
(2.35)
donde <j)(y) es una funcin de y, ya que integramos con respecto de x.
Derivando (2.35) parcialmente con respecto a y se obtiene
pero como deseamos que tambin se satisfaga (2.34), igualando (2.36) con (2.34), se sigue
que
Luego
Integrando ahora ambos lados respecto a y obtenemos
con C\ una constante arbitraria.
Sustituimos <f>(y) en (2.35) y se tiene que
Finalmente, igualamos f(x,y) con una constante k para obtener la siguiente solucin de
(2.32), definida implcitamente
Renombrando c = k c\, resulta la solucin
EJEMPLO 2. Resolver
(2.37)
Solucin. Escribamos la ecuacin en su forma diferencial,
Esta ecuacin diferencial es exacta ya que
luego, existe una funcin / tal que
Integrando respecto a x
es decir
donde (f>(y) es una funcin que depende nicamente de y.
Derivamos parcialmente a (2.38) respecto a y
pero sabemos que
por lo que
De esta ecuacin resulta
con c\ una constante.
Luego, sustituimos cf>{y) en (2.38) y se tiene que
Finalmente, tomando en cuenta que f(x,y) = k da la solucin implcita, obtenemos
EJEMPLO 3. Resolver
Solucin. Esta ecuacin en su forma diferencial nos queda de la siguiente forma
(2.39)
en donde
y como
Se tiene que nuestra ecuacin a resolver es una ecuacin diferencial exacta, por lo que
existe una funcin / tal que
Luego
Pero entonces
Sustituyendo (/>(y) se tiene que
La solucin est dada por
de donde
Nota: Para resolver las ecuaciones exactas anteriores primero integramos con respecto
a x en la igualdad pero se puede proceder en forma anloga si en lugar
de esto integramos con respecto a y en
en los tres ejemplos siguientes.
Ilustraremos dicho procedimiento
EJEMPLO 4. Resuelva
Solucin. Escribamos (2.40) en su forma diferencial
En este caso
y dado que
tenemos que (2.40) es una ecuacin diferencial exacta, por lo que existe / tal que
Integrando esta ecuacin respecto a y obtenemos
Derivando parcialmente con respecto a x resulta
Recordando que e igualando con la expresin anterior, tenemos que
Sustituyendo (f>(x) tenemos
La solucin est dada por
o bien
EJEMPLO 5. Resolver
y como
se trata de una ecuacin diferencial exacta. Buscamos una funcin / tal que
Por otro lado, se requiere que
Igualando las dos expresiones anteriores se tiene que
Sustituyendo </>(#) obtenemos
Luego, la solucin viene dada por
Aplicando la condicin y(0) = 6 en la ltima ecuacin tenemos que
As, concluimos que nuestra solucin particular est definida implcitamente mediante la
ecuacin
EJEMPLO 6. Resolver
Solucin. Como
tenemos que (2.42) es una ecuacin diferencial exacta y por lo tanto existe una funcin
/ tal que
de donde
Entonces
Pero de donde
Luego
o bien, la solucin y est definida implcitamente en la ecuacin
Como y est sujeta a la condicin y(0) e, se tiene que
As, y est definida implcitamente en la ecuacin
EJEMPLO 7. Resolver
Solucin. Se tiene que
por lo cual (2.43) no es una ecuacin diferencial exacta.
62 Captulo 2, Ecuaciones Diferenciales de Primer Orden
EJERCICIOS 2.3
De las siguientes ecuaciones resuelva aqullas que sean exactas.
2.4 Factores Integrantes
Definicin 2.4.1 Si la ecuacin diferencial
no es exacta, pero existe una funcin n{x,y), tal que al multiplicar (2.44) P
or
M
x
>/)> l>
a
ecuacin resultante
es exacta, entonces se dice que /i(x, y) es un factor integrante de la ecuacin diferencial
(2.44)-
(2.44)
2.4. Factores Integrantes 63
Debemos observar que la solucin de (2.45) es la solucin de (2.44) y que en general no
es fcil encontrar un factor integrante para una ecuacin no exacta de la forma (2.44).
Sin embargo, si M(x,y) y N(x,y) cumplen ciertas condiciones entonces los factores in-
tegrantes son conocidos. Veamos algunos casos.
CASO I. Factor integrante dependiente de x. Suponga que
es una funcin que depende nicamente de x, la cual denotaremos por g{x). Entonces,
un factor integrante para la ecuacin dada es
CASO II. Factor integrante dependiente de y. Si se tiene que
es una funcin de y nicamente, denotada por h(y), entonces
es un factor integrante para la ecuacin diferencial (2.44).
CASO III.Factores de integracin de la forma Si existen m y n tales que
entonces
es un factor integrante para (2.44).
CASO IV. Si existen funciones P(x) y Q(y) que satisfacen
entonces un factor integrante para (2.44) es
Obsrvese que el Caso IV incluye a los Casos I, II y III si tomamos
respectivamente.
Tambin queremos advertir que al aplicar las frmulas anteriores estamos interesados
en obtener solamente un factor integrante, por lo cual despus de calcular las integrales
indefinidas implicadas en dichas expresiones basta considerar un valor fijo de la constante
de integracin; cero por simplicidad. Continuaremos con esta prctica ms adelante, al
emplear las frmulas (2.60), (4.15) y (4.85).
EJEMPLO 1. Resolver
Solucin. En este caso
Como tenemos que (2.46) no es una ecuacin diferencial exacta. Buscaremos
un tactor integrante para (2.46) investigando si M y N cumplen con las condiciones
mencionadas en el Caso I.
La expresin en (2.47) no es una funcin exclusivamente de x. Por lo que investigaremos
si M y TV son funciones que cumplen con la condicin mencionada en el Caso II.
La expresin en (2.48) si es una funcin exclusivamente de y, luego un factor integrante
es de la forma
As
Ya que se conoce un factor integrante de (2.46), multiplicamos dicha ecuacin por y
2
y
procedemos a resolver la ecuacin diferencial resultante, cuya solucin es igual a la de
(2.46). Tenemos
Ahora se tiene en (2.49) que
Luego, (2.49) ya es una ecuacin diferencial exacta, cuya solucin est definida implcitamente
en la ecuacin
EJEMPLO 2. Resolver
(2.50)
De modo que (2.50) no es una ecuacin diferencial exacta. Veamos si M y N cumplen
con la condicin del Caso I.
Como (2.51) es exclusivamente funcin de x, un factor integrante es
Multiplicamos (2.50) por /(x) y obtenemos la ecuacin
la cual es una ecuacin diferencial exacta que tiene la misma solucin que (2.50), definida
implcitamente en la ecuacin
EJEMPLO 3. Resolver
Solucin. Ahora
(2.52)
Entonces (2.52) no es una ecuacin exacta. Hallaremos un factor integrante. Primero,
veamos si M y N cumplen con las condiciones mencionadas en el Caso I.
Vemos que (2.53) es una funcin constante y que se puede considerar como funcin de x
o bien de y. En este caso nos interesa considerarla como funcin nicamente de x. Luego
Multiplicamos a (2.52) por lix) y se obtiene
Se puede observar que (2.54) si es una ecuacin exacta, cuya solucin es
(2.54)
no es funcin exclusivamente de x. Por otra parte
no es una funcin exclusivamente de y, por lo que tampoco podemos aplicar el Caso II.
Busquemos un factor integrante de la forma
el cual se puede construir slo si existen constantes m y n, tales que
Esto nos lleva al sistema
cuya solucin es m = 1 y n = 2. De modo que /x(x, y) = xy
2
es un factor integrante para
(2.55). Por lo que al resolver la ecuacin
obtenemos la solucin de (2.55). Su solucin y est definida implcitamente en la ecuacin
EJEMPLO 5. Resolver
Solucin. Como
observamos que Luego la ecuacin no es exacta. Busquemos un factor
integrante. Se tiene que
no es funcin exclusivamente de x, por lo que no podemos calcular un factor integrante
dependiente solamente de x. En forma similar vemos que
no es funcin exclusivamente de y, por lo que el Caso II no puede aplicarse. Vamos ahora
a buscar un factor integrante de la forma /i(x, y) = x
m
y
n
. Para ello buscamos constantes
m y n tales que
De la ltima igualdad se sigue que m = n = 1. As que un factor integrante es
Si multiplicamos la ecuacin diferencial dada por fi(x,y) , obtenemos la ecuacin exacta
cuya solucin est definida implcitamente por la ecuacin
EJEMPLO 6. Resolver
Solucin. La ecuacin no es exacta, ya que
(2.57)
Procedemos a determinar un factor integrante. Primero vemos que la ecuacin no cae
dentro del Caso I, puesto que
no es funcin solo de x. Adems
no es funcin solo de y, y el Caso II falla. Trataremos ahora de hallar un factor integrante
de la forma n{x,y) = x
m
y
n
. Tenemos que
Observamos en la ltima igualdad, que no existen ra, n que la satisfagan. Por ltimo, de
acuerdo con el Caso IV, buscamos un factor integrante de la forma
Es necesario encontrar funciones P{x) y Q{y) tales que
Es decir
En consecuencia Q(y) = 1, P(x) y un factor integrante es
As que multiplicamos la ecuacin (2.57) por este /(x,y) y obtenemos
la cual es una ecuacin diferencial exacta, cuya solucin es la misma que la de la ecuacin
original y est definida implcitamente en la ecuacin
EJERCI CI OS 2.4
Mediante un factor integrante adecuado resuelva las siguientes ecuaciones diferenciales.
2.5 Ecuaciones Diferenciales Lineales
La ecuacin diferencial lineal de primer orden, tiene la forma general
Dividiendo entre ai(x), resulta la forma mas til
Es fcil verificar que la ecuacin diferencial (2.59) tiene como factor integrante a la
funcin
Si multiplicamos la ecuacin (2.59) por (E), se sigue que
(2.60)
2.5. Ecuaciones Diferenciales Lineales 71
e integrando ambos miembros de esta ltima igualdad, obtenemos la solucin general
de la ecuacin.
EJEMPLO 1. Resolver
Solucin. La ecuacin es lineal. Escribamos primero la ecuacin diferencial como
(2.61)
Un factor integrante es
Multiplicando la ecuacin diferencial por (x + 1)
2
, podemos escribir la ecuacin como
o equivalentemente
e integrando
Por lo tanto, la solucin general es
EJEMPLO 2. Resolver
Solucin. La ecuacin es lineal. La escribimos en la forma
(2.62)
De modo que
La solucin general, viene dada por
La condicin inicial y(0) = 1 da c = 1. Por consiguiente la solucin del problema de
valor inicial es
EJEMPLO 3. Resolver
Solucin. Escribimos la ecuacin diferencial en la forma
Multiplicamos por
As que
La solucin general es
EJEMPLO 4. Resolver
(2.64)
2.5. Ecuaciones Diferenciales Lineales 73
Solucin. La ecuacin diferencial (2.64) no es separable, ni homognea, ni exacta, ni
lineal en la variable y. Sin embargo considerando los recprocos, tenemos
o bien
La ltima ecuacin es lineal en x. El factor de integracin correspondiente es
Luego
EJERCICIOS 2.5
Diga si la ecuacin diferencial dada es lineal en y o en x. En caso de serlo determine su
solucin general.
2.6 Ecuacin de Bernoulli
A una ecuacin diferencial de la forma
(2.65)
con n un nmero real, se le llama ecuacin de Bernoulli.
Si n 0 o n 1, (2.65) es una ecuacin diferencial lineal. Adems si n = 1, la
ecuacin se puede resolver mediante separacin de variables. As que nos concentramos
en el caso en que n ^ 0,1.
El mtodo para resolver una ecuacin de Bernoulli consiste en transformarla en una
ecuacin diferencial lineal mediante un cambio de variable, veamos.
Dividiendo ambos lados de (2.65) por y
n
, resulta
Sea
entonces
por lo cual
Sustituyendo (2.67) y (2.68) en la ecuacin diferencial (2.66) obtenemos
75
(2.69)
(2.70)
(2.71)
2.6. Ecuacin de Bernoulli
que es una ecuacin diferencial lineal.
EJEMPLO 1. Resolver
Solucin. Dividiendo la ecuacin (2.69) por y
2
, resulta
Sea
Sustituyendo en (2.70)
Resolviendo la ecuacin diferencial lineal tenemos
y recordando que w
de donde
EJEMPLO 2. Resolver
Solucin. Veamos si es una ecuacin de Bernoulli.
As, efectivamente se trata de una ecuacin de Bernoulli. Dividiendo por y
3
, se sigue que
Sea w Entonces
Resolviendo la ecuacin diferencial lineal se obtiene
EJEMPLO 3. Resolver
Solucin. Ntese que
luego la ecuacin (2.72) no es de Bernoulli en la variable y, pero si la escribimos como
tenemos que
2.6. Ecuacin de Bernoulli 77
la cual es una ecuacin diferencial de Bernoulli en x. Dividiendo por x
2
, resulta
Sea w Sustituyendo en (2.73) resulta
(2.73)
y resolviendo la ecuacin diferencial lineal en w obtenemos
Ya que w se tiene
de donde
EJERCI CI OS 2.6
Resuelva las siguientes ecuaciones diferenciales de Bernoulli.
2.7 Miscelnea de Ecuaciones Diferenciales
No existe un mtodo general para resolver ecuaciones diferenciales. Por consiguiente,
para resolver una ecuacin diferencial recomendamos primero investigar si es de alguno
de los tipos estudiados en las secciones anteriores: separable, homognea, exacta, con
un factor integrante, lineal, etctera, y posteriormente aplicar el mtodo de solucin
correspondiente.
Algunas ecuaciones diferenciales pueden resolverse por varios de los mtodos vistos
en las secciones anteriores.
EJEMPLO. Resolver
(2.74)
Solucin. Veamos primero si la ecuacin (2.74) es de variables separables.
As, (2.74) result de variables separables. Integrando y reduciendo, obtenemos
Por otra parte, si en (2.74) denotamos
entonces
Por lo tanto, la ecuacin (2.74) tambin es exacta. Luego, existe una funcin / tal
que A partir de estas relaciones encontramos que
2.7. Miscelnea de Ecuaciones Diferenciales
De donde
79
La solucin est dada por f(x,y) = C2, con c< una constante, o equivalentemente
Finalmente investiguemos si la ecuacin (2.74) es lineal. Tenemos
De aqu, es claro que no es lineal, pero s es de Bernoulli. Resolvindola como tal,
encontramos
En conclusin, la ecuacin (2.74) se pudo resolver por tres mtodos distintos y por
supuesto el lector desarrollar su propia estrategia al resolver una ecuacin diferencial.
EJERCICIOS 2.7
Resuelva las siguientes ecuaciones diferenciales.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
2.7. Miscelnea de Ecuaciones Diferenciales 81
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
61,
2.7. Miscelnea de Ecuaciones Diferenciales 83
Encuentre el valor de n para el cual ecuacin:
exacta y resulvala para ese valor de n.
Resuelva:
Resuelva:
Resolver la siguiente ecuacin diferencial
Una curva parte desde el origen por el primer cuadrante. El rea bajo la curva
desde (0,0) hasta (x, y) es un tercio del rea del rectngulo que tiene a esos puntos
como vrtices opuestos. Hallar la ecuacin de esa curva.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78. Sea t la tangente a una curva C en un punto P. Sea F el punto sobre el eje x tal
que FP es perpendicular al eje x y sea T el punto de interseccin de t y el eje x.
Encontrar la ecuacin de la familia de curvas C las cuales tienen la propiedad de
que la longitud TF es igual a la suma de la abscisa y de la ordenada de P.
79. Sea t la tangente a una curva C en un punto P. Encontrar la ecuacin de la familia
de curvas C las cuales tienen la propiedad de que la distancia del origen a t es igual
a la abscisa de P.
80. Sea t la tangente a una curva C en el punto P y sea F el punto del eje x tal que
PF es perpendicular a dicho eje. Encontrar la ecuacin de la familia de curvas C
las cuales tienen la propiedad de que la distancia de F a t es constante.
2.8 Cambios de Variables Diversos
El cambio de variable es una idea genrica en Matemticas, la cual como se ha visto
en las secciones anteriores, tambin puede aplicarse al tratar de encontrar soluciones de
ecuaciones diferenciales. Sin embargo, deseamos enfatizar que no hay un mtodo general
para proponer cambios de variable. Veamos algunos ejemplos adicionales.
EJEMPLO 1. Resolver
Solucin. Sea Entonces
(2.75)
Sustituyendo en (2.75)
As, la ecuacin se ha transformado en la ecuacin diferencial lineal
cuya solucin es
Usando que u = e
2y
, encontramos la solucin de la ecuacin original en la forma
2.8. Cambios de Variables Diversos 85
EJEMPLO 2. Resolver
Solucin. Sea u = \ny. Tenemos que
y como y
Sustituyendo en (2.76) resulta
Luego, hemos obtenido una ecuacin diferencial lineal con solucin
Por lo tanto
de donde
EJEMPLO 3. Resolver
(2.77)
Solucin. Sea u Entonces
Note que la ecuacin diferencial que acabamos de obtener es separable. Resolvindola
llegamos a que
En consecuencia
EJERCI CI OS 2.8
Resuelva las siguientes ecuaciones diferenciales
1.
2.
3.
4.
5.
6.
7.
8.
9.
Haciendo el cambio de variable z = y/x
n
, esto es, y = zx
n
y escogiendo un valor
adecuado de n, demuestre que las ecuaciones diferenciales siguientes pueden trans-
formarse en ecuaciones de variables separables y resulvalas:
Resuelva
Resolver las ecuaciones
(Sugerencia: haga v = y')
2.9. Ecuacin de Ricatti 87
2.9 Ecuacin de Ricatti
La ecuacin
se llama ecuacin de Ricatti. El mtodo de solucin que proponemos para resolver
(2.78) supone que se conoce, o bien puede calcularse, una solucin particular de la
ecuacin. Sea y\ una solucin particular de (2.78), entonces es fcil mostrar que la
relacin
define una familia de soluciones de la ecuacin de Ricatti, donde v satisface cierta ecuacin
diferencial lineal, como veremos a continuacin.
Derivando (2.79) con respecto a x obtenemos
Sustituyendo esta expresin junto con (2.79) en (2.78) tenemos
o equivalentemente
Ahora bien, ya que que y\ es solucin de la ecuacin diferencial, se cumple que
lo cual reduce (2.80) a
Es decir, v satisface la ecuacin
que efectivamente es una ecuacin diferencial lineal.
EJEMPLO . Resolver
Solucin. Es claro que (2.82) es una Ecuacin de Ricatti con p(x) 1, q(x) =
- 1/ x, r(x) = 1/x
2
. Observando que y\ 1/x es una solucin particular de (2.82),
en este caso la ecuacin (2.81) para v toma la forma
esto es
Un factor integrante de (2.83) es
(2.83)
Multiplicando (2.83) por fi(x) y resolviendo igual que antes, encontramos que
donde c = 2ci. Por lo tanto, de acuerdo con (2.79)
EJERCICIOS 2.9
Resuelva las siguientes ecuaciones diferenciales de Ricatti, utilizando la solucin particu-
lar y\ indicada.
2.10. Ecuacin de Clairaut 89
2.10 Ecuacin de Clairaut
La ecuacin diferencial de la forma
(2.84)
(2.85)
con / una funcin arbitraria, se llama ecuacin de Clairaut.
Mostraremos que las rectas
con c un nmero real, constituyen una familia de soluciones de la ecuacin de Clairaut
y que adems dicha ecuacin posee la siguiente solucin definida en forma paramtrica
por las ecuaciones
(2.86)
En primer lugar, derivando (2.85) con respecto a x se sigue que
y sustituyendo en (2.84) vemos que sta se reduce a la identidad
lo cual muestra que (2.85) es una solucin de (2.84) para todo c. Obsrvese que la
expresin (2.85) de esta familia de soluciones se obtiene trivialmente poniendo y
1
c en
el lado derecho de (2.84).
Por otra parte, de las ecuaciones (2.86) tenemos que
Sustituyendo y en la ecuacin diferencial de Clairaut llegamos a la identidad
De modo que efectivamente (2.86) definen otra solucin de (2.84).
A la solucin paramtrica se le llama solucin singular y requiere que /"() ^ 0.
Note adems que esta solucin no se puede obtener de la familia de rectas (2.85).
EJEMPLO. Resolver
(2.87)
Solucin. La familia de soluciones es
Identificando f(y') = 1 lny', la solucin singular est dada por
y como
resulta
de donde, y = 2 + \nx es la solucin singular.
EJERCICIOS 2.10
con a una constante.
Captulo 3
Aplicaciones de las Ecuaciones
Diferenciales de Primer Orden
3.1 Trayectorias Ortogonales
Se dice que dos curvas son ortogonales si se intersectan y en los puntos de corte sus
rectas tangentes son perpendiculares entre s.
Si todas las curvas de una familia de curvas F(x,y,C\) 0 son ortogonales a todas
las curvas de otra familia G(x\ y, C2) = 0, entonces se dice que las familias son, cada una,
trayectorias ortogonales de la otra.
Una aplicacin elemental de las trayectorias ortogonales es la siguiente. Se tiene
un imn y se han esparcido limaduras de hierro alrededor de l. Ver figura 3.1. Las
Figura 3.1: Lneas equipotenciales.
lneas punteadas (las limaduras) son las lneas de fuerza. Las lneas continuas son las
trayectorias ortogonales y se llaman lneas equipotenciales (lneas de igual potencial).
91
92 Captulo 3. Aplicaciones de las Ecuaciones Diferenciales de Primer Orden
EJEMPLO 1. Encontrar las trayectorias ortogonales de la familia de parbolas y = ex
2
.
Solucin. De la ecuacin de la familia dada y ex
2
se sigue que
Sustituyendo
obtenemos
Luego, las trayectorias ortogonales deben cumplir
Resolviendo la ecuacin diferencial encontramos
As, las trayectorias ortogonales de las parbolas con vrtice en el origen y cuyo eje
es el eje y, son elipses con centro en el origen y eje mayor en el eje x.
EJEMPLO 2. Determine el miembro de la familia de trayectorias ortogonales de
que pasa por (0,4).
Solucin. Se tiene que
Derivando esta expresin obtenemos
de donde
3.2. Mecnica 93
Luego, la familia de trayectorias ortogonales debe satisfacer la ecuacin diferencial
La solucin general de dicha ecuacin es
Finalmente, de la condicin inicial y(0) = 4 resulta c = 4. Por lo tanto la curva buscada
es
3.2 Mecnica
EJEMPLO 1. Al abrir su paracaidas un paracaidista est cayendo con una velocidad
de 176 ft/s, si la fuerza de resistencia del aire es Wv
2
/256 Ib, donde W es el peso total
del hombre y del parcaidas y v la velocidad con que va cayendo, hallar la velocidad en
cualquier tiempo despus de abierto el paracaidas.
Solucin. El diagrama de cuerpo libre es
Figura 3.2: Diagrama de cuerpo libre
De la segunda ley de Newton, usando que W mg y g = 32 ft/s
2
, se sigue que
o equivalentemente
(3.1)
La ecuacin diferencial (3.1) es separable, resolvindola y usando la condicin inicial
v(0) = 176, obtenemos
de donde se observa que para t muy grande v(t) tiende al valor de 16 ft/s. Esta es la
llamada velocidad terminal y es la velocidad con que el cuerpo llega al suelo.
EJEMPLO 2. Pensar en un viaje interplanetario antes de mediados del siglo XX era
ubicarse en el terreno de la ficcin, pero hoy es una realidad. Consideremos un viaje a la
Luna. Con qu velocidad debe salir una nave de la Tierra para poder llegar a la Luna?
Solucin. La figura ?? ilustra la situacin.
Figura 3.3: Viaje a la Luna.
Denotamos por Mr, M, m, r y d a la masa de la Tierra, la masa de la Luna, la masa
de la nave, la distancia del centro de la Tierra a la nave, y la distancia de la Luna a la
nave, respectivamente. Vase la figura 3.3.
De la ley de la gravitacin universal de Newton tenemos que sobre la nave actan dos
fuerzas
Si ignoramos la influencia de la Luna y dems planetas distintos a la Tierra, as como
otras fuerzas de resistencia entonces
3.2. Mecnica 95
es decir
con K una constante positiva. El valor de K puede expresarse en trminos de
R y g como sigue. Cuando r
As
Por otra parte
En consecuencia
(3.2)
Resolviendo la ecuacin diferencial (3.2) mediante separacin de variables obtenemos
Ahora, si v = i>o para r = R entonces
de donde
y
Note que si vi 2gR < 0 entonces existe un valor de r tal que v ser igual a cero, lo
cual implicara que la velocidad v cambiara de positiva a negativa y la nave volvera a
la Tierra.
Por lo tanto, para que la nave escape de la Tierra debemos pedir que
o equivalentemente
De aqu que la velocidad mnima llamada velocidad de escape es
Tomando en cuenta que R = 3960 millas y g encontramos millas
millas
3,3 Aplicaciones a problemas relacionados con cre-
cimiento y decrecimiento
3.3.1 Desintegracin Radiactiva.
Ley de desintegracin radiactiva. La velocidad de desintegracin de una sustancia
radiactiva en un instante dado es proporcional a la cantidad de sustancia presente en ese
instante.
La vida media de una sustancia radiactiva es el tiempo necesario para que se desin-
tegren la mitad de los tomos de una cantidad inicial de dicha sustancia.
EJEMPLO 1. La velocidad con que se desintegran ncleos radiactivos es proporcional
al nmero de ncleos que estn presentes en una muestra dada. La mitad del nmero
original de ncleos radiactivos ha experimentado la desintegracin en un perodo de 1500
aos.
a)Qu porcentaje de ncleos radiactivos originales continuarn despus de 4500 aos?
b)En cuntos aos quedar solamente un dcimo del nmero original de ncleos radiac-
tivos?
Solucin. Sea x(t) la cantidad de ncleos radiactivos presente despus de t aos y sea
XQ el nmero original de ncleos radiactivos. Entonces x(0)
la velocidad con la que se desintegran los ncleos al tiempo t.
As , este problema queda formulado por la siguiente ecuacin diferencial
dnde k es la constante de proporcionalidad, junto con las condiciones
La solucin de la ecuacin (3.3) es ya conocida
Usando la condicin inicial x(0) = XQ encontramos que
a) Para calcular el porcentaje de ncleos radiactivos originales despus de 4500 aos,
determinamos x(4500). Considerando que z(1500) = x
o
/2 obtenemos
3.3.1. Desintegracin Radiactiva 97
Sustituyendo k en (3.4) resulta
Luego
lo cual nos dice que despus de 4500 tenemos un 12.5% de x
0
.
b) Para determinar en cuntos aos quedar solamente un dcimo del nmero original
de ncleos, es necesario hallar el valor de t tal que x(t) = XQ/10, es decir
EJEMPLO 2. Se sabe que cierto material radiactivo se desintegra a una razn propor-
cional a la cantidad presente. Si inicialmente hay 50 miligramos de material y despus
de dos horas se observa que el material ha perdido el 10% de su masa original, encuentre
a) Una expresin para la masa de material restante en un momento t.
b)Cuntos miligramos del material quedan despus de cuatro horas?
c) Cul es la vida media de este material?
Solucin. Sea x(t) la masa del material restante despus de cierto tiempo t. Como al
cabo de dos horas el material se ha desintegrado el 10% de su masa original, es decir
dx
el 10% de 50 mg que son 5 mg, tenemos que x(2) = 45 mg. Igual que antes, es la
CLL
velocidad con que se desintegra el material radiactivo.
As este problema queda formulado con la siguiente ecuacin diferencial y sus condi-
ciones
con k una constante de proporcionalidad, y las condiciones
Sabemos que la solucin general de (3.5) es
(3.7)
Empleando la condicin inicial (3.6), resulta
por lo cual
y
Por otra parte, de (3.7) tenemos que
b) El nmero de miligramos del material despus de cuatro 4 horas es
c) Para calcular la vida media, determinamos el valor de t para el cual x(t)
Es decir,
Por lo tanto la vida media de este material es de 13 horas.
Mtodo del Carbono 14 Este mtodo se debe al qumico Willard Libby cuyo
descubrimiento le vali el Premio Nobel de Qumica en 1960. La teora se basa en lo
siguiente. La atmsfera terrestre es continuamente bombardeada por rayos csmicos, los
cuales producen neutrones libres que se combinan con el nitrgeno de la atmsfera para
3.3.1. Desintegracin Radiactiva 99
producir el istopo C-14 (Carbono 14 o bien radiocarbono). Este C-14 se combina con el
bixido de carbono presente en la atmsfera, el cual es absorbido por las plantas y stas
a su vez son alimento para los animales. As es como se incorpora el radiocarbono a los
tejidos de seres vivos.
El cociente de la cantidad de C-14 y la cantidad de carbono ordinario presentes en
la atmsfera es constante, y en consecuencia la proporcin de istopo presente en todos
los organismos vivos es la misma que en la atmsfera. Cuando un organismo muere,
la velocidad de incorporacin de radiocarbono a l se hace nula y entonces comienza el
proceso de desintegracin radiactiva del C-14, que se encontraba presente en el momento
de su muerte. As comparando la proporcin de C-14 que hay en un fsil con la proporcin
constante encontrada en la atmsfera es posible obtener una estimacin razonable de su
edad.
EJEMPLO 3. Se ha encontrado que un hueso antiguo contiene | de la cantidad original
de C-14 de un hueso al tiempo actual. Cul es la antigedad del fsil?
Solucin. Sea x(i) la cantidad de C-14 presente en el hueso al tiempo t y sea x
0
la
cantidad de C-14 cuando se form la muestra, es decir x(0) = x
Q
. La vida media del C-14
es de 5,568 aos, por lo cual x(5568) Adems es la velocidad de desintegracin
radiactiva del C-14.
Determinaremos la edad del fsil al encontrar el valor de t para el cual x(t)
Para eso, partimos de que
cuya solucin es
Considerando que x(5568) = XQ/2, obtenemos
y as
Buscamos el valor de t para el cual
Tenemos que
As, el fsil tiene una antigedad de 16705 aos.
EJEMPLO 4. En 1950 se hicieron excavaciones en Nipur (Babilonia), en las cuales se
encontraron muestras de carbn que reportaron 4.09 desintegraciones por minuto y por
gramo. Una muestra actual report 6.68 desintegraciones por minuto y por gramo. Se
sabe que la primer muestra se form en la poca del reinado de Hammurabi. Con estos
datos, determine hace cuanto tiempo Hammurabi rein en Babilonia.
Solucin. Sea x(t) la cantidad de C-14 presente en el tiempo t. Entonces
velocidad de desintegracin del C-14 al tiempo y
es la
Sabemos por la ley de decaimiento radiactivo que el modelo a seguir es
(3.8)
(3.9)
(3.10)
(3.11)
Adems
Como se vio en el problema anterior
Sustituyendo (3.10) en (3.9), se tiene
Considerando (3.8) en (3.11), resulta
donde k es una constante de proporcionalidad.
EJEMPLO 1. Una barra metlica a una temperatura de 100F se pone en un cuarto
a una temperatura constante de 0F. Despus de 20 minutos la temperatura de la barra
es 50F.
a) Cunto tiempo tardar la barra para llegar a una temperatura de 25F?
b) Cul ser la temperatura de la barra despus de 10 minutos?
Solucin. Sea T(t) la temperatura de la barra al tiempo t, luego JT(O) = 100^ y
TY20) = 50F. La temperatura del medio ambiente, T
A
, es T
A
= 0F. Ntese que es
dt
la velocidad a que se enfra la barra.
Por la ley de enfriamiento de Newton se tiene que
3.3.2. Problemas de Enfriamiento 101
Ahora bien, para determinar hace cuanto tiempo Hammurabi rein en Babilonia, ten-
dremos que calcular para que valor de , se cumple que 4.09
Aproximadamente hace 3941 aos que Hammurabi rein en Babilonia.
3.3.2 Problemas de Enfriamiento
Ley de enfriamiento de Newton. En un cuerpo que se est enfriando la tasa de
cambio de la temperatura T(t) con respecto al tiempo t es proporcional a la diferencia
entre la temperatura del cuerpo T(t) y la temperatura T
A
del medio que lo rodea. Esto
es
y como T
A
= 0, este problema queda formulado con la siguiente ecuacin diferencial y
sus condiciones
La solucin general de la ecuacin diferencial es ya conocida
Como T(0) = 100 se tiene que
Usando adems que T(20) = 50 resulta
(3.12)
Sustituyendo k en (3.12) tenemos que
a) El tiempo necesario para que la temperatura de la barra sea de 25^ se obtiene
resolviendo la ecuacin T() = 25, esto es
de donde
As que la barra tardar 40 minutos en alcanzar una temperatura de 25F.
b) La temperatura de la barra despus de 10 minutos es igual a
es decir, ser aproximadamente de 71F.
EJEMPLO 2. Un cuerpo a una temperatura desconocida se pone en un refrigerador a
una temperatura constante de I
o
F. Si despus de 20 minutos la temperatura del cuerpo
es de 40F y despus de 40 minutos la temperatura del cuerpo es de 20F, hallar la
temperatura inicial de ste.
Solucin. Denotemos nuevamente con T(t) a la temperatura del cuerpo en un instante
dado. As 7(20) = 40F, T(40) = 20F
;
es la velocidad con que se enfra el cuerpo.
Ahora la temperatura constante del medio ambiente es T& = lF.
Por la ley de enfriamiento de Newton, este problema se formula de la siguiente forma
103
(3.13)
(3.14)
(3.15)
(3.16)
3.3.2. Problemas de Enfriamiento
La solucin general de la ecuacin diferencial es
Para obtener c y k utilizamos las condiciones dadas, como siempre.
de donde
y
Aplicando logaritmo natural en (3.14) y (3.15) se obtiene
De aqu que
o bien
Sustituyendo (3.16) en (3.14) resulta
Usando los valores de c y k en (3.13), obtenemos que
Luego
La temperatura inicial del cuerpo era de 81F.
3.3.3 Modelos de Poblacin
Sea x(t) el nmero de individuos en el tiempo t. La ley de Malthus de crecimiento de
poblaciones dice que la razn de cambio de la poblacin es proporcional al nmero de
individuos en ese tiempo, es decir
Este modelo lineal para crecimiento de poblaciones, son satisfactorios siempre que la
poblacin no sea demasiado grande o bien que no se aplique a un futuro distante.
Cuando la poblacin es demasiado grande, este modelo no puede ser exacto, ya que
no refleja el hecho de que los individuos compiten entre s por el limitado espacio vital,
por recursos naturales, etc. As pues, hay que agregar un trmino de competicin para
que el crecimiento de la poblacin est representado en forma ms realista. Una eleccin
adecuada del trmino competitivo es bx
2
, llamada ley logstica (Verhulst, en 1837):
Ahora bien, en general la constante b es muy pequea comparada con a, de tal modo que
si x no es demasiado grande, entonces el trmino bx
2
es insignificante comparado con
ax. Sin embargo, si x es grande entonces el trmino bx
2
debe tomarse en cuenta ya que
disminuye la tasa de crecimiento.
EJEMPLO 1. En un cultivo de bacterias se tenan x nmero de familias. Despus de
una hora se observaron en el cultivo 1000 familias de la bacteria y despus de cuatro
horas, 3000 familias. Encontrar la expresin para el nmero de familias de la bacteria
presentes en el cultivo al tiempo t y el nmero de familias de la bacteria que haba
originalmente en el cultivo.
Solucin. Sea x(t) el nmero de familias de la bacteria que hay en t horas. De ah que
es la velocidad a la que crece el cultivo de bacterias.
Por la ley malthusiana este problema se formula de la siguiente manera
cuya solucin es ya conocida
y considerando las condiciones se tiene que
3.3.3. Modelos de Poblacin 105
es la expresin que nos d el nmero de familias presentes en un momento t.
Observamos que el nmero de familias que haba originalmente en el cultivo es
familias.
EJEMPLO 2. La poblacin x(t) de una cierta ciudad satisface la ley logstica
donde el tiempo t se mide en aos. Suponiendo que la poblacin de esta ciudad es 100,000
en 1980, determine:
a) La poblacin como una funcin del tiempo t.
b) La poblacin en el ao 2000.
c) El ao en que se duplicar la poblacin de 1980.
d) El comportamiento de la poblacin cuando t > oo.
Solucin.
a) Debemos resolver el problema de valor inicial
Separando variables, tenemos que
y descomponiendo en fracciones parciales el miembro izquierdo de esta ecuacin, encon-
tramos que
Al integrar ambos lados, resulta
al despejar x llegamos a que
(3.17)
Empleando la condicin inicial :r(1980) = obtenemos el valor de
Sustituyendo el valor de c
2
en (3.17) y simplificando se tiene que
b) La poblacin en el ao 2000 es
es decir en el ao 2000 habr aproximadamente 119,500 habitantes.
c) Para encontrar el ao en que se duplicar la poblacin de 1980 buscamos el valor de
tal que
Para el ao del 2061 tendremos duplicada la poblacin de 1980.
d) Tenemos que
Luego, en el transcurso de los aos la poblacin de esta ciudad se estabilizar en un
milln de habitantes.
EJEMPLO 3. Este es un modelo para la propagacin de una infeccin o un rumor
en una poblacin fija. Supngase que un estudiante portador de un virus de gripe,
regresa a un campus universitario, aislado, que tiene 1000 estudiantes. Supongamos que
la rapidez con que el virus se propaga, es proporcional no slo al nmero de estudiantes
contagiados, sino tambin, al nmero de estudiantes no contagiados. Determinar el
3.3.3. Modelos de Poblacin 107
nmero de estudiantes contagiados despus de 6 das, si adems se observa que despus
de 4 das ya eran 50 los contagiados.
Solucin. Denotemos con x(t) al nmero de estudiantes contagiados en t das. Entonces
expresa el nmero de estudiantes no contagiados y
es la velocidad con la que aumenta el nmero de estudiantes contagiados. Por hiptesis
es proporcional a
Este problema queda formulado as
Podemos observar que
es la ecuacin logstica con a = 1000/c y b = k. Separamos variables en (3.18) y por
fracciones parciales se tiene que
Integrando en ambos lados, obtenemos
y simplificando, se tiene
de donde
Como x(0) = 1 tenemos que c<i = 1/999 y sustituyendo el valor de c
2
en (3.19), x(t)
queda de la forma
o bien
Adems x(4) = 50, por lo cual
De esta expresin despejamos k y obtenemos que
As , sustituyendo el valor k en (3.20), tenemos que x(i) queda al fin de la forma
El nmero de estudiantes contagiados despus de 6 das est dado por
es decir, 276 estudiantes han sido contagiados.
3.4 Mezclas
Vamos a considerar ahora los problemas relacionados con mezclas, en los cuales se supone
que una sustancia S fluye hacia una mezcla en un recipiente, con una cierta rapidez, y la
mezcla se mantiene uniforme mediante agitacin. Adems, la mezcla uniforme sale del
recipiente y pasa a otro. Nos interesa determinar la cantidad de la sustancia 5 presente
en la mezcla para el tiempo t.
Si denotamos por A(t) la cantidad de S al tiempo , entonces la derivada es la
razn de cambio de A con respecto a t. Si R\ indica la razn, rapidez o tasa con la que S
entra a la mezcla y R
2
representa la razn con la que sale, tenemos la ecuacin diferencial
lineal bsica
de la cual determinaremos la cantidad A{t) de S en el tiempo t. A continuacin pre-
sentaremos algunos ejemplos.
EJ EMPLO 1. Un gran tanque est parcialmente lleno con 200 gal de agua en las cuales
se disuelven 20 Ib de sal. Una salmuera que contiene 2 Ib de sal por galn, se bombea al
tanque con una rapidez de 6 gal/min y la mezcla bien agitada sale a la misma tasa.
a) Halle el nmero de libras de sal en el tanque en cualquier tiempo.
b) Cunta sal est presente despus de 30 min?
3.4. Mezclas 109
c) Cunta sal estar presente despus de un tiempo largo?
Solucin. Denotemos con A(t) el nmero de libras de sal en el tanque despus de t
minutos. Entonces mide la tasa de cambio de A(t) con respecto al tiempo.
Por conservacin de masa, tenemos que
(3.21)
donde Ri y i?2 son la rapidez con que entra y sale la sal del tanque, respectivamente.
Sean G\ y G2 el gasto volumtrico de las soluciones de entrada y salida al tanque y
Ci, C2 sus concentraciones de sal. Entonces
En consecuencia, la ecuacin (3.21) se reduce a
o equivalentemente
la cual resolvemos sujeta a la condicin inicial A(0) = 20.
a) La solucin a este problema de valor inicial es
que nos da la cantidad de sal al tiempo i (en minutos),
b) Despus de 30 minutos la cantidad de sal es
c) Despus de un tiempo largo, esto es, cuanto t tiende a infinito, vemos que A se aproxima
al valor de 400 Ib.
EJEMPLO 2. Suponga ahora que en el ejemplo anterior la solucin adecuadamente
mezclada se bombea hacia afuera a una tasa de 4 gal/min. Determine A(t).
Solucin. El volumen V{t) de la solucin en el tanque vara a una razn de
luego
as que
Por consiguiente tenemos ahora
o bien
junto con la condicin inicial A(0) 20.
Resolviendo
EJEMPLO 3. Un tanque contiene inicialmente 60 galones de agua pura. Entra al
tanque, a una tasa de 2 gal/min, salmuera que contiene 1 Ib de sal por galn, y la
solucin (perfectamente mezclada) sale de l a razn de 3 gal/min. Obtenga el nmero
de libras A{t) de sal que hay en el tanque en un instante cualquiera. Cunto demorar
el tanque en vaciarse? Cul es la mxima cantidad de sal que llega a tener el tanque?
Solucin. Continuaremos usando la notacin introducida en los ejemplos anteriores. En
este caso tenemos
Por lo tanto, la ecuacin diferencial es
es decir
La condicin inicial es
(3.22)
(3.23)
3.4. Mezclas 111
La solucin del problema de valor inicial (3.22)-(3.23) viene dada por
El tanque se vaca despus de 60 minutos. Por otro lado
de modo que A'(i) == 0 si y slo si
concluye que la cantidad mxima de sal que llega a tener el tanque es
Como
EJEMPLO 4. Una cierta presa, en su mxima capacidad, contiene 1,000 millones de m
3
de agua. En un instante dado, estando llena la presa, tiene una masa de 2 toneladas de
contaminantes, distribuida en forma homognea. Suponga que en temporada de lluvias
entra agua a la presa a razn de 10 millones de m
3
por da, con una masa de contaminantes
de 0.09% toneladas por milln de m
3
de agua y sale con la misma rapidez. Determine
la cantidad de contaminantes en la presa en cualquier instante. En cunto tiempo se
reducir la contaminacin total de la presa a 1.2 toneladas?
Solucin. Denotemos con A{t) el nmero de toneladas de contaminantes despus de t
das.
En este caso, tenemos la ecuacin diferencial
junto con la condicin inicial ^4(0) = 2. La solucin est dada por
Buscamos ahora el valor de t para el cual A(i) = 1.2, es decir
de donde se obtiene el valor de t 129.9 das.
EJEMPLO 5. Un tanque contiene inicialmente 100 di de agua, en el cual se disuelven
80 kg de sal. Se introduce en el tanque agua pura a velocidad de 4 dl/min y la mezcla,
conservada homognea mediante agitacin, sale a la misma velocidad y va a parar a un
segundo tanque que contiene al principio 100 di de agua pura. Agitando se mantiene
homognea la mezcla que sale de este segundo tanque a la misma velocidad ya citada.
Hallar la cantidad de sal en el segundo tanque al cabo de 1 hora.
Solucin. Denotemos por A\(t) y A
2
(t) la cantidad de sal en los tanques uno y dos,
respectivamente. Para el primer tanque, tenemos que
y resolviendo este problema de valor inicial, obtenemos
Para el segundo tanque, la concentracin de la solucin de entrada est dada por
Luego
o equivalentemente
La solucin de esta ecuacin diferencial, junto con la condicin inicial ^( O) = 0 es
Por lo tanto la cantidad de sal en el segundo tanque despus de una hora es
3.5 Circuitos Elctricos LR y RC en Serie
En Mecnica se tiene como base fundamental las leyes de Newton, de manera anloga,
en electricidad se cuenta con las leyes de Kirchhoff que describen el comportamiento de
los circuitos elctricos. En particular estamos interesados en aplicar la Segunda Ley de
Kirchhoff, que enunciaremos mas adelante.
Para el estudio de los circuitos LR y RC en serie haremos uso de las Tablas 1 y 2.
3.5. Circuitos Elctricos 113
Tabla 1
Cantidad
Voltaje, fem o potencial
Resistencia
Inductancia
Capacitancia
Corriente
Carga Elctrica
Smbolo
E
R
L
C
i
q
Unidad
Volt (V)
Ohm (0)
Henry (H)
Farad (F)
Amper (A)
Coulomb (C)
Tabla 2
Elemento
Resistencia
Inductor
Condensador
Cada de Potencial
E = Ri
Se acostumbra indicar los diferentes elementos de un circuito como se ilustra en la
figura 3.4.
Figura 3.4: Elementos de un circuito
El siguiente es un enunciado de la Segunda Ley de Kirchhoff.
Segunda Ley de Kirchhoff. La suma algebraica de todas las cadas de potencial en
cualquier camino cerrado de un circuito elctrico es igual a cero.
SMBOLO ELEMENTO
Generador o batera
Resistencia
Inductor o bobina
Condensador
Llave o interruptor
Convencin. La corriente fluye del lado positivo (+) de la batera o generador a travs
del circuito hacia el lado negativo (-).
3.5.1 Circuito LR en Serie
Consideremos el circuito elctrico que se muestra en la figura 3.5.
vw
R
Figura 3.5: Circuito LR en serie
Aplicando la segunda Ley de Kirchhoff a este circuito, la suma de las cadas de poten-
cial a travs del inductor Ldi/dt y de la resistencia Ri, es igual a la fuerza electromotriz
(fem) E(t) aplicada al circuito, es decir
di
=:
dt
(3.24)
Como se observa la ecuacin diferencial que describe el comportamiento de la corriente
elctrica a travs del circuito es una ecuacin diferencial lineal y puede resolverse con el
mtodo descrito anteriormente.
EJEMPLO 1. Resuelva la ecuacin (3.24) si E(t) = E
o
, donde L, R y E
o
son constantes.
Solucin. En este caso la ecuacin diferencial para el circuito LR en serie, toma la forma
di R. _E
0
Resolviendo se obtiene
donde c es una constante arbitraria que depende de la condicin inicial i(0) = o-
Ntese que independientemente del valor de
0
, cuando t tiende a infinito el segundo
trmino de la solucin se aproxima a cero. A un trmino como este se le llama usualmente
trmino transitorio o corriente transitoria. Al trmino (o trminos) restante se le llama
3.5.1. Circuito LR en Serie 115
parte estacionaria de la solucin. En este caso E
o
/R es la corriente estacionaria; esto es,
cuando t tiende a infinito, la corriente se comporta como si estuviese gobernada por la
ley de Ohm (i = E
o
/R).
EJEMPLO 2. Un generador con una fem de 50 V se conecta en serie con una resistencia
de 6 7 y un inductor de 2 henrys. Si el interruptor K se cierra a i = 0, determine la
corriente para todo t. Ver figura 3.6
Solucin. La ecuacin diferencial del circuito es
o equivalentemente
sujeta a la condicin inicial i(0) = 0.
^ + 6 = 50
dt
^ + 3 = 25
dt
E(t) = 50 V Ct
L = 2H
^v^.
\
Figura 3.6: Circuito del ejemplo 2
Resolviendo la ecuacin diferencial se obtiene
i(t) +
25
y empleando la condicin inicial, encontramos que c =
j
.
Por lo tanto
f
EJEMPLO 3. Determine i(t) para el circuito elctrico del problema anterior si el
generador de 50 V se reemplaza por otro con una fem de E(t) = 10sen7 (figura 3.7).
Solucin. Ahora la ecuacin diferencial se reduce a
di
+ 3z = 5 sen 7
dt
cuya solucin general est dada por
i(t) = ( 3s en7 - 7cos7) +ce~
3t
.
58
De la condicin inicial se sigue que c = , as que
58
5 35
i(t) = (3sen7 - 7cos7) + e~
3
.
58 58
L = 2H
E(t)=10sen7tV ( i
= 6 Q
\K
3.5.2 Circuito RC en Serie
Estudiaremos ahora el circuito elctrico de la figura 3.8.
Procediendo en forma anloga a nuestra discucin anterior, aplicamos la segunda Ley
de Kirchhoff y los resultados de la Tabla 2, al circuito RC en serie. Obtenemos
(3.25)
(3.26)
Pero la corriente i y la carga q estn relacionadas por
dq
por lo cual (3.25) se transforma en la ecuacin diferencial lineal
3.5.2. Circuito RC en Serie 117
E(t)Q
Figura 3.8: Circuito RC en serie
EJEMPLO 4. Una batera cuya fem est dada por E(t) = 200e~
5
se conecta en serie
con una resistencia de 20 l y un condensador de 0.01 F. Suponiendo que q(0) = 0
encuentre la carga y la corriente en cualquier tiempo. Muestre que la carga alcanza un
mximo, calcule su valor y halle el valor de t para el cual se alcanza.
Solucin. La ecuacin diferencial para la carga elctrica es
20-^ + lOOq = 200e"
5
.
dt
Resolviendo (3.27) sujeta a la condicin inicial q(0) = 0, obtenemos
q(t) = 10te~
5t
A partir de aqu y usando (3.26), se sigue que
(3.27)
Finalmente, empleando el criterio de la segunda derivada encontramos que
q-max q(-z) = 0.74 coulombs.
5
Puede el lector interpretar fsicamente el comportamiento de q(t)l
EJEMPLO 5. Una resistencia de R Q vara con el tiempo t (en segundos) de acuerdo
a R = 1 + O.Ol. Se conecta en serie con un condensador de 0.1 F y un generador con
una fem de 100 V. La carga inicial en el condensador es de 5 coulombs. Encuentre
a) La carga y la corriente como una funcin del tiempo.
b) La carga mxima terica.
Solucin. El circuito se muestra en la figura 3.9.
AAAA
= 1+0.01tQ
= 100v
Q)
C = 0. 1F
Aplicando la segunda ley de Kirchhoff al circuito obtenemos
(1 + 0.01t)i + lOq = 100.
Sustituyendo (3.26) en (3.28) se sigue que
dq lOq 100
dt 1 + O.Ol 1 + O.Ol
a) La solucin general de la ecuacin diferencial (3.29) es
q(t) = 10 + c(l + O.Ol)"
1000
.
Empleando la condicin inicial #(0) = 5, encontremos que c = 5, por lo cual
q(t) = 10 - 5 ( 1 + O.Ol)"
1000
y consecuentemente
(3.28)
(3.29)
= ^ = 50(1 + 0.01)-
1001
.
dt
dq
b) Ya que > 0 para todo > 0, q(t) es una funcin creciente. De manera que la carga
mxima terica est dada por
Qmax = lim q(t) = 10 coulombs.
Ejercicios del Captulo 3 119
EJERCI CI OS 3
1. Encuentre las trayectorias ortogonales para la familia de curvas dada. Trace la
grfica de algunas curvas de la familia y de las trayectorias ortogonales.
a) y ce~
x
b) x
2
+ y
2
= 2cx
2. Determine el miembro de la familia de trayectorias ortogonales a la familia de curvas
dada, que pase por el punto indicado
a) y = csenx, (0>2).
3. El uranio se descompone a una velocidad proporcional a la cantidad presente. Si
inicialmente hay 10 g y despus de 2 horas se ha perdido el 5% de su masa original,
hallar
a) La cantidad restante de uranio como funcin del tiempo.
b) La cantidad de uranio despus de 5 horas.
4. Cierto material radiactivo se desintegra con una rapidez proporcional a la cantidad
existente en cada instante. En una prueba realizada con 60 mg de este material, se
observ que despus de 3 horas, solamente permaneca el 80% de la masa original.
Hallar
a) La cantidad restante de masa en cualquier instante.
b) Qu cantidad de material hay despus de 5 horas?
c) Cuanto tiempo debe transcurrir para que la cantidad de material sea un
cuarto de la cantidad inicial?
5. Se ha observado en el laboratorio que el radio se desintegra a una rapidez pro-
porcional a la cantidad de radio presente. Su vida media es de 1600 aos. Qu
porcentaje desaparecer en un ao?
6. En un cultivo de levadura la rapidez de cambio es proporcional a la cantidad exis-
tente. Si la cantidad de cultivo se duplica en 4 horas, qu cantidad puede esperarse
al cabo de 12 horas?
7. Un cultivo tiene inicialmente una cantidad AT
0
de bacterias. Para t = 1 hora, en
nmero de bacterias medido es |AT
0
. Si la rapidez de multiplicacin es proporcional
al nmero de bacterias presentes, determine el tiempo necesario para que el nmero
de bacterias se triplique.
8. En cierto zoolgico se ha observado que la cantidad de animales aumenta propor-
cionalmente al nmero actual de dichos animales. Si despus de 5 aos su nmero
se ha duplicado y despus de siete aos el nmero de animales es 576, hallar el
nmero de animales con que se contaba el da de la inauguracin.
9. Supngase que la poblacin P de bacterias en un cultivo al tiempo , cambia a una
razn directamente proporcional a P
2
P. Si inicialmente hay 1000 bacterias y
despus de 5 horas la poblacin se redujo a 100 bacterias, determine:
a) La poblacin como funcin del tiempo.
b) La poblacin despus de un tiempo grande.
10. Al apagar un motor su temperatura es de 98C y el medio en que se encuentra
se conserva a 21C Si despus de 10 minutos el motor se ha enfriado a 88C,
encuentre:
a) La temperatura del motor como funcin del tiempo.
b) El instante en el cual su temperatura es de 35C
11. Un cuerpo a una temperatura de 50F se coloca al aire libre donde la temperatura
es de 100F. Si despus de 4 minutos la temperatura del cuerpo es de 60F, cunto
tiempo transcurrir para que la temperatura del cuerpo sea de 75F? Cul ser
su temperatura despus de 20 minutos?
12. Un cuerpo a una temperatura desconocida se coloca en un cuarto que se mantiene
a una temperatura constante de 30F. Si despus de 10 minutos la temperatura del
cuerpo es de 0F y despus de 20 minutos la temperatura del cuerpo es de 15F,
Cul es la temperatura inicial (desconocida) del cuerpo?
13. Un tanque contiene 100 litros (/) de una solucin que consta de 100 kg de sal
disueltos en agua. Se bombea agua pura hacia el tanque a razn de 5 l/s y la
mezcla, que se mantiene uniforme mediante agitacin, se extrae a la misma razn.
Cunto tiempo pasar antes de que queden solamente 10 kg de sal en el tanque?
14. Un tanque de 500 galones contiene inicialmente 300 galones de solucin salina en
la que se han disuelto 50 libras de sal. Se agrega solucin salina que contiene 3
libras de sal por galn con una rapidez de 4 gal/min. Determine cunta sal hay en
el tanque en el momento que ste se desborda.
15. Un tanque tiene 60 galones de agua pura. Una solucin con 3 Ib de sal por galn
entra a 2 gal/min y la mezcla bien agitada sale a 2.5 gal/min.
a) Halle el nmero de libras de sal que hay en el tanque en cualquier tiempo t.
b) Encuentre la concentracin de sal en el tanque cuando contenga 30 gal de agua
salada.
Ejercicios del Captulo 3 121
16. El lago Ontario tiene un volumen de 1636 km
3
y una concentracin inicial de
contaminantes del 0.25%. Hay un ingreso anual de 209 km
3
de agua con una
concentracin de contaminentes del 0.05% y un derrame anual de igual cantidad,
bien mezclada en el lago. Cunto tiempo pasar para que la concentracin de
contaminantes en el estanque se reduzca al 0.10% ?
17. Suponga que un cuarto contiene 32 m
3
de aire, originalmente libres de monxido de
carbono. En el instante t = 0 se empieza a introducir al cuarto humo de cigarrillo,
con un contenido del 4% de monxido de carbono, con una rapidez de 0.002 m
3
/min
y se deja salir la mezcla bien circulada, con la misma rapidez.
a) Encuentre una expresin para la concentracin x(t) de monxido de carbono
en el cuarto, en cualquier instante t > 0.
b) Para un ser humano, quedar expuesto a una concentracin de monxido de
carbono tan baja como 0.00012 puede ser nocivo. Encuentre el tiempo en el
cual se alcanza esta concentracin.
18. Un tanque de 500 gal contiene inicialmente 100 gal de agua, en la cual se han
disuelto 50 Ib de sal. Comenzando en t = 0, una salmuera cuya concentracin es
de 2 Ib de sal por galn entra al tanque a razn de 5 gal/s. La mezla se mantiene
uniforme mediante agitacin, y estando bien agitada sale del tanque con una rapidez
de 3 gal/s. Qu cantidad de sal contendr el tanque cuando est lleno de salmuera?
19. Un tanque A contiene 100 litros de salmuera, que se obtuvo al disolver 40 kg de sal
en agua. Se introduce en este tanque una salmuera, cuya concentracin es de 3 kg//,
a una rapidez de 2 //min. La mezcla se conserva homognea, sale con la misma
rapidez y va a parar a un segundo tanque B que contiene al principio 100 litros
de salmuera a una concentracin de 0.1 kg//. Agitando se mantiene homognea la
mezcla en el tanque B y sale de ste con una rapidez de 1 //min. Hallar la cantidad
de sal en cada uno de los tanques en cualquier instante.
20. A un circuito en serie, en el cual la inductancia es de 0.1 H y la resistencia es de
50 , se le aplica una tensin de 30 V. Determine la corriente i(t) si z(0) = 0. Cul
ser el valor de la corriente despus de un tiempo largo?
21. A un circuito en serie, en el cual la resistencia es de 200 Q y la capacitancia es de
10~
4
F, se le aplica una tensin de 100 V. Si q(0) = 0, calcule la carga q(t) en el
capacitor y obtenga la corriente i(t).
22. Un inductor de L henrys vara con el tiempo t (en segundos) de acuerdo a L
0.05 + O.OOl. Se conecta en serie con un generador cuya fem es de 40 V y una
resistencia de 10 Q. Si la corriente i es cero inicialmente encuentre i(t) para todo
t > 0. Cul es la corriente mxima terica?
23. Una resistencia de 20 l y un inductor de 5 H se conectan en serie en un circuito
elctrico en el cual hay un flujo de corriente de 20 A al tiempo t = 0. Encuentre la
corriente para t > 0 si la fem es cero para t > 0.
24. Un condensador de 5 x 10~
3
F est en serie con una resistencia de 25 O y una fem
de 50cos6 volts, t > 0. El interruptor se cierra en t = 0. Si la carga inicial en el
condensador es cero, determine la carga y la corriente en cualquier tiempo.
25. Una resistencia de 20 Q se conecta en serie con un condensador de 0.01 F y una
fem en volts dada por 40e~
3
+ 20e~
6
. Si q(0) 0, muestre que la carga mxima
en el condensador es de 0.25 coulombs.
26. Un circuito consiste de una resistencia constante de R ohms en serie con una fem
constante de E volts y una inductancia constante de L henrys. Si la corriente inicial
es cero, muestre que la corriente crece a la mitad de su valor terico mximo en
Ll n2
Captulo 4
Ecuaciones Diferenciales Lineales de
Segundo Orden
4.1 Conceptos Bsicos
Definicin 4.1.1 Se dice que las funciones f y g son linealmente dependientes (Id.) en
el intervalo (a, 6) si existen constantes C\,c
2
no ambas cero tales que
cif(x) + c
2
g{x) = 0
para todo x (a, b).
Definicin 4.1.2 Decimos que las funciones f y g son linealmente independientes (l.i.)
en el intervalo (a, b) si no son Id. , es decir si las nicas constantes para las cuales
cif(x) + c
2
g(x) = 0
para todo x en el intervalo (a, b), son c\ = c
2
= 0.
Es fcil entender las definiciones anteriores, por ejemplo, si / y g son l.d en un intervalo
/ , entonces existen constantes Cx,C2, no ambas nulas tales que
cif(x) + c
2
g(x) = 0
para todo x G / . Sin prdida de generalidad podemos suponer que c\ ^ 0, de modo que
f(x) = -g(x)
Esto es, si dos funciones son l.d. en un intervalo / , entonces una es simplemente un
mltiplo constante de la otra, para todo x en I. Recprocamente si para alguna constante
c
2
se tiene que f(x) = c
2
g(x), x 6 / , entonces
123
124 Captulo 4. Ecuaciones Diferenciales Lineales de Segundo Orden
para todo x en / . Por lo tanto, las funciones / y g son l.d. en / puesto que al menos una
de las constantes es diferente de cero (c\ 1).
Concluimos entonces que dos funciones son l.i. en un intervalo I cuando ninguna es
un mltiplo constante de la otra en I.
EJEMPLO 1. Muestre que las siguientes parejas de funciones son l.d. en R.
a) f(x) = x, g(x) = -2x.
b) f(x) = e\ g{x) = e
x
.
c) f(
x
) = x
2
, g{x) = -V2x
2
.
Solucin. En cada inciso, las parejas de funciones son claramente l.d. en R ya que una
es un mltiplo escalar de la otra. En efecto, tenemos que
a) g(x) = - 2/ ( x) , o equivalentemente 2/(x) + g(x) = 0 para todo x E i
b) g(x) = | / ( x) , \f(x) + (-l)g(x) = 0 para todo xGR.
c) g(x) = - V2/ ( x) , y/2f(x) + g(x) = 0 para todo x G R.
EJEMPLO 2. Compruebe que las funciones f(x) = x y g(x) = x
2
son l.i. en R.
Solucin. Supngase por contradiccin que / y g son l.d. en R, es decir que existen
constantes c\ y c< no simultneamente nulas tales que
c\X + C2X
2
= 0, x G
Si derivamos (4.1) con respecto a x obtenemos
C\ + 2c2X = 0, x G I
(4.1)
(4.2)
Luego, el sistema de ecuaciones lineales (4.1)-(4.2) tiene una solucin no nula {c\ ^.0
c<i ^ 0). De lgebra sabemos que para un sistema homogneo, sto es posible solamente
si el determinante del sistema es cero. As que
x x
1 2x
= 2x
2
- x
2
= x
2
para todo x en R.
Por consiguiente, la suposicin de que las funciones f(x) x y g{x) x
2
son l.d. en
R nos llev a concluir que x
2
0 para todo x en R, lo cual claramente es falso. Esto nos
muestra que f y g son linealmente independientes en R.
Podemos generalizar las ideas expuestas en el ejemplo 2. Primero hacemos la definicin:
Definicin 4.1.3 Sean f y g funciones derivables en el intervalo I. A la funcin definida
por el determinante
f(x) g(x)
f'(x) g\x)
= f{x)g\x) - f(x)g{x),
4.1. Conceptos Bsicos 125
se le llama el wronskiano de f(x) y g{x) y se denota por W(f(x),g(x)) o W(f,g). El
valor de W(f,g) en el punto x se indicar por W(f,g)(x) o simplemente W(x).
Para determinar si dos funciones son l.i. podemos emplear el siguiente teorema.
Teorema 4.1.1 Sean f y g funciones derivables en un intervalo I. Si el wronskiano
W(/? 9)
es
diferente de cero en por lo menos un punto XQ del intervalo I, entonces f y g
son linealmente independientes en I.
Demostracin. Haremos la demostracin por contradiccin. Supngase que W(f, g) (XQ) ^
0 para algn XQ fijo en el intervalo / y que f y g son linealmente dependientes en / . En-
tonces existen constantes c\ y c<i no ambas cero tales que
c
1
f(x) + c
2
g(x) = 0,
para todo x en L Derivando (4.3) resulta
(4.3)
x) = 0. (4.4)
Luego, el sistema de ecuaciones (4.3)-(4.4) tiene una solucin diferente de la trivial para
cada x en / , por lo cual
W(f,g)(x) =
g{x)
Sf{x)
= 0,
para todo x en / . Esto contradice la hiptesis de que W(f,g)(xo) 7^ 0. Por lo tanto se
concluye que f y g son linealmente independientes en / .
EJEMPLO 3. Demuestre que las funciones f(x) = e
x
y g(x) = e~
x
son l.i. en R.
Solucin. Ya que
W(e
x
,e~
x
) =
e
x
e~
x
= -e
x
e~
x
- e
x
e~
x
= - 2 0,
para todo x en R, del Teorema 4.1.1 se concluye que f y g son l.i. en R.
Finalizaremos est seccin recordando algunas nociones muy elementales sobre los
nmeros complejos.
Un nmero complejo z puede representarse en la forma z a + bi donde a, b G R
se denominan parte real y parte imaginaria de z, respectivamente y el smbolo i denota
al nmero imaginario puro que tiene la propiedad de que i
2
= 1.
La igualdad i
2
= 1 nos permitir obtener races de nmeros negativos, puesto que
i == y/^1. Por ejemplo
Para sumar o restar dos nmeros complejos, simplemente se suman o se restan las
partes reales o imaginarias correspondientes. Para multiplicar dos nmeros complejos, se
aplica la propiedad distributiva y el hecho de que i
2
= 1.
La posibilidad de calcular races cuadradas de nmeros negativos permite afirmar
ahora que toda ecuacin cuadrtica, ax
2
+ bx + c = 0 con a, 6, c, nmeros reales, tiene
dos races que pueden ser reales y distintas, reales e iguales o complejas.
Tambin aplicaremos la identidad (Frmula de Euler)
e

con 9 un nmero real.
EJERCICIOS 4.1
Pruebe que los pares de funciones dados son l.i. en R.
1. / (x) = x, g(x) = xe
x
2. / (x) = e
ri
*, g(x) = e
r2
*; n ^ r
2
3. f(x) senx, g(x) = cosa;
4. f(x) = senx, p(x) = sen2x
5. f[x) = senx, ^(x) = x
6. / (x) = senx, g(x) = e
x
7. / (x) = e
x
eos 2x, g'(x) = e
x
sen 2x
8. f(x) = e
x
eos 2x, p(x) = e
2x
eos 2x
9. f(x) = 1,
ff
(x) = x
10. / (x) = x
m
, p(x) = x
n
; m T n
4.2 Solucin de Ecuaciones Diferenciales Lineales de
Segundo Orden
Definicin 4.2.1 Una ecuacin diferencial de la forma
) ^ + a
o
(x)y = f(x) (4.5)
donde ao(x),ai(x),a2(x) y f(x) son funciones de x nicamente, se llama ecuacin dife-
rencial lineal de segundo orden.
4.2. Solucin de la Ecuacin Homognea 127
Para este tipo de ecuaciones, el siguiente teorema garantiza la existencia de una solucin
nica.
Teorema 4.2.1 Sean ao(x),ai(x),d2{x),f(x), funciones continuas en un intervalo I y
supngase que a2(x) ^ 0 para toda x G / . Si Xo es cualquier punto en I, entonces existe
una y sola una solucin y{x) de la ecuacin diferencial
que satisface las condiciones iniciales
y'(x
0
) = 2/i,
donde 2/o>2/i son constantes arbitrarias.
Definicin 4.2.2 La ecuacin diferencial
d v dv
a
*)- +
a
i (
x
) +
a
(
x
)y = (
4
-
6
)
que se obtiene de (4-5) haciendo f(x) 0, se llama ecuacin homognea, reducida
o complementaria.
Si la funcin f(x) no es idnticamente nula, entonces la ecuacin (4.5) recibe el nombre
de ecuacin no homognea .
Para hallar la solucin de la ecuacin (4.5) es necesario resolver la ecuacin homognea
asociada (4.6). Por esta razn veremos primero algunos resultados concernientes a la
solucin de esta ltima.
4.2.1 Solucin de la Ecuacin Homognea
Teorema 4.2.2 Si y\(x) y 2/2(2) son soluciones de la ecuacin diferencial (4-6) entonces
la combinacin lineal y C\y\{x) + 02^/2(2) tambin es solucin de (4-6), donde C\ y C2
son nmeros reales o complejos cualesquiera.
El teorema anterior nos dice que si 2/1,2/2 son soluciones de (4.6), entonces es posible
elaborar un nmero infinito de soluciones de dicha ecuacin, de la forma cij/i(x)+C2i/2(^),
con Ci,C2 constantes. Una pregunta natural sera si esta infinidad de soluciones incluye
a todas las soluciones posibles de (4.6). La respuesta a esta pregunta es si, pero siempre
y cuando y\ y y
2
sean l.i. . Hacemos la siguiente definicin.
Definicin 4.2-3 Si y\ y 2/2 son dos soluciones de la ecuacin (4-6), que adems son
linealmente independientes en un intervalo I, entonces decimos que 2/1 y 2/2 constituyen
un conjunto fundamental de soluciones de (4-6) en I.
Por consiguiente nos interesa saber cuando dos soluciones de la ecuacin diferencial
homognea (4.6) son l.i. en algn intervalo. El siguiente teorema nos da una condicin
necesaria y suficiente para la independencia lineal de soluciones (comprese con el Teo-
rema 4.1.1).
Teorema 4.2.3 Supongamos que ao(x), a\{x), a
2
{x) son funciones continuas en un in-
tervalo I y que a
2
(x) ^ 0 para toda x G / . Sean 2/1,2/2 dos soluciones de la ecuacin (4-6)
en I. Entonces, y'1,2/2
s
n l-i- (forman un conjunto fundamental de soluciones) en I si
y solo si W(yi,y
2
)(XQ) 7^ 0, para algn XQ G / .
A continuacin enunciamos el resultado principal de esta seccin.
Teorema 4.2.4 Supongamos que ao(x), a\(x), a
2
(x) son funciones continuas en un in-
tervalo I y que a
2
(x) ^ 0 para toda x G / . Entonces la ecuacin diferencial homognea
(i6)
x
d
2
y dy
tiene dos soluciones 2/i(x), 2/2(^)7 que son linealmente independientes en I. Adems,
para cualquier otra solucin y = (f>(x) de (4-6) en I se pueden encontrar constantes c\ y
c
2
tales que
) , xel. (4.7)
En vista de (4.7), la solucin general de la ecuacin (4.6) se define como
V = c i 2/ 1(2; ) + c
2
y
2
( x ) , x e l ,
donde 2/1,2/2 son soluciones l.i. en / de (4.6) y Ci,C2 son constantes arbitrarias.
4.2.2 Solucin de la Ecuacin no Homognea
Estamos ya en posicin de explicar cmo determinar la solucin general de la ecuacin
diferencial lineal de segundo orden no homognea (4.5). En primer lugar enunciamos la
siguiente propiedad del conjunto de sus soluciones.
Teorema 4.2.5 La diferencia de dos soluciones cualesquiera de la ecuacin diferencial
no homognea (4-5) es una solucin de la correspondiente ecuacin homognea (4-6).
Como una consecuencia inmediata de este teorema, a continuacin damos el resultado
principal, que nos dice cmo determinar la solucin general de (4.5).
Teorema 4.2.6 Dada una solucin y
p
(x) de la ecuacin diferencial (4-5) en un intervalo
I, entonces para cualquier solucin y = (f>(x) de esta ecuacin, existen constantes C\,c
2
tales que
<t>{x) = ciyi(x) + c
2
y2(x) + y
p
(x), xel, (4.8)
donde 2/1 y 2/2 son soluciones l.i. en I de la ecuacin homognea correspondiente (4-6).
4.3. Mtodo de Reduccin de Orden 129
De acuerdo con (4.8), la solucin general de la ecuacin diferencial lineal no ho-
mognea (4.5) se define como
y(x) = CI T/ I ( X) + c
2
y
2
(x) + y
p
(x), xel, (4.9)
Llamaremos a y
p
una solucin particular de (4.5). A la solucin general de la ecuacin
homognea (4.6) se le denomina solucin complementaria y la denotaremos por y
c
(x).
Es decir
y
c
(x) = c
1
y
1
(x) + c
2
y2(x).
En consecuencia, podemos escribir la solucin general (4.9) de la ecuacin no homognea
(4.5) en la forma
y(x)=y
c
(x)+y
p
(x). (4.10)
Posponemos hasta la seccin 4.6 la discucin de cmo determinar y
p
.
4.3 Mtodo de Reduccin de Orden
Dada una solucin y\(x) de la ecuacin diferencial de segundo orden
y
ff
+p{x)y
/
+ q(x)y = 0, (4.11)
puede determinarse una segunda solucin y
2
(x) que sea linealmente independiente con
2/i(x), de la forma v(x)y\(x), para cierta funcin v(x) distinta de una constante.
Sea y(x) yi(x)v(x), entonces
y' =
y" =
Sustituyendo las expresiones anteriores para y, y' y y" en (4.11) y simplificando resulta
{y
t
v
ff
+ 2y'
1
v' + y
/
;)+p(x)(y
1
v
/
+ y[v) + q{x)v = 0,
v(y" +P(x)y
/
1
+ q(x)y
l
)+v
f
(2y
/
l
+p(x)y
1
)+y
l
v
/f
= 0.
Y como y\ es una solucin de (4.11), el primer trmino en el lado izquierdo de la igualdad
anterior es igual a cero. As que
yiv" + (2y[ + p(x)
yi
)v' = 0
o bien
v" + (p(x) + ^ - V = 0. (4.12)
V vi J
Luego, para que la funcin yi(x)v(x) sea una solucin de la ecuacin diferencial (4.11),
v(x) debe satisfacer la ecuacin diferencial de segundo orden (4.12). Ntese que haciendo
la sustitucin u(x) = v'{x) entonces u'{x) = v"{x) y (4.12) se reduce a la ecuacin
u'+ (p(x) + ^ ) u = 0, ( 4. 13)
V Vi
la cual es ahora de primer orden para la funcin incgnita u. Es por esta razn que
al mtodo que estamos desarrollando para calcular y
2
se le conoce como Mtodo de
Reduccin de Orden.
La ecuacin (4.13) es lineal en u y tambin de variables separables. Separando varia-
bles tenemos
du ( ,
x
2y[\
e integrando y simplificando, obtenemos
ln u p(x)dx 2 ln yi + lii c
\nu lnq/f
2
/ p(x)dx,
donde c es una constante arbitraria. Aplicando exponencial a ambos lados de la ltima
igualdad encontramos que
Por consiguiente
v(x) = c - ^ e
J y{{x)
Tomando c = 1 tenemos el siguiente resultado.
Teorema 4.3.1 Si yi(x) es solucin de la ecuacin diferencial
y"(x) + p(x)y'(x) + q{x)y(x) = 0, (4.14)
entonces una segunda solucin 2/2(2) de (4-14).linealmente independiente con y\{x) es
r
e
~ Jp(
x
)
dx
=
Vl
(x) j dx. (4.15)
EJEMPLO 1. Dado que y\(x) = x~
2
es solucin de la ecuacin diferencial
x
2
y" - 7xy' - 20y = 0, (4.16)
encuentre su solucin general en el intervalo (0,00).
Solucin. Verifiquemos que y\(x) es solucin de la ecuacin diferencial (4.16). Tenemos
que
y[(x) = - 2z"
3
, y'i{x) - 6x~
4
.
4.3. Mtodo de Reduccin de Orden 131
x
2
6x"
4
- 7x(-2x~
3
) - 20z"
2
= 0
6x~
2
+ 14x~
2
- 20z~
2
- 0
0 = 0.
As, efectivamente yi es una solucin de (4.16).
Ahora utilizaremos el resultado (4.15) del teorema anterior para determinar una se-
gunda solucin de la ecuacin diferencial, l.i. con y\.
Primero, reescribimos (4.16) en la forma
_!,_,-,
7
de aqu que en este caso p(x) = y entonces
x
-J
P
(x)dx -j-i
d
,
e
7 1 n
/ 2( \
dx =
/ ~7Tr
dx =
/
J y{{x) J (x~
2
)
2
J x~
- dx - x
ll
dx -
~ J x~* J 12*
Por lo tanto
T
i2
r
i o
Note que una segunda solucin l.i. con y\{x) es simplemente 2/2(2) = x
10
. De modo
que la solucin general en (0,00) de la ecuacin diferencial (4.16) es
y{x) = Cix~
2
+ c
2
x
10
.
EJEMPLO 2. Encuentre la solucin general en (0,00) de la ecuacin diferencial
x
2
y" + xy
f
+ y = 0, (4.17)
si yi(x) = coslnx, es una solucin de la ecuacin.
Solucin. Nuevamente emplearemos (4.15) para obtener una segunda solucin y
2
de
(4.17). En este caso p(x) = - , por lo cual
x
x)dx r
e
-
l
*
d
*
dx
I 0/ \
dx =
/
J y\(x) J
cos
2
mx
1
= / ~dx = tanlnx.
eos
2
ln x
En consecuencia
y
2
{x) (coslnx)(tanlnx) = sen lnx.
De donde la solucin general en (0, oc) de (4.17) es
y(x) c\ eos ln x + C2 sen ln x.
EJERCI CI OS 4.3
Verifique si la funcin y\ indicada es una solucin de la ecuacin diferencial dada. En
caso de serlo determine la solucin general de la ecuacin.
1. i / " - 9 / = 0, ?/i =e
3
*
2. y" + 9y = 0, Vi = eos 3x
3. x
2
y" - 2xy
f
+ 2y = 0, y
x
= x
4. x
2
y" + 2xy
f
- 6y = 0,
Vl
= x
2
5. x
3
y" + x
2
y' + xy = 0, y\ = sen (lnx)
6. xV
7
- 4xy
7
+ 6y = 0, y
x
= x
2
7. x V - 4:n/ + 4y = 0, ?/i = x
8. (2x + 1 ) ^ - 4(x + I)?/
7
+ 4j/ = 0, yi = x + 1
9- / + -y
7
22/ = 0,
10. x
2
y
;/
+ 3x2/ = 0, i/! - 1
11. x
2
y" + xy
1
4y 0, yi = ^
2
12. ( 1 - X V - 2 X 2 /
/
+ 22/ = O,
Vl
=x
13. x
2
/ + xy
7
+ (x
2
- -} y = 0, y
x
= x"
1/ 2
eos x
14. x V + xy
;
= 0,
yi
= 1
15. x
2
y" xy
r
+ y = 0, yi = x ln x
16. (4cotx)y" + (4 senx)y
r
- y = 0, yi = sonx
4.3. Ecuaciones Diferenciales de Segundo Orden con Coecientes Constantes 133
4.4 Ecuaciones Diferenciales Lineales de Segundo Or-
den Homogneas con Coeficientes Constantes
En esta seccin estudiaremos la ecuacin diferencial de la forma
ay" + by' + cy = 0, (4-18)
donde a, 6, c G R y a ^ 0.
Para este tipo de ecuaciones proponemos una solucin de la forma
y(x) = e
rx
.
Entonces
y\x) = re
r x
, y"(x) = r
2
e
rx
y sustituyendo en la ecuacin (4.18) resulta
ar
2
e
rx
+ bre
rx
+ ce
rx
= 0
e
rx
(ar
2
+ br + c) = 0.
Luego, si r es una raz de la ecuacin
ar
2
+ br + c = 0, (4-19)
llamada ecuacin auxiliar o ecuacin caracterstica, la funcin y = e
rx
es una
solucin de (4.18). Debemos considerar tres casos, segn sean las races caractersticas:
reales y distintas, reales e iguales o complejas.
CASO 1. Si ri y r-2 son races reales y distintas entonces y\{x) e
TlX
y 2/2(2) = e
r2X
son
dos soluciones linealmente independientes de la ecuacin diferencial (4.18), de donde su
solucin general es
y(x) = c
x
e
rix
+ c
2
e
r2X
.
CASO 2. Si las races son reales e iguales entonces r\ r
2
= = r. As que, una
2a
solucin de (4.18) es y\(x) = e
rx
.
Podemos encontrar una segunda solucin y
2
linealmente independiente con y\ emple-
n-iido la frmula (4.15), del mtodo de reduccin de orden estudiado en la seccin anterior.
Tenemos
2/2(2) - y
x
(x)
-dx
O 1(
o
zt
xe~
b
2a
^
e~
i*
Luego, para obtener una segunda solucin linealmente independiente con j/i(x), basta
con multiplicar yi{x) por x.
La solucin general de (4.18) en este caso es
y(x) = cie
rx
+ c
2
xe
rx
.
o bien
y(z) = (c
1
+ c
2
x)e
rx
.
CASO 3. Supongamos finalmente que las races son complejas y denotmoslas por
r\ a + ip; r
2
= a i/3.
Entonces dos soluciones l.i. de la ecuacin diferencial son
yi{x) =
a
n*
9
y
2
(x) = *-W*.
Sin embargo, estamos interesados en encontrar soluciones con valores reales. Tenemos
que y(x) = kie^
a+i
^
x
+ k
2
e(
a
~~
i3>)x
, es solucin de la ecuacin diferencial, para cualesquier
constantes ki y k
2
. Pero, de la frmula de Euler se sigue que
y(x) = k
x
e
ax
^
x
+ k
2
e
ax
~
il3x
e
ax
[ki(cos3x + sen/3x) + k
2
{cosf3x i sen/?x)]
= e
ax
[(ki + k
2
) cos/fa + (ki - k
2
)i sen/fo].
Ahora bien, si en particular tomamos k\ = k
2
.= - , obtenemos que la funcin
yi(x) = e
ai;
(cos fix + 0) = e
ax
cosâ;
es una solucin de la ecuacin de variable real con valores reales.
Anlogamente, si fci = -; k
2
= - resulta la funcin
y2
(x) = e
ax
[o + (-
%
- - ^\ i sen(3x\ = e
ax
sen(3x,
que tambin es una solucin de la ecuacin que toma valores reales. Por lo tanto, dos
soluciones de la ecuacin diferencial linealmente independientes con valores reales son
() fa; 2/2(2) e
ax
sen(3x.
La solucin general de (4.18) ahora es
y(x) e
ax
(ci cos/3x + C2sen/3x).
EJEMPLO 1. Resolver
y" - y' - 6y = 0.
Solucin. La ecuacin tiene soluciones de la forma y = e
rx
donde r es solucin de la
ecuacin auxiliar
r
2
-
r
- 6 = 0
( r - 3) ( r + 2) = 0,
de modo que las races caractersticas son r\ = 3, r
2
= - 2. Por lo tanto, dos soluciones
linealmente independientes de la ecuacin diferencial son y\{x) = e
3x
y 2/2(2) = e~
2x
. De
donde la solucin general es
y(x) =
3 2
EJEMPLO 2. Resolver
y" - oy = 0.
Solucin. La ecuacin auxiliar es
r
2
- 5 = 0.
Las races caractersticas son r\ \ / 5, r
2
= y/E. de donde, dos soluciones l.i. de
la ecuacin son
yi(x) = c ^ , y
2
(x) = e-^
x
.
Por lo tanto, la solucin general de la ecuacin diferencial es
EJEMPLO 3. Resolver
y" + Yly
1
+ 36?/ = 0.
36 = 0
(r + 6)
2
= 0
y las races caractersticas son
n = r
2
= - 6 .
Por lo tanto, dos soluciones linealmente independientes son
De donde la solucin general es
y(x) = cie~
6x
+ c
2
xe~
6x
= (c
x
+ c
2
z)e~
6x
.
EJEMPLO 4. Resolver
y" - 16j/ + 64y = 0.
r
2
- 16r + 64 = 0.
Sus races son r\ = r
2
= 8. Por lo tanto, dos soluciones linealmente independientes son
yi
(x) = e
8x
, y
2
(x) = xe
8x
.
As, la solucin general de la ecuacin diferencial es
y(x) = (ci + c
2
x)e
8x
.
EJEMPLO 5. Resolver
y" + 2y' + 17y = 0.
r
2
+ 2r + 17 = 0.
Las races estn dadas por
- 2
v
/
4 ^ 6 8 - 2 ^ 6 4 - 28 , , , .
n,
2
= -
2
= = -
T
- = -i4z
o bien
r
1
= 1 + 4
7
r2 == 1 4.
Por lo tanto, dos soluciones linealmente independientes son
yi(x) = e~
x
cos4x, V2{x)
de donde, la solucin general es
y(x) e~
x
{c\ cos4x + C2
EJEMPLO 6. Resolver
y" + y
f
+ y = 0.
r
2
+ r + 1 = 0,
cuyas races son
1 \ / 3 . 1
2 2 '
2
2
De donde, dos soluciones linealmente independientes son
1
2 2 '
2
2 2
/ \ - x V^ , v _i
x
y / %
yi{x) = e
2
eosx, y
2
{x) = e * s enx.
Por lo tanto, la solucin general es
y{x) = e
2
I ci eos x + c
2
sen x I .
y 2 2 y
EJEMPLO 7. Resuelva la ecuacin diferencial
y" + y = o,
sujeta a las condiciones iniciales
y(0) = l, i/(0) = - 1 .
Solucin. La ecuacin caracterstica es r
2
+1 y sus races son r = . Luego, la solucin
general de la ecuacin diferencial es
y(x) = Ci cosx +
Ahora, la primera condicin inicial implica que
1 =y(0) = ci, ci = 1.
Por otra parte, tenemos que
y'(x) = Cisenx + C2COSX
y usando la segunda condicin inicial
_l =y
/
( 0) = c
2
, c
2
= - l .
Por lo tanto, la solucin que satisface las condiciones iniciales es
y(x) = cosx senx.
EJEMPLO 8. Resolver el problema de valores iniciales
y"-y'-2y = 0; 2/(0) = 1, j/(0) = 4.
Solucin. La ecuacin caracterstica es
r
2
- r - 2 = 0
( r - 2) ( r + l) = 0,
de modo que las races caractersticas son n = 2 y r
2
= 1. La solucin general es
entonces
y(x) = c
x
e
2x
+ c
2
e~
x
.
De las condiciones iniciales, se sigue que
1 = y(0) Ce
0
+ C2e = c\ + C2
4 = y
;
(o) = 2c\e C2e = 2ci C2.
Resolviendo el sistema de ecuaciones
C\ + C2 = 1,
2ci - c
2
= 4,
obtenemos
5 2
ci = - , c
2
= r.
Por lo tanto, la solucin que satisface las condiciones iniciales es
7/ T* I P f*
3 3
4.5. Ecuaciones Diferenciales Lineales de Orden n 139
EJERCI CI OS 4.4
Resuelva las siguientes ecuaciones diferenciales.
1. y" + 5y = 0
2. y" + 9y'-10y = 0
3. y" - Ay' + 12y = 0
4. 4y" + ly' - 2y = 0
5. y" + 6y' + 9y = 0
6. 2y" - Uy' + 3y = 0
7. y" - 18y' + 81y = 0
8. Ay" + y' + 2y = 0
9. y" - 6y' + 5y = 0 con y(0) = 3, j/(0) = 11
10. y" + Ay' + 5y = 0 con y(0) = 1, y'(0) = 0
11. y" + 4y = 0 con y(0) = 1, y'(0) - - 2
12. y" - 2y' + y = 0 con y(0) = - 1 , y'(0) = 3
4.5 Ecuaciones Diferenciales Lineales de Orden n
Presentaremos ahora la generalizacin de los resultados enunciados en la secciones 4.1 y
4.2.
Definicin 4.5.1 Se dice que un conjunto de funciones fi(x), f
2
(x),..., /
n
(x) es lineal-
mente dependiente (l.d.) en un intervalo I si existen constantes c\, c
2
,.. , c
n
no todas
cero tales que
c\f\{x) + c
2
/
2
(z) + + Cnf
n
(x) = 0,
para todo x en I. Mientras que fi(x), ^( a: ) , . . . , f
n
(x) son linealmente independien-
tes (l.i.) en I si no son l.d. en I, es decir, la igualdad
Clfl(x) + C
2
/
2
(z) + + Cnf
n
(x) = 0,
para todo x G / implica que C\ = 0, c
2
= 0, . . . , c
n
= 0.
Definicin 4.5.2 Sean / i , / 2, . ., f
n
, funciones que tienen al menos n 1 derivadas en
un intervalo abierto I. Para x en / , el wronskiano de dichas funciones se define como
el determinante
/
/ \ r / \ r / \
l\X) J2\
X
) ' ' ' Jn\X)
Definicin 4.5.3 Se dice que una ecuacin diferencial de orden n es lineal si tiene la
forma
) + a
o
(x)y(x) = g(x), (4.20)
donde las funciones a
o
(x), ai ( x) , . . . , a
n
(x) y g(x) dependen solamente de la variable x.
Como ant es , l a ecuaci n (4. 20) es no homognea si g(x) ^ 0 y haci endo g(x) = 0
obt enemos l a ecuaci n homognea r educi da o compl ement ar i a
an(x)yW(x) + On-i^y*"-
1
^) + ... + a
x
{x)\/{x) + a
o
(x)y(x) = 0. (4.21)
Teorema 4.5.1 Sean yx^yi, , y
n
un conjunto fundamental de soluciones de la ecuacin
homognea (4-21) en un intervalo I, esto es,y\,y2,...,y
n
son li. en I o equivalentemente
W(yi,y2,.. ,/n)(ô) T 0 P
ara
algn x
0
6 / . Entonces la solucin general de (4-21) est
dada por
y(x) = ciyi(x) + c
2
y
2
(z) + . . . + c
n
y
n
(x), x e J,
donde ci, C2,..., c
n
son constantes arbitrarias.
Teorema 4.5.2 Sea y
p
una solucin dada de la ecuacin diferencial no homognea (4-20)
en el intervalo I y sea
y
c
(x) = cij/i(x) + c
2
y2(x) + . . . + c
n
y
n
(x)
la solucin general de la ecuacin homognea asociada (4-21) en el intervalo I . Entonces
la solucin general de (4-20) es
y(x) = CI/I(X) + c
2
y2x) + . . . + Cny
n
(x) + y
p
(x) = y
c
(x) + y
p
(x).
4.5.1 Ecuaciones Diferenciales Lineales de Orden n Homogneas
con Coeficientes Constantes
En la seccin 4.4 estudiamos como resolver una ecuacin diferencial lineal homognea de
segundo orden con coeficientes constantes. Ahora, de manera ms general, considerare-
mos la ecuacin diferencial de orden n
)
+ a
n
.
iy
(n
-
l)
+ - + a
2
y
n
+ a
lV
' + a
o
y = 0, (4.22)
donde ao, a i , . . . , a
n
son nmeros reales.
Es fcil ver que una funcin exponencial de la forma
y = e
rx
,
es solucin de (4.22) si y slo si r es una raz de la ecuacin
a
n
r
n
+ dn^r
71
-
1
+ + a
2
r
2
+ a
x
r + a
0
= 0. (4.23)
A (4.23) se le llama ecuacin auxiliar o ecuacin caracterstica de la ecuacin
diferencial (4.22).
No podemos hacer un anlisis general de las races de la ecuacin auxiliar (4.23)
semejante al que hicimos para las ecuaciones de orden 2, ya que pueden aparecer muchas
combinaciones si n es mayor que dos. Por ejemplo una ecuacin de quinto grado puede
tener cinco races reales diferentes, tres races reales diferentes y dos races complejas,
una raz real y cuatro complejas, cinco races reales e iguales, cinco races reales con tres
de ellas iguales, y as sucesivamente. En su lugar mencionaremos los siguientes tres casos.
CASO 1. Si todas las races de (4.23), r i , r 2, . . . , r
n
, son reales y distintas entonces la
solucin general de (4.22) es
CASO 2. Si r\ es una raz de multiplicidad k de (4.23), es decir k races son iguales a
r
l 5
entonces correspondiendo a esta raz se tienen las siguientes k soluciones l.i. de (4.22)
y la solucin general de (4.22) debe contener la combinacin lineal
de
rix
+ c
2
xe
rix
+ c
3
x
2
e
rix
+ + c
k
^x
k
~
l
e
riX
.
CASO 3. Cuando T\ a + i(5 es una raz compleja de multiplicidad k de (4.23), su
conjugado r
2
= a ifi es tambin raz de multiplicidad k. En este caso, la solucin general
de la ecuacin diferencial (4.22) debe contener una combinacin lineal de las siguientes
2k soluciones l.i.
e
ax
eos j3x, xe
ax
eos (3x, x
2
e
ax
eos (3x,..., x
k
~
l
e
ax
eos /3x,
e
ax
sen (3x, xe
ax
sen (3x, x
2
e
ax
sen /3x,..., x
k
~
1
e
ax
sen (3x.
EJEMPLO 1. Resolver
y(
4
>-
1
/ '
/
- 7i /
/ /
+
I
/ + 6y = 0. (4.24)
Solucin. La ecuacin auxiliar de (4.24) es
r
4_
r
3_
?r
2
+ r + 6 = 0
( r - l ) ( r + l)(r + 2) ( r - 3) = 0,
cuyas races son
ri = l, r
2
= - l , r
3
= - 2, r
4
= 3.
Luego, cuatro soluciones linealmente independientes de la ecuacin son
yi(x) = e
x
, y
2
(x) = e~
x
, y
3
(a;) - e"
2x
, y
4
(a:) = e
3x
y la solucin general de (4.24) es
y(x) = cie
x
+ c
2
e~
x
+ c
3
e~
2x
+ c
4
e
3a
\
EJEMPLO 2. Resolver
2,(6) _
8/
(5)
+ 17?/
(4)
+ 6y
/// _
44j/
//
+
gy
+ 32y = 0
(425)
Solucin. La ecuacin auxiliar de (4.25) es
r
6
- 8r
5
+ 17r
4
+ 6r
3
- 44r
2
+ 8r + 32 = 0
( r - 2)
3
( r - 4) ( r + l)
2
= 0.
De la ltima expresin es claro que las races de la ecuacin auxiliar, con sus respectivas
multiplicidades son
T\ = 2, de multiplicidad tres;
?2 = 4, de multiplidad uno (raz simple);
r
3
= 1, de multiplicidad dos (raz doble).
Tomando en cuenta cada raz y sus multiplicidades respectivas, tenemos las siguientes
soluciones l.i. de (4.25)
e
2x
, xe
2x
, x
2
e
2x
correspondientes a la raz r\ = 2,
e
4x
correspondiente a la raz r<i 4,
e~
x
,xe~
x
correspondientes a la raz r = 1.
Por consiguiente la solucin general de (4.25) es
y = (ci + c
2
x + c
3
x
2
)e
2x
.+ c
4
e
4x
+ (c
5
+ c
6
x)e~
x
.
EJEMPLO 3. Resolver
y(*)
+
$y +
16/ =
0. (4.26)
r
4
+ 8r
2
+ 16 = 0
(r
2
+ 4)
2
=
o.
Entonces, las races caractersticas son en este caso r\ 2% y r<2 = 2i de multiplicidad
dos, cada una de ellas. Correspondientemente tenemos las cuatro soluciones l.i.
2/i(x) = cos2x, 1/2(^) sen2x, y$(x) xcos2x, ?A(X) = x sen2x.
Por lo tanto la solucin general de (4.25) es
y(x) = C\ eos 2x + C2sen 2x + x(cs eos 2x + C4sen 2x).
EJEMPLO 4. Resolver
y(*)-2y'" + 2y"-2y
f
+ y = 0. (4.27)
3
+ 2r
2
- 2r + l = 0
( r - l )
2
( r
2
+ l) = 0,
cuyas races son
ri = 1, de multiplicidad dos;
r2 = , de multiplicidad uno;
r3 = , de multiplicidad uno.
De donde, la solucin general de (4.27) es
y(x) = C\e
x
+ c
2
xe
x
+ c
3
sen x + c
4
eos x.
EJEMPLO 5. Resolver
y
(8) _
1Qy
(7)
+ 46l /
(6) _
88r
4
+ 146r
3
- 350r
2
+ 98r + 343 - 0 = 0
(r + 1)
2
(r
2
- 4r + 7)
3
- 0
y sus races son
ri = 1, de multiplicidad dos;
r
2
= 2 + \/3, de multiplicidad tres;
r3 = 2 \/3, de multiplicidad tres.
De modo que, la solucin general de (4.28) es
y(x) = (ci 4- c
2
x)e~
x
+ e
2a:
(c3 eos VSx + c
4
sen \/3x)+
xe
2x
(c5 eos v3x + cesen v3x) + x
2
e
2x
(c7 eos A/3X + cssen y/3x).
EJEMPLO 6. Cul es la solucin general de una ecuacin diferencial cuya ecuacin
auxiliar tiene races: 2, 1,0, 0,3 5, 2,0,, 3 5z?
Solucin. La ecuacin auxiliar es de grado 10 y sus races, con sus respectivas multipli-
cidades son
0, de multiplicidad tres;
1, de multiplicidad uno;
2, de multiplicidad dos;
3 + 5, de multiplicidad dos;
3 5z, de multiplicidad dos.
Por lo tanto la solucin general tiene la forma
y = ci + c
2
x + c
3
x
2
+ ce~
x
+ (c
5
+ c
6
x)e
2:c
+
xe
3x
(cgsen5x + CIQCOS5X).
EJERCICIOS 4.5
i. y"
f
-y = o
2. yW - 162/ = 0
3. y
/
" + 3y" + 3
2
/ +
2
/ = O
4.
y
W + 1 3 / + 36y = 0
4.6. Mtodo de Coecientes Indeterminados: Superposicin 145
5. y& + 32y" + 256y = O
6. y'" - 3y" - 10y' = O
7. y& + 8y'" + 16y' = O
8. y'" + y" + y'-3y = 0
9. j/
5
> - y^ + 2y'" + 22y" - 35y' + 7by = O
10. y'" + 3y" + 2y' = 0 con y(0) = 1; y'(0) = 2; y"(0) = - 1
11. y"' - 6y" + 12y' - 8y = 0 con j/(0) = - 2 ; y'(0) - 1; y"(0) = 2
12. ?/
4
) + 5y" + Ay = 0 con y(0) = 3; y'(0) = 0; y"(0) = 0; y'"(0) = - 3
4.6 Mtodo de Coeficientes Indeterminados: Enfoque
de Superposicin
Este mtodo nos permite encontrar una solucin particular y
p
(x) para las ecuaciones
diferenciales lineales de segundo orden de la forma
i_ h~ -L cu CI(T\ (A 2Q^1
donde a, 6, c son constantes y
9\
x
) '
b
k
x
k
+ b
k
_
x
x
k l
+
e
ax
senbx
cosbx
f b\x + be funcin
funcin
funcin
funcin
polinomial
exponencial
seno
coseno
El mtodo es aplicable tambin cuando la funcin g{x) en (4.29) consiste de una suma
y productos finitos de funciones polinomiales, exponenciales, seno y coseno. Asimismo,
pueden considerarse ecuaciones diferenciales no homogneas con coeficientes constantes
de orden superior.
El enfoque del mtodo de coeficientes indeterminados que presentamos en esta seccin
se basa esencialmente en tres principios u observaciones que la prctica de derivacin de
funciones nos ha enseado.
1. Cuando derivamos un polinomio, el grado de ste disminuye en uno.
Sip(x) = b
k
x
k
+b
k
_
1
x
k
~
1
-\ \-biX+b
Q
entonces g'{x) = kb
k
x
k
'
1
+ (k-l)b
k
îx
k
~~
2
' +
- - - + b\. Evidentemente si derivamos dos veces p, su grado disminuye en dos.
2. Al derivar una funcin exponencial, la funcin "casi no cambia". Si g(x) = e
ax
entonces g'(x) ae
ax
ag(x). La derivada es casi la funcin g (salvo por la
constante multiplicativa a).
3- Si derivamos g{x) = senmx pasamos al coseno: g'{x) = mcosmx.
Si derivamos g{x) = cosmx pasamos al seno: g'{x) = msenmx.
Si derivamos dos veces g{x) = senmx regresamos casi ag(x), g"(x)
=
~~
m2
senmx.
Si derivamos dos veces g(x) = cosmx regresamos casi a g(x), g"{x) = m
2
cosmx.
Luego, es razonable pensar que una solucin particular de (4.29) tendr la misma
forma que #(x), excepto cuando g es una solucin de la ecuacin homognea.
En esencia, el mtodo consiste en proponer una solucin particular de (4.29) que
contenga uno o ms coeficientes desconocidos. Entonces sustituimos esta solucin pro-
puesta en la ecuacin diferencial y escogemos los coeficientes de tal manera que la funcin
efectivamente satisfaga la ecuacin.
A continuacin discutiremos algunos casos para hallar una solucin particular de
(4.29), dependiendo de la forma de g{x).
CASO I. g(x) = P
n
(x) = a
n
x
n
+ aîx
71
"
1
+ + a
x
x + a
0
.
En este caso la ecuacin diferencial (4.29) toma la forma
d xi dxi
a\ + bj- + cy = a
n
x
n
+ aîx
71
"
1
+ + a
x
x + a
0
. (4.30)
Proponemos una solucin particular de la forma
y
p
(x) = A
n
x
n
+ A
n
-ix
n
~
l
+ + A
2
x
2
+ A
x
x + A
o
.
Sustituyendo y
p
, y'
p
y y^ en (4.30) resulta
a[n{n - l)A
n
x
n
"
2
+ + 2A
2
] + b[nA
n
x
n
-
1
+ + A
x
}+
c(A
n
x
n
+ An-xx
71
"
1
+ + Ao) = a
n
x
n
+ a
n
-
X
x
n
~
x
+ + ao,
o equivalentemente
cA
n
x
n
+(cA
n
^
1
+nbA
n
)x
n
^
l
+' +(cA
0
+bA
l
+2aA
2
) = a
n
x
n
+a
n
_i x
n
~
1
+- -+a
0
, (4.31)
y comparando coeficientes obtenemos el sistema de ecuaciones
cA
n
= a
n
nbA
n
+ cA
n
_i = a
n
_i
2a A
2
+ bAi + cA
0
= a
Q
.
Si c ^ O de la primera ecuacin determinamos A
n
y de las restantes los dems coefi-
cientes.
Si c 0 pero b ^ 0, el polinomio en el miembro izquierdo de (4.31) es de grado n 1
y dicha ecuacin no puede satisfacerse. As que si c 0 proponemos
y
p
(x) = x(A
n
x
n
+ A
n
-ix
n
~
l
H h A
2
x
2
+ A
x
x + A)),
y procedemos como antes para determinar A
n
, A
n
_i , . . . , AQ. Ntese adems que si c = 0
una constante es solucin de la ecuacin diferencial homognea.
Si tanto b = 0 como c = 0 (1 y x son soluciones de la homognea), se propone
y
p
(x) = x
2
(A
n
x
n
+ An^x
71
-
1
+ > + A
2
x
2
+ A
x
x + A
o
) ,
aunque ahora la ecuacin diferencial puede integrarse directamente.
CASO II. g(x) = e
ax
P
n
(x), donde P
n
(x) es un polinomio de grado n.
Tenemos ahora la euacin
P
^ y = e
a
*P
n
(x). (4.32) b^+cy eP
n
(x).
dx
Son posibles los siguientes subcasos.
a) a no es una raz de la ecuacin auxiliar ar
2
+ br + c = 0.
En este caso, es preciso hallar una solucin particular de la forma
y
p
(x) = (A
n
x
n
+ An-xx"'
1
.*.. - + A
0
)e
ax
- Q
n
(x)e
ax
.
En efecto, introduciendo y
p
, y'
v
y y^ en (4.32) y dividiendo por e
ax
se sigue que
aQ'^x) + (2aa + b)Q'
n
{x) + (aa
2
+ ba + c)Q
n
(x) = F
n
(x). (4.33)
Ya que grad (Q
n
(z)) = n, grad(Q^(x)) = n - l y grad(Q^(x)) = n - 2, los
polinomios en ambos miembros de (4.33) son de grado n. Igualando los coeficientes
de las mismas potencias de x se obtiene un sistema de n+1 ecuaciones que determina
los valores de A
n
, A
n
_ i , . . . , A
o
.
b) a es una raz simple de la ecuacin auxiliar.
En este caso aa
2
+ ba + c = 0, de modo que en lado izquierdo de (4.33) se tiene un
polinomio de grado n - l y dicha igualdad no puede satisfacerse sin importar cuales
sean los valores de A
n
, A
n
_i , . . . , A
o
. Debido a esto, ahora buscamos una solucin
particular en la forma de un polinomio de grado n + 1 sin trmino independiente
(pues ste se anula durante la derivacin) por e
ax
. As hacemos
y
p
{x) = xQ
n
(x)e
ax
.
c) a es una raz doble de la ecuacin auxiliar.
Entonces ao
2
+ bo + c = 0, y 2aa + b = 0 (a b/2a es la raz doble), de manera
que el lado izquierdo de (4.33) se reduce a aQ'^x) y para satisfacer la igualdad
se requiere buscar una solucin particular de la forma de producto de e
ax
por un
polinomio de grado n + 2. Los trminos independiente y lineal se anulan al derivar
dos veces, por lo que pueden omitirse en la forma de y
p
. Por consiguiente en este
caso proponemos
y
p
(x) = x
2
Q
n
(x)e
ax
.
CASO III. g{x) = P(x)e
ax
eos Px + Q(x)e
ax
sen px, donde P(x) y Q(x) son polinomios.
Podemos examinar este caso en forma anloga al caso II, usando que
e
i(5x
+ e
-i(3x
e
i(3x _
e
-if3x
eos px , sen px
:
,
2 2z
por lo cual
P{X) + ~Q(x)
Y considerando de manera independiente las partes real e imaginaria, podemos hallar
soluciones que no contengan nmeros complejos de la siguiente forma:
a) Si a + i(3 no es raz de la ecuacin auxiliar, buscamos una solucin particular de
la forma
y
p
(x) = u(x)e
ax
eos (5x + v(x)e
ax
sen Px, (4.34)
donde u(x) y v{x) son polinomios cuyo grado es igual al mayor de los grados de
P(x) y Q(x).
b) Si a + ifi es raz de la ecuacin auxiliar, hacemos
y
v
[x) = x[u(x)e
ax
eos (3x + v(x)e
ax
senfix], (4.35)
con u(x) y v(x) como antes.
Un caso particular es cuando g(x) tiene la forma
g(x) = a eos fix + b sen px, a, 6 e R,
es decir P{x) y Q{x) son de grado cero y a = 0. Entonces los resultados anteriores se
reducen a los siguientes:
a.l) Si fii no es una raz de la ecuacin auxiliar, buscamos una solucin de la forma
y
p
(x) A eos fix + B sen 3x
b.l) Si (5% es una raz de la ecuacin auxiliar, proponemos
y
p
(x) x{Acos3x + Bsenfix).
Finalmente enfatizamos que las formas propuestas (4.34) y (4.35), para la solucin
particular, tambin son vlidas cuando P(x) = 0 o Q{x) 0 y en el caso particular
cuando a 0 o b = 0.
EJEMPLO 1. Resolver
y
11
+ 3T/' + 2y = 3x
2
- x + 1. (4.36)
Solucin. De acuerdo con (4.10), la solucin general de (4.36) tiene la forma y = y
c
+ y
P
,
donde y
c
es la solucin general de la ecuacin homognea
y" + 3i/ + 2y = 0, (4.37)
y y
p
es una solucin particular de (4.36). La ecuacin auxiliar de (4.37) es r
2
+ 3r + 2 = 0,
cuyas races son n = I y r 2 = 2. Luego
y
c
(x) = c
x
e~
x
+ c
2
e~
2x
.
Por otra parte, proponemos una solucin particular de la forma
y
p
(x) = Ax
2
+ Bx + C,
ya que el lado derecho de (4.36) es un polinomio de grado 2 y 0 no es raz caracterstica.
Tenenemos que y'
p
= B + 2Ax, y' = 2Ay sustituyendo en (4.36), resulta
2A + 3(2Ax + B) + 2{Ax
2
+ Bx + C) = 3x
2
- x + 1
Comparando coeficientes en la ltima igualdad obtenemos el sistema de ecuaciones linea-
les
2A = 3
6A + 2B = - 1
2A + W + 2C = 1,
del cual se sigue que
A-\, B.^. C-l
As
y la
que
solucin general es
Vi
y(x) = ce
,(x) = - x
2
-
3 4- c
2
e
5x4
3
+
2
X
13
13
+
y
EJEMPLO 2. Resolver
/ + %' + 4y = e
3x
. (4.38)
Solucin. En este caso la ecuacin caracterstica es r
2
+ 4r + 4 = 0 y tiene las races
T\ = r<i 2, por lo que
y
c
(x) = c
x
e + c
2
xe
es la solucin general de la ecuacin homognea y" + 4y' + 4y = 0.
Buscamos ahora una solucin particular de la forma
y
p
(x) = Ae**.
Al derivar y sustituir en (4.38) tenemos
9Ae
3x
+ 12Ae
3x
+ Ae
3x
= e
3x
25Ae
3x
= e
3x
,
por lo cual A = ,
VM=y-
y{x) = c
x
e-
2x
+ c
2
xe~
2x
+ e
3x
.
Zo
EJEMPLO 3. Resolver
y
//
+ j / = cos2x. (4.39)
Solucin. La ecuacin caracterstica r
2
+ r = 0 tiene por races r\ = 0, r
2
= 1, por lo
que
y(x) = ci + c
2
e~
x
.
En este caso se propone
y
v
A eos 2x + B sen 2x.
4.6. Mtodo de Coeficientes Indeterminados: Superposicin 151
Sustituyendo y
p
,y'
p
y y
p
en (4.39) encontramos que
(4^4 eos 2x 45 sen 2x) + (2A sen 2x + 25 eos 2x) = eos 2x,
por lo cual
-4A + 2B = 1,
-2A-4B = 0.
Resolviendo el sistema de ecuaciones obtenemos
A - -
1
- B - l
A
~ 5'
fl
"l0-
En consecuencia
y
p
(x) = - - eos 2x + sen 2x
o 1U
y(x) c\ + c
2
e~
x
- eos 2x + sen 2x.
5 1(J
EJEMPLO 4. Resolver
y" + y' + y = xe
x
. (4.40)
Solucin. La ecuacin diferencial homognea tiene la ecuacin auxiliar r
2
+ r + 1 = 0,
cuyas races son r\ \ 2 * Entonces
-x/2 / \/5 \/3 \
?/
c
(a;) = e ' I Ci eos -r-aJ + C2 sen - r - x I .
\ 2 2 /
Una solucin particular de (4.40) tiene que ser un producto de funciones, un polinomio
de grado uno por la funcin exponencial, es decir
y
p
= (Ax + B)e
x
.
Se tiene que
y'
p
= {Ax + B)e
x
+ Ae
x
y'; = (Ax + B)e
x
+ Ae
x
+ Ae
x
= (Ax + B)e
x
+ 2Ae
x
.
(Ax + B)e
x
+ 2Ae
x
+ (Ax + B)e
x
+ Ae
x
+ (Ax + B)e
x
= xe
x
Z(Ax + B)e
x
+ ZAe
x
= xe
x
ZAxe
x
+ (3A + 3B)e
x
= xe
x
.
Por consiguiente
ZA = 1
3A + 3B = 0,
de donde
As
y
p
(x) = - ( x - l)e
x
\ -x/2 ( V^ A/ 3 \ 1
/
,
x
y{x) = e ' \ c\ eos x + c
2
sen -x + ~(x l)e .
y Z Z J o
EJEMPLO 5. Resolver
y" -y = e
x
senx. (4.41)
Solucin. La ecuacin auxiliar de la homognea es r
2
1 = 0. Luego
y
c
(x) = 0^ + c
2
e~
x
.
Ahora proponemos
y
p
= e
x
(Acosx + Bsenx).
Se tiene que
y^ e
x
(Asenx + Bcosx) + e
x
(Acosx + Bsenx),
2/p = e
x
(-ylcosx - Bsenx) + 2e
x
(-Asenx + Bcosx) + e
x
(Acosx + Bsenx),
y
f
p = 2e
x
(-Asenx +Bcosx).
Sustituyendo en (4.41) hallamos que
2e
x
(Asenx + Bcosx) e
x
(Acosx + Bsenx) = e
x
senx
(A + 2B)cosx + (-2A B)senx = senx,
de donde
-A + 2B = 0,
-2A-B = 1.
2 1
La solucin del sistema es A = --,B = -. Por lo tanto
5 5
j/
p
(x) = e
x
( - - c o s x - - s enx) ,
y(x) = cie
x
+ c
2
e~
x
e
x
(2cosx + senx).
EJEMPLO 6. Resolver
y" - 9y = 81x
2
+ 7. (4.42)
Solucin. La solucin general de la ecuacin homognea es
y
c
(x) = ci + c
2
e
9x
.
Como se observa, una constante arbitraria es solucin de la ecuacin homognea, o equi-
valentemente 0 es una raz simple de la ecuacin auxiliar, y de acuerdo con la discusin
del caso I una funcin del tipo Ax
2
+ Bx + C no satisface la ecuacin no homognea.
Ms bien, debemos proponer una solucin particular de la forma
y
p
= x(Ax
2
+ Bx + C).
Tenemos que
y
p
= Ax
3
+ Bx
2
+ Cx,
y'
p
=
QAx + 2B - 9(3Ar
2
+ 2Bx + C) = 81x
2
+ 7
-27Ax
2
+ (6A - 18B)x + 2B - 9C - 81x
2
+ 7.
As que
-27A = 81,
6A - 185 = 0,
2B-9C = 7,
de donde A = 3, B = 1, C = - 1 y entonces
En consecuencia la solucin general de (4.42) es
y(x) = ci + c
2
e
9x
- 3x
3
- x
2
- x.
EJEMPLO 7. Resolver
4y" + 4y' + y = e-
x/2
. (4.43)
Solucin. Primeramente ntese que la ecuacin auxiliar de la ecuacin homognea es
(2r + I)
2
- 0,
de modo que su solucin general es
y
c
(x) = (ci + c
2
x)e~
x/2
.
Como 1/2 es una raz doble de la ecuacin auxiliar, y en consecuencia Ae~
x
l
2
\ Axe~
x
^
2
son soluciones de la ecuacin homognea para cualquier constante A, proponemos una
solucin particular de la forma
y
P
=
segn se vio en el caso II-(c). Luego
y'p
=
j
x e x
'
i JT\- O T-/9 r* A T/ 9 r\ A T/9
^ - ~
/ o
" -^
+
2Ae-
x
'
2
) + {~x
2
e-
x
/
2
+ 2Axe-
x
/
2
) + Ax
2
e-
x
>
2
= e~
x
'
2
(Ax
2
- 8Ax + 8A - 2Ax
2
+ 8Ax + Ax
2
)e~
x/2
= e~
x
'
2
8Ae~
x/2
= e~
x
'
2
8A = 1.
Por lo tanto A = - ,
8
y
P
-
8
X
EJEMPLO 8. Resolver
y
n
+ 25y = 10 sen 5x. (4.44)
Solucin. Las races de la ecuacin auxiliar son 5, as que
y
c
{x) = C\ eos hx + C2 sen 5x.
y ahora proponemos una solucin particular de la forma
y
p
x(Acos5x + 3sen5x),
como se estableci en el caso Ill-(b.l). Se tiene que
y
p
= x(5^4 sen 5x + 55 eos 5a;) + (Acosbx + 5sen5x),
y'p = x(25 A eos 5x 25B sen 5x) + (-5A sen 5x + 5B eos 5x) +
(5^4 sen 5x + 5B eos 5x)
= x(25Acos5x 255 sen 5x) + 2(5^4 sen 5x + 5Bcos5x).
x(25A eos 5x 25B sen 5x) + 2(5A sen 5x + 5B eos 5x)+
25(Acos5x + Bsen5x) = 10sen5x
10^1 sen 5x + 10B eos 5x = 10sen5x.
Por consiguiente A = 1, JB = 0,
y
p
= x eos 5x
y(x) = c\ eos 5x + C2 sen 5x x eos 5x.
A continuacin enunciaremos una propiedad de las ecuaciones diferenciales lineales
que nos permitir resolver ecuaciones no homogneas, cuyo lado derecho es ms general.
Teorema 4.6.1 (Principio de Superposicin.) Si y
Pl
es una solucin particular de
de la ecuacin diferencial
y" + p(x)y' + q{x)y = #i(z),
y y
P2
es una solucin particular de la ecuacin
entonces y
Pl
+ y
P2
es una solucin particular de la ecuacin
y" + p(x)y' + q{x)y = g
x
(x) + g
2
(x).
EJEMPLO 9. Resolver
/ + 6y' + 8y = 3e
5x
+ x
2
- 1. (4.45)
Solucin. Claramente la solucin general de la ecuacin homognea es
y
c
(x) C\e~
2x
+ C2e~
Ax
.
Una solucin particular de
y" + 6y +8y = 3e
5x
, (4.46)
tiene la forma
y
Pl
= Ae
5x
.
Derivando y sustituyendo en (4.46) se obtiene
30Ae
5x
+ 8Ae
5x
= 3e
5x
QZAe
bx
- 3e
bx
,
de donde A ~ y entonces
21
J
Por otra parte, una solucin particular de
y" + 6y' + 8y = x
2
- l , (4.47)
es de la forma
y
P2
= Ax
2
+ Bx + C,
que al derivar y sustituir en (4.47) nos conduce a
2A + 6{2Ax + B) + 8(Ax
2
+ Bx + C) = x
2
- 1
8Ax
2
+ (12A + 8B)x + (2A + 6B + 8C) =x
2
-l.
Igualando los coeficientes de las respectivas potencias de x se sigue
8A = 1,
12A + 8B = 0,
2A + 6B + 8C = - 1 ,
por lo cual hallamos que
A-
1
5 - - 1 C- -
_ 1
2
3 1
y
*-8
X
~~ 16
X
~64'
Por el principio de superposicin, una solucin particular de (4.45) es y
p
y
Pl
+ y
P2
, es
decir
1
5x 1 2 3 1
=
2
C +
8 ^ " 1 6 * "64-
Por lo tanto la solucin general de (4.45) es
21 8 16 64
EJEMPLO 10. Resolver
y" - 2y' -3y = xe
3x
+ sen x. (4.48)
Solucin. La ecuacin auxiliar de la ecuacin diferencial homognea es r
2
2r 3 = 0,
cuyas races son r\ = 3, T2 = 1, por lo cual
Una solucin particular de
y" - 2y' - 3y = ze
3
* (4.49)
tiene la forma
2/
Pl
= x(Ax + )e
3x
,
puesto que 3 es una raz caracterstica simple (ver Caso II-(b)). Luego
y
pi
= (Ax
2
+ Bx)e
3x
,
y'
p
\ = 9{Ax
2
+ Bx)e
3x
+ 6(2Ax + B)e
3x
+ 2Ae
3x
,
y sustituyendo en (4.49) se sigue que
9(Ax
2
+ Bx)e
3x
+ 6(2Ax + B)e
3x
+ 2Ae
3x
~6(Ax
2
+ Bx)e
3x
- 2(2Ax + B)e
3x
- 3{Ax
2
+ Bx)e
3x
= xe
3x
4(2Ax + B)e
3x
+ 2Ae
3x
= xe
3x
(8Ax + 2A + 4B)e
3x
= xe
3x
Q A rp _|_ O A _J_ A D rp
KjTXJb 1^ ali. "I '-X.JLJ Jb ^
de donde
8A = 1,
2A + AB = 0.
Resolviendo el sistema hallamos que A= -,B = y entonces
8 l o
1^2 _ J_
8 16
Por otro lado, una solucin particular de
y" - 2y - 3y = senx, (4.50)
es de la forma
y
P2
= A eos x + B sen x.
Derivando
y'
P2
^4senx + JScosx,
y
f
p
2
Acosx Bsenx,
y sustituyendo en (4.50) se obtiene
A eos x B sen x 2(A sen x + Bcosx) 3(A cosx + B sen x) = sen x
( - 4 A- 25) eos x + ( 2A- 45) sen x = senx,
de donde
i,B=-i,
-4>1-2B - 0,
2, 4- 45 - 1,
y
P2
= eos x - - sen x.
1U 5
Por el principio de superposicin, una solucin particular de (4.48) es
= ^-x - xj e + eos x - - sen x.
/ \ 3x -x
1
2
1
\ 3x
1
1
y(x) = cie
d:E
+ c
2
e * + -x^ - x e
dx
+ cosx - -
\ 8 lo / 10 5
senx.
4.7 Mtodo de Coeficientes Indeterminados: Enfoque
del Operador Anulador
4.7.1 Operadores Diferenciales
d
n
y
La derivada de orden n , tambin se denota por D
n
y. De modo que la ecuacin
dx
n
diferencial
d
n
y d
n
~
1
y dy
a + a + +a
4.7. Mtodo de Coecientes Indeterminados: Operador Anulador 159
puede escribirse en la forma
a
n
D
n
y + a
n
^D
n
~
l
y + ... + a
x
Dy + a
o
y = g(x)
o bien
(a
n
D
n
+ a
n
^D
n
-
1
+ ... + a
x
D + a
o
)y = g(x). (4.51)
La expresin
P(D) = a
n
D
n
+ dn-xD
71
-
1
+ ... +
ai
D + a
0
(4.52)
se llama Operador Diferencial de Orden n. Ntese que P(D) es un polinomio en el
smbolo D y que es posible expresar la ecuacin (4.51) en la forma compacta
Si los coeficientes a
o
,ai, ...,a
n
en (4.52) son nmeros reales, entonces P(D) tiene las
siguientes propiedades.
1. P(D) se puede factorizar como un producto de operadores diferenciales de orden 1
y operadores diferenciales de orden 2 que son irreducibles, con coeficientes reales.
2. Los factores de P(D) pueden conmutarse.
3. P(D)(f + g) = P(D)f + P(D)g, para cualesquier funciones / y g que sean al menos
n veces derivables.
EJEMPLO 1. Factorice si es posible.
a) P(D) = D
2
+ 5D + 6
b) P(D) = D
2
+ D + 1
c) P(D) = D
3
+ 4D
2
+ 3D
d) P(D) = D
S
+ 4D
e) P(D) = D*- 8D
2
+ 16
Solucin.
b) Es un operador cuadrtico irreducible.
c) P(D) = D(D
2
+ 4D + 3) = D{D + 1)(D + 3).
d) P(D) = D(D
2
+ A).
e) P(D) = (>
2
- 4)
2
= [(D + 2){D - 2)]
2
= (D + 2)
2
(D - 2)
2
.
Definicin 4.7.1 Sea y = f(x) una funcin que tiene al menos n derivadas. Si
(a
n
D
n
+ On-i/T
1
"
1
+ ...aiD + ao)/(x) = 0
entonces decimos que el operador diferencial P(D) = a
n
D
n
+ a
n
^\D
n
~
l
+ ... + aiD + a
0
anula a f.
Observe que la definicin 4.7.1 es equivalente a decir que la funcin / es solucin de
la ecuacin diferencial P(D)y 0.
Al operador diferencial del tipo (4.52) con a
n
\ (mnico) y de menor orden posible
que anula a / le llamaremos el operador anulador de / .
Los siguientes resultados nos dicen como encontrar el operador anulador para algunas
funciones y estn basados en nuestros conocimientos sobre las soluciones de las ecuaciones
diferenciales con coeficientes constantes.
Proposicin 4.7.1 El operador diferencial D
n
anula a cada una de las funciones
Ms generalmente, la funcin polinomial f(x) akX
k
+ ak-\X
k
~
l
+ ... + a\X + a$ es
anulada por el operador D
n
con n > k + 1.
EJEMPLO 2. Encuentre el operador anulador P\(D) de la funcin dada.
a) f(x) = x
3
b) f(x) = x
5
- 2x + 1
c) f(x) = x
4
(l + 2x-3x
2
)
Solucin.
b)
c) P
1
(D) = D
7
.
Proposicin 4.7.2 El operador diferencial (D a)
k
anula a cada una de las funciones
ax ai 2 ax fc-1 ax
EJEMPLO 3. Hallar el operador anulador P\{D) de la funcin indicada.
a) f(
x
) = e
7x
b) f(x) = 3e"
4x
+ 2xe~
4x
c) h(x) = e
2x
+ 5xe~
3x
d) g(x) = x- 7xe
6x
e) f(x) = (2 - e*f
Solucin.
a) Tomando a = 7, k = 1 en la proposicin 4.7.2, tenemos que P\{D) = D 7.
b) Ahora, a = - 4 y A; = 2 de modo que P\{D) = (D + 4)
2
.
c) Se tiene que (D - 2)e
2x
= 0 y (D + 3)
2
5ze-
3:c
= 0. Entonces el producto Pi{D) =
(D -2)(D + 3)
2
es el anulador de h. En efecto
Pi(D)h(x) = (D-2)(D + Z)
2
(e
2x
+ 5xe-
3x
)
= (D + 3)
2
(D - 2)e
2x
+ 5(D - 2){D + 3)
2
xe'
3x
= (D + 3)
2
(0) + 5( D- 2) ( D- 2) ( 0)
= 0.
d) Para esta funcin aplicamos las proposiciones 4.7.1 y 4.7.2. En primer lugar D
2
x = 0
y por otra parte (D - 6)
2
7xe
6x
= 0. Por lo tanto P
X
(D) = D
2
(D - 6)
2
.
e) La funcin / dada no es de los tipos mencionados en las proposiciones anteriores.
Sin embargo, f(x) = 4 - 4e
x
+ e
2x
y ya que L>4 = 0, (D - l)4e* = 0, (D - 2)e
2x
= 0
concluimos que Pi(D) = D(D - 1)(D - 2).
Proposicin 4.7.3 El operador diferencial [D
2
2aD + (a
2
+ {3
2
)]
k
anula a cada una
de las funciones
x, xe
ax
cos/3x, x
2
e
ax
cos/3x, ...,x
k
~
1
e
ax
cos/3x,
x, xe
ax
sen fix, x
2
e
ax
sen fix, ...,x
k
~
l
e
ax
senf3x.
En particular, si en la proposicin 4.7.3, consideramos a = 0yfc = l, se sigue que
(JD
2
+/ 3
2
)COS/ ?X = 0,
(D
2
+ (3
2
)sen(3x = 0.
Ejemplo 4. Obtenga el operador anulador P\{D) de la funcin dada
a) f(x) = 2e
x
cos3x e
x
sen3x
b) g(x) = 2e
x
eos 3x + xe
x
eos 3x e
x
sen 3x
c) h(x) = 4 + 3x - 5 eos 2x
d) f(x) = 8x
3
- senx + 10cos5x
e) g(x) = x
2
e~
x
senx 6e
2x
cos3x
Solucin.
a) En la proposicin 4.7.3 usamos los valores de a = 1, /3 = 3 y k = I para encontrar
h) Igual que en el inciso anterior, a = 1 y fi = 3, pero ahora k = 2. Entonces
(D
2
- 2D + 10)
2
.
c) Aplicando la proposicin 4.7.1 se tiene que D
2
(A + 3x) = 0. Mientras que, de la
proposicin 4.7.3 con a = 0 y /? = 2, se sigue que (D
2
+ 4)(5cos2x) = 0. Por lo tanto
d) De las proposiciones 4.7.1 y 4.7.3 tenemos que
D
4
(Sx
3
) = 0, (D
2
+ l )senz = 0, (D
2
+ 25)10 eos 5x = 0.
Por consiguiente P
X
(D) = D
A
(D
2
+ 1)(D
2
+ 25).
e) Aplicamos la proposicin 4.7.3 primero on a = 1, / ? = l y f c = 3 y luego con
a = 2, /? = 3 y k = 1, para obtener
- 0
- 0,
de manera que el operador anulador en este caso viene dado por
P
X
(D) = (D
2
+ 2D + 2f{D
2
- 4D + 13).
4.7.2 Mtodo de los Coeficientes Indeterminados
Consideremos la ecuacin diferencial no homognea
P(D)y = g(x), (4.53)
donde P(D) = a
n
D
n
+ a
n
-iD
n
~
1
+ ... + aiD+ao es un operador diferencial con coeficientes
constantes y g(x) es
a) un polinomio en x,
b) una funcin exponencial e
ax
,
c) sen/fc, cos/3x
d) sumas finitas y productos de las funciones mencionadas en (a), (b) y (c).
En este caso siempre es posible encontrar el operador anulador P\(D) de g(x). Apli-
cando P\(D) a (4.53) resulta
P
1
(D)P(D)y = P
1
(D)g(x) = 0. (4.54)
Resolviendo la ecuacin diferencial homognea (4.54) es posible descubrir la forma de
una solucin particular y
p
de la ecuacin diferencial no homognea (4.53). Los siguientes
ejemplos esclarecern el procedimiento a seguir.
EJEMPLO 1. Resolver
y - 2y' + y = x
2
+ 4x. (4.55)
Solucin. Primero resolvemos la ecuacin homognea asociada
y" -2y' + y = 0.
De la ecuacin caracterstica r
2
2r +1 = (r I)
2
, se obtiene la solucin complementaria
y
c
= (ci +c
2
x)e
x
.
Ahora bien, el operador anulador de la funcin g(x) = x
2
+ 4x que figura en el lado
derecho de la ecuacin diferencial (4.55) es P\{D) = J9
3
; es decir
D
3
{x
2
+ 4x) = 0.
Aplicando P\(D) D
3
a (4.55), obtenemos la ecuacin homognea
D\D
2
-2D + l)y - D\x
2
+ 4x) = 0. (4.56)
La ecuacin caracterstica de (4.56) es
r
3
( r
2
- 2 r + l) = 0
o bien
r
3
( r - l )
2
= 0
y por lo tanto su solucin general es
y=(k
x
+ k
2
x)e
x
+ k
3
+ k
A
x + k
5
x
2
, (4.57)
donde fc
1?
k
2
,..., k$ son constantes arbitrarias.
Es posible argumentar que toda solucin de (4.56) es tambin una solucin de (4.55).
Dado que la solucin complementaria y
c
= {c\ +C2x)e
x
aparece como parte de la solucin
de (4.56), los trminos restantes en (4.57) deben ser la estructura bsica de y
p
y
p
= A + Bx + Cx
2
. (4.58)
Para simplificar la notacin hemos reemplazado Aft, k y k$ por A, B y C respectiva-
mente. Finalmente, calcularemos los valores A, B y C de modo que (4.58) sea en efecto
una solucin de la ecuacin diferencial (4.55).
Tenemos que
y'
p
= B + 2Cx,
y; = 2C.
Sustituyendo estas expresiones en (4.55) resulta
2C - 2(B + 2Cx) + (A + Bx + Cx
2
) - x
2
+ 4x.
Equivalentemente
Cx
2
+ (B
2
Igualando coeficientes en la ltima identidad, se obtiene el sistema de ecuaciones
lineales
A-2B + 2C = 0,
B-AC = 4,
C = 1.
Resolviendo el sistema, resulta A = 14, B = 8 y C 1. En consecuencia
y
p
= 14 + 8x + x
2
.
La solucin general de (4.55) es y = y
c
+ y
p
, o sea
y
= (
Cl
+ c
2
x)e
x
+ 14 + Sx + x
2
.
EJEMPLO 2. Resolver
y" - 2y' -8y= 16e
2x
- 21
e
-
3 x
. (4.59)
Solucin. La ecuacin caracterstica de la ecuacin diferencial homognea asociada a
(4.59) es r
2
- 2r - 8 = (r - 4)(r + 2) = 0, as que
y
c
= c
x
e~
2x
+ c
2
e
4x
.
Note que la ecuacin diferencial (4.59) puede escribirse como
(D - 4)(D + 2)j/ = 16e
2x
- 21e~
3x
. (4.60)
El operador anulador de la funcin g(x) = 16e
2x
21e~
3x
es
y al aplicarlo a ambos miembros de (4.60) resulta
(D - 2){D + 3){D - 4)(D + 2)y = 0. (4.61)
( r - 2) ( r + 3) ( r - 4) ( r + 2) =0.
Por consiguiente
y = fc
ie
2x
+ k
2
e'
3x
+ fc
3
e
4x
+ k
4
e~
2x
.
y
p
= Ae
2x
+ Be~
3x
.
Sustituyendo yp,y'
p
,y'p en (4.59) y simplificando, obtenemos
-8Ae
2x
+ 7Be'
3x
= 16e
2x
- 21e"
3x
.
Igualando coeficientes, se sigue que
-8A = 16,
IB = - 21.
Encontramos A = 2, B = 3 y por lo tanto
y y
=
B &
La solucin general de (4.59) es entonces
y
=
Cl
e~
2x
+ c
2
x
- 2e
2x
- 3e"
31
.
EJEMPLO 3. Resuelve
y" - 2y +3y = e~
x
eos x. (4.62)
Solucin. Podemos escribir (4.62) como
(D
2
-2D + 3)y = e~
x
eos x. (4.63)
La ecuacin caracterstica de la ecuacin diferencial homognea asociada es r
2
2r+3 = 0,
cuyas races son 1 y/2i. Luego
y
c
(x) = e
x
(ci eos y/2x + C2sen y/2x).
El operador anulador de la funcin en el lado derecho de (4.62) es
Aplicando P\(D) a (4.63) obtenemos
(D
2
+ 2D + 2)(D
2
-2D + 3)y = 0. (4.64)
La ecuacin caracterstica de (4.64) es (r
2
+ 2r + 2)(r
2
2r + 3) = 0 y las races son
1 \/2, 1 i de multiplicidad uno, por lo cual
y = e
x
(ki eos \2x + fc
2
sen\/2x) + e~
x
(k^cosx + fc
4
senx).
Luego
y
p
= e~
x
(Acosx + Bsenx).
Tenemos que
y'
v
= -e~
x
(^4cosx + 5senx) + e~
x
(-Asenx + Bcosx),
Sustituyendo estas expresiones en (4.62) y simplificando, resulta
e~
x
[(5A- AB)cosx + (4A +5B)senx] = e~
x
Comparando coeficientes, se obtiene el sistema de ecuaciones
5A-AB = 1,
AA + 5B = 0,
5 4
de donde A , B y en consecuencia
41 41
y
p
(x) = t~
x
( cos x- senx).
La solucin general de (4.62) es
y e
x
(ci eos y/2x + C2sen V2x) + e~
x
( eos x Tjsen x).
EJEMPLO 4. Resuelve
y" + 3y
f
- lOy = x(e
2x
+ 1). (4.65)
Solucin. Escribimos la ecuacin diferencial como
(D
2
+ 3D- 10)y = x(e
2x
+ 1)
o bien
(D + 5)(D - 2)y = x{e
2x
+ 1). (4.66)
Se ve que
y
c
(x) = C!e-
5x
+ c
2
e
2x
.
Ahora bien, sabemos que D
2
x = 0, (D 2)
2
{xe
2x
) = 0. Por consiguiente el operador
anulador de la funcin g(x) = x(e
2x
+ 1) = xe
2x
+ x, es P\(D) = D
2
(D 2)
2
. Aplicando
P\(D) a ambos miembros de (4.66) se obtiene
D
2
(D + 5)(D-2)
3
y = 0. (4.67)
La ecuacin caracterstica de (4.67) es r
2
(r + 5)(r 2)
3
= 0 y sus races son 0, - 5 y 2 de
multiplicidades dos, uno y tres, respectivamente. Luego
y = ki + k
2
x + k
3
e'
5x
+ (fc
4
+ k$x + k
6
x
2
)e
2x
.
Excluyendo en esta expresin la combinacin lineal de trminos correspondientes a y
c
, se
llega a la forma de y
p
y
p
= A + Bx + (Cx + Ex
2
)e
2x
.
Se tiene que
y
'
p
= B + (C
Sustituyendo en (4.65) y simplificando, resulta
(-1L4 + 3JB) - lOBx + (7C+2E)e
2x
+ UExe
2x
= xe
2x
+ x.
Esto implica que
-10A + 3B = 0,
-10J3 = 1,
7C + 2E = 0,
UE - 1.
3 1 1 1
As que A = , B = , C = , E y la solucin general de (4.65) es
y =
Cl
e + c
2
e ,
+
( ,
EJEMPLO 5. Resuelve
y" - 2j/' + 5y = e
x
eos 2x + e
x
sen 2x. (4.68)
Solucin. La ecuacin caracterstica de la ecuacin diferencial homognea asociada es
r
2
- 2r + 5 = 0 y sus races son 1 2z, de modo que
y
c
(x) e
x
(c\ eos 2x + C2sen 2x).
Ntese entonces que el operador anulador del lado derecho es precisamente P$D) =
D
2
2D + 5. En consecuencia
(D
2
-2D + 5)(D
2
-2D + 5)y = 0
o bien
(D
2
-2D + 5)
2
y = 0.
Las races de la ecuacin caracterstica de esta ltima ecuacin diferencial homognea
son 1 2z, de multiplicidad dos. De sto deducimos que
y e
x
(ki eos 2x + ^sen 2x) + xe
x
(ks eos 2x + fc
4
sen 2x).
Sustituyendo
y
p
xe
x
(A eos 2x + Bsen 2x)
en (4.68) y simplificando obtenemos
e
x
(4tB eos 2x - A Asen 2x) = e* eos 2x + e
x
sen 2x.
Igualando coeficientes resultan las ecuaciones
-4A = 1,
4B = 1,
de donde se sigue que A y J3 = - . Por lo tanto la solucin general de (4.68) es
X
y = e
x
(ci eos 2x + C2sen 2x) He
x
(sen 2x eos 2a;).
EJEMPLO 6. Resuelve
y" + 2y' + y = -Ze~
x
+ 8xe~
x
+ 1. (4.69)
Solucin. Escribimos la ecuacin diferencial (4.69) en la forma
(D + ify = -3e~
x
+ Sxe~
x
+ 1. (4.70)
Luego
y
c
(x) = (ci + c
2
x)e~
x
.
Se tiene que (D+l)
2
(-3e-
x
+8xe-
x
) = 0 y D{) = 0. Por consiguiente P\(D) = D{D+$
2
es el operador anulador de la funcin g(x) = -3e~
x
+ 8xe~
x
+ 1. Aplicamos P\(D) a
(4.70) para obtener.
D(D + l)
4
y = 0. (4.71)
La ecuacin caracterstica de (4.71) es r(r + I)
4
= 0 y sus races son 0 y -1 de multipli-
cidades uno y cuatro, respectivamente. Por consiguiente
y = ki + (fc
2
+ h% + hx
2
+ k
5
x
3
)e~
x
.
As , una solucin particular de (4.69) tiene la forma
Sustituyendo y
p
en (4.69), resulta
A + 2Be~
x
+ 6Cxe~
x
= -3e~
x
+ 8xe'
x
+ 1,
O A
de donde AI, B= y C = - , d e modo que la solucin general de (4.69) es
z o
y
= (
Cl
+
C2
x)e-
X
+ (^x
3
- \x
2
)e~
x
+ 1.
EJEMPLO 7. Resuelve
y
" + 4y = 3x eos 2z + sen 2x. (4.72)
Solucin. La ecuacin (4.72) puede escribirse como
(D
2
+ 4)j/ = 3x eos 2x + sen 2x. (4.73)
Es claro que
y
c
(x) = Ci eos 2x + C2 sen 2x.
Adems P\(D) = (D
2
+4)
2
es el operador anulador de la funcin g(x) = 3x eos 2x+ sen 2x.
Aplicamos P\{D) a (4.73) para obtener
(>
2
+ 4)
3
j/ = 0. (4.74)
Las races de la ecuacin caracterstica de (4.74) son 2z de multiplicidad tres, por lo
que
y = ki eos 2x + A/2 sen 2x + x(fc
3
eos 2x + A:
4
sen 2x) + x
2
(fc
5
eos 2x + k$ sen 2x).
Luego, buscamos una solucin particular de la forma
y
p
= x(^4cos2x + Bsen2x) + x
2
(Ccos2x + Esen2x).
Tenemos que
y'
v
= (A eos 2x + B sen 2x) + 2x{-A sen 2x + B eos 2x) + 2x(C eos 2x + E sen 2x)
+2x
2
(-C sen 2x + eos 2x),
/p = 4(-i4 sen 2x + S eos 2x) - 4x(Aeos 2x + B sen 2x) + 2(Ceos 2x + E sen 2x)
+8x(-Csen2x + 'cos2x) - 4x
2
(Ccos2x + Esen2x).
Sustituyendo en (4.72), resulta
(4B + 2C) eos 2x + (2E - AA) sen 2x - 8Cx sen 2x + 8Ex eos 2x = 3x eos 2x + sen 2x.
Igualando coeficientes, se obtiene el sistema de ecuaciones
4B + 2C = 0,
2E-AA = 1,
- 8C - 0,
8E = 3,
1 3
de donde se sigue que A --, B = 0, C = 0, E -. Por lo tanto la solucin general
16 8
de (4.72) es
1 3
y C\ eos 2x + C2 sen 2x -x eos 2x + -x
2
sen 2x.
16 8
EJEMPLO 8. Determine la forma de una solucin particular de
y" + 61/ + 13y = xe"
3x
sen 2x + x
2
e~
2x
sen 3x + 6. (4.75)
Solucin. La ecuacin caracterstica de la ecuacin diferencial homognea asociada a
(4.75) es r
2
+ 6r + 13 = 0, cuyas races son 3 2%. En consecuencia
y
c
= e~
3x
(ci eos 2x + C2 sen 2x).
Ahora bien, tenemos que (D
2
+6D+13)
2
(xe~
3a;
sen2x) = 0, (D
2
+4
J
D+13)
3
(x
2
e-
2x
sen3x)
0 y D(6) = 0, por lo cual el operador anulador de g(x) = xe"
3x
sen 2x + x
2
e~
2x
sen 3x + 6
es
Pi(D) = {D
2
+ 6D + 13)
2
(D
2
+ 4D + 13)
3
D.
Al aplicar P\(D) a (4.75) resulta
D(D
2
+ 4L> + 13)
3
(D
2
+ 6D + 13)
3
y = 0. (4.76)
Puesto que las races de la ecuacin caracterstica de (4.76) son 3 2z, 2 3z y 0 de
multiplicidad tres, tres y uno, respectivamente. Concluimos que
y = e~
3x
(ki eos 2x + k
2
sen 2x) + xe~
3x
(k
3
eos 2x + fc
4
sen 2x)
+ e~
2x
(fc
7
cos3x + A:
8
sen3x)
x
2
e~
2x
(k
u
cos3x + fci2
Por lo tanto se puede encontrar una solucin particular de (4.75) que tiene la forma
y
p
= xe~
3x
(A eos 2x + B sen 2x) + x
2
e~~
3x
(C eos 2x + E sen 2x)
+ e~
2x
(F eos 3x + G sen 3x) + xe~
2x
(H eos 3x + / sen 3x)
+ x
2
e~
2x
( J eos 3x + /f sen 3x) + L.
EJEMPLO 9. Determinar la forma de una solucin particular de
y"'
- 5y" + y
1
- 5y = x
2
- 3 + x
3
e
bx
- 3x sen x. (4.77)
Solucin. La ecuacin caracterstica de la ecuacin homognea asociada es
r
3
- 5r
2
+ r - 5 = 0
o bien
(r
2
+ l)(r - 5) - 0.
De modo que
y
c
C\ eos x + C2 sen x + C3e
5x
.
La ecuacin diferencial (4.77) puede escribirse como
(D
2
+ 1)(D -5)y = x
2
-3 + x
3
e
5x
- 3xsenx. (4.78)
El operador anulador de la funcin en el lado derecho de (4.77) es
Aplicando P\(D) a (4.78) obtenemos
D
3
(D
2
+ 1)
3
(D - bfy = 0. (4.79)
4.7. Mtodo de Coeficientes Indeterminados: Operador Anulador 171
Ahora, la ecuacin caracterstica de (4.79) es r
z
(r
2
+ l)
3
(r 5)
5
= 0 y sus races son
0, i y 5 de multiplicidad tres, tres y cinco, respectivamente. Por consiguiente
y = ki + k
2
x + k
3
x
2
+ (k
4
+ k
5
x + k
e
x
2
+ k
7
x
3
+ (kg + kwx + k
u
x
2
) eos x + (ki2 + k
i3
x + k
u
x
2
) sen x,
y una solucin particular de (4.77) es de la forma
y
p
= A + Bx + Cx
2
+ (Ex
+ (Ix + Jx
2
)cosx + (Kx + Lx
2
)senx.
EJERCICIOS 4.7
En los problemas 1-6 factorice el operador diferencial dado, cuando sea posible
1. 16>
2
- 9
2. D
2
-llD + 24
3. D
3
+ 12D
2
+ 36D
4. D
3
- 5D
2
+ 2D + 8
5. >
4
6. D
4
En los problemas 7-15 encuentre el operador anulador de la funcin dada
7. 7 + 8 z - 5 z
2
8. x
9
(l + Ax)
9. 6 - l i e
4 1
10. 9x + 5xe
101
11. 3 + x - 6x
2
- eos 9x + 2 sen9x
12. 6x + eos 2x + 21 sen7 x
13. e~
x
+ 2xe
x
- x
2
e
x
14. 5 + e
3z
sen4x
15. e~
2x
eos 3a; e
6x
sen3x + e~
2x
sen3x
En los problemas 16-49 resuelva la ecuacin diferencial dada, por el mtodo de los
coeficientes indeterminados
16. y" + 8y' + 7y = 32e
x
- 27e
2x
17. y" - 2y' + y 10 eos x + 8 sen x
18. y" - 2y' + Yly = 289x
2
+ 9
19. y" + 8y' = 64x
20. y" + Ay
1
+ 4y = 8e~
2x
21. y" + 4y' + 3y = 4
e
-
3 x
+ 18x + 15
22. y" + 10y' + 25y = (2 + 6x)e~
5x
23. y" + 3y' = 36xe-
3x
- 9e-
3x
+ 7
24. y" + 5y' + 6y = 10(1 - z)e~
2x
25. y" + 2y' + 2y = l + x
26. y" + y' + y = (x + x
2
)e
x
27. y" + 4y' -2y = 8sen2x
28. y" + y = 8 sen x + 4 eos x + 1
29. y" + y = 4x eos x
30. y" + y = 2 sen x sen 2x
31. y" + y' = eos
2
x + e
x
+ x
2
32. y" - 4y' + 5y = e
2x
( senx + 2 eos x)
33. y" 4y' + 4y = 4x + sen x + sen 2x
34. y" + 2y' + y = 1 + 2 eos x + eos 2x sen 2x
35. y" + 6y' + 9y = 6xe"
3x
+ 9 + 50 sen x
36. y" + 4y' + 4y = 2x + 6
37. y" - y' - 12y = 14e
4x
38. y" - 2y' - 3y = 4e
x
- 9
4.8. Mtodo de Variacin de Parmetros 173
39. y"-y= 12xV - 10
40. y" - 2y' + 5y = 3e
x
sen x
41. y" + 25y = 20sen5x
42. y" + 2y' + 5y = 50x sen x + lOOx eos x
43. y"' + y" = l
44. y<
5)
- y^ = 2xe
x
+ 24
45. y"
1
+ y" = 12x
2
46. y'" + 2y" = 24x
2
+ 24x + 4
47. y'" - 3y" + 3y' - y = e
x
- x + 16
48. y<
4)
- 2y'" + y" = e
x
+ 1
49. 16y
(4)
- y = 16e
En los problemas 50-54 determine la forma de una solucin particular de la ecuacin
diferencial dada
50. y" - y = e
x
(2 + 3x eos 2x)
51. y" - 3y' = 6 + xe
3x
+ x
2
e"
x
senx
52. y" + 4y = senxsen2x
53. y'" y' 4xe
x
+ e
x
eos x sen 2x
54. y
(4)
+ 2y'" - 3y" = x
3
+ 1 + x
2
sen x + xe
x
4.8 Mtodo de Variacin de Parmetros
Este es un mtodo general para determinar una solucin particular de una ecuacin
diferencial lineal. Sin prdida de generalidad, consideremos la ecuacin lineal de segundo
orden escrita como
y" + p(x)y' + q(x)y = g(x). (4.80)
El mtodo consiste en buscar una solucin de la forma
y
p
(x) = i(x)yi(x) + it
2
(x)y
2
(x), (4-81)
donde y\ y y
2
son dos soluciones l.i. de la ecuacin diferencial homognea asociada, y
ti,i2 son dos funciones a determinar de modo que (4.81) sea una solucin de (4.80) y
satisfagan una condicin arbitraria, pero seleccionada de tal forma que se simplifiquen
los clculos.
Derivando (4.81) se tiene que
y'
p
= u
x
y
x
+ uiy[ + u
2
y
2
+ u
2
y
2
[ + u
2
y
f
2
) + (u[yi + u'
2
y
2
).
Podemos simplificar esta expresin, imponiendo a u\ y u
2
la condicin de que
u\yi + u
2
y
2
= 0.
En tal caso
y'
p
= uiy[ + u
2
y'
2
y por consiguiente
y'p
Sustituyendo las expresiones de y
p:
y
f
p
^ y
f
en (4.80), y usando el hecho de que y\ y y
2
son soluciones de la ecuacin diferencial homognea, resulta
u[y[ + my" + u
2
y
2
+ u
2
y
2
+ p(uiy[ + u
2
y
2
) + q(u
l
y
1
+ u
2
y
2
) = g{x)
+ vy'
2
+ 92/2) + u[y[ + u
2
y
2
= g(x)
As, buscamos una solucin particular de la forma (4.81), con ii,i2 funciones que
satisfacen las ecuaciones
u'
iyi
+u'
2
y
2
1= 0, (4.82)
= 9{?)' (4-83)
Es fcil resolver el sistema de ecuaciones (4.82)-(4.83) para las incgnitas u[ y u
2
,
empleando la regla de Cramer. Obtenemos
V2{x)g(x) yi(x)g(x)
'
( }
donde W(x) denota al wronskiano W(yi,y2)(x). Finalmente, integrando las expresiones
(4.84) resulta
Sustituyendo (4.85) en (4.81) se obtiene la solucin particular deseada.
EJEMPLO 1. Resolver
y" - 2y' + 2y = e
x
sec x. (4.86)
Solucin. Deteminamos primero la solucin general de la ecuacin homognea asociada
a (4.86), a saber
y" - 2y' + 2y = 0. (4.87)
r
2
- 2r + 2 = 0,
y sus races son ri = 1 + i y r^ 1 i- En consecuencia
Ve = e
x
(ci cosx + C2senx).
Denotemos por
Luego
= e
x
eos x, ?/2 = e
x
sen x.
e
x
eos x
= e
2x
.
(eos x sen x)e
x
(eos x + sen x)e
x
Buscamos una solucin particular de (4.86) de la forma y
p
i
1
y
1
+ i2?/2, donde las
funciones u\ y t2 se calculan utilizando las ecuaciones (4.85 ). Se tiene que
dx = I tanxdx = ln | cosx|,
f (e
x
cosx)(e
x
secx) r
1
) = J 2 = J dx = X.
Luego
j /
p
= (ln | eos x |) (e* eos x) + xe
x
sen x
y por lo tanto, la solucin general de la ecuacin diferencial no homognea (4.86) est
dada por
y = y
c
+ y
P
= e
x
(cicosx + c
2
senx) + (ln | eosx\)(e
x
cosx) + xe
x
senx.
EJEMPLO 2. Resolver
p
-2x
x > 0. (4.88)
Solucin. Puesto que la ecuacin caracterstica es
r
2
+ 4r + 4 - (r + 2)
2
= 0,
se sigue que
y
c
(x) = k
x
e~
2x
+ k
2
xe~
2x
.
Sean y^ = e~
2x
y y
2
= xe~
2x
. Entonces
e
"
2x
= e
"
4
*
-2e~
2x
(l-2x)e-
2x
Usando las expresiones (4.85) obtenemos
u$x) = dx dx lnx
J e~
4x
J x
y
u
2
(x) = / dx = / zdx .
v ;
J e~
Ax
J x
2
x
As que
y
p
= ~(lnx)e-
2x
-(-)xe-
2x
x
= (k
l
+ k
2
x)e-
2x
-(lnx)e-
2x
-e-
2x
.
Ya que fci es una constante arbitraria, ntese que podemos escribir la solucin de
(4.88) simplemente como
?/ = (ci + c
2
x)e~
2x
- (\nx)e~
2x
,
siendo c\ k\ 1 y c
2
= k
2
constantes arbitrarias.
EJEMPLO 3. Para la ecuacin diferencial
xy" + 2y' -xy = 2e
2x
. (4.89)
a) Compruebe que las funciones y\ x~
1
e
x
,y
2
= x~
l
e~
x
, forman un conjunto fun-
damental de soluciones en el intervalo (0, oo) de la ecuacin diferencial homognea
correspondiente.
b) Obtenga la solucin general de la ecuacin no homognea dada.
Solucin.
a) Un clculo directo muestra que
y[ = {x-
1
- x'
2
)e
x
, y'i = (2x~
3
- 2a;"
2
+ x'ê*
y
>
2
= (_a;-i _ x-
2
)e~
x
, y'
2
' = (2x~
3
+ 2x~
2
+ x-
l
)e~
x
y sustituyendo en la ecuacin diferencial, encontramos que
'i + 2y[ - xy
x
= (2x~
2
- 2x~
l
+ $ + 2{x'
1
- x'
2
)e
x
- e
x
= 0
= (2x~
2
2x~
x
l)e~
x
- 2{x~
x
- e"
z
= 0.
Adems
W(
yi
,y
2
) =
(x-
x~
x
e
x
1
x~
2
)e
x
(-X
x~
l
e~
x
~
x
- x~
7
)e
= -2x
- 2
o,
para todo x G (0, oo). Lo cual demuestra que yi y y
2
son dos soluciones l.i. en (0, oo) de
la ecuacin diferencial homognea.
b) Obtenemos ahora una solucin particular de (4.89). Escribimos primero la ecuacin
diferencial en la forma
x x
Consecuentemente
z -dx = / e
x
dx = e
x
,
_ 2
/
T * i
u
2
{x) = j
(x~
1
e
:r
)(2x-
1
e
2x
)
-2x~
2
dx = - (
J
-
y
P
=
La solucin general de (4.89) es
y =
2
3
:
+ c
2
x~'
e-
x
+ -x~
l
e
2x
.
O
EJERCI CI OS 4.8
En los problemas 1 a 12, obtenga una solucin particular de la ecuacin diferencial
dada, utilizando el mtodo de variacin de parmetros. Escriba la solucin general de la
ecuacin diferencial.
1.
2.
3.
4.
5.
y"-
y"-
y"-
y"-
y"-
\-y =
\-y =
h3y'
h6y'
= tan x
- cscx
+ 3, =
1
1 + e
x
: xe~
x
e
-3x
X
2
6. y" - Ay
f
+ % = e
2x
arctan x
7. y" - 3y' + 2y = e
x
sen2x
8. y" + 5y
f
+ 6y = sene
x
9. y" -6y' + 9y = e
3x
lnx
10. y" + y xsenx
12. 9 / - 12?/ + Ay - -
En los problemas 13 a 15 determine la solucin general de la ecuacin diferencial no
homognea dada, usando que la funcin y\ indicada es una solucin de la ecuacin ho-
mognea correspondiente.
13. x
2
y" -2y = x
2
lnx, y
1
= x
2
14. x
2
y" + xy
f
- y = - , y
x
= x
x
15. x
2
y" + xy' + y = tan ln x, y\ eos ln x
Captulo 5
Aplicaciones de las Ecuaciones
Diferenciales Lineales de Segundo
Orden
5.1 Movimiento Armnico Simple
Supngase que un cuerpo de masa m est sujeto al extremo de un resorte flexible (de
peso despreciable), suspendido de un soporte rgido.
Cuando el peso est en reposo, describimos su posicin como la posicin de equilibrio.
Si el cuerpo se desplaza hacia abajo una cierta distancia y luego se suelta, estar bajo
un movimiento vibratorio alrededor de la posicin de equilibrio (ver figura 5.1). Nuestro
propsito es estudiar el movimiento del cuerpo, conocido como movimiento armnico
simple, en el cual se ignora cualquier fuerza de friccin con el medio que lo rodea.
Resorte libre
Posicin de equilibrio
mg - ks = 0
Cuerpo en movimiento
Figura 5.1: Sistema masa-resorte
179
180 Captulo 5. Aplicaciones de las Ecuaciones Diferenciales Lineales de Segundo Orden
En este caso, las nicas fuerzas que actan son:
Una fuerza de restitucin, /
r
, opuesta a la direccin del alargamiento y pro-
porcional a su magnitud (Ley de Hooke). En trminos simples f
r
= kd, donde k es
una constante de proporcionalidad y d la magnitud del alargamiento.
El peso del cuerpo, dado por W = mg.
Adoptaremos la siguiente convencin. Todas las cantidades (desplazamiento, ve-
locidad, aceleracin y fuerza), medidas hacia abajo desde la posicin de equilibrio se
considerarn como positivas. Las que se miden hacia arriba, son negativas.
En la posicin de equilibrio
mg ks = 0.
Ahora, al desplazar el cuerpo de esta posicin en una magnitud x y soltarla, de la Segunda
Ley de Newton se sigue que
d
2
x , , x
m2~ = mg k\s + x)
= mg ks kx,
y usando la condicin de equilibrio, resulta
d
2
x ,
N
m-r = kx. (5-1)
dt
El signo negativo indica que la fuerza de restitucin del resorte acta en direccin opuesta
a la del movimiento.
Podemos escribir la ecuacin (5.1) en la forma
d
2
x k
+ -x = 0,
dt
m
o bien
g+^-0, (5.2)
donde u
2
= k/m.
La ecuacin (5.2) es la ecuacin diferencial del movimiento armnico simple
o movimiento vibratorio no amortiguado.
Hay dos condiciones iniciales asociadas con (5.2), a saber
x(0) = xo, x
f
(0) = vo,
que representan el deplazamiento y velocidad iniciales, respectivamente. Por ejemplo, si
Xo < 0 y ô > 0 entonces el movimiento se inicia en un punto que est |xo| unidades
arriba de la posicin de equilibrio y con una velocidad inicial dirigida hacia abajo. Si
5.1. Movimiento Armnico Simple 181
xo > O y v
0
= O, la masa est inicialmente en reposo a XQ unidades abajo de la posicin
de equilibrio.
La ecuacin auxiliar de (5.2) es
r
2
+ J
2
= 0,
cuyas races son imaginarias puras
T\ = O>Z, T2 = CO.
En consecuencia, la solucin general de la ecuacin diferencial (5.2) es
x(t) = c\ eos or + C2 sear, (5.3)
donde c\ y c<i son constantes que dependen de x
0
y ô-
Ntese que independientemente de los valores de c\ y c
2
, la ecuacin del movimiento
armnico simple (5.3), define una funcin peridica de periodo T = 2TT/UJ y describe un
movimiento ideal en el que el cuerpo se mueve alternadamente hacia arriba y hacia abajo
de la posicin de equilibrio, infinitas veces.
El periodo T es el tiempo necesario para que se complete un ciclo y su recproco
/ = 1/T se llama la frecuencia. El desplazamiento mximo del cuerpo, medido desde
la posicin de equilibrio, se llama la amplitud.
EJEMPLO 1. Se encontr experimentalmente que un peso de 4 Ib estira un resorte 6
pulgadas. Si el peso se suelta desde la posicin de equilibrio con una velocidad dirigida
hacia abajo de 4 pulg/s, determine:
a) La ecuacin diferencial y condiciones iniciales que describen el movimieno.
b) La ecuacin del movimiento.
c) La posicin, velocidad y aceleracin del peso 2 segundos despus.
d) El periodo, la frecuencia y la grfica de la solucin.
Solucin. Ya que 6 pulgadas equivalen a 1/2 ft, de la Ley de Hooke tenemos que
4 = (Jfe
de donde
Adems m = W/g = 4/32 = 1/8 slug.
a) Luego, de (5.2), la ecuacin diferencial que describe el movimiento es
cPx 8
o bien
sujeta a las condiciones iniciales
dt
2
+ 64x = 0, (5.4)
z(o) = o, Ao) = \-
b) La ecuacin auxiliar de (5.4) es r
2
+ 64 = 0, cuyas races son r = 8z. En consecuencia
la solucin general de (5.4) viene dada por
x(t) = c\ eos 8 + c
2
sen 8.
La condicin incial x(0) = 0 implica que c\ = 0, mientras que x
f
(0) = 1/3 conduce a
8c
2
= 1/3. De modo que c\ = 0, c
2
= 1/24 y la solucin requerida es
x{t) sen8.
Z4
Figura 5.2: Solucin del ejemplo 1
c) La posicin, velocidad y aceleracin del peso 2 segundos despus, estn dadas, respec-
tivamente por
x
{2) = sen 16 = -0.011996,
i' (2) = - eos 16 = -0.31922,
x"{2) = - ^ sen 16 = 0.76774,

lo cual indica que el cuerpo se encuentra a 0.011996 ft arriba de la posicin de equilibrio
movindose hacia arriba.
d) El periodo y la frecuencia son
2 7 T 7 T 4
1 ;
Claramente, la amplitud es de 1/24 ft. La solucin muestra que una vez que el sistema se
pone en movimiento, permanece en tal estado con la masa desplazndose alternadamente
1/24 ft hacia cada lado de la posicin de equilibrio x = 0. La grfica se muestra en la
figura 5.2.
EJEMPLO 2. Suponga que en el ejemplo anterior la masa se desplaza 3 pulgadas por
debajo de la posicin de equilibrio y luego se le da una velocidad hacia abajo de 4 pulg/s.
Determine la ecuacin de movimiento.
Solucin. Como antes
x(t) = C\ eos 8 + C2 sen 8,
pero ahora, las condiciones iniciales son
x(0) = \, x'(0) = i
La condicin x(0) = 1/4 exige de inmediato que c\ = 1/4, en tanto que usando x'(0) =
1/3 se obtiene c^ 1/24. As que, la solucin es
x(t) = - eos 8* + sen 8. (5.5)
Cuando c\ ^ 0 y c< ^ 0 en (5.3), como en el ejemplo 2, la amplitud real A de las
oscilaciones libres no se obtiene en forma inmediata. Para este caso hacemos uso del
siguiente resultado.
Proposicin 5.1.1 (Forma Alternativa de la Solucin)
Si x(t) = Cicosut + C2senu;, con C\ ^ 0 y C2 ^ 0, es conveniente escribir la solucin
x(t) en la forma ms simple
x(t) = Asen(ut + <f)), (5.6)
donde la amplitud A est dada por
A = v/c? + 4 (5.7)
y (f) es un ngulo de fase definido por las relaciones
sen</>=~- y cos(/>=-j. (5.8)
Jx /x
Tomando en cuenta las expresiones (5.8), junto con el hecho de que el rango de la
funcin f(x) = arctan z es el intervalo ( , ), concluimos que el valor de </> puede
calcularse simplemente como sigue
arctan ^ si C2 > 0,
<t> = { (5.9)
arctan ^ + TT si c<i < 0.
EJ EMPLO 3. Reescriba la solucin (5.5) del ejemplo 2 en la forma alternativa (5.6) y
determine el primer valor de t para el cual la masa pasa por la posicin de equilibrio en
direccin hacia abajo.
Solucin. En el ejemplo 2
de modo que, usando (5.7)
x(t) = - eos 8 + sen 8,
A) V24/ 24
Adems c
2
= 1/24 > 0, por lo cual de (5.9), se sigue que
1/4
(f) = arctan - - = arctan 6 = 1.4056 rad.
Por consiguiente
x(t) = ^ sen (8 + 1.4056).
O.
Los instantes en los cuales el cuerpo pasa por la posicin de equilibrio, estn deter-
minados por la condicin
x(t) = 0,
es decir
sen (8 + 1.4056) = 0.
Las soluciones de esta ecuacin vienen dadas por
8t +1.4056 = 7ITT,
con n un nmero entero. Despejando t y recordando que representa una cantidad positiva
(el tiempo), obtenemos la sucesin de valores
Luego, el tiempo requerido es t
2
= 0.6097 segundos (vese la grfica 5.3).
EJEMPLO 4. Una fuerza de 9 Ib estira un resorte 3 pulgadas. Un cuerpo que pesa 24
Ib se sujeta al resorte y se suelta desde un punto que est 3 pulgadas abajo de la posicin
de equilibrio con una velocidad dirigida hacia arriba de 36 pulg/s .
a) Determine la ecuacin del movimiento x(t).
b) En qu instante pasa el cuerpo por la posicin de equilibrio en direccin hacia
arriba por tercera vez?
c) En qu instantes est el cuerpo 3 pulgadas abajo de la posicin de equilibrio?
Solucin. Primero debemos observar que es necesario convertir a ft las longitudes ex-
presadas en pulgadas, usando la equivalencia
1 ft = 12 pulgadas.
a) Por la Ley de Hooke, se sigue que el valor de la constante del resorte k es
9 i, 16
Adems, m 24/32 = 3/4 slug. De modo que la ecuacin diferencial del movimiento es
cPx
dt
2
+ 48x = 0.
En este caso, las condiciones iniciales son
1
x(0) = -, x'(0) = - 3 .
(5.10)
(5.11)
Al resolver el problema de valores iniciales (5.10)-(5.11), obtenemos
i
x(t) = - c os 4>/ 3 -
b) Escribimos la solucin (5.12) en la forma alternativa. Tenemos que
y como C2 < 0
Luego
1 3 __ 1
16
+
16 ~ 2'
/ 1 \ 7T 5
= arctan = +TT = h TI* = -TT.
V v 3 / 6 6
= - s en
2
La grfica de la ecuacin del movimiento se muestra en la figura 5.4.
(5.12)
Los instantes en los cuales el cuerpo pasa por la posicin de equilibrio vienen dados
por las soluciones de la ecuacin x(t) 0, es decir
1 / 5 \
- sen ( 4\ / 3 +- ; 7r ) - 0 .
2 V 6 /
De aqu obtenemos la sucesin de valores de t
717T | ?T
1, 2, o , . . .
El tiempo pedido es claramente
5
= 1.8894 segundos.
c) Ahora, debemos hallar los valores de t para los cuales x(t) 1/4, esto es
sen
Notemos primero que la ecuacin sen# = 1/2 tiene como soluciones todos los nmeros
9 de la forma | + 2nir y |TT + 2ri7r, con n un nmero entero. Luego, el cuerpo est 3
pulgadas abajo de la posicin de equilibrio en los instantes
(i) _ 2_^
+ 2 m r
P ^ = ^ = 0,1,2,3,
Obsrvese que en los tiempos tty el cuerpo se mueve hacia abajo de la posicin de
equilibrio, mientras que en los tiempos ffl lo hace hacia arriba.
EJERCICIOS 5.1
1. Una masa de 1/2 kg est suspendida de un resorte cuya constante es de 18 N/m.
a) Si el cuerpo en reposo se suelta desde un punto que est a 0.1 m abajo de la
posicin de equilibrio, determine la ecuacin del movimiento.
b) Cul es el periodo del movimiento?
2. Una fuerza de 10 N estira un resorte 0.125 m. Despus, al extremo libre de ese
resorte se fija una masa de 5 kg.
a) Encuentre la ecuacin del movimiento si la masa se suelta desde un punto
que est a 0.4 m arriba de la posicin de equilibrio con una velocidad dirigida
hacia abajo de 1.2 m/s.
b) Escriba la ecuacin del movimiento en su forma alternativa.
c) Cuntas oscilaciones completas realiza el cuerpo durante un intervalo de 8TT
segundos?
3. Cuando se sujeta una masa de 100 kg al extremo de un gran resorte, ste se estira
0.98 m. Se quita esta masa y se reemplaza por una de 40 kg, la cual se suelta desde
un punto que est 0.6 m debajo de la posicin de equilibrio, con una velocidad
dirigida hacia arriba de 4 m/s.
a) Determine la ecuacin del movimiento.
c) Obtenga los instantes en los cuales el cuerpo pasa por la posicin de equilibrio.
d) Grafique la ecuacin del movimiento.
4. Un cuerpo de 2 kg se suspende de un resorte de constante 162 N/m.
a) Encuentre la ecuacin del movimiento si la masa se suelta desde un punto a
0.1 m sobre la posicin de equilibrio con una velocidad dirigida hacia arriba
de 1.2 m/s.
c) Grafique la ecuacin del movimiento.
d) Obtenga los instantes en los cuales el cuerpo pasa por la posicin de equilibrio
movindose hacia arriba.
e) En qu posicin se encuentra el cuerpo para t = ?r/8, TT/9, TT/3 ?
f) Calcule la velocidad de la masa para los tiempos del inciso anterior y diga en
que direccin se est moviendo?
5. Al sujetar un peso de 48 Ib a un resorte, ste se alarga 6 pulgadas y luego permanece
en reposo. El cuerpo se desplaza 3 pulgadas por debajo de la posicin de equilibrio
y se suelta.
b) Cul es el periodo del movimiento?
c) Cuntas oscilaciones completas realiza el cuerpo durante un intervalo de 8TT
segundos?
d) Obtenga los instantes en los cuales el cuerpo pasa por la posicin de equilibrio.
6. Suponga ahora que en el ejercicio 5 el peso se suelta desde un punto que se encuentra
3 pulgadas por debajo de la posicin de equilibrio y con una velocidad dirigida hacia
abajo de 4 ft/s .
a) Obtenga la ecuacin del movimiento.
c) Determine los instantes en los cuales el cuerpo pasa por la posicin de equi-
librio.
d) Grafique la ecuacin del movimiento.
7. Encuentre la posicin para la cual un peso sujeto a un movimiento armnico simple
alcanza su velocidad mxima. Cunto tiempo transcurre entre dos mximos o
mnimos consecutivos?
5.2. Movimiento Vibratorio Amortiguado 189
8. Interprete como un movimiento armnico simple los siguientes problemas de valores
iniciales.
ld
2
x
a) - ^ 2 + x = 0; z(0) - 2, z'(0) = - 4
1 d?x
b) ^ + 25x = 0; s(0) = - 0. 1, x(0) = 3
9. Un peso de 25 Ib estira un resorte 6 pulgadas. El resorte est suspendido de un
techo y se encuentra en reposo. Posteriormente el peso se desplaza 4 pulgadas por
debajo de la posicin de equilibrio y se suelta con una velocidad de 2 ft/s, dirigida
hacia arriba.
a) Obtenga la ecuacin del movimiento.
c) En qu instantes el peso se encuentra 5/24 ft abajo de la posicin de equili-
brio?
10. Un cuerpo que pesa 20 libras sujeto al extremo de un resorte lo estira 0.32 ft. El
peso se desplaza 6 pulgadas hacia abajo de la posicin de equilibrio y desde ah se
le comunica una velocidad dirigida hacia arriba de 5 ft/s.
b) En qu instante pasa el cuerpo por la posicin de equilibrio en direccin hacia
arriba por tercera vez? Qu velocidad lleva?
c) En qu instantes est el cuerpo 1/3 ft abajo de la posicin de equilibrio?
d) En qu instantes alcanza el cuerpo sus desplazamientos extremos hacia uno
u otro lado de la posicin de equilibrio?
5.2 Movimiento Vibratorio Amortiguado
En la seccin anterior se supuso que no actan fuerzas retardadoras sobre la masa en
movimiento, lo cual no es cierto a menos que se encuentre suspendida en un vaco perfecto.
Vamos a considerar ahora el efecto de la resistencia del medio sobre la masa. Su-
pongamos que sobre el cuerpo acta una fuerza amortiguadora, dada por un mltiplo
dx
constante de la velocidad .
dt
De la segunda ley de Newton, en ausencia de fuerzas externas, se sigue que
d
2
x . Ax
k 0
donde /? es una constante de amortiguacin positiva y el signo se debe a que la fuerza
amortiguadora acta en direccin opuesta al movimiento. Obtenemos as la ecuacin
diferencial del movimiento vibratorio amortiguado libre
d
2
x (3 dx k
m
o bien
dt
2
m dt
dt
2
con 2A = /3/m y cu
2
k/m.
La ecuacin auxiliar de (5.13) es
(5.13)
r
2
+ 2Ar + u
2
= 0,
cuyas races estn dadas por
r
2
= - A - VA
2
- u
2
n
= _A + VA
2
- u
2
,
Dependiendo del valor de A
2
tu
2
, distinguimos los tres casos siguientes.
(5.14)
CASO I. Movimiento Sobre-Amortiguado. Si A
2
u
2
> 0, las races (5.14) son
reales y distintas, y en consecuencia la solucin general de (5.13) es
~
xt
x
{t) - e
que representa un movimiento suave y no oscilatorio.
Algunas grficas posibles de (5.15) se muestran en la figura 5.5.
(5.15)
Figura 5.5: Movimiento sobreamortiguado
CASO II. Movimiento Crticamente Amortiguado. Si A
2
L
2
general de (5.13) es
x{t) = e~
xt
{c
l
+c
2
t),
puesto que T\ = r
2
= A.
0 la solucin
(5.16)
X
Figura 5.6: Movimiento crticamente amortiguado
En esta situacin, una pequea disminucin de la fuerza de amortiguamiento pro-
ducira un movimiento oscilatorio. Algunas grficas posibles de la solucin (5.16) se
muestran en la figura 5.6.
Un examen de las derivadas de las soluciones (5.15) y (5.16), de los casos I y II
respectivamente, permite ver que estas funciones pueden tener a lo ms un mximo
relativo o un mnimo relativo para t > 0, por lo que el cuerpo puede pasar a lo ms una
vez por la posicin de equilibrio.
CASO III. Movimiento Subamortiguado. Si A
2
cu
2
< 0, las races (5.14) son
complejas y se pueden escribir como
n = -A + VC
2
- A
2
, r
2
= -A - \ / o;
2
- A
2
,
de modo que la solucin general de (5.13) es en este caso
x(t) = e"
A
(ci eos v ^
1
A
2
1 + c
2
sen Va;
2
- X
2
t). (5.17)
Ahora, el movimiento es oscilatorio, pero la amplitud de las oscilaciones tiende a cero
cuando t tiende a infinito.
Ntese que, en analoga a lo que hicimos en el movimiento armnico simple, la solucin
(5.17) puede expresarse en forma compacta, de acuerdo a la siguiente proposicin.
Proposicin 5.2.1 (Forma Alternativa de la Ecuacin del Movimiento Sub-
amortiguado) Cualquier funcin de la forma
x(t) = e"
A
(ci eos VUJ
2
- X
2
t + c
2
sen Veo
2
- A
2
),
con C\ T 0 y c
2
^ 0, puede escribirse como
donde
x(t) = Ae~
xt
sen {Vu;
2
-X
2
t + 0),
A= Jc\ + 4
(5.18)
(5.19)
y el ngulo de fase </ > es tal que
Es decir
sen 0 =
arctan
y
arctan f- +
, c
2
eos (p .
si c
2
> 0,
si c
2
< 0.
(5.20)
En la forma alternativa (5.18) el coeficiente Ae
xt
se denomina la amplitud amor-
tiguada de las soluciones, 2TT/^/UJ
2
A
2
es el cuasiperiodo y y/u
2
A
2
/2TT es la cuasifre-
cuencia. El cuasiperiodo es el intervalo de tiempo transcurrido entre dos mximos suce-
sivos de x(), tambin es igual al doble de tiempo entre dos ceros sucesivos de la solucin.
Figura 5.7: Movimiento subamortiguado
Para representar grficamente la solucin (5.18), es til tomar en cuenta las siguientes
observaciones. Ver figura 5.7. En primer lugar, las intersecciones con el eje t se obtienen
resolviendo la ecuacin x(t) = 0, esto es
- 0,
de donde
X
2
t + () =
727T,
T7T
n = 0,1,2,
Por otra parte, la grfica de x(t) es tangente a las curvas exponenciales x\(t) =
Ae~
xt
,X2{t) Ae~
xt
en los valores de t, tales que
( v ^
2
- A
2
+ 0) = 1.
Resolviendo esta ecuacin encontramos las soluciones t = t%, dadas por
_ (2fc + 1)TT/2 - 4>
con k en N.
EJEMPLO 1. Se encontr experimentalmente que un cuerpo de 4 Ib estira un resorte 6
pulgadas. El medio ofrece una resistencia al movimiento del cuerpo numricamente igual
a 2.5 veces la velocidad instantnea. Encuentre la ecuacin del movimiento si el peso se
desplaza 4 pulgadas por debajo de la posicin de equilibrio y se suelta.
Solucin. La ecuacin diferencial del movimiento es
d
2
x_ dx
32dP " " * *

^Tt
o equivalentemente
^ + 2 0 ^ + 64x = 0. (5.21)
dt
dt
Las condiciones iniciales son
x(0) = ^, x'(0) = 0.
La ecuacin auxiliar de (5.21) es r
2
+ 20r + 64 = 0 y sus races son r\ = 4,r
2
= 16,
de modo que
La condicin x(0) = 1/3 implica que
ci +c
2
= i , (5.22)
en tanto que x'(0) = 0 nos lleva a
-4ci - 16c
2
= 0. (5.23)
Resolviendo el sistema (5.22)-(5.23) obtenemos los valores
4 1
Cl =
9'
2 =
"9"
0. 2 0. 4 0. 6 0. 8 1 1. 2
Por consiguiente
*<*> - 9
e
" " 9
e
~
W
-
( 5
'
2 4 )
La grfica de la solucin (5.24) se da en la figura 5.8. Como se observa no ocurren
oscilaciones ya que el peso tiene tanto amortiguamiento que slo retorna gradualmente a
la posicin de equilibrio sin pasar por esta. Se trata de un movimiento sobreamortiguado.
EJEMPLO 2. Resuelva nuevamente el ejemplo 1, suponiendo ahora que /? = 2.
Solucin. En este caso la ecuacin diferencial del movimiento es
4 d
2
x dx
32d2~
=
~
x
~~
o bien
d
2
x dx
_
+ 1 6
_ = 0. (5.25)
La ecuacin caracterstica de (5.25) es r
2
+ 16r + 64 = 0, cuyas races son r\ = r
2
= 8,
por lo cual
x
(t) =de~
8t
+ c
2
te~
st
.
Las condiciones iniciales x(0) 1/3, x'(0) = 0 conducen al sistema de ecuaciones lineales
1
Cl =
3'
-8ci +c
2
= 0.
Por consiguiente
xt) = - e"
8
+ -6"
8
= ~e"
8
(l + 8). (5.26)

La grfica de la solucin (5.26) se muestra con lnea continua en la figura 5.9, donde
con lnea interrumpida aparece tambin la solucin del ejemplo 1, a fin de hacer una
comparacin. Aunque ahora es un movimiento crticamente amortiguado, vemos que las
grficas son muy similares.
1/3
0.3
0. 2
0. 2 0. 4 0. 6 0.
Figura 5.9: Soluciones de los ejemplos 1 y 2
Si modificamos las condiciones iniciales el tipo de movimiento se conserva, pero algu-
nas de sus caractersticas si pueden cambiar.
EJEMPLO 3. Si en el ejemplo 2 se cambian las condiciones iniciales a x(0) = 0 y
z'(0) = 1/3, determine x(t).
~
st
Solucin. La solucin general sigue siendo
x(t) = cie"
8
+ c
2
te
Sin embargo, las condiciones iniciales nos conducen a los valores de c\ = 0 y c
2
= 1/3.
En consecuencia
x(t) = i e"
8
. (5.27)
Se tiene que
A*) = ^ - 8e-
8t
) = i e-
8 t
( l - 8),
de dnde se observa que x(t) alcanza un desplazamiento mximo en t = 1/8 segundos
igual a x
max
= x(l/8) = 0.01533 ft. Ver figura 5.10.
EJEMPLO 4. Tomando en cuenta que (5 = 1 repita el ejemplo 1. Escriba la solucin
en la forma alternativa (5.18). Determine los instantes en los que el cuerpo pasa por la
posicin de equilibrio y realice la grfica de la ecuacin del movimiento.
4 d
2
x dx
32dt* * ~ di '
0. 0025
esto es
0. 125 0. 25 0. 375 0. 5 0. 625
d
2
x dx
+ 8
_
+ 6 4 l = 0
.
dt
2
dt
Ya que las races caractersticas son r ^ = 4 4\/3i, su solucin general est dada por
x(t) =e"
4
(
De las condiciones iniciales x(0) = 1/3 y x'(0) = 0 se sigue que
-4ci
1
" 3'
= 0.
De modo que c< = \/3/9 y
x
(t) = i
e
-
4
(3cos4\/3 + V3 sen4V3t). (5.28)
Ahora bien, usando las expresiones (5.19) y (5.20) obtenemos
1/3 7T
i
iH^-S*
= arctan
V3/9 3"
Luego, escribimos la solucin (5.28) en la forma
2
z() = ^V
/
3e-
4t
sen(4v
/
3+0. (5.29)
Por otra parte, las soluciones de la ecuacin x(t) = 0 se encuentran directamente de
(5.29) y estn dadas por
4\ / 3* + 77 = n7r, n N.
As que, el cuerpo pasa por la posicin de equilibrio en los instantes
- -
3
-
1
) * . n = l,2,...
Finalmente, la grfica de x(t) es tangente a las curvas exponenciales x\(t) 2>/3e
4
/9
y %2(t) ~ 2\/3e~
4
/9 en los valores de t dados por la condicin
sen
es decir
o bien
(Gn + lO
2 4 ^
n<EN,
1 1 / \ J X , J <> . . .
La grfica de (5.29) se muestra en la figura 5.11.
x
x
( t )
( t )
EJ EMPLO 5. Despus de que un cuerpo que pesa 10 Ib se sujeta a un resorte de 5 ft
de largo, el resorte mide 7 ft. Se quita el cuerpo de 10 Ib y se le reemplaza por uno de 8
Ib. El sistema completo se coloca en un medio que ofrece una resistencia numricamente
igual a la velocidad instantnea.
a) Obtenga la ecuacin del movimiento si el peso se suelta desde un punto que se
encuentra 1/2 ft abajo de la posicin de equilibrio con una velocidad dirigida hacia
abajo de 1 ft/s.
b) Encuentre los instantes en los cuales el cuerpo pasa por la posicin de equilibrio
en direccin hacia abajo.
c) Grafique la ecuacin del movimiento.
Solucin.
a) De la ley de Hooke es claro que
10 Ib _ Ib
*~T I" ~
5
a"
Adems
8 Ib 1
m=
32ft 7?
=
4
y P 1, as que la ecuacin diferencial del movimiento es
d x dx
-^ + 4 + 20x = 0. (5.30)
dt
dt
Resolviendo (5.30) sujeta a las condiciones iniciales x(0) = 1/2 ft y x
f
(0) = 1 ft/s,
obtenemos
x{t) = ~e"
2
(cos4 + sen4). (5.31)
b) Escibimos primero la solucin (5.31) en la forma alternativa. Tenemos que
2
2,
= arctan 1 = .
4
De modo que
por lo cual x(t) = 0 si y slo si
7T
4 H = nn, n N.
4
Por lo tanto, los valores
7T 7T
t
n
= n- - ,
4 16
con n un entero positivo par, son los instantes en los que el cuerpo pasa por la posicin
de equilibrio movindose hacia abajo.
c)La grfica se muestra en la figura 5.12.
- 0 . 3
EJERCICIOS 5.2
1. Un peso de 2 Ib est sujeto a un resorte el cual tiene una constante de elasticidad
de 4 lb/ft. El peso se suelta desde un punto que se encuentra 6 pulgadas abajo
de la posicin de equilibrio con una velocidad dirigida hacia abajo de 2 ft/s, en
un medio que presenta una resistencia al movimiento numricamente igual a la
velocidad instantnea. Determine:
a) La ecuacin del movimiento.
b) Los instantes en los cuales el cuerpo pasa por la posicin de equilibrio.
c) El desplazamiento extremo del peso.
d) La grfica de la ecuacin del movimiento.
2. Supngase que de un resorte se suspende una masa de 1 slug. Para los valores
de las constantes y condiciones iniciales dadas a continuacin, determine en cada
inciso: la ecuacin del movimiento en su forma alternativa, los instantes cuando el
cuerpo pasa por la posicin de equilibrio en direccin hacia arriba y la grfica del
movimiento.
a) k = 3 lb/ft, 0 = 6; x(0) = 4 pulg, x'(0) = - 2 ft/s.
b) k = 2 lb/ft, (3 = 2] x(0) = - 3 pulg, z'(0) = - 1 ft/s.
3. Un peso de 16 Ib estira un resorte 4 ft, el sistema completo se sumerge en un medio
viscoso que opone una fuerza de amortiguacin numricamente igual a tres veces
la velocidad instantnea. Determine:
a) La ecuacin del movimiento si el peso se suelta desde un punto que est 1 ft
arriba de la posicin de equilibrio con una velocidad dirigida hacia abajo de 7
ft/s.
b) El instante en que cruza por la posicin de equilibrio.
c) La grfica de la ecuacin del movimiento.
4. Un peso de 8 Ib estira un resorte 2 ft y el sistema se encuentra en un medio que ofrece
una resistencia numricamente igual a /? veces la velocidad instantnea, donde p es
una constante positiva. Determinar los valoes de fi para que el movimiento sea:
a) Sobreamortiguado.
b) Crticamente amortiguado.
c) Subamortiguado.
5. Una masa de 1/2 slug estira un resorte 4 ft y el medio que rodea al sistema masa-
resorte ofrece una resistencia al movimiento numricamente igual a 4.5 veces la
velocidad instantnea. El peso se suelta 6 pulgadas abajo de la posicin de equilibrio
con una velocidad dirigida hacia arriba de VQ ft/s. Cmo debe ser vo para que la
masa pase por la posicin de equilibrio ?
6. Dar una posible interpretacin fsica del siguiente problema de valores iniciales
4 d
2
dx 1
+ 2 +
x(0) = - 1 , 2/(0) - - 1
7. Indique si las siguientes grficas corresponden a un movimiento amortiguado, en
cuyo caso clasifquelo.
(b)
5.3. Movimiento Vibratorio Forzado 201
5-3 Movimiento Vibratorio Forzado
En las dos secciones anteriores estudiamos el problema de un resorte donde slo se con-
sideraron las fuerzas restauradora y amortiguadora. Veremos ahora casos dnde actan
otras fuerzas externas que varan con el tiempo. Dichas fuerzas pueden ocurrir, por ejem-
plo, cuando el soporte que sostiene al resorte se mueve verticalmente de cierta manera
dada, tal como en un movimiento peridico o cuando al peso se le da un pequeo empuje
cada vez que alcanza la posicin ms baja.
Denotemos con f(t) la fuerza exterior que acta sobre la masa. De la segunda ley de
Newton, la ecuacin diferencial del movimiento es
^ (5-32)
o bien
^ + 2A^+a;
2
x = F(), (5.33)
donde 2A = (3/m,uj
2
= k/m y F(t) = f(t)/m.
Para resolver la ecuacin no homognea (5.33) podemos emplear el mtodo de los
coeficientes indeterminados o el de variacin de parmetros, segn sea ms conveniente.
EJEMPLO 1. Un resorte vertical con constante de 6 lb/ft tiene suspendida una masa
de 1/2 slug. Se aplica una fuerza externa dada por f(t) = 40sen2, > 0. Supngase que
acta una fuerza amortiguadora numricamente igual a dos veces la velocidad instantnea
y que inicialmente el cuerpo est en reposo en su posicin de equilibrio. Determine la
posicin del cuerpo en cualquier tiempo t > 0.
Solucin. Con los valores de k = 6 lb/ft, m = 1/2 slug y (3 = 2, la ecuacin diferencial
de movimiento resultante es
^
+ 4
^
+
l
2 x =
80sen2. (5.34)
dt
dt
La solucin complementaria de (5.34) es
x
c
(t) = e'
2
\c
x
cos2\2t + c
2
sen2v
/
2).
Usando el mtodo de los coeficientes indeterminados proponemos una solucin particular
de (5.34) de la forma
x
p
(t) = A eos 2 + B sen 2.
En tal caso
x'
p
(t) = -
x'(t) = -4Acos2t-4Bsen2t.
Sustituyendo en (5.34), se sigue que
(8 A + 8B) eos 2 + (8B - 8A) sen 2 = 80 sen 2t.
El sistema de ecuaciones resultante
8.4 + 85 = 0,
-8,4 + 85 = 80,
conduce a los valores A = 5 y 5 = 5. As que
x(t) = e~
2t
(c
1
cos2\/2 + C2sen2v
/
2) + 5(sen2 - cos2).
Empleando las condiciones iniciales x(0) = 0 y x
f
(0) = 0 encontramos que
c
2
= 0. Por lo tanto
x(t) = 5e~
= 5y
Obsrvese que en el ejemplo anterior la solucin complementaria
x
c
(t) = 5e~
2
co
tiene la propiedad de que
lim x
c
(i) = 0,
-+OO
por lo cual se dice que x
c
(t) es un trmino transitorio o una solucin transitoria.
As para valores grandes de , x(t) se aproxima a x
p
(t). A x
p
(t) se le llama solucin
estacionaria o de estado permanente. Ver figura 5.13.
Solucin estacionara
Solucin completa
De hecho, si en la ecuacin diferencial (5.33) ponemos F(t) = Fosenat o F(t) =
F
o
eos at, donde F
o
y a son constantes, entonces su solucin general consiste en la suma
de dos trminos: trmino transitorio ms trmino estacionario.
5.3.1 Resonancia
Estudiaremos la ecuacin (5.33) en el caso especial en que F(t) = Fo sen ai , i > 0, donde
F
o
y a son constantes positivas. La ecuacin diferencial bsica es
d x dx
-Tp + 2A + u
2
x = F
o
sen ai , (5.35)
donde 2A = (3/m y u
2
= k/m.
Supondremos que (3 es suficientemente pequeo de modo que el amortiguamiento es
menor que el crtico. En otras palabras, consideraremos que A < u. Luego, la solucin
complementaria de (5.35) tiene la forma
x
c
(t) = e
xt
(c
1
eos yjuj
2
- XH + c
2
sen y/u
2
- A
2
),
con C\ y C2 constantes arbitrarias, que dependen de las condiciones iniciales, o equivalen-
temente
x
c
(t) = Ae~
xt
sen {y/u
2
- XH + 0),
donde A ye
2
+ c
2
, sen<f> = c\/A y eos4> =
Ahora determinaremos una solucin particular de (5.35), utilizando el mtodo de los
coeficientes indeterminados. Sea
x
p
(t) = B eos at + C sen at.
Entonces
x'
p
(t) = aB sen at + aC eos ai ,
Zp(t) = - a
2
B eos at - a
2
C sen ai .
Sustituyendo x
p
(t),x'
p
(t) y x^() en (5.35) se obtiene
(-a
2
B + 2XaC + u
2
B) eos ai + [{u
2
- a
2
)C - 2aAJ5] sen at = F
o
sen ai .
Igualando los coeficientes en la ltima igualdad, resulta el sistema de ecuaciones
(u
2
- a
2
)B + 2XaC = 0,
-2aAB + (u
2
- a
2
)C = F
o
,
cuya solucin es
2aAF
0
B =
(o;
2
- a
2
)F
0
Consecuentemente
2a\F
0
(
2
- a
2
)F
0
Xp{t) = - cosa+ sen ai.
(LJ
2
a
2
)
2
+ 4a
2
X
2
(u
2
a
2
)
2
+ 4a
2
A
2
Podemos escribir x
p
(t) en la forma
x
p
(t) = Asen (at + 9),
donde
A =
2a\F
0
( u;
2
- a
2
)
2
+ 4a
2
A
2
_
2
(o;
2
-a
2
)F
0
]
2
(o;
2
-a
2
)
2
+ 4a
2
A
2
J
N
es decir, simplificando
* ,
F
^(u
2
- a
2
)
2
+ 4a
2
A
2
El ngulo 6 est determinado por las relaciones
2aA
sen 6 = T=
^ 2 _
Q
2)2
+ 4 a
2
A
2 ' ^(^2 _
Q
2)2
+ 4 Q
2
A
2
As que
p
xJt) = , = sen (at + 6).
f/ 2
- a
2
)
2
+ 4a
2
A
2
Obsrvese que la solucin completa es la suma de dos trminos
El primero
x
c
(t) = Ae~
xt
sen (Vu
2
- XH + 0),
representa la oscilacin amortiguada, que sera todo el movimieno del sistema si la fuerza
externa F(t) no actuara. El segundo trmino
xJt) = . sen (ai + 0),
V
/
( a;
2
- a
2
)
2
+ 4a
2
A
2
que resulta de la presencia de la fuerza externa, representa un movimiento armnico
simple de periodo 2n/a y amplitud
- a
2
)
2
+ 4a
2
A
2
Para F
0
,C y A fijos, la amplitud es funcin de a. Consideremos la funcin g(a) en el
intervalo (0,oo). Se tiene que
. 2F
0
a(
2
- a
2
- 2A
2
)
Luego, g'(a) = 0 si y slo si a = a
0
= 0 o a = ai = \ / ^
2
- 2A
2
. Se puede verificar
fcilmente que la amplitud de las oscilaciones alcanza un valor mximo cuando
a = ai - Va;
2
- 2A
2
.
El valor mximo de la amplitud es
2 2
-
J
Definicin 5.3.1 Cuando la frecuencia de la fuerza exterior es

ai _ VC
2
- 2A
2
2TT ~ 2^ '
5e dzce que el sistema est en resonancia.
En un sistema en resonancia (a = ai ), la amplitud de la oscilacin vara inversamente
con la constante de amortiguamiento. De hecho, se observa que
lim g(ai) lim .
n
= oo
l i m a
x
= C.
Diremos que hay resonancia pura si /? = 0. En tal caso a = ai = a;, o sea que la
frecuencia de la fuerza externa es igual a la frecuencia natural del sistema.
Finalmente obsrvese que la resonancia puede ocurrir solamente si
u
2
> 2A
2
,
k
2
- >
2
m
esto es
p < V2km. (5.36)
EJ EMPLO 2. Un peso de 4 Ib se suspende de un resorte cuya constante es de k 8
lb/ft. Suponga que una fuerza externa dada por f(t) 4cos8 se aplica al resorte y que
no hay amortiguamiento. Describa el movimiento que resulta si se asume que inicialmente
el peso est en la posicin de equilibrio y que su velocidad inicial es cero.
Equivalentemente
d?x
- + 64x = 32cos8. (5.37)
dt
2 v }
La solucin complementaria de (5.37) es
x
c
(t) = c\ eos 8 + C2 sen 8.
x
p
(t) = (A eos 8 + B sen 8).
Sustituyendo en (5.37) se encuentra que A = 0 y J5 = 2. As que
x () = ci eos 8 + c
2
sen 8 + 2 sen 8.
De las condiciones iniciales x(0) = 0 y x'(0) = 0 encontramos que c\ = 0 y c
2
= 0. Por
consiguiente la ecuacin del movimiento es
x(t) = 2tsen8t.
Su grfica se muestra en la figura 5.14.
Se observa que en este caso hay resonancia pura en vista de que (3 = 0 y la frecuencia
de la fuerza externa aplicada es igual a la frecuencia natural del sistema no amortiguado.
x=2t
x=-2t
EJERCICIOS 5.3
1. Una masa de 1 slug se encuentra suspendida de un resorte de constante de elasti-
cidad igual a 4 lb/ft y el sistema est inmerso en un medio que opone una fuerza
de resistencia numricamente igual a 5 veces la velocidad instantnea. Si la masa
se suelta 6 pulgadas arriba de la posicin de equilibrio con una velocidad dirigida
hacia abajo de 4 ft/s, obtenga:
a) La ecuacin del movimiento, si acta una fuerza externa sobre la masa dada
por f(t) = 20cos2 + 10sen2
b) Las grficas de la solucin transitoria y de la solucin estacionaria utilizando
el mismo sistema de ejes coordenados.
2. Un peso de 32 Ib se sujeta a un resorte de constante de elasticidad igual a 5 lb/ft. El
peso y el resorte se sumergen en un medio que ofrece una resistencia numricamente
igual a 6 veces la velocidad instantnea. El movimiento se inicia en un punto que
se encuentra a 4 pulgadas abajo de la posicin de equilibrio y partiendo del reposo.
Determine:
a) La ecuacin del movimiento si sobre el peso se aplica una fuerza externa igual
a f(t) = e-*.
b) La grfica de la ecuacin del movimiento.
3. Un resorte tiene una constante de elasticidad igual a 1 lb/ft. Un peso de 8 Ib se
suspende de un extremo del resorte y el sistema se coloca en un medio que ofrece
una resistencia numricamente igual a la velocidad instantnea. Si el peso se suelta
en reposo, 4 pulgadas sobre la posicin de equilibrio y sobre l acta una fuerza
externa /() = 25sen4, obtenga la ecuacin del movimiento y su grfica.
4. Un peso de 3.2 Ib estira un resorte 6.4 ft. Si el peso se suelta 3 pulgadas abajo
de la posicin de equilibrio con una velocidad dirigida hacia abajo de 6 ft/s y el
medio en que est el sistema masa-resorte ofrece una fuerza de amortiguamiento
numricamente igual a la quinta parte de la velocidad instantnea, determine la
ecuacin del movimiento si adems se aplica al peso una fuerza externa dada por
f(t) = e~
cos2. Grafique la solucin obtenida.

5. Resuelva el ejercicio 1 en ausencia de la fuerza de resistencia.
6. Resuelva el ejercicio 4 en ausencia de la fuerza de resistencia.
7. Un resorte sujeto a un soporte tiene suspendida una masa de 2 kg y la constante
de elasticidad del resorte es de 4 N/m. El sistema est en reposo cuando el soporte
empieza a oscilar de acuerdo a la expresin h(t) = 2cos3. Determine:
a) La ecuacin diferencial del movimiento si el sistema completo est inmerso en
un medio que opone una fuerza de resistencia numricamente igual a 6 veces
la velocidad instantnea.
b) La ecuacin del movimiento (tome en cuenta que el peso est en reposo en la
posicin de equilibrio cuando el soporte empieza a oscilar).
8. Resuelva el ejercicio 7 en ausencia de amortiguamiento.
5.4 Circuito LRC en Serie
Ahora aplicaremos la teora antes vista para determinar la carga q(t) y la corriente i(t)
en un circuito como el mostrado en la figura 5.15, en el que se conectan un inductor o
bobina de L henrys, una resistencia de R ohms, un condensador o capacitor de C farads
y un generador de voltaje cuya fuerza electromotriz est dada por una funcin E(t) volts.
De la segunda ley de Kirchhoff se tiene
di
dt
dq
dt'
1
di
~dt
T
(5.38)
dt C
Ya que
d
2
q
5.4. Circuito LRC en Serie 209
E(t) ( 1
Figura 5.15: Circuito LRC
sustituyendo en (5.38) resulta la ecuacin diferencial para la carga elctrica en el con-
densador
Not tambin que si primero derivamos con respecto a t en (5.38) obtenemos
d
2
i di 1 dq dE
d^
+
dt
+
C~dt
=Z
~dt'
y si luego sustituimos las expresiones (5.39), esto nos conduce a la ecuacin diferencial
de la corriente elctrica
T
d
2
i di 1 . dE
L + R +
Cabe adems destacar la similitud entre las ecuaciones (5.40) y (5.32), lo cual permite
resolver un problema de movimiento vibratorio en base al anlisis del correspondiente
circuito elctrico y viceversa, identificando
la carga q con la posicin x,
la inductancia L con la masa m,
la resistencia R con la constante de amortiguamiento /?,
el recproco de la capacitancia \/C con la constante del resorte k,
la fuerza electromotriz E(t) con la fuerza externa /(.) y
la corriente elctrica i = dq/dt con la velocidad v = dx/dt.
Es claro entonces que podemos aplicar todos los resultados de la seccin anterior al
estudio de un circuito LRC en serie.
EJ EMPLO 1. Un circuito en serie consta de un inductor de 0.25 H, una resistencia de
40 , un capacitor de 4 x 10~
4
F y una fuerza electromotriz dada por E(t) = 5 sen 100
V. Si la corriente inicial y la carga inicial en el capacitor son ambas cero, determine la
carga en el capacitor y la corriente elctrica del circuito para cualquier tiempo t > 0.
Solucin. Sustituyendo los valores de L = 0.25 H, R = 40 Q, C = 4 x 10"
4
F y
E(t) 5 sen 100 V en la ecuacin diferencial (5.40) obtenemos
o bien
4 4 + 160^ + 10000? = 20 sen 100. (5.42)
dt
dt
La ecuacin auxiliar de (5.42) es r
2
+ 160r + 10000 = 0, cuyas races son r
x
= - 80 4- 60
y
r
2
= - 80 - 60. Luego
q
c
(t) = e~
80t
(ci eos 60 + c
2
sen 60).
Adicionalmente, empleando el mtodo de coeficienes indeterminados encontramos que
una solucin particular de (5.42) es
En consecuencia, la solucin general de (5.42) es
q(t) = e"
80
(c! eos 60 4- c
2
sen 60) - eos 100.
800
De las condiciones iniciales q(0) = 0 y q'(0) = 0 se sigue que
Cl
-
85
-80ci + 60c
2
= 0,
respectivamente. A partir de estas ecuaciones encontramos que
1 1
V 1
800 " 600
Por consiguiente, la carga en el capacitor es
5.4. Circuito LRC en Serie 211
y la corriente elctrica viene dada por
5 1
i(t) = -e~
80
sen60 + -senlOO.
EJERCICIOS 5.4
1. Un inductor de 1 H, una resistencia de 2 2, un condensador de 0.2 F y un generador
con una fuerza electromotriz dada por E(t) = 35 Volts se conectan en serie. Si
la corriente inicial es cero y la carga inicial en el condensador es de 1 Coulomb,
determine la carga y la corriente para todo tiempo t > 0.
2. Se conecta un circuito en serie con un inductor de 0.5 H, una resistencia de 6 7,
un condensador de 0.02 F y una fuente de voltaje alterno dado por 24senl0.
Determine la carga y la corriente al tiempo , si la carga en el condensador y la
corriente en el circuito son cero al tiempo t ~ 0.
3. Un inductor de 4 H, una resistencia de 20 , un capacitor de 0.008 F y un generador
con una fuerza electromotriz dada por E(t) = 500 Volts se conectan en serie. Si
inicialmente la carga y la corriente son ambas cero, obtenga.
a) La carga y la corriente para todo tiempo.
b) La carga y la corriente despus de un tiempo largo.
4. Se conectan en serie un inductor de 1 H, una resistencia de 2 , un capacitor de 0.5
F y una fuente de voltaje alterno dado por E(t) = 20cos2 V. Si la carga inicial
almacenada en el capacitor es de 1 Colulomb y la corriente inicial es igual a cero
Amper, encuentre:
a) La carga que contiene el capacitor en el tiempo t > 0.
b) La corriente de estado estacionario y exprsela en la forma alternativa.
5. Un inductor de 0.4 H, un condensador de 0.001 F y un generador con una fuerza
electromotriz de 20 V se conectan en serie. Si en t = 0 la carga y la corriente son
cero, determine:
a) La carga y la corriente para todo tiempo.
b) Los valores mximos de la carga y la corriente.
6. Un circuito en serie contiene un inductor de 1 H,un capacitor de 10~
4
F y una
tensin E(t) = 100sen50 Volts. Inicialmente la carga y la corriente son nulas,
encuentre:
a) La carga y la corriente en un instante cualquiera.
b) Los instantes en los cuales la carga del capacitor es cero.
7. Un circuito en serie contiene un inductor de 5/3 H,una resistencia de 10 fi, un
condensador de 1/30 F y una bateria de 300 Volts. Si en t = 0 la carga y la
corriente son cero, determine:
a) La carga y la corriente en un instante cualquiera.
b) La carga mxima en el condensador.
c) Manteniendo fijos los valores de la inductancia y la capacitancia (ignorando la
bateria) para qu valores de la resistencia el circuito llega a estar: sobreamor-
tiguado, crticamente amortiguado o subamortiguado?
8. Suponga que en un circuito LRC en serie L = 1 H y C 1/3 F.
a) Determine los valores de R para los cuales el circuito est subamortiguado.
b) Para qu valores de R se puede producir resonancia?
5.5 Otras Aplicaciones
Existen una gran cantidad de aplicaciones de las ecuaciones diferenciales ordinarias de
segundo orden. En esta seccin presentaremos dos ejemplos ms para ilustrar las posi-
bilidades.
5.5.1 Vigas Horizontales
El problema consiste en determinar la flexin de una viga rectangular sometida a una
carga. Inicialmente la viga es recta y su eje central coincide con el eje x, como se muestra
en la figura 5.16. Posteriormente, dicho eje se ha desplazado debido a la accin de la
carga (ver figura 5.17). Lo que se desea es obtener la ecuacin de la curva punteada,
llamada curva elstica, que nos da la deformacin de la viga.
Por simplicidad consideraremos la curva elstica y un punto P(x, y) sobre ella. De
los cursos de Fsica se sabe que el momento M en el punto P es la suma algebraica de
los momentos de las fuerzas externas que actan sobre el segmento de la curva. Aqu
supondremos que las fuerzas hacia arriba dan momentos positivos y las fuerzas hacia
abajo dan momentos negativos. El momento est dado por
donde E es el mdulo de elasticidad de la viga e / es el momento de inercia. Luego, si
queremos conocer la ecuacin de la curva elstica debemos resolver la ecuacin diferencial.
%
= M
-
( 5
-
4 3 )
5.5. Otras Aplicaciones
YA
213
y
Figura 5.16: Viga horizontal
Figura 5.17: Aplicacin de una carga a una viga
Veamos un caso concreto.
EJEMPLO. Una viga de 8 m de longitud est apoyada en dos columnas verticales. Si
la viga tiene una carga uniforme de 500 kg por metro de longitud y una carga al centro
de 5000 kg, cul es la curva elstica de la viga?
Solucin. En la figura 5.18, las fuerzas que actan sobre OP son
1) Una fuerza aplicada e n O a x metros de P, dirigida hacia arriba e igual a la
carga total, es decir -(5000 + 8 500).
2) Una fuerza de 500x dirigida hacia abajo que se supone concentrada en el punto
medio de OP.
Figura 5.18: Viga del ejemplo
Q
As el momento flexionante (flector) en P es.
M = Fidi - F
2
d
2
= \ (5000 + 8 500)x - 500x (?)
= 4500z-250z
2
,
y la ecuacin diferencial (5.43), en este caso, tiene la forma
= 4500z - 250x
2
.
dx
2
(5.44)
Podemos resolver (5.44) integrando directamente. Integrando una vez resulta
dx
y volviendo a integrar obtenemos
2225
3
125
Ely = x

c
x
x + c
2
.
En O, x = y = 0 de modo que c
2
= 0. En Q, x = 8, y = 0, por lo cual c\ = 36800. Por
lo tanto
es la curva elstica de la viga.
6
5.5. Otras Aplicaciones 215
5.5.2 El Pndulo Simple
Un pndulo simple consiste en una partcula de masa m suspendida de una cuerda (o un
hilo inelstico) de largo I y de masa despreciable. Suponiendo que la cuerda est siempre
tensa, que las oscilaciones son en un plano vertical y que las nicas fuerzas que actan
son el peso de la partcula y la tensin en la cuerda, deseamos hallar la ecuacin del
movimiento.
mg
Figura 5.19: Pndulo Simple
Sean 9 y s como en la figura 5.19. Se tiene que s = 19, de donde
Descomponiendo el peso mg en dos componentes, una en la direccin de la tangente a
la trayectoria y la otra perpendicular a sta, vemos que la componente perpendicular
se compensa por la tensin. La magnitud de la componente tangencial es mg sen 9. Ver
figura 5.20.
Luego, de la segunda ley de Newton se sigue que
,<P0
ma = ^"772"
=
mg sen 9.
Obtenemos as la ecuacin diferencial
o equivalentemente
d
2
9
(5.45)
Figura 5.20: Pndulo Simple. Componentes del peso.
La ecuacin (5.45) no es lineal y no puede resolverse en trminos de funciones ele-
mentales. Sin embargo, para ngulos pequeos sen# 9; que sustituyendo en (5.45) nos
conduce a la ecuacin diferencial lineal
(5.46)
La solucin general de (5.46) es claramente
6(t) = ci eos J-t + c
2
sen J-t,
V V
que corresponde a un movimiento armnico simple con perodo
T =
217
Conclusiones
El material que hemos expuesto permitir que el lector resuelva algunas ecuaciones
diferenciales ordinarias. No hemos pretendido presentar todas las tcnicas de solucin
conocidas o posibles. En lugar de ello, animamos ahora al estudiante a explorar las
opciones que ofrecen las diversas calculadoras o paquetes computacionales existentes,
tiles para este propsito.
Por supuesto, el usuario de tales recursos siempre deber proceder con sumo cuidado,
desde el almacenamiento de los datos hasta la interpretacin de la solucin que obtenga.
Por ejemplo, para resolver ecuaciones diferenciales ordinarias de primer orden con la
calculadora TI-92 la instruccin general es
deSolve(ecuaein
)
variable independiente,variable dependiente),
la cual conduce a la solucin general.
EJEMPLO 1. Resolver
d,x
Solucin. En este caso debemos teclear:
deSolve(x*y * 2y '=y * 3-x * 3,x,y),
lo que da por resultado la expresin
donde @9 denota una constante. En general con @n, n = 1,2,... ,25, la calculadora
indica constantes.
Tambin pueden resolverse problemas de valor inicial empleando la calculadora. Ver
pginas 120, 121 en [14]. Asimismo, la instruccin para resolver ecuaciones diferenciales
de orden dos es anloga.
218 Conclusiones
EJEMPLO 2. Resolver
y" + 4j/ + Ay = 10x
3
e~
2x
.
Solucin. Tecleamos ahora:
deSolve(y"+4y'+4y=10x~3 e' (-2x),x,y),
obtenindose
x
5
y
= [+@ 10 x+@ll]e-
2x
Obsrvese que (@ 10 x+@ll)e~
2
'
x
es la solucin complementaria.
Para este ejemplo es claro que la calculadora es una gran herramienta que ahorra
tiempo y esfuerzo en el trabajo necesario para calcular la solucin. Hasta el momento
slo se pueden resolver mediante ella ecuaciones diferenciales de orden 1 y 2.
Con el paquete computacional Mathematica, la solucin de las ecuaciones de los ejem-
plos anteriores se obtiene mediante las instrucciones:
DSolve[x*(y[x]) *2*y'[x]==(y[z]) *3-x*3,yfx],x]
DSolve[y"[x]+4 *y '[x]+4 *y[x]==10*x" 3*Exp[-2x],y[x],x],
respectivamente.
Las posibilidades al utilizar Mathematica son mayores, ya que adems pueden resol-
verse ecuaciones diferenciales y problemas con valores iniciales de orden superior. Vase
[13].
En relacin a las aplicaciones queremos destacar lo siguiente. Para un estudio intere-
sante y completo de las curvas de persecucin se recomienda [1]. El mtodo de datacin
empleando Carbono 14 est en el libro del ganador del Premio Nobel de Qumica [7]. Una
vasta coleccin de modelos de poblaciones se encuentra en [11]. El anlisis de circuitos
elctricos y la deformacin de vigas se pueden estudiar en [12] y [6], respectivamente.
Finalmente, el lector interesado en consultar otro texto para un primer curso de
ecuaciones diferenciales puede ver [10], que contiene una buena cantidad de referencias
histricas de ecuaciones diferenciales y de matemticas en general, as como otros ejem-
plos atractivos de aplicaciones. Por otra parte en los textos [2], [3], [4], [5] y [9] encontrar
una exposicin completa y con diferentes grados de profundidad de la teora de ecuaciones
diferenciales.
219
Respuestas a los problemas
Ejercicios 1.2, Pgina 33
1. Primer orden, y si es solucin.
2. Primer orden, y no es solucin.
7. Segundo orden, y si es solucin.
8. Tercer orden, y si es solucin.
9. Segundo orden, y no es solucin.
10. Segundo orden, y si es solucin.
1. x = y/cy/i- 1
2 3 3
2. ^ + 2y + l l
3. 0 = arccot (sen t + c)
5. y = .(
6. (y - 1)6 = e
x
- e~
x
+ c
7. y-2 ln(y + l)=x-5 \n(x + 2) + c
220 Respuestas a los problemas
8. x =
9. r = 4
1-t
2
1 + t
2
10. \{y - I)
3
+ -Vx
2
+ 4 + 21n(> + Vx
2
+ 4) = c
Z
1. y ~ x ln
2. y\2x
2
-
3. sen - =
X
A . > J.XX 0>0
5. x
3
+ y
3
6. xy + y
2
7. y - 2z H
8. ln(2x +
a;
- y
2
) = ex
2
ex
(x
2
-4- v
2
)
3
/
2
= cxy
= 2x
3
h7 = c(x + y4
3y + 2) = 2y -
9. y = x arctanln x + 1)
- I )
4
- x + c
x x
10. sen + t an = cy
2
y y
1. x
A
y
2
x
3
y = c
2. y =
3. j/ =
(x-l)
4. No es exacta.
5. x sen y y eos x + ln xy = c
6. ysent + t
2
e
y
+ 2y = c
Respuestas a los problemas 221
7. ln xy + e
xy
= e
Q-x
2
9. xy + e
y/x
= c
10. No es exacta.
1. /J,(X) x~~
2
, 4x
2
2y + xy
2
ex = 0
2
c
4 ^
5
/ \ 2
c
4 ^
2. /i(x) = x , y = arceos
ox
3. /j,(y) e~
y
, xe~
y
+ \nx 3y = c
4. /i(y) - y
2
, 2: rV + x
2
+ c = 0
5. fJ>(y) = y~
1
, x\nxy + x
2
-y
2
= c
6. /i(x) = x
3
, x
6
y
2
- x
4
j/
4
= c
7. ^{x)x
2
, xsenxy + x
3
y
3
= c
8. /x(i/) = y
2
, xlnxy -y
3
= c
9. //(z) = x~
3
, 2y = x
2
(c + eos 2y)
10. /i(y) = eos"
3
y,
x
= (y + c) eos
2
y
_ x
3
- 3x + c
'
y
~ 3 (
2
)
3.
y
= _I ( i +
x
2
) + c ( 1
4. ?/ = x
3
senx
5. x = ^y"
2
e
2/
(2y
2
- 2y
222
v c
6. x = y\ny- - + -
2 y
8. y = cscx
9.
x
= -e
y
+ ce
10. y - l -
1.
2.
y--
X
2
+
V c
y
2
x
'
X
ex
Q / i
4. y = (x + ex
2
)"
3
5. y = \/3 + ce-
3x2
6. x
2
= "
7. y =
2 + C/
5
1
cVT- x
2
- 7
8. y = ^ce-
8
*
2
+ 2
9. y(l + 21nx + cx
2
10. x(2 - y
2
+ ce'
2
'
2
/
= 4
= 1
1. y
2
+ 2j/(l - x) - 4x
2
+ 12x + c = 0
2. (l + e
5
' )
2
(l+ sen2x)
5
= c
3. y - v
/
TTce
3
A
4. y
7
+ 7e
y
senx = c
, x ln(\/l + x
2
+ x) + c
5. y'
1
= + ,
6. e
x
sen y + 2y eos x = c
7. y = x aresen (c ln x)
x
2
e
x
e
x
1
9. y= r - -l nrr + c
x x^ 2
10. x
2
4- 3y
2
+ 4xy -5x -y = c
11. ex e os - = 1
x
12. (l + cosx)(l + e
2/
) = c
13. (xe
x
e
x
) sen y + y eos ye^ = c
y
1 ex
15. (
v
/
y +l )
2
+ 21n(
A
/ y- 1) + ( </ - I)
2
+ 21n( v^+1) = c
16. xe
y
/
x
ye
y
x
+ x lnx + ex = 0
1?
"
y =
(i
+ e
-*)(i
+ e
*)
a;
2
18. y = eos 3a; + ca;
2
19. 4x = 6v + 3 + ce
25
'
21. y = xln(lncx)
22. 3x
4
+ 4x
3
+ 6x
2
y
2
= c
23. (x +1/)"
1
= ~- e
3 x
+ ce
x
. No se puede resolver con los mtodos estudiados hasta
ahora.
24. 2,= .)
25. xy = y eos y + seny + c
26. ex y + 3x ln x 4x ln y = O
27. x + 2-y/y + 3x 2 + 8 ln(4 -^/y + 3a; 2) = c. No se puede resolver con los mtodos
estudiados hasta ahora.
28. x
3
= -3e2
y2
1 e-*
y2
dy + ce**
2
29. j / = c l n | - - 2 )
\V J
30. (6x
2
y + y
2
)e
x
= c
7/
31. ln(x
2
+ y
2
) - lnx = arctan - + c
x
32. x Iny + e
x
+ sen y = c
33. xV
Ai
( _
T
2 _ 1 34. y = tan
\ x J
35. y
3
= - 3z
2
+ ex
3
/
2
36. y
2
+ xy + ex
3
= 0
37. 3 z V + y
4
= c
X
38. y = xsenx H
39. x
2
eos y + x sen y = c
40. 2x
2
+ cx
2
y
3
-3y = 0
41. x
2
y + x + y
3
+ cy = 0
42. y = ce y
43. x'
1
= -y + 2 +
44. y =
3(x - x
3
)
45. y/x
2
+ y
2
+ y = ex
2
46. y
4
+ y
6
= c
47. sen x eos x sec y + c 0
6v 3 x 1 2 \
4 8 . l n [ ( 3 y - 7 )
2
- ( 3 * - 2 ) ( 3 y - 7 ) + ( 3 x - 2 )
2
] - 2 v
/
3 a r c t a n ^ ^ J
4 9 . y + 2 1 n | y - 1| - z - 5 1 n | z - 3 | = c
5 0 . (
x
y - 2xy
2
+ x
2
-2x + 2y
2
+ 2)
e
2arct an3/
= c
51. l n| l + v\ + 2^/v 2arctan^/; =
52. y = arcsen (c sen )
54. y
6
x sen x y = 1
\l-xyj
56. 4a;
2
y
4
- x
4
y
2
= c
57. y x 3 ln(cc + y 1) = c
58. y = i(z
4
-l)
59. ^/xy + xy
3
/
2
- lO^y + 4 = 0
60. y = v^
3
+ ex
2
61. l 2p
62. 2z + e*" - y
2
= c
63. y ex
64. y = -
3cosrr 6 sena; 6 eos a; c
x x
2
x
3
x
2
65. y=\ sena; + - +
66. e(v) +cy + l = 0
67. y sen - =
68. a:
2
+ 2x
2
y + y
2
+ 2y = c
226 Respuestas a Jos problemas
69. y
2/ 3
= - x
5
+ ex
2
%7
70. y
72. y
1
73. y
7A ,
=
V
/(x3 +
C
)
2
-
1
3
- 6y + 3e*
2
=
ln(c x
2
)
2x
6-x
2
16
3
75.
76.
77. y = ex
2
, con c > 0
78. y = cey
79. x
2
+ y
2
= ex
80. y = ce
x/fc2
+ /c
2
, con k la distancia de F a .
1. tan(x + y) sec(x + y) = x + c
2. x y ln(c x) 5 = 0
3 v-
2
x
- y
x
_ 3x ~
Ce o
4. x + y ln(c eos x) = 0
5. y = t an(x + c) x
6. x = tan(x j / + 1) + c
7. a) ln y + xy = ln x + c
b) eos r H- sen - = ex
x
8. a) y c\ l nx + e<
i \ x
2
e\
b) y = + - + c
2
o x
9. a) y
2
= (x + ci)
2
+ c
2
b) y
e 2

lX C2
-i
1 + ce
x
Ejercidos 2.10, Pgina 90
2
1. y = ex x , la solucin singular es y = -x
3
'
2
2. y ex tan c, la solucin singular es y = tan( arcsec y^) + arcsec

3. y = ex + lnc, la solucin singular es y = lnx 1
4. y = ex ac c
2
, la solucin singular es y = ( )
Ejercicios 3, Pgina 119
1. a) y
2
= 2x + c
b) x
2
+ y
2
= cy
2. a) y
2
= 21ncosx+ 4
b) * = # * *
3. a) y() = 10e-
0256t
b) 8.8 g
4. a) z() = 60e-
007438<
b) 41.37 mg
c) 18.64 horas
5. 0.0433% de la cantidad inicial.
6. 8yo, donde yo denota la cantidad inicial.
7. 2.7 horas
8. 218 animales
9. a) P(t) =
1O
oo - 999e--
0018t
b) 1
10. a) T() - 21 + 77e-
a0139113
b) 122.54 minutos
11. 12.425 minutos. T(20) = 83.61F
12. -30F
13. 46 segundos
14. 650 Ib
15. a) A(t) = | (120 -t)- 180 ( l - - M
b) 2.8125 Ib/gal
16. 10.85 aos.
17. a) x{t) = 0.04(1 _
e
-(-25xio-)
t)
b) 48 minutos.
18. 986.58 Ib
19. A(t) = 300 - 260e-
/5
B(t) = 3(100
+
,
) +
39000+ 260t
e
_
t / 50
_ 68000
100 +1 100 +
20. i(t) = 0.6(1 - e-
500
), 0.6 Amperes
a o 2 t ) 1

=
4
4
22. i(t) = 4 -
( 1 + ao2t)10
ooo

^ =
4
Amperes
23. i(t) = 20e = 20e~
4t
O A A
24- g() = - sen 6 + - eos 6 - - e "
8
'
25 25 25
/ x 18 24 32 _
8t
*
=
25
C
~ 25
S6 n +
25
6
25. g() = e~
3
- e-
6t
, g^ = ?(ln2/3) = 0.25 coulombs
nf
. . . E _
t
.(L\n2\ E
1. y(x) = c
x
e
3x
+ c
2
e~
3x
2. y[x) C\ eos 3x + C2 sen 3x
3. y{x) = Cix + C2X
2
4. ?/(:r) - cix
2
+ 4
4
5. y(x) = c\ sen lnx + C2coslno;
6. y(x) = C1X
2
+ c
2
x
3
7. y(x) = Cix + c
2
x
4
8.

/(x) = ci(a: + l)
2
9. y(x) = c
x
x + c
2
x~
l
10. y(x) =
2
11. y(x) =
12. j/(x) = cix + c
2
f I ln
\2
1 x
13. y(x) = CizT
1
/
2
eos x + c
2
x"
1
^
2
sen x
14. y(x) = c\ + C2 ln x
15. y(x) = cixlnx + c
2
x
16. 2/1 (x) no es solucin de la ecuacin diferencial.
1- y = ci eos V5 x + C2 sen \/5 x
2. y = exe"
10
* + c
2
e
x
3. y = cie
2x
cos2V2x +
4. y =
Cl
e-
2x
+ c
2
e
x
'
4
5. y = cie~
3x
+ c
2
xe~
3x
6 y
_
Cie
(14-2V43)x/4
7. y =
Cl
e
9a;
+ c
2
xe
9x
8. y - de"
1
/
8
eos ^ x + c
2
e~
x
l
%
sen
o
9. y = e
x
+ 2e
5x
10. y = e~
2x
cosx + 2e~
2x
senx
11. y = cos2x sen2x
12. 2 / -
1- y = Cie
x
+ c
2
e"
x/ 2
eos -x + c
3
e~
x/2
sen x
2. / = exe"
21
+ c
2
e
21
+ c
3
eos 2a; + c
4
sen2x
3. y = (ci + c
2
x + csxê"
1
4. y = Ci eos 2x + c
2
sen 2x + c
3
eos 3a; + c
4
sen 3x
5. y = ci eos 4a; + c
2
sen 4a; + c
3
x eos 4a; + c
4
a; sen 4a;
6- y = ci + c
2
e
5x
+ cge"
21
7. y = ci + c
2
eos 2x + c
3
sen 2x + c
4
x eos 2x + c
5
x sen 2a;
8. y = ce
1
+ C2e~
x
eos vâ; + c
3
e~
x
sen y/2x
9. y = cie~
3z
+ c
2
e
x
eos 2x + c
3
e
x
sen2a; + c
4
a;e
I
eos 2a; + c
5
xe
x
sen2a;
l0
-
V =
I
+
\
e
~
2X
~
3e
~
X
11. y =-2e
2x
+ 5xe
2x
- 5x
2
e
2x
12. y = 4 eos x sen x eos 2a; + - sen 2a;
1. (4>-3)(4> + 3)
2. (D-8)(D-3)
3. D(D + 6)
2
4. (D + l ) ( Z?- 2) ( D- 4)
5.
6.
7. >
3
8. D
n
9. >(>-4)
10. >
2
(> - 10)
2
11. D
3
(D
2
+ 81)
12.
13. (>
14. D(D
2
-6D + 25)
15. (>
2
+ W + 13)(>
2
- 12D + 45)
16. y = c
x
e~
lx
+ c
2
e~
x
+ 2e
x
- e
2x
17. y = (ci + C2x)e
x
+ 4cosx 5 sen x
18. y (ci eos 42; + c
2
sen42;)e
x
1 + 42; + 172;
2
19. y = ci + c
2
e~
81
- x + 4x
2
20. y = (ci + c
2
2;)e-
21
+ 42; V
2x
21. y = cie~
3x
+ c
2
e~
x
- 2xe~
3x
+ 6x - 3
22. y = (ci + c
2
a;)e-
51
+ (o;
2
+ z
3
)e-
51
23. y = ci + c
2
e~
31
+ \x - (x + 6x
2
)e~
3x
24. y = ce"
31
+ c
2
e-
2
* + (202; - 5x
2
)e-
21
25. y = e~
x
{c\ cosx + C2senx) -\x
i
-x/2 i ^ y/3 \ (I o 1 1
26. y = e
x /
I c
x
eos a: + c
2
s en x 1 + I - i
2
- - ^ -
27. y =
Cl
e~
( N/ 5+2) x
+ c
2
e
(
^-
2 ) x
- -^(12sen2a; + 16eos 2x)
28. y c\ eos x + C2 sen x + 1 + x(2 sen x 4 eos x)
29. y = C\ eos x + C2 sen x + x
2
senx + x eos x - sen x
30. y C\ eos x + C2 sen x + - x sen x + - eos 3x
2 8
31. y = ci + c
2
e~
a:
+ - x
3
- x
2
+ - x + - e
x
+ sen 2x - eos 2x
O Z A Zi\j 1U
32. y = (ci eos x + c
2
sen x - x eos x + x sen x)e
2x
33. y = (
Cl
+ c
2
x)e
2x
+ x + 1 H (4cosx + 3senx) + -cos2x
25 8
34. y (ci + c
2
x)e~
x
+ 1 + senx + (cos2x + 7sen2x)
o
35. y = (ci + c
2
x)e~
3x
+ 1 + x
3
e~
3x
+ 4senx - 3cosx
36. y = (ci + c
2
x)e"
2x
+ - x + 1
37. y = ce"
3
* + c
2
e
4x
+ 2xe
4x
38. y = de"* + c
2
e
3x
- e
x
+ 3
39. y = de~
x
+ c
2
e
x
+ (2x
3
- 3x
2
+ 3x)e
x
+ 10
40. y = e
x
(ci eos 2x + c
2
sen 2x) + e
x
sen x
41. y c\ eos 5x + c
2
sen 5x 2x eos 5x
42. y = e~
x
(cicos2x + C2sen2x) + (20x 9)senx + (15x 13)cosx
43. y = ci + c
2
x + c
3
e"
x
+ - x
2
i
44. y = ci + c
2
x + c
3
x
2
+ c
4
x
3
+ c
5
e
x
+ (x
2
8x)e
x
x
4
45. y = ci + c
2
x + c
3
e~
x
+ 12x
2
- 4x
3
+ x
4
46. y = ci + c
2
x + c
s
e~
2x
+ x
2
+ x
4
47. y = (c
1
+ c
2
x + c
3
x
2
)e
x
+ - x V + x - 13
48. y = ci + c
2
x + (c
3
+ c
A
x)e
x
+ -x
2
e
x
+ - x
2
49. y = cie
x/ 2
+ c
2
e"
x/ 2
+ c
3
eos - + c
4
sen - + 2xe
x/2
50. y
p
= >lxe
x
+ e
x
( 5 eos 2x + Csen 2x) + xe
x
(Eeos 2x + Fsen2x)
51. y
p
= Ax + (Bx + Cz
2
)e
3
* + (? + Fx + Gx
2
)e"
x
cosx + (H + Ix + Jx
2
)e~
x
senx
52. y
p
A eos x + B sen x + C eos 3x + E sen 3x
53. t/p = (A + Bx)xe
x
+ (Ccosx + Esenx)e
x
+ Feos2x + Gsen2x
54. y
p
- (A + Bx + Cx
2
)cosx + (E + Fx + Gx
2
)senx + ix
2
+ /x
3
+ Jx
4
+
(Lx
2
1. y{x) Cicosx + c
2
senx - (cosx)ln | secx + tanx|
2. y(x) c\ eos x + c
2
sen x x eos x + (senx) ln | sen x|
3. y(x) = c
ie
~
x
+ c
2
e-
3x
+ \[e~
2x
+ {e~
x
- e~
3x
) ln(l + e
x
)]
4. y(x) = c
ie
~
x
+ c
2
e~
2x
+ (x
2
- x) e~
x
5. y(x) = (ci + c
2
x)e~
3x
- e~
3x
lnx
6. y(x) = c
x
e
2x
+ c
2
xe
2x
+ ~e
2x
[x
2
arctanx - arctanx - xln(l + x
2
)]
i
7. y(x) = cie
x
+ c
2
e
2x
+ e
x
eos 2x e
x
sen 2x
8. y(x) = c^-
2
* + c
2
e~*
x
- 2e"
3l
cose
x
- e n s e n e *
9. y(x) = (
Cl
+ c
2
x)e
3x
+(^\nx-^j x
2
e
3x
10. y(x) = Ci eos x + C2 sen x + -x sen x - x
2
eos x
11. y(x) = (ci + c
2
x)e
2x
- -e
2x
\n(x
2
+ 1) + xe
2x
arctanx
12. y(x) =
13. y(x) =
14. y(x) =
15. y (x) =
J_
xarcsenx)
+ C2X~
l
Hx
2
ln
2
x x
2
ln x
6 9
+ c
2
x -
lnx
zx
eos ln x + C2 sen ln x (eos ln x) ln(sec ln x + tan ln x)
1. a) x(t) = 0.1 cos
b) T = TT/3
2. a) x(t) = -0. 4 eos 4 + 0.3 sen 4
b) x(t) = 0.5sen(4-0.9273)
c) 16 oscilaciones.
3. a) x(t) = 0.6 eos 5 - 0.8 sen 5
b) x() - sen (5 + 2,4981)
T17T
c)t
n
= - 0.4996, n = 1,2,3,.
5
d)
Figura 5.21: Grfica de x() = sen (5 + 2.4981).
1 2
4. a) x{t) = cos9 sen9
1U lo
b) x(t) = - sen (9 + 3.7851)
6
c)
235
Figura 5.22: Grfica de x(t) = \ sen (9 + 3.7851)
d) t
n
= nn/9 - 0.4206, n = 3, 5, 7, . . .
e) X(TT/18) = -0./333, x(n/9) = 0.1, x(n/3) = 0.1
f) x'(7r/18) = 0.9, x'{ir/9) = 1.2, x'(ir/3) = 1.2 , en direccin hacia arriba en
los tres casos.
5. a) x(t) = \cos%t
b) T = TT/4
c) 32 oscilaciones
d) t
n
= (2n + 1)TT/16, n = 0,1, 2, . . .
6. a) x(t) = - cos8 + - sen8
4 2
b)
x
( ) = ^ sen (8 + 0.4636)
c) t
n
= mr/8 - 0.05795, n G N
d)
Figura 5.23: Grfica de x(t) = & sen (8 + 0.4636)
7. En la posicin de equilibrio. 2TT/UJ.-
8. a) m = 1/4 slug, k = 1 lb/ft. El movimiento inicia en un punto que est 2 ft
abajo de la posicin de equilibrio y con una velocidad inicial de 4 ft/s dirigida
hacia arriba.
b) m 1/25 slug, A; = 25 lb/ft. El movimiento inicia en un punto que est 0.1 ft
arriba de la posicin de equilibrio y con una velocidad inicial de 4 ft/s dirigida
hacia abajo.
9. a) x{t) | eos 8 \ sen 8
b)
x
{t) = sen (8t + 2.2143)
c)
n
= n7r/4 - 0.2113, n = 1, 2, 3, . . . y i
n
= nyr/4 + 0.0505, rz = 0,1, 2, . . .
10. a) x() = | eos 10 - \ sen lO
b) En t = 1.34 segundos, con una velocidad de 7.071 ft/s.
c) t
n
= nn/5 - 0.1865, n. = 1, 2, 3, . . . y
n
= n?r/5 + 0.02945, n = 0, 1, 2, . . .
d) t ; = (2n + 1) - 0.23562, n = 1, 2, 3, . . .
,-st
b) Nunca pasa por la posicin de equilibrio
c) El desplazamiento mximo es de fe"
1
/
3
ft.
d)
Figura 5.24: Grfica de x(t) = ( + 6)
, - 8t
237
\/l3 _ 7T
2. a) x(t) = -e~
M
sen (2 + 2.5536),
n
= n 1.2768, n > l , e impar.
2 . 5 3 3 . 5
Figura 5.25: Grfica de x(t) = ^ e ~
3 t
s e n (2 + 2.5536)
ftyct
b) x() = - c " * sen ( + 3.3390),
n
= rnr - 3.3390, n > 2, y par.
-o . i
- 0 . 2
- o . 3
5 6
Figura 5.26: Grfica de x(t) = ^ e " * sen ( + 3.3390)
3- a) x(t) = \ \e
- 4
- O . 2
- O . 4
- O . 6
- O . 8
Figura 5.27: Grfica de x(t) = e~
2
4. a) (3 > 2
b) p = 2
c) O < P < 2
5. ô < - 4 - \ < v
0
< 0
6. m = 1/8 slug, k = 1/4 lb/ft. El movimiento inicia en un punto que est 1 ft
arriba de la posicin de equilibrio y con una velocidad inicial de 1 ft/s dirigida
hacia arriba. El medio ofrece una resistencia al movimiento numricamente igual a
dos veces la velocidad instantnea.
7. a) Movimiento sobreamortiguado o crticamente amortiguado. b) No.
1 2
1. a) x(t) = -~e~
4t
+ -
6 3
b)
- eos 2 + 2 sen 2
Solucin Transitoria
Solucin Estacionaria
Figura 5.28: Soluciones transitoria y estacionaria del problema 1
c)
- 1
- 2
Figura 5.29: Grfica de x(t) = -he~
4t
+ fe~* - cos2 + 2sen2
-
&)X[t)
48 48 4
b)
1 2 3 4 5 6 t
Figura 5.30: Grfica de x(t) = ge"* - e~
5
+ Je
11 58
3.
x
(t) = (-=- + ^- - 4cos4 - 3sen4
Figura 5.31: Grfica de x(t) = ( + f i 4 eos 4 3 sen 4
4.
1 25 5
= e~*(- eos 2 H sen 2) + -e~
sen 2
4 8 2
3 4 5
Figura 5.32: Grfica de x(t) = e~*(J eos 2 + f sen 2) + fe"* sen 2<
1 13 1
5. x{t) = cos2+ sen2 + 5(sen2 cos2)
4 Z
Figura 5.33: Grfica de x(t) = -\ eos 2 + ^
s e n 2
+ 5( sen 2t-\ eos 2)
6. x() = -
3
fa ,
u
^
- eos v 5 H sen
+ (eos 2 - 2 sen 2t)e~
l
Figura 5.34: Grfica de x(t) = -feos \/Et + ^ ^ sen y/bt + (eos 2 - 2 sen 2t)e
, - t
d x dx
7. a) La ecuacin diferencial es m- = k(x h)P, que en este caso se reduce a
(JLL LLL
d x dx
+3 + 2x = 4cos3; sujeta a las condiciones iniciales x(0) = 0, ^(O) = 0.
o o o o
b) x() = -;r-e~
+ e~
2t
- - eos 3 + - sen 3
c)
Figura 5.35: Grfica de x() = - e "
+ f cos3 + sen 3
8.
x
(t) = -(cos\/2-cos3)
241
Figura 5.36: Grfica de x(t) = |(cos y/2t cos3)
1. q(t) = 7- e- ' (
i(t) = 12e~~
sen2
/2 3 \ 2
2. q(t) = e~
u
- eos 8 + sen 8 ) - - eos lO
\5 1 0 / 5
i(t) = 4 sen 10* - 5e~
6
sen 8
3. a) q[t) = -e~i' (4cos5
b) t() = 25-i
c) 4 Coul y 0 A.
+4
4. a) g() = e"
4
(3cos- 5sen) - 2cos2 + 4sen2
b) i
e
(t) = v^s e n (2 + 1.1071)
5. a) g() = -0.02eos50 +0.02, i(t) = senhOt
b) 0.04 Coul y 1 A.
6. a) q(t) = - - sen 100 + sen 50
150 75
2 2
b) () = eos 100< + - eos 50
O O
c) t
n
= , n = 1, 2, 3, . . .
50
7. a) q(t) = -10e"
3
(cos3 + sen3) + 10
b) i(t) = 60e~
3i
sen3
c) q
m
ax = g(7r/3) = 10.432 Coul
d) Est: sobreamortiguado si R > \/200, crticamente amortiguado si R v
/
200
y subamortiguado si 0 < R < \/200.
8. a) 0 < R < 2VS
b) 0 < R < y/6
Bibliografa
[1] Bornhart A., Curves of Pursuit, Scripta Mathematica, Vol 23 (1957), 49-66.
[2] Birkhoff, G. y Rota, G. C, Ordinary Differential Equations, 3a. edicin, Wiley, New
York, 1978.
[3] Braun, M., Differential equations and their applications, 4a. edicin, Springer-Verlag,
New York, 1993.
[4] Coddington, E. A. y Levinson, N., Theory of Ordinary Differential Equations, Me
Graw-Hill, New York, 1955.
[5] Elsgoltz, L., Ecuaciones Diferenciales y Clculo Variacional, 2a. edicin, MIR,
Mosc, 1977.
[6] Hibbeler, R., Anlisis Estructural, 3a. edicin, Prentice Hall Hispanoamericana,
Mxico, 1997.
[7] Libby, Willard, Radiocarbon Dating, 2a. edicin, University of Chicago Press, 1955.
[8] Marsden, Jerrold E. y Tromba, Anthony J. , Clculo Vectorial 3a. edicin, Addisson-
Wesloy, 1991.
[9] Pontryagun, L. S., Ordinary differential equations, Addisson-Wesley, 1962.
[10] Siminons, F., Ecuaciones Diferenciales con Aplicaciones y Notas Histricas, 2a.
edicin, Me Graw-Hill, Mxico, 1995.
[11] Skmk, W., Matemticas Aplicadas a la Ecologa, CINVESTAV, Mxico, 1987.
[12] Svoboda, D., Circuitos Elctricos, 3a. edicin, Alfa Omega, Mxico, 2000.
[13] Wolfram, S., Mathematica, Addison-Wesley, 1991.
[14] Mdulo TI-92 plus, Suplemento para el manual del usuario de la TI-92. 1998, Texas
Instruments.
243
ndice
Prologo 7
Introduccin 9
1.1 Ecuaciones Diferenciales y Modelos Matemticos 9
1.2 Conceptos Bsicos de Ecuaciones Diferenciales 22
1.2.1 Mtodo de Isoclinas 28
Ecuaciones Diferenciales de Primer Orden 35
2.1 Ecuaciones Diferenciales de Variables Separables 35
2.2 Ecuaciones Diferenciales Homogneas 44
2.3 Ecuaciones Diferenciales Exactas 53
2.4 Factores Integrantes 62
2.5 Ecuaciones Diferenciales Lineales 70
2.6 Ecuacin de Bernoulli 74
2.7 Miscelnea de Ecuaciones Diferenciales 78
2.8 Cambios de Variables Diversos 84
2.9 Ecuacin de Ricatti 87
2.10 Ecuacin de Clairaut . . 89
Aplicaciones de las Ecuaciones Diferenciales de Primer Orden 91
3.1 Trayectorias Ortogonales 91
3.2 Mecnica 93
3.3 Aplicaciones a problemas relacionados con crecimiento y decrecimiento . 96
3.3.1 Desintegracin Radiactiva 96
3.3.2 Problemas de Enfriamiento 101
3.3.3 Modelos de Poblacin 104
3.4 Mezclas 108
3.5 Circuitos Elctricos LR y RC en Serie 112
3.5.1 Circuito LR en Serie 114
3.5.2 Circuito RC en Serie 116
Ecuaciones Diferenciales Lineales de Segundo Orden 123
4.1 Conceptos Bsicos . . . 123
4.2 Solucin de Ecuaciones Diferenciales Lineales de Segundo Orden 126
4.2.1 Solucin de la Ecuacin Homognea 127
4.2.2 Solucin de la Ecuacin no Homognea 128
4.3 Mtodo de Reduccin de Orden 129
4.4 Ecuaciones Diferenciales Lineales de Segundo Orden Homogneas con Coe-
ficientes Constantes 133
4.5 Ecuaciones Diferenciales Lineales de Orden n 139
4.5.1 Ecuaciones Diferenciales Lineales de Orden n Homogneas con Coe-
ficientes Constantes 140
4.6 Mtodo de Coeficientes Indeterminados: Enfoque de Superposicin . . . . 145
4.7 Mtodo de Coeficientes Indeterminados: Enfoque del Operador Anulador 158
4.7.1 Operadores Diferenciales 158
4.7.2 Mtodo de los Coeficientes Indeterminados 162
4.8 Mtodo de Variacin de Parmetros 173
Aplicaciones de las Ecuaciones Diferenciales Lineales
de Segundo Orden 179
5.1 Movimiento Armnico Simple 179
5.2 Movimiento Vibratorio Amortiguado 189
5.3 Movimiento Vibratorio Forzado 201
5.3.1 Resonancia 203
5.4 Circuito LRC en Serie 208
5.5 Otras Aplicaciones 212
5.5.1 Vigas Horizontales 212
5.5.2 El Pndulo Simple 215
Conclusiones 217
Bibliografa 243
Ecuaciones diferenciales. Tcnicas de
solucin y aplicaciones se termin de
imprimir en el mes de abril de 2004, en
los talleres de AGES, en la Ciudad de
Mxico. Se utilizaron los tipos Palati-
no y Carleton. Los interiores estn im-
presos en papel kromos ahuesado de
90 g y la portada en Multiart de 250
g. Se tiraron 1,000 ejemplares. El cui-
dado de la edicin estuvo a cargo de
Silvia Guzmn Bofill y los autores.
DAVID ELIZARRARAZ MARTNEZ
Licenciado en Matemticas, obtuvo el
grado de Maestro en Matemticas con
la tesis "Anlisis comparativo de dis-
tintos mtodos para la resolucin nu-
mrica de problemas de Stefan" y es
doctor en Ciencias por la Universidad
Autnoma Metropolitana, Iztapalapa.
El ttulo de su tesis de doctorado es
"Mtodos algebraicos para resolver
ecuaciones funcionales lineales". Des-
de 1993 ejerce como profesor titular
del departamento de Ciencias Bsicas
en la UAM, unidad Azcapotzalco; don-
de es miembro del rea de Anlisis
Matemtico y sus Aplicaciones. Ha
impartido el curso de Ecuaciones Di-
ferenciales en varias ocasiones y ha
publicado algunos artculos en revis-
tas especializadas. Sus investigaciones
se centran en ecuaciones diferenciales,
Mtodos operacionales y Transfor-
madas integrales.
ISBN 970-31-0230-1
Este libro est diseado para un curso trimestral de ecuaciones diferen-
ciales ordinarias. Presentamos los teoremas y tcnicas de solucin que
consideramos bsicos en un estudio introductorio de esta importante disci-
plina de las Matemticas. Aunque no hemos puesto nfasis en las demos-
traciones, proporcionamos una buena cantidad de ejercicios resueltos, de
modo que un estudiante de Ingeniera podra obtener, mediante su an-
lisis, un nivel satisfactorio en los diferentes mtodos de solucin de ecua-
ciones diferenciales y sus aplicaciones elementales ms comunes. Los
diferentes temas se exponen en forma clara y sencilla para su inmediata
comprensin. En las primeras secciones los desarrollos se hacen de ma-
nera exhaustiva. Ms adelante, lo que ya es conocido no se desarrolla com-
pletamente, sino que se dejan al lector los detalles que en ese momento
ya est en capacidad de realizar. En consecuencia, el texto se puede utili-
zar para un curso tradicional o bien, para un curso en el Sistema de Apren-
dizaje Individualizado, en el que el alumno estudia por su propia cuenta.
3O AOS
...transformando el dilogo por la razn

Ecuaciones Diferenciales Tecnicas de Solucion y Ecuaciones

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Ecuaciones Diferenciales Tecnicas de Solucion y Ecuaciones

Cargado por

Información del documento

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

Ecuaciones Diferenciales Tecnicas de Solucion y Ecuaciones

Cargado por

Copyright:

Formatos disponibles

ECUACIONES DIFERENCIALES

TCNICAS DE SOLUCIN Y APLICACIONES

x"{2) = - ^ sen 16 = 0.76774,

Definicin 5.3.1 Cuando la frecuencia de la fuerza exterior es

cos2. Grafique la solucin obtenida.

2. y ex tan c, la solucin singular es y = tan( arcsec y^) + arcsec

También podría gustarte