Módulo Ecuaciones

Introducción a las Ecuaciones
Diferenciales Ordinarias
(NOTAS DE CLASE)
Documento de trabajo
Para los estudiantes de la Universidad del Magdalena
Compilado por el docente

Darwin Dacier Peña González
UNIVERSIDAD DEL MAGDALENA

FACULTAD DE INGENIERÍAS
Santa Marta
2007
cuaciones Diferenciales 2
Darwin D. Peña G.
PRESENTACIÓN
Las Ecuaciones Diferenciales se encuentran inmersas en una rama de las matemáticas

(Análisis Matemático) desde hace 300 años. Estas constituyen el objetivo natural del
cálculo elemental y la parte de la matemática que nos permite comprender los fenómenos
físicos.
Por esta razón, es una de las asignaturas de gran importancia en la formación de todo
ingeniero, pues su rigor matemático es decir, el papel de la demostración es importante para
el entendimiento y la convicción razonable de los temas planteados.
Las siguientes notas de clase no pretende suplir ningún libro texto del mercado, solo
intenta facilitar el desarrollo de la clase, proporcionándole al estudiante de la Universidad
del Magdalena un material de apoyo con ejercicios resueltos y planteados que le permitan
un mejor y rápido desarrollo de los conceptos desarrollados en clase.
Mostraremos de manera sencilla desde los orígenes de las Ecuaciones Diferenciales hasta la
transformada de Laplace y los sistemas de ecuaciones diferenciales, pasando por los
algunos métodos de solución para las ecuaciones de primer orden y de orden superior, todo
esto complementado por un buen número de ejercicios resueltos y planteados.
3 cuaciones Diferenciales
Darwin D. Peña G.
Capítulo I
Darwin D. Peña G.
1. ORIGEN DE LAS ECUACIONES DIFERENCIALES
Las ecuaciones diferenciales fueron inicialmente tratadas por Newton para estudiar el
movimiento planetario, luego fue progresando a medida que se afianzo en la ciencia
natural, especialmente en la física con problemas importantes como, la ley del movimiento
de Newton, las ecuaciones de Euler para la hidrodinámica, la ecuación de calor por Fourier,
etc.
Actualmente las ecuaciones diferenciales no solo se utilizan en el campo de la física, sino
también en la ingeniería, de la química, economía, agronomía, etc., de ahí que su estudio
sea indispensable para la especulación de toda ciencia natural. Consideremos ahora tres
maneras de formular las ecuaciones diferenciales.
1.1 Problemas geométricos
Ejemplo 1.1 se define una curva con la condición de que cada uno de sus puntos (x,y) su
dy
pendiente es igual al doble de la suma de las coordenadas del punto. Expresar la
dx
condición mediante una ecuación diferencial.
dy
Solución: La ecuación que representa la condición es  2( x  y )
dx
Ejemplo 1.2 hallar la ecuación diferencial de un espejo curvado tal que la luz de una fuente
situada en el origen se refleje en èl como un haz de rayos paralelos al eje x.
Solución:
Figura 1
De la ley de reflexión tenemos que    y por la geometría del problema tenemos que
   y     180    180 es decir       2 .
y 2 tan  dy
Puesto que tan   y tan   tan 2  , además conocemos que  tan  ,
x 1  tan 
2
dx
Darwin D. Peña G.
con lo que
dy
2 2
y
 dx , despejando dy se tiene 2 dy  y  y  dy 
 
x  dy 
2
dx dx x x  dx 
1  
 dx 
2
x  2x  x 4x 2  4 y 2
2
 2  
 y   4  2 
 dy  x dy dy y   dy y y2
   2  1  0    
 dx  y dx dx 2 dx 2
dy  2 x  2 x  y dy  x  x  y
2 2 2 2
  
dx 2y dx y
1.2 problemas físicos
Ejemplo 1.3 Una masa m de peso W se suspende del extremo de una varilla de longitud
constante l. Suponiendo que el movimiento se realiza en un plano vertical, se trata de
determinar al ángulo de desplazamiento  , medido con respecto a la vertical, en función
del tiempo t.
Figura 2
Solución:
Consideremos   0 a la derecha de OP y  < 0 a la izquierda de OP. Recordemos que el
arco s de un círculo de radio l se relaciona con el ángulo del centro  por la fórmula
s  l
d 2s d 2
Por lo tanto, la aceleración angular es a  2  l 2 . Por la segunda ley de Newton,
dt dt
d 
2
F  ma  ml 2
dt
Darwin D. Peña G.
En la figura (2) vemos que la componente tangencial de la fuerza debida al peso W es mg

sin . Si no se tiene en cuenta que la varilla y se igualan las dos expresiones de la fuerza
d 2 d 2 g
tangencial se obtiene ml 2  mg sin   2  sin   0
dt dt l
Ejemplo 1.4 halle la variación de la atmósfera con la altura en la atmósfera terrestre

suponiendo que la temperatura es uniforme en todas sus partes,( lo cual no es cierto).
dp
p  po  gy    g
dy (1)
para hallar la densidad sabemos que pV  Rnt ,luego:

pM dp pM
 Reemplazando en la ecuación (1),  g. Despejando variables e
Rt dy Rt
dp 2 2 Mg p  Mg
integrando
1 p  
1 Rt
dy  ln  2   
 p1  Rt
( y 2  y1 )
p1  y1 Son condiciones iniciales, por tanto p1  p o y y1  0

Mg
 y2
p 2  p o e Rt
La anterior formula, nos sirve para hallar la presión en cualquier parte del mundo
Ejemplo 1.5 Considere el circuito simple conectado en serie que se muestra en la figura
(3) y que consta de un inductor, un resistor y un capacitor. La segunda ley de Kirchoff dice
que la suma de las caídas de voltaje a través de cada uno de los componentes del circuito es
igual a la tensión E(t) aplicada. Si llamamos q(t) a la carga del capacitor en un instante
dy
cualquiera, entonces la corriente i(t) esta dada por i  . Encontrar la ecuación diferencial
dc
Figura 3
di d 2q
Solución: Sabemos que las caídas de voltaje son: En un inductor = L  L 2
dt dt
Darwin D. Peña G.
1 dq
En un capacitor = q ; En un resistor = iR  R . En donde L, C y R son constantes
C dt
llamadas inductancia, capacitancia y resistencia, respectivamente. Para determinar q(t)
d 2q dq 1
debemos, por lo tanto, L 2  R  q  E (t )
dt dt C
1.3 Primitivas
Ejemplo 1.6 Obtener la ecuación diferencial asociada con la primitiva
y  x 2   x  
Solución: Como hay tres constantes arbitrarias, se consideran las cuatros ecuaciones
dy d2y d3y
y  x 2  x     2x    2  2  3  0
dx dx dx
La última ecuación está libre de constantes arbitrarias y es del orden adecuado, luego está
es la ecuación pedida.
Ejemplo 1.7 Obtener la ecuación diferencial asociada con la primitiva

y   cos x   sin x , Siendo  ,  constantes arbitrarias y  una constante fija
Solución: Derivando tenemos

dy
  sin x   cos x
dx
d2y
2
 2 cos x  2 sin x   2 (cosx  sinx)  - 2 y
dx
d2y
Luego la ecuación diferencial pedida es 2
 2 y  0
dx
Darwin D. Peña G.
EJERCICIOS 1
1.1 ¿Cuál es la ecuación diferencial de la velocidad v de un cuerpo de masa m que cae

verticalmente a través de un medio que opone una resistencia proporcional al cuadrado
de la velocidad instantánea?
1.2 ¿Un medicamento se inyecta en el flujo sanguíneo de un paciente con una intensidad
constante de r gramos /segundos. Simultáneamente, la sustancia se elimina con una
rapidez proporcional a la cantidad de sustancia x(t) presente en cada instante. Halle la
ecuación diferencial que rige la cantidad x(t).
1.3 Halle la ecuación diferencial de la familia de circunferencias de radio r cuyos centros

están en el eje x.
1.4 Hallar la ecuación diferencial de la familia de parábolas cuyos focos están en el origen
y cuyos ejes están sobre el eje x.
1.5 Obtenga la ecuación diferencial de la familia de curvas. x 2  y 2  1
1.6 Obtenga la ecuación diferencial de la familia de curvas. y  e x cos x  e x sin x
1.7 Hallar la ecuación diferencial de todas las circunferencias del plano.

Sugerencia: utilice x 2  y 2  2ax  2by  c  0 .
Darwin D. Peña G.
2. DEFINICIONES BÁSICAS
Definición 2.1 Si una ecuación contiene las derivadas o diferenciales de una o más
variables dependientes con respecto a una o más variables independientes, se dice que es
una ecuación diferencial.
2.1 Clasificación de las ecuaciones diferenciales según su tipo
Si una ecuación contiene sólo derivadas ordinarias de una o más variables dependientes con
respecto a una sola variable independiente, entonces se dice que es una ecuación diferencial
ordinaria (EDO).
Ejemplo 2.1
 dy
 dx  5 y  1

( x  y )dx  4 ydy  0

 du dv Son ecuaciones diferenciales ordinarias.
 dx  dx  x
 2
 d y  2 dy  6 y  0
 dx 2 dx
Una ecuación que contiene las derivadas parciales de una o más variables dependientes de
dos o más variables independientes se llama ecuación parcial.
Ejemplo 2.2
 u v
 y  x  0

 u u Son ecuaciones diferenciales Parciales.
x y u
 x y
  2u  2u u
 2  2 2
 x t t
2.2 Clasificación de las ecuaciones diferenciales según el orden
El orden de las más alta derivada en una ecuación diferencial se llama orden de la ecuación
Ejemplo 2.3
3
d2y  dy 
 5   4 y  x
 dx 
2 Es una ecuación diferencial de segundo orden.
dx
Darwin D. Peña G.
Ejemplo 2.4
 4u  2u u
 2 2 Es una ecuación diferencial parcial de cuarto orden.
x 4
t t
2.3 Clasificación de las ecuaciones diferenciales según la linealidad o no linealidad.
Se dice que una ecuación diferencial es lineal si tiene la forma

dny d n 1 y dy
a n ( x) n  a n 1 ( x) n 1  ...  a1 ( x)  a 0 ( x) y  g ( x)
dx dx dx
Es decir, la variable dependiente y junto con todas sus derivadas son de primer grado esto
es, la potencia de cada término en y es 1 y, cada coeficiente depende sólo de la variable
independiente x.
Ejemplo 2.5

 xdy  ydx  0

 y ' '2 y ' y  0 Son ecuaciones diferenciales lineales ordinarias.
 d3y 2
x 3 2 d y dy
 x  3x  5 y  e x
 dx 3 dx 2
dx
Ejemplo 2.6
 yy ' '2 y '  x
 3
d y No son ecuaciones diferenciales lineales ordinarias. ¿Por qué?
 3  y 2
 0
 dx
Definición 2.2 Se dice que una función f cualquiera, definida en algún intervalo I, es
solución de una ecuación diferencial en el intervalo, si sustituyéndola se reduce a una
identidad. En otras palabras, si y = f (x) tiene por lo menos n derivadas es solución de la
 dy dny
ecuación diferencial ordinarias F  x, y, ,..., n   0. Es decir,
 dx dx 
 
F x, f ( x), f ' ( x),..., f n ( x)  0 , Para toda x  I.
Ejemplo 2.7 La función y = xex es una solución de la ecuación diferencial. y' '2 y' y  0
En    x  
Solución: Derivando 2 veces tenemos

y '  xe x  e x
y ' '  xe x  2e x
Reemplazando en la ecuación diferencial se tiene xe x  2e x  2( xe x  e x )  xe x  0
Darwin D. Peña G.
Ejemplo 2.8 Verificar si x = 4t +1 (2), y = t2 – 2 (3) es una solución para la ecuación

diferencial
2
 dy 
4   y  2 (4)
 dx 
dy dy dt
Solución: 
Aplicando la regla de la cadena es decir, (5)
dx dt dx
dy x 1 x 1 dt 1
De (3) tenemos  2t , y de (2) t     
dx 4 4 4 dx 4
Reemplazando en (5) obtenemos

t2 
2 2 2
dy dy dt 1 1  dy  1  t2  dy 
  (2t )   t      t    4   4 
dx dt dx 4 2  dx  2  4  dx  4
2 2
 dy   dy 
4   t 2  y  2  4   y  2
 dx   dx 
Obs. 2.1 Al resolver una ecuación de n – eximo orden F (x, y, y’,…., y(n)), en donde y(n)
dny
significa , esperamos obtener y una familia n – paramétrica de soluciones
dx n
G (x, y, c1,….cn) = 0.
Obs. 2.2 Una solución de una ecuación diferencial que no contiene parámetros arbitrarios
se llama solución particular.
Obs. 2.3 A veces una ecuación diferencial tiene una solución que no puede obtenerse dando
valores específicos a los parámetros en una familia de soluciones. A tal solución se le llama
solución singular.
Darwin D. Peña G.
EJERCICIOS 2
2.1 Diga si las ecuaciones diferenciales dadas son lineales o no lineales. Indique el orden
de cada ecuación.
1. xdy  dx  0
2. x 2 dy ( y  xy  xe x )dx  0
3. y iv  2 y '3  0
2
dy d2y
4.  1   2 
dx  dx 
d3y d2y dy
5. 3
 2  3x  5 y  e x
dx dx dx
6. (sin x) y ' '(cos x) y  1
2.2 Verifique que la función indicada es una solución de la ecuación diferencial dada.
1. x 2 y' ' xy '2 y  0 ; y  x cos(ln x), x  0
d3y 2
2 d y dy
2. x 3 3
 2 x 2
 x  y  12 x ; y  c1  c2 ln x  4 x 2 , x0
dx dx dx
3. y   25  y 2 ; y  5 tan(5x)
1
4. y  2 xy   y( y ) 2 ; y 2  c1 ( x  c1 )
4
 
5. x y  xy  2 y  x ln x ;
2
y  c1 x cos(ln x)  c2 xsen(ln x)  x ln x
dp ac1e at
6.  p(a  bp) ; p
dt 1  bc1e at
2.3 Encuentre valores de m tales que y= xm sea solución de las ecuaciones diferenciales.
1. x 2 y  y  0
2. x 2 y  6 xy '4 y  0
3. y   5 y   6 y  0
Darwin D. Peña G.
3. ECUACIONES DIFERENCIALES DE PRIMER ORDEN
El problema de las ecuaciones diferenciales consiste esencialmente en encontrar la

primitiva que dio origen a la ecuación. En otras palabras, resolver una ecuación diferencial
de orden n es, en realidad, hallar una relación entre las variables contenido n constantes
arbitrarias, que, junto con las derivadas obtenidas de ella, satisfaga la ecuación diferencial.
Ejemplo 3.1
Ecuación diferencial Primitiva
3
d y y  x   x  
2
0
dx 3
d3y d2y dy y  e 3 x  e 2 x  e x
3
 6 2  11  6 y  0
dx dx dx
2
 dy  (x   )2  y 2  r 2
y    y2  r2
2
 dx 
Las condiciones que ha de cumplir una ecuación diferencial para poder ser resuelta se dan
en los teoremas de existencias. Por ejemplo una ecuación diferencial de la forma y’= g(x, y)
g
en la que g (x,) es continua y uniforme en una región R de puntos (x, y), existe y es
y
continua en todos los puntos R, Admite infinitas soluciones f (x, y,C) = 0 (siendo C una
constante arbitraria), tales que por cada punto de R pasa una y sola una curva de la familia
f (x, y,C) = 0
Variables separables
Definición 3.1 Se dice que una ecuación diferencial de la forma

dy g ( x)
   h( y )dy   g ( x)dx Es separable o que tiene variables separables.
dx h( y )
Ejemplo 3.2 Resolver (1 + x )dy – ydx = 0
Solución:
dy dx
(1  x)dy  ydx  0  (1  x)dy  ydx   y

1 x
 ln y  ln 1  x  c1
ln 1 x  c1 ln 1 x
ye  ye e c1  y  / 1  x)e c1  y  c(1  x)
dy
Ejemplo 3.3 Resolver  y 2  4, sujeta a y(0) = 0
dx
Darwin D. Peña G.
dy dy
Solución :  y2  4  2  dx
dx y 4
Utilizando fracciones parciales en el lado izquierdo tenemos
1 A B
   A( y  2)  B( y  2)  1  ( A  B) y  2 A  2 B  1
y 4 y2 y2
2
A B  0
2 A  2B  1
Luego A  1 ; B   1
4 4
Así que:
1 1  1 1 
   
1 4  4  4 dy  dx  1 ln y  2  1 ln y  2  x  c
 4 
y 4 y2 y2  y 2 y  2
2 1
4 4
 
 
 y 2 y2
ln    4 x  c 2   ce 4 x  y  2  ( y  2)ce 4 x  y  ce 4 x y  2  2ce 4 x 
 y  2  y  2
1  ce 4 x
y2
1  ce 4 x
Ahora al sustituir x  0, y  0 se tiene

1 c
02  c  1  0  c  1
1 c
1  e4x
Luego y2
1  e4x
dy
Una ecuación diferencial de la forma  f (x  by  c), b  0, puede reducirse siempre a
dx
una ecuación de variables separables por medio de la sustitución u  x  by  c.
dy 1
Ejemplo 3.4 Resolver 
dx x  y  1
du dy
Solución: Sea u = x + y + 1, derivando obtenemos 1 con lo que
dx dx
du 1 u u 1 1  1 
1   du  dx  du  dx  1  du  dx
dx u 1 u 1 u  1 u 
Integrando obtenemos:
u  ln 1  u  x  c1  x  y  1  ln x  y  2  x  c1  x  y  2  ce y
Darwin D. Peña G.
Aplicaciones a la química
Existe una variada gama de aplicaciones de las ecuaciones diferenciales de primer orden a procesos
químicos. Entre las aplicaciones más comunes se tienen los problemas de mezclas, es decir, los
problemas que genéricamente se conocen como problemas de salmuera y los relativos a las reacciones
químicas, o sea, problemas derivados de la aplicación de la ley masa-acción.
Problema resuelto 1.
Se tiene un depósito con un volumen: Vi, de una solución que tiene una concentración: c1, Se desea
cambiar la concentración de la solución mediante la apertura de dos válvulas tales que por una de ellas
entre solución de concentración conocida: ce a una rala determinada: A y otra por la que la mezcla
homogénea escape a una rata: B . Se desea determinar el volumen y la concentración en todo instante
posterior a la apertura de las válvulas.
Solución: Las figuras 4 y 5 muestran al sistema en los instantes: t = 0 y t > O. Se toma como referencia
el momento en que se abren las válvulas.
Fig.4 fig. 5
Es claro que tanto el volumen como la cantidad de soluto están variando con respecto al tiempo. La
q(t )
concentración en todo instante se define como: c(t )  . Siendo q(t) la cantidad de soluto en todo
V (t )
instante. La variación del volumen con respecto al tiempo debe ser igual a la rata de entrada menos la
rata de salida, es decir, el volumen en todo instante esta regido por el problema de valor inicial:
dV (t )
 A B V (0)  Vi . Sí las ratas de entrada y salida son constantes, la solución del problema
dt
es: V (t )  ( A  B)t  Vi
En cuanto a la cantidad de soluto, la variación con respecto al tiempo viene dada por la rata de entrada
de soluto: Ace menos la rala de salida: Bc(t}, es decir, la cantidad de soluto está regida por el problema
dq(t ) B
de valor inicial:  q(t )  Ac e q(0)  q i
dt V (t )
La ecuación diferencial para la cantidad de soluto es lineal y se resuelve con las técnicas previamente
estudiadas. Particularmente, si se tiene un depósito con 400 gramos de salmuera que tiene una
concentración de 0.1 Kilogramos de sal por litro al que entra agua pura a una rata constante de 12 Litros
por minuto y la mezcla homogénea sale a una rata de 8 Litros por minuto, el volumen y la cantidad de
sal en todo instante se calculan de la siguiente manera:
Darwin D. Peña G.
dV (t )
a) Para el volumen: 4 V (0)  400. . El volumen en todo instante viene dado por:
dt
V (t )  4t  400
b) Para la cantidad de sal, puesto que la cantidad inicial es q1  400(0.1)  40, , el problema de valor
dq(t ) 2
inicial asociado es:  q(t )  0 q(0)  40
dt t  100
La solución del problema, esto es, la cantidad de sal en todo instante, vendrá dada por:
q(t )  400000(t  100) 2
En cuanto a la concentración en lodo instante, se tiene: c(t) = 100000(t + 100)-3. Las figuras 6 y 7
muestran las graficas de la cantidad de sal y de la concentración en todo instante.
El análisis del problema nos leva a decir que el tiempo necesario para que la concentración de la
solución en el tanque baje a 0.01 Kilogramos de sal por Litro es de, aproximadamente, una hora y 55
minutos.
Figura 6 figura 7
Problema resuelto 2
Antes del mediodía se introduce una torta en un horno precalentado a 300 grados centígrados. A los
doce meridianos la temperatura de la torta es de 200 grados y a la una de la tarde es de 250 grados. Sí la
temperatura inicial de la torta es de 30 grados, determine la temperatura de la torta en lodo instante y
calcule la hora aproximada en que se introdujo al homo.
Solución: Supongamos que el lapso de tiempo transcurrido entre el instante en que se introdujo la torta
y las doce del día es: t1 . Del enunciado del problema se sigue que: T (t1 )  200 T (t1  1)  250 .
Ahora bien, la temperatura en todo instante viene dada por:
t t
T (t )  300  (30  300)e T  300  270e T
A partir de las dos condiciones dadas resulta el sistema de ecuaciones:
t1

200  300  270e T
t1 1
250  300  270e T
El sistema se puede escribir en la forma:

Darwin D. Peña G.
t1

270e T
 100
t 1
1
270e T
 50
1
Dividiendo las ecuaciones, se tiene: e T  2
Por tanto, la constante de tiempo del sistema es: T  1.443 horas . En cuanto a: t1 se encuentra que:
t1  T ln(2.7)  1.443 . Se concluye que la torta se introdujo al horno aproximadamente a las once horas
y 34 minutos. La figura 8 muestra la gráfica de la temperatura en todo instante.
T(t)
Figura 8
De la figura se concluye que la torta se puede sacar del homo al cabo de cinco horas a partir del instante
en que se introdujo, es decir, a las tres y media de la tarde. Evidentemente se trata de una especie de
torta poco común.
Una partícula de masa constante: m es atraída al origen con una fuerza proporcional a la distancia.
Determine la posición y la velocidad de la partícula en lodo instante si se suelta desde un punto que
dista: a metros del origen.
Solución: La figura 9 muestra la situación planteada.
Figura 9
dv
El problema de valor inicial asociado a la velocidad en todo instante es: m  kx v(a)  0
dt
dv
La regla de la cadena permite escribir la ecuación diferencial en la forma: mv  kx  0
dx
La ecuación diferencial es de variables separables y su solución general viene dada por:

k
v 2  C  x 2 . Puede verse que el coeficiente de x2 tiene unidades de frecuencia al cuadrado. Por tanto,
m
hacemos:  2  k , con lo que la solución general se puede escribir en la forma: v 2  C   2 x 2 .
m
Darwin D. Peña G.
Aplicando la condición inicial: v(a)  0 , se obtiene que la velocidad de la partícula en todo instante es:
v(t )   a 2  x 2 . En cuanto a la posición de la partícula en todo instante, tenemos el problema de
dx
valor inicial:  dt x(0)  a . Integrando y evaluando la condición inicial se encuentra
a  x2
2
que la posición de la partícula en todo instante es: x(t )  a cos(t ) . La granea de la posición es una
señal cosenoidal de frecuencia:  . Se trata de un movimiento armónico simple. En ausencia de
rozamiento, la partícula describe un movimiento oscilatorio de amplitud: a alrededor del origen.
EJERCICIOS 3
3.1 Resuelva las ecuaciones diferenciales dadas sujetas a la condición inicial que se indica.
1. senx(e  y  1)dx  (1  cos x)dy, y ( 0)  0
2. (1  x 4 )dy  x(1  4 y 2 )dx  0, y (1)  0
 
3.
d `x
 
 4 x2 1 , x   1
dy 4
5
4. y '2 y  1, y (0) 
2
3
5. 1  y 2 dx  1  x 2 dy, y (0) 
2
3.2 Resolver
dy
1.  ( x  y  1) 2
dx
dy 1  x  y
2. 
dx x y
dy
3.  1  e y  x 5
dx
dy
4.  2  y  2x  3
dx
dy
5.  ( x  y  1) 2  2 6.( x  y) y   2 x  2 y  3
dx
dy xy  3x  y  3
7. 
dx xy  2 x  4 y  8
8.x 3 dx  yx 3 y 2  1dy  y y 2  1dy  0
Darwin D. Peña G.
4. ECUACIONES HOMOGÉNEAS
Definición 4.1 Una función f (x, y) se llama homogénea de grado n si

f (tx, ty )  t n f ( x, y)
Ejemplo 4.1 f ( x, y)  x 4  x 3 y es homogénea de grado 4, ya que
 
f (tx, ty )  (tx ) 4  (tx ) 3 (ty )  t 4 x 4  x 3 y  t 4 f ( x, y)
Ejemplo 4.2 f ( x, y)  x 2  sin x cos y no es homogénea, ya que

f (tx, ty )  t 2 x 2  sin(tx ) cos(ty )  t n f ( x, y)
Una ecuación de la forma M ( x, y)dx  N ( x, y)dy  0 (*), en donde M, N tienen el mismo grado de
homogeneidad, puede reducirse a una ecuación de variables separables usando cualquiera de las
situaciones y  ux o bien x  vy, en donde u y v son nuevas variables dependientes. En particular si
elegimos y  ux, entonces dy  udx  xdu. Por tanto en la ecuación (*) se tiene que:
M ( x, ux)dx  N ( x, y)udx  xdu   0
Ahora por la homogeneidad de M y N se es posible escribir
x n M (1, u )dx  x n N (1, u )udx  xdu   0

M (1, u )dx  N (1, u )udx  N (1, u ) xdu  0
M (1, u )  N (1, u )u dx  xN (1, u )du  0
M (1, u )  N (1, u )u dx   xN (1, u )du
dx N (1, u )du

x M (1, u )  N (1, u )u 
dx N (1, u )du
 0
x M (1, u )  N (1, u )u 
dy y 2  x 2
Ejemplo 4.3 Resolver 
dx 2 xy
Solución:
   
2 xydy  y 2  x 2 dx  2 xydy  y 2  x 2 dx  0  hacemos y  ux
  
2ux 2 (udx  xdu )  x 2 u 2  x 2 dx  0  2ux 2 (udx  xdu )  x 2 u 2  1 dx  0 
   
2u (udx  xdu )  u  1 dx  0  2u  u  1 dx  2uxdu  0  u  1 dx  2uxdu
2 2 2
 2

dx 2udu
  2 Integrando tenemos que y 2  K x  x 2  y   K x  x 2
x 
u 1 
Darwin D. Peña G.
4.1 Forma alternativa
Una ecuación diferencial de primer orden y’=f(x,y) se llama homogénea si f(x,y) puede expresarse como
g(y/x), donde g es una función de una variables. Una ecuación diferencial homogénea y’=g(y/x) se
transforma en una ecuación de variables separables mediante el cambio de variable: v = y/x.
Ejemplo 4.4 Utilizando la forma alternativa en el ejemplo anterior tenemos

Solución
Si dividimos el numerador y el denominador del miembro derecho de la ecuación por x 2, tendremos:
2
 y
  1
dy  x 
 . Si ahora hacemos el cambio de variables u = y/x tenemos:
dx y
2
x
2
 y
  1
v2 1 v2 1 v 2  1  2v 2
 v' x  v   
y dy x dv
v   y  vx    x v  
x dx y 2v dx 2v 2v
2
x
dv v 1
2
2v dx 1
 x  2 dv    ln(v 2  1)  ln x  C  2  K1 x 
dx 2v v 1 x (v  1)
1 1
 (v 2  1)   v2  1  y2  K x  x2  y   K x  x2
K1 x K1 x
En cada punto: P(x,y) de una curva del plano, el ángulo formado por la tangente y la ordenada es
bisecado por la recta que une al punto con el origen. Halle la ecuación de la curva sabiendo que pasa
por el punto: (1,2)
Solución:
La figura (10) ilustra gráficamente la situación. Podemos relacionar el ángulo  con las variables del
problema así:
Figura 10
Darwin D. Peña G.
 
  2  
2 2
  
De las relaciones anteriores se sigue que:   2  . Puesto que: tan x     cot( x), se tiene que:
2  2
dy 1  tan 2 ( )
tan( )  p    cot(2 )  
dx 2 tan( )
dy y 2  x 2
Puesto que: tan(a) = y/x, tenemos la ecuación diferencial del problema, así;  La
dx 2 xy
ecuación diferencial hallada es homogénea y su solución general es la familia de circunferencias: x2 + y2
= Cx. En el punto (1,2), la circunferencia es: x2 + y2 = 5x. La figura (11) muestra
la gráfica de la función y la propiedad expresada en el enunciado del problema. El estudiante puede
comprobar que la ecuación de la recta tangente a la curva en el punto (1,2) viene dada por:
3
y  ( x  1)  2
4
Figura 11.
Darwin D. Peña G.
EJERCICIOS 4
4.1 Determine si las siguientes funciones son homogéneas, si lo es indique el grado de homogeneidad
x3 y  x2 y 2
1. f ( x , y ) 
x  8y
 x2 
2. f ( x, y )  cos 
 x  y 

3. f ( x, y )  x 1  y 1 
2
x
4. f ( x , y ) 
y2  x4  y4
4.2 Resuelva las siguientes ecuaciones diferenciales
dy y x dy x  3 y
1.   8. 
dx x y dx 3 x  y
dy y x 2 x 3  xy 2
2.   1 9 .8 x
dy
 8y 
dx x y 2 dx 64 x
2 x
dx 10.2 xdy  2 ydx  x 2  4 y 2 dx
3. y  x  4 ye y
dy
x2
dy 11 .( y  x ) dx   dy
4. 2 x 2  3 xy  y 2 ; sujeta a y (1)  2 y
dx
dy yx 2  2 xy 2  y 3
5.  ydx  ( x  xy )dy  0 12. 
dx 2 y 3  xy 2  x 3
6.( y 2  xy )dx  x 2 dy  0
dx x y
13.   1
dy y x
dy 2 xy  y 2
7.  sujeta a y (1)  2
dx 2 xy  x 2 y
dy
14.x  ye x  x
dx
Darwin D. Peña G.
Una ecuación diferencial de la forma
dy  ax  by  c 
 f  , aB - bA  o,
dx  Ax  By  C 
Siempre puede ser reducida a una ecuación homogénea mediante las sustituciones x = u + h, y
= v + k. en los problemas 15 y 16 use las sustituciones indicadas para reducir la ecuación a una con
coeficientes homogéneos. Resuélvala.
dy x - y - 3
 ; x  u  2, y  v - 1
15. dx x  y - 1
dy x  y - 6
16.  ; x  u  3; y  v  3
dx x-y
17. suponga que M (x, y) dx + N(x, y) dy = 0 es una ecuación homogénea. Demuestre que la sustitución
x = vy reduce la ecuación a una de variables separables.
18. suponga que M(x, y) dx + N(x, y) dy = 0 es una ecuación homogénea. Pruebe que las sustituciones
x= rcos  , y  rsen reducen la ecuación a una de variables separables.
19. suponga que M (x, y) dx + N (x, y) dy = 0 es una ecuación homogénea. Demuestre que la ecuación
tiene la forma alternativa dy/dx = G(x/y).
Darwin D. Peña G.
5. ECUACIONES DIFERENCIALES EXACTAS Y REDUCIBLES A EXACTAS
La condición necesaria y suficiente de que
M ( x, y)dx  N ( x, y)dy  0 (6)

M N
Sea exacta es  . Si (6) es exacta existe una función f para la cual
y x
f f
M ( x, y)dx  N ( x, y )dy  dx  dy (7)
x y
5.1 Método de solución
Dada la ecuación M ( x, y)dx  N ( x, y)dy  0

M N
Primero demuestre que  (8)
y x
Supongamos una solución f(x,y) = c
f
Por (7) tenemos que:  M ( x, y )  f  M ( x, y)x . Integramos con respecto a x
x
f ( x, y)   M ( x, y)dx  g ( y) (9) . Donde g(y) es la constante de integración. Derivando (9) con
f
 N ( x, y )
respecto a y por (7) tenemos que y
f  

y y  M ( x, y )dx  g ' ( y )  N ( x, y )  g ' ( y )  N ( x, y )   M ( x, y )dx (10)
y
integrando
     
 g ' ( y)    N ( x, y)  y  M ( x, y)dx dy  g ( y)    N ( x, y)  y  M ( x, y)dx dy (11)
  
sustituyendo (11) en (9) f ( x, y )   M ( x, y )dx    N ( x, y )   M ( x, y )dx  dy
 y 
  
luego la soluciòn f ( x, y )  c es  M ( x, y )dx    N ( x, y )   M ( x, y )dx  dy  c
 y 
Ejemplo 5.1 Resolver 2 xydx  ( x 2  1)dy  0
M N M N
Solución: Verifiquemos que    2x 
y x y x
Darwin D. Peña G.
f
Utilizando  2 xy , E integrando respecto a x, mientras que y se mantiene constante se tiene
x
 f   2xyx  f ( x, y)  x y  g ( y) . Derivando parcialmente con respecto a y
2
f
 x 2  g ' ( y )  x 2  1  x 2  g ' ( y )  x 2  1  g ' ( y )  1 integrando se tiene
y
g ( y)   y con lo que f ( x, y )  x 2 y  y
Luego la solución es f ( x, y)  c  x 2 y  y , es decir x2 y  y  c
5.2 Factores integrantes
Si (6) no es exacta se buscará un factor integrante que nos permita convertir (6) en exacta
M N

y x
, Función solo de x, entonces e 
f ( x ) dx
Si f ( x)  . Es un factor integrante para (6)
N
N M

x y
, Función solo de y, entonces e 
g ( y ) dy
Si g ( y )  .Es un factor integrante para (6)
M
Ejemplo 5.2  
Resolver x 2  y 2  x dx  xydy  0
M N
Solución:  2 y, y La ecuación no es exacta. Sin embargo
y x
M N

y X 2 y  y 1
dx
 f ( x ) dx x
f ( x)    e e  e ln x  x
N xy x
Es un factor integrante. Luego

     
x 2  y 2  x dx  xydy  0  x x 2  y 2  x dx  xxydy  0  x 3  xy 2  x 2 dx  x 2 ydy  0 Luego es
exacta y existe una función f(x,y) para la que
f f
 x 3  xy 2  x 2 y  x2 y
x y
f
Tomando  x 3  xy 2  x 2
x
x4 x2 y 2 x3
integrando se tiene f ( x, y )    g ( y)
4 2 3
Derivando parcialmente con respecto a y
Darwin D. Peña G.
f
 x 2 y  g ' ( y )  x 2 y  x 2 y  g ' ( y )  x 2 y  g ' ( y )  0 integrando
y
x4 x2 y2 x3
g ( y )  c1 luego f ( x, y )    c
4 2 3
x4 x2 y 2 x3
   c o bien 3 x 4  4 x 3  6 x 2 y 2  c
4 2 3
EJERCICIOS 5
5.1 Resolver las siguientes ecuaciones exactas

1. (2 x  3 y  4)dx  (3 x  4 y  5)dy  0
 1  1
2.  4 x 3 y 3  dx   3 x 4 y 2  dy  0
 x  y
 
3.  x  y  dx  2 xy  x 2  1 dy  0
2
sujeta a y (1)  1
4. ( seny  ysenx )dx  (cos x  x cos y  y )dy  0
5. (3 x 2  3xy 2 )dx  (3 x 2 y  3 y 2  2 y )dy  0
6. ( x  y cos x)dx  (sin x)dy  0
7. (e 2 y  y cos xy )dx  (2 xe 2 y  x cos xy  2 y )  0
8. (2 xy 2  3)dx  (2 x 2 y  4)dy  0
 1  dx
9.  x 2 y 3  2 
 x3 y 2  0
 1  9 x  dy
10. 5 y  2 x y´2 y  0
11. tan x  senx seny dx  cos x cos y dy  0
12. 3x cos 3x  sen3x  3dx  2 y  5dy  0

13. 1  2 x 2  2 y 
dy
dx
 4 x 3  4 xy
  
14. 2 ysenx cos x  y  2 y 2 e xy dx  x  sen 2 x  4 xye xy dy
2 2

   
15. 4 x 3 y  15x 2  y dx  x 4  3 y 2  x dy  0
1 1 y   y x 
16.   2  2 
2 
dx  
 ye  dy  0,
2 
y1  1
 x x x  y   x 2
 y 

17. x  y  dx  2 xy  x  1 dy  0
2 2

   
18. e x  y dx  2  x  ye y dy  0, y0  1
19. 4 y  2 x  5dx  6 y  4 x  1dy  0, y 1  2
5.2 Encuentre el factor integrante para las siguientes ecuaciones y resolverlas

Darwin D. Peña G.
 
1. 6 xydx  4 y  9 x 2 dy  0
2. (12 x y  24 x y )dx  (9 x 4  32 x 3  32 x 3 y  4 y )dy  0
3 2 2
3. ( x  y )dx  ( x  y )dy  0
4. (3 y 3  xy )dx  ( x 2  6 xy 2 )dy  0
5. (2 y  3xy 2 )dx  ( x  2 x 2 y )dy  0
6. ( xy  2 y 2 )dx  ( x 2  3xy )dy  0
7. (3x 2 y 2 )dx  4( x 3 y  3)dy  0
8. (1  y 2 )dx  ( x  x 2 )dy
6. ECUACIONES LINEALES
La forma general de la ecuación lineal de orden n es :

dny d n 1 y dy
a n ( x) n  a n 1 ( x) n 1  ...  a1 ( x)  a 0 ( x) y  g ( x) (12)
dx dx dx
Veamos el caso n = 1 es decir

dy
a1 ( x)  a 0 ( x) y  g ( x) dividiendo por a1 ( x)
dx
a1 ( x)dy a 0 ( x) g ( x)
 y es decir
a1 ( x)dx a1 ( x) a1 ( x)
dy
 P( x) y  f ( x), (13) con P( x), f ( x) continuas
dx
e
P ( x ) dx
Solucionar (13) nos obliga a utilizar un factor integrante que seria . Multiplicando (13) por el
factor integrante nos queda
e  P ( x )e  y  f ( x )e 
P ( x ) dx dy P ( x ) dx P ( x ) dx
(14)
dx
El término de la izquierda es la derivada del producto del factor integrante por la variable dependiente
por tanto (3) nos queda
d   P ( x ) dx 
 d e 
p ( x ) dx 
y   f ( x )e  y   f ( x )e 
P ( x ) dx P ( x ) dx
 e  integrando
dx    
 P ( x ) dx   P ( x ) dx dx 
e y   f ( x )e  dx  c  y  ce  e 
 P ( x ) dx

P ( x ) dx P ( x ) dx
 f ( x ) e 
dy
Ejemplo 6.1 Resolver x  y  2x sujeta a y(1)  0
dx
1
dy 1  x dx
Solución:  y2 el factor integrante es e  e ln x  x
dx x
Darwin D. Peña G.
Por consiguiente
dy 1 d c
x  x y  2 x  ( xy )  2 x  xy  x 2  c  y  x 
dx x dx x
c
Aplicando las condiciones iniciales y (1) = 0 0  1   c  1
1
1
Luego yx , 0 x
x
Ejemplo 6.2 Hallar una solución continua que satisfaga
dy 1, 0  x  1
 y  f ( x) en donde f ( x)  
dx 0, x  1
Y condición inicial y(0) = 0
Solución: Resolvemos el problema en dos partes

dy
Para  y  1, 0  x  1 . El factor integrante es ex luego,
dx
d x
dx
 
e y  e x  e x y  e x  c1  y  1  c1e  x Utilizando las condiciones iniciales tenemos
o  1  c1  c  1  y  1  e  x , 0  x 1
Para
dy
dx
 y  0, x 1
dx
 
d x
e y  0  e x y  c2  y  c2 e  x
 x
1  e , 0  x  1
Por consiguiente y   x

c 2 e , x 1
Como es continua, por lo visto en el calculo diferencial

1 1
 x
1  e , 0  x 1
lim y( x)  y (1)  1  e  c2e  c2  e  1 Luego y
x 1  x
(e  1)e , x 1
Es continua pero no diferenciable en x = 1
Problemas de crecimiento y decrecimiento
Son innumerables las aplicaciones que involucran el crecimiento o decrecimiento de una variable física,
económica o de cualquier otra índole. Si x(t) es una variable cualquiera, la tasa de variación de dicha
d
variable con respecto al tiempo, viene dada por: x(t ) .
dt
Cualquiera que sea la situación o problema, nos debemos preguntar: De qué depende la tasa de variación
de la variable? Se han desarrollado algunos modelos simples para problemas tales como: crecimiento de
una población, crecimiento de dinero en el tiempo, desintegración de materiales radioactivos, etcétera.
Ilustraremos algunos de dichos modelos.
Darwin D. Peña G.
Problema resuelto 5
Una población tiene, en un instante determinado, una cantidad inicial de habitantes: Xi. Se desea
determinar la cantidad de habitantes en cualquier instante suponiendo que la tasa de variación es
directamente proporcional al número de habitantes en todo instante.
Solución:
Sea x(f) la población en el instante: t > 0. La tasa de variación de la oblación, según el modelo,
dx(t )
es:  kx(t ) . Siendo k una constante de proporcionalidad. En consecuencia, se debe resolver el
dt
dx
problema de valor inicial:  kx x(0)  X i . La ecuación diferencial es de variables separables y su
dt
solución general viene dada por: x(t) = Cek'. Aplicando la condición inicial, se tiene: x(t )  X i e kt . El
valor de la constante se determina a partir de una condición del problema. Por ejemplo, si la población
se incrementa un 50% en un período de tiempo, se tiene:
x(l)=1.5Xi y por tanto: k = ln(1.5). En consecuencia, la población en todo instante viene dada por:
x(t )  X i e t ln(1.5)  X i (1.5) t
El modelo usado no es bueno ya que implicaría que la población crecería indefinidamente, es decir, no
daría cuenta de las limitaciones propias de la población, tales como espacio y alimentos.
En un instante determinado se tiene una muestra de radio: Xi. Determina la cantidad de sustancia en todo
instante sabiendo que el tiempo de vida media del elemento es 1700 años.
Solución:
El tiempo de vida media es el que se requiere para que una sustancia se desintegre en un 50%. Si x(t) es
la cantidad de sustancia remanente en el instante: t, se tiene:
dx
 kx x(0)  X i
dt
La solución del problema de valor inicial dado es:
x(t )  X i e  kt
ln( 2)
Ahora, puesto que: x(l700) = 0.5Xi , se tiene que la constante k viene dada por; k  .
1700
En consecuencia, la cantidad remanente en el instante: t viene dada por:
t

x(t )  X i (2) 1700
. Al cabo de 10 años, la cantidad remanente será
aproximadamente: x(10)  0.996 X i , . Es decir, se ha desintegrado el 0.4% de la muestra.
En el crecimiento de una población surgen circunstancias que impiden que su número exceda de cierto
máximo: M. En consecuencia, la tasa de variación de la población es directamente proporcional al
número de habitantes en todo instante y a la diferencia entre el máximo y la población instantánea.
Determine la población en todo instante sabiendo que inicialmente tenía N habitantes.
Solución: El modelo sugerido para la tasa de variación de la población nos conduce a la ecuación
 kx(t )M  x(t )
dx(t )
diferencial:
dt
Darwin D. Peña G.
La ecuación diferencial recibe el nombre de ecuación logística y es la más usada para estimativos de
crecimiento de poblaciones. Con base en lo estudiado previamente, la ecuación diferencial es de
variables separables, pero también es una ecuación de Bernoulli. Escribimos convenientemente la
ecuación, así:
dx 1 1 1 
 kdt . Descomponiendo en fracciones parciales, resulta:  dx  kdt .
x( M  x) M  x M  x 
Integrando, resulta: ln(x) - ln(M - x) = kt + C
x
La expresión anterior se puede escribir en la forma:  Ae kt .La constante A se determine con la
M x
N
condición inicial, así: A = . Finalmente, la población en todo instante viene dada por:
M N
MN
x(t )  . La constante k se determina a partir de un dato del problema, por ejemplo, sí
N  ( M  N )e  kt
la población se duplica al cabo de diez años, la constante k viene dada por:
 2( M  N ) 
k  0.1 ln   . La gráfica de la figura (12) corresponde a la población en todo instante con los
 M  2N 
siguientes datos: N=5 M=20 x(10)=10
Figura 12
Un país tiene en circulación papel moneda por valor de mil millones de pesos. Los saldos en los bancos
suman cinco millones diarios. El gobierno decide introducir una nueva moneda cambiando todos los
billetes antiguos que lleguen a los bancos por otros nuevos. Determine el tiempo aproximados para que
el papel moneda circulante sea dinero nuevo en un 90%.
Solución: Sea x(f) la cantidad de dinero nuevo en un instante t cualquiera y y(t) = 1000 - x(t) la
cantidad remanente de dinero viejo. La tasa de variación de dinero nuevo, es decir, la velocidad con que
se incrementa el dinero nuevo debe ser proporcional a la cantidad de dinero viejo presente en todo
instante. Matemáticamente, se tiene que resolver el problema de valor inicial:
 k 1000  x(t )
dx(t )
x(0)  0
dt
Darwin D. Peña G.
La ecuación diferencial es de variables separables y también es lineal. Resolviendo como lineal, se tiene:
dx
 kx  1000k . El factor integrante es:   e kt y por tanto, la solución general es:
dt
x(t )  Ce kt  1000. Aplicando la condición inicial resulta que la cantidad de dinero nuevo en todo
 
instante viene dada por; x(t )  1000 1  e  kt . Para calcular k partimos del hecho de que el primer día se
introducen cinco millones de moneda nueva, con lo que resulta:
 k
 k k
5  1000 1  e  1  e  0.005  e  0.995  k  ln(1.005)  k  0.005
 
Así las cosas, la cantidad de dinero nuevo en todo instante es: x(t )  1000 1  e 0.005(t ) . Por tanto, el
tiempo en días para que el dinero se renueve en un 90% se calcula de la siguiente
 
manera: 900  1000 1  e 0.005(t )  1  e 0.005(t )  0.9  e 0.005(t )  0.1  e 0.005(t )  10  t  460
EJERCICIOS 6
6.1 Resolver
1. x 3
dy
dx
 
 2  3x 2 y  x 3
2. y ln ydx  ( x  ln y )dy  0
dy
3. ( x  2)  y  2 ( x  2) 2
dx
di
4.  6i  10 sen2t
dt
5. ydx  ( xy  x  3 y )dy  0
6. dr  (2r cot   sen2 )d  0
7. y (1  y 2 )dx  2(1  2 xy 2 )dy
8. xdy  2 ydx  ( x  2)e x dx
dy
9.  y  2  2 x
dx
10. (.2  y 2 )dx  ( xy  2 y  y 3 )dy  0
En los problemas 11-16, resuelva la ecuación diferencial dada sujeta a la condición inicial que se indica.
dy
11.  5 y  20, y(0)  2
dx
12. y  2 y  x(e3 x  e2 x ), y(0)  2
di
13. L  Ri  E; L, R, y E son constantes; i(0)  i0
dt
dy
14. y  x  2 y 2 , y(1)  5
dx
Darwin D. Peña G.
15. y  (tan x) y  cos 2 x, y(0)  1

dQ
16.  5 x 4Q, Q(0)  7
dx
dy
17.  y  sinh x
dx
18. ydx  ( x  2 xy 2  2 y)dy  0
19. ydx  ( ye y  2 x)dy
dr
20.  r sec   cos 
d
dP
21.  stP  P  4t  2
dt
dy
22. ( x  2) 2  5  8 y  4 xy
dx
dy
23. ( x 2  1)  2 y  ( x  1) 2
dx
24. xy  (1  x) y  e  x sin 2 x
25. cos 2 x sin xdy  ( y cos 3 x  1)dx  0
26. (1  cos x)dy  (2 y sin x  tan x)dx  0
6.2 Obtenga una solución continua
dy 1, 0 x3
1.  2 y  f ( x) en donde f ( x)   , y (0)  0
dx 0, x3
 x, 0  x 1

2. 1  x 2  dy  2 xy  f ( x) en donde f ( x)   , y (0)  0
dx  x x 1
Darwin D. Peña G.
7. ECUACIONES REDUCIBLES A LINEALES
7.1 ECUACIÓN DE BERNOULLI
Existen ecuaciones no lineales que se pueden reducir a una ecuación lineal, por ejemplo la ecuación de
dy
Bernoulli  P( x) y  f ( x) y n (15)
dx
En donde n es un número real cualquiera. Para n ≠ 0 y n ≠ 1, la sustitución w = y1-n lleva a una
ecuación lineal de la siguiente manera. Derivamos
dw dy y n dw dy
 (1  n) y  n o bien  Luego sustiyendoen (15)
dx dx (1  n) dx dx
y n dw dw (1  n) (1  n)
 P( x) y  f ( x) y n o bien  P( x) y  f ( x) y n es decir
(1  n) dx dx y n
yn
dw
 (1  n) P( x) w  (1  n) f ( x) (16)
dx
Que es la ecuación que necesitamos resolver.
dy 1
Ejemplo 7.1 Resolver  y  xy 2 ´
dx x
Solución
1
Podemos identificar que P( x)  , f ( x)  x y que n = 2 luego en (16) se tiene
x
dw 1 dw 1
 (1  2) w  (1  2) x es decir  w  x
dx x dx x
El factor integrante de esta ecuación es x-1 con lo que
d 1
dx
 
x w  1 Integrando obtenemos x 1 w   x  c  w   x 2  cx
1 1 1
Como w  y 1 , y con lo que   x 2  cx  y 
w y  x  cx
2
7.2 ECUACIÓN DE RICATTI
Es otra ecuación no lineal que puede ser reducida a una ecuación lineal. Es decir, realizaremos el
siguiente procedimiento:
Ecuaciòn de Ricatti  Ecuaciòn de Bernoulli  Ecuaciòn lineal
dy
La ecuación es  P ( x )  Q( x ) y  R ( x ) y 2 (17)
dx
Para resolverla suponemos una solución particular conocida y1 de tal forma que y = y1+ u es una
solución de (17), con esto reducimos la ecuación de Bernoulli, veámoslo. Si y es una solución entonces
dy dy1 du
 
dx dx dx
Sustituyendo den (1) obtenemos lo siguiente
Darwin D. Peña G.
dy1 du
  P( x)  Q( x)( y1  u )  R( x)( y1  u ) 2
dx dx
dy1 du
  P ( x)  Q ( x) y1  Q( x)u  Rxy 2 1  2 Rxy 1u  Rxu 2
dx dx
 P ( x)  Q ( x) y1  Rxy 12   Q ( x)u  2 R ( x) y1u  R ( x)u 2
dy1 du

dx dx   
dy 1
dx
dy1 du dy1
   Q ( x)u  2 R ( x) y1u  R ( x)u 2
dx dx dx
du
 Q ( x)u  2 R ( x) y1u  R ( x)u 2
dx
du
 (Q ( x)  2 R ( x) y1 )u  R ( x)u 2 es una ecuaciòn de Bernoulli
dx
Y haciendo la sustitución w = u-1obtenemos la ecuación lineal
 Q( x)  2 R( x) y1 w   R( x)
du
(18)
dx
Resolviendo (18) Hallamos el factor integrante.
P( x)  Q( x)  2 R( x) y1  e P ( x )  e Q ( x )  2 R ( x ) y1
 Q( x)  2 R( x) y1 we Q ( x )  2 R ( x ) y1   R( x) e Q ( x )  2 R ( x ) y1
du Q ( x )  2 R ( x ) y1
e
dx
d
( we Q ( x )  2 R ( x ) y1 )   R( x) e Q ( x )  2 R ( x ) y1  we Q ( x )  2 R ( x ) y1    R( x) e Q ( x ) 2 R ( x ) y1 dx  c
dx
  R( x) e
Q ( x )  2 R ( x ) y1
dx  c e Q ( x )  2 R ( x ) y1
w , w  u 1  u 
e Q ( x )  2 R ( x ) y1   R( x) e
Q ( x )  2 R ( x ) y1
dx  c
dy
Ejemplo 7.2  2  y  y 2
Resolver con y1  2
dx
Solución: De la ecuación podemos identificar que P(x) = -2, Q(x) = -1, R(x) =1
  1  2(1)(2) w  1 
du du
Sustituyendo en (18) obtenemos  3w  1
dx dx
El factor integrante seria e3x entonces, e 3 x
du
dx
 e 3 x 3w  e 3 x 
d 3x
dx
e w  e 3 x Integrando  
1 1
e 3 x w  3e 3 x  c  w  3  ce 3 x , Pero w = u-1  3  ce 3 x  u  Luego la
u  3  ce 3 x
1 2ce 3 x  5
solución y = y1 + u es: y  2  y 
 3  ce 3 x ce 3 x  3
Darwin D. Peña G.
EJERCICIOS 7
7.1 Resolver las siguientes ecuaciones de Bernoulli
1. x 2
dy
dx
 y 2  xy 
2. 3 1  x 2
dy
dx

 2 xy y 3  1  
1 dy 3 dy
3. y 2
 y 2  1, y (0)  4 4. 2 xy  y 2  2x 2
dx dx
dy 1 dy y
5.  y  y2 y( )  1 6.  (4  x 2 y 3 )
dx 2 dx x
dx x dy 1
7. 2   x 3 cos y  0 8. x 2  2 xy  3 y 4 y (1) 
dy y dx 2
7.2 Resolver las siguientes ecuaciones de Ricatti

dy
dx
1
1.  2 x 2  y  2 y 2 ,
x
y1  x
dy

2.  e 2 x  1  2e x y  y 2 ,
dx
 y1  e x
7.3. Sin resolverlas clasifique cada una de las ecuaciones siguientes como: separable, homogénea,
exacta, lineal, de Bernoulli, de Ricatti o de Clairaut.
dy 1 dy x  y
a).  b). 
dx y  x dx x
2
 dy  dy dy 1
c).    2 y  2 x d). 
 dx  dx dx xx  y 
dy y  y
2
dy
e).  f).  4  5y  y2
dx x 2  x dx
x
g). ydx  y  xy dy 2

 h). x
dy
dx
 ye y  x
i). xyy´ y 2  2 x j). 2 xyy´ y 2  2 x 2
dx  3  Lnx dy

 2y 
k). ydx  xdy  0 l).  x 2  2
 x 
dy x y y dy
  1  e2x  y  0
3 2
m). n). 2
dx y x x dx
o). y  xy´ y´3 p). y´5 y 2  3x 4  2 xy
2
Darwin D. Peña G.
8. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE PRIMER ORDEN
Aplicaciones al vaciado de recipientes 9.

Consideremos un recipiente que tiene la forma de un sólido de revolución como el que se muestra en la
figura 13. A partir de un instante determinado se abre una válvula colocada en la parte inferior del
recipiente, con lo que empieza a vaciarse. Se pide calcular el tiempo de vaciado del recipiente.
Figura 13.
Solución:
Supongamos que el orificio correspondiente a la válvula tiene un área: a que es mucho menor que las
dimensiones del recipiente.
En todo instante, el volumen que sale por unidad de tiempo es proporcional al área del orificio y a la
velocidad con que sale el líquido, esto es, la tasa de variación del volumen por unidad de tiempo viene
dV (t )
dada por:  kav(t )
dt
En la expresión anterior se tiene que la velocidad en todo instante es: v(t) y el volumen en todo instante
es: V(t). La constante depende básicamente de la forma del orificio. Ahora bien, la rapidez con que sale
el líquido se determina aplicando la ley de la conservación de la energía para una gota de líquido, así:
Supongamos que en el instante en que una gota de líquido de masa: m se encuentra en el nivel: y, en ese
momento su energía potencial medida respecto al orificio viene dada por: Ep = mgy. Al nivel del
1
orificio, la energía de la gota es cinética, así: Ec = mv 2 . De lo anterior se sigue que la velocidad con
2
que sale es:
v(t )  2 gy . De otro lado, la variación del volumen con respecto al tiempo determina tomando un
dV (t ) dy
infinitésimo de volumen, así:  A( y ) . Con todo lo anterior resulta el problema de valor inicial
dt dt
siguiente que relaciona la altura del nivel del líquido en todo instante:
dy 1
A( y)  ka 2 g y 2 y(0)  H
dt
En la ecuación diferencial obtenida, A{y} es la sección transversal en todo instante y H es el nivel
inicial del líquido. En algunos casos el área del orificio no es constante y su magnitud puede depender
de la altura del nivel del líquido, esto es, el área del orificio puede ser controlada.
Darwin D. Peña G.
Problema resuelto 10
Un recipiente cónico con la base en la parte superior tiene un orificio de área constante: a en la parte
inferior. Calcule el tiempo de vaciado.
Solución: La figura 14 muestra la situación en un instante cualquiera: t. De la geometría del problema
se sigue que la sección tiene un área: A  x 2 . De otro lado, a partir de la figura 15, por semejanza de
triángulos, se tiene la relación: Hx = Ry. Con lo anterior, la sección tiene un área que depende de la
R 2 y 2
altura del nivel del líquido en todo instante, así: A( y )  .Por simplicidad, definimos la
H2
R2
constante: B   kH 2 2 g . Con lo anterior resulta el problema de valor inicial de variables
a
3
separables: By 2 dy  dt  0 y(0)  H . El tiempo de vaciado resultante de resolver el problema es:
2R2
H
tv  . Particularmente, sí las dimensiones del recipiente son: H = 2metros, R = 1metros, y,
5ka 2g
además: a = 5cm2 k = 1 g = 9.8 \m/seg2, el tiempo de vaciado será: t v  13.4 minutos.
Figura 14 figura 15
Un depósito en forma de canal, tal como lo ilustra la figura 7 tiene un orificio de área a practicado en el
fondo. Calcule el tiempo de vaciado.
Figura 16
Darwin D. Peña G.
Solución: Teniendo en cuenta la figura 16, el diferencial de volumen viene dado por: dV  2Lxdy . De
otro lado, puesto que el punto: P{x,v) está sobre la circunferencia de radio: R, se tiene que:
( x  R) 2  y 2  R 2 . Con base en lo anterior, el diferencial de volumen es una función de la variable
dependiente: y. Así: dV  2L 2 Ry  y 2 dy . En consecuencia, resulta el problema de valor inicial:
dy 1
2 L 2 Ry  y 2  ka 2 g y 2 y(0)  R
dt
2l
Por simplicidad hacemos: B  , con lo que la ecuación diferencial a resolver es:
ka 2 g
 B 2R  y dy  dt . Por integración directa, la solución general de la ecuación diferencial es:
2 3
t  C  B(2 R  y ) 2 Aplicando la condición inicial: y(0) = R, se encuentra que el tiempo de vaciado
3
4 RL R
es: t v 
3ka 2 g

Determine la velocidad en todo instante de una partícula que se deja caer desde una distancia muy
grande de la superficie cíe la tierra. Suponga que no hay rozamiento.
Solución: La fuerza con que la tierra atrae a la partícula obedece a la Ley de la gravitación universal,
esto es, si denotamos por y (t ) a la distancia en todo instante entre la partícula y el centro de la tierra, tal
como lo ilustra la figura 11, la segunda Ley de Newton queda en la forma:
dv kmM
m   2 . Donde M es la masa de la tierra. La constante: k se determina con la condición de
dt y
que la fuera en el momento de tocar la superficie es igual al peso de la partícula, esto es:
kmM gR 2 dv gR 2
 mg k  . En consecuencia, la ecuación diferencial del movimiento es:  
R2 M dt y2
Figura 17 figura 18
Aplicando la regla de la cadena, la ecuación diferencial se puede escribir en la forma:

dv gR 2
v  2  vdv   gR 2 y 2 dy . Integrando, se obtiene la solución general: v 2  2 gR 2 y 1  C La
dy y
constante de integración se encuentra con base en la condición inicial, así: supongamos que el cuerpo se
Darwin D. Peña G.
dejó carca desde una distancia: R de la superficie, siendo R el radio de la tierra; en tal caso, la condición
inicial es: v(2R)  0 . La constante de integración es:
C   gR En consecuencia, la velocidad en todo instante es: v(t )  2 gR 2 y 1  gR
La velocidad con que la partícula llega a la superficie viene dada por: v(r )  gR

Un hombre y su paracaídas están cayendo con una velocidad: Vi en el instante en que se abre el
paracaídas. Sí el aire opone una resistencia proporcional al cuadrado de la velocidad instantánea,
encuentre una expresión para la velocidad en todo instante.
Solución: Tomando como referencia el punto en el que se abre el paracaídas, como lo ilustra la figura
dv
14 la segunda Ley de Newton establece que: m  mg  kv . La ecuación se puede escribir en la
dt
m dv mg mg
forma:   v 2 . Definimos la velocidad límite como: V L  . En consecuencia, el
k dt k k
dv g
problema de valor inicial asociado a la velocidad es:  2 dt v(0)  Vi . Integrando, se
VL  v
2 2
VL
1 V  v  g
obtiene: ln  L  tc
2V L  V L  v  V L2
Evaluando la constante y despejando, se encuentra que la velocidad en todo instante viene dada por:
2 gt

V L  Vi  Vi  V L e VL
v(t )  V L 2 gt
. Claramente se observa que: VL  lim v(t )
 t 
V L  Vi   VL  Vi e VL
A manera de ejemplo, sí e! hombre y su paracaídas tiene una masa de 120 Kilogramos y la constante de
fricción es numéricamente igual a 240, la velocidad límite será: VL  5mts / seg . La figura 19 muestra
la velocidad en todo instante cuando la velocidad inicial es: Vi = 50 mts/seg. .
Figura 19
Puede verse que la velocidad límite se alcanza un segundo después de abrir el paracaídas. Se deja como
ejercicio al estudiante que determine la posición en todo instante.
Darwin D. Peña G.
Capítulo II
Darwin D. Peña G.
9. ECUACIONES DIFERENCIALES LINEALES DE ORDEN SUPERIOR
dny d n 1 y
a n ( x) n
 a n 1 ( x ) x 1
dy
   a1 ( x)  a 0 ( x) y  g ( x) 19
dx dx dx
sujeta : y ( x0 )  y 0 , y ' ( x0 )  y ' 0 ,...., y n 1 ( x0 )  y 0n 1
d2y
Para el caso n  2 tenemos : a 2 ( x) 2
dy
 a1 ( x)  a 0 ( x) y  g ( x) 20
dx dx
y ( x0 )  y 0 , y ' ( x0 )  y ' 0
Teorema: (Existencia y unicidad) Sean an x , an1 x ,....a1 x , a0 x  y g(x) continuas en un intervalo

I y sea an x   0 para todo x en este intervalo. Si x  x0 es cualquier punto de este intervalo,
entonces existe una solución y(x) del problema de valor inicial (19) y esa solución es única.
9.1 Dependencia e independencia lineal de funciones
Definición 9.1 si f 1 , f 2 ,... f n son funciones entonces f 1 , f 2 ,... f n son linealmente independientes si
existen constantes c1 , c 2 ...c n todas nulas tal que c1 f 1  c 2 f 2  ...c n f n  0  c1  c 2  ....c n  0 .
Definición 9.2 f 1 , f 2 ,... f n son linealmente dependientes si existen c1 ,c 2 ...c n , no todas nulas tales
que: c1 f 1  c 2 f 2  c n f n  0  no todas c1 , c 2 c n son cero.
Ejemplo 9.1 Verificar si f 1 ( x)  sen2 x, f 2 ( x)  senx. cos x son linealmente independientes o

linealmente dependientes.
Solución
c1 f1  c2 f 2  0  c1sen2 x  c2 senx. cos x  0
c1 2 senx. cos x  c2 senx. cos x  0  c1  1; c2  2
2senx. cos x  2senx. cos x  0
f 1 ( x)  e x , f 2 ( x)  e 2 x , f 3  ( x)  e 3 x
Ejemplo 9.2
c1 e x  c 2 e 2 x  c3 e 3 x  0
¿Como podemos hallar estas constantes para determinar la dependencia o independencia lineal?
9.2 Wronkiasno
Definición 9.3 sea f 1 , f 2 , , f n , con n-1 derivadas entonces el wronkiasno se define:

Darwin D. Peña G.
 f1 f2  fn 
f' f '2  f n 
w   f1 , f 2 , f n   
1
    
 n 1 n 1

 f 1 f2 f nn 1 
Si w = 0 entonces f 1 , f 2 ,  f n son linealmente dependientes

Si w  0 entonces f 1 , f 2 ,  f n son linealmente independientes
Ejemplos 9.3
1.
 sen 2 x senx. cos x 
wsen 2 x, senx. cos x    
2 cos x sen x  cos 2 x 
2
 
 sen 2 x  sen 2 x  cos 2 x  2 cos 2 xsenx cos x 
 
 2 senx. cos x cos 2 x  sen 2 x  2 sen. cos x cos 2 x  sen 2 x  0  
2.
e x e 2x e 3x 
 

w e x , e 2 x , e 3x   e x 2e 2 x 3e 3 x 
e x 4e 2 x 9e 3 x 

  
 e x 18e 5 x  12e 5 x  e x 9e 5 x  4e 5 x  e x 3e 5 x  2e 5 x   
 6e 6x
 5e 6x
e 6x
 2e 6x
0
x
x
3. e , e , e
4x
   x   , son linealmente independiente. Ver la figura 20.
Figura 20
0  x   , son linealmente independiente. Ver la figura 21.

1
4. x 2
, x2
Darwin D. Peña G.
Figura 21
9.3 Elaboración de una segunda solución a partir de una conocida.

Reducción de orden (Caso general).
dny d n 1 y
a n ( x) n  a n 1 ( x) x 1
dy
   a1 ( x)  a 0 ( x) y  0 22
dx dx dx
Para el caso n  2 tenemos
d2y dy
a 2 ( x) 2  a1 ( x)  a 0 ( x) y  0 23
dx dx
y ( x0 )  y 0 , y ' ( x0 )  y ' 0
Supongamos que de (23) conocemos una solución y1 (x) entonces podemos encontrar y2 (x) de la
siguiente manera
a 2 ( x) y' ' a1 ( x) y ' a0 ( x) y

   0  y' ' P( x) y'Q( x) y  0 (24)
a 2 x  a 2 ( x) a 2 ( x)
e 
 P ( x ) dx
y 2  y1  dx
Luego la segunda solución esta dada por: y12
Demostración: Suponemos que y1 es solución de (24) en algún intervalo I y que y1  0 , para todo x
en el intervalo. Si definimos y  u ( x) y1 tenemos
y   u ( x) y1  y1u ( x)
y   u ( x) y1  2u ( x)  y1u ( x)1
Reemplazando en (21)
Darwin D. Peña G.
u ( x) y1  2u ( x)  y1u ( x)1  P ( x)u ( x) y1  y1u ( x)  Q( x)u ( x) y1   0

u ( x) y1  P ( x) y1  Q( x) y1   2u ( x) y1  u ( x) y1  P( x)u ( x) y  0
2 y1 u ( x)
  P( x)  0
y1 u ( x)
2 y1 u ( x)
 y1   u ( x)    P( x)dx
 
ln y1 u ( x)    P( x)dx
2
e 
 P ( x ) dx
e   P ( x ) dx y2
u ( x)   2
dx; u ( x)   y 2  y1  dx
yi y1 y12
e 
 P ( x ) dx
y  y1  dx
Queda como ejercicio para el estudiantes mostrar que 2 y12 , es solución de (21) y
además que y1, y2 son linealmente independientes.
Ejemplo 9.4: Resolver. Hallar la segunda solución de: (1  x) y' ' xy ' y  0 y1  x
Solución:
e 
 P ( x ) dx
(1  x)e  x
y2  x dx  y 2  x  x2 queda como ejercicio.
y12
Ejemplo 9.5: a) Sea y1 y y 2 dos soluciones de a2 ( x) y  a1 ( x) y  a0 ( x) y  0 muestre que:

'' '
dw
a2 ( x )  a1 ( x ) w  0
dx
Solución:
a2 ( x) y1  a1 ( x) y1  a0 ( x ) y1  0
'' '
a2 ( x) y2  a1 ( x ) y 2  a0 ( x) y 2  0
'' '
Multiplicando una ecuación por y 2 y la otra por  y1 y restándolas se obtiene lo siguiente:

a2 ( x) y2 y1  a1 ( x) y2 y1  a0 ( x) y2 y1  0
'' '
 Resolviendo el sistema,

 a 2 ( x ) y y
1 2
''
 a1 ( x ) y y
1 2
'
 a0 ( x ) y y
1 2  0
 a2 ( x)( y2 y1  y1 y2 )  a1 ( x)( y2 y1  y2 y1 )  0
'' '' ' '
y1 y2
w  w  y1 y2'  y2 y1'  w  ( y2 y1'  y1 y2' )
y1' ' y2'
Darwin D. Peña G.
dw
 y1 y 2''  y 2' y1'  ( y 2 y1''  y1' y 2' )
dx
dw dw
 y1 y 2''  y 2 y1''   ( y 2 y1''  y1 y 2'' )
dx dx
Remplazando en la ecuación se tiene:
dw dw
 a2 ( x )  a1 ( x) w  0  a2 ( x)  a1 ( x) w  0
dx dx
b) Deduzca la formula de Abel
a ( x)
  a12 ( x ) dx
w  ce
Tenemos que:
dw dw a ( x)
a2 ( x)  a1 ( x) w  0   1 dx
dx w a2 ( x)
dw a ( x) a ( x)
 w
  1
a2 ( x )
dx  ln w  c    1
a2 ( x )
dx
a1 ( x )
 dxc
a ( x)
ln w    1 dx  c  e ln w  e a2 ( x )
a2 ( x )
a1 ( x ) a1 ( x )
 dx  dx
c
we a2 ( x )
e  w  ke a2 ( x )
9.4 Conjunto fundamental de soluciones
Definición 9.4 Se llama conjunto fundamental de soluciones en un intervalo I a cualquier conjunto

,f1,f2,…,fn de n soluciones linealmente independientes de la ecuación diferencial homogénea de orden n
(22) en algún intervalo I.
Definición 9.5 Si ,f1,f2,…,fn es un conjunto fundamental de soluciones de la ecuación diferencial

homogénea de orden n (22) en algun intervalo I y si F es cualquier solución de (22) entonces es
posible encontrar c1 , c 2 , c3 ,..., c n tales que F = c1 f 1  c 2 f 2  ...c n f n
Darwin D. Peña G.
Ejemplo 9.5:
d3y d2y dy
3
 6 2
 11  6y  0
dx dx dx
y1  e x
y2  e 2 x
y3  e 3 x ;
Son solución de la ecuación diferencia l
1) e x  6e x  11e x  6e x  0
2) 8e 2 x  24e 2 x  22e 2 x  6e 2 x  0
3) 27e 3 x  54e 3 x  33e 3 x  6e 3 x  0
Veamos que la combinación lineal también es solución:
y  c1 e x  c 2 e 2 x  c 3 e 3 x
y '  c1 e x  2c 2 e 2 x  9c 3 e 3 x
y ' '  c1 e x  8c 2 e 2 x  27c 3 e 3 x
Definición
     
c1 e x  8c 2 e 2 x  27c 3 e 3 x  6 c1 e x  4c 2 e 2 x  9c 3 e 3 x  11 c1 e x  2c 2 e 2 x  3c 3 e 3 x  6 c1 e x  c 2 e 2 x  c 3 e 3 x  0
 6c1 e x  24c 2 e 2 x  54c 3 e 3 x  11c1 e x  22c 2 e x  22c 2 e 2 x  33c 2 e 3 x  6c1 e x  6c 2 e 2 x  6c 3 e 3 x  0
00
9.6 Sean f1,f2,…,fk solución de (22) en algún intervalo I y sea fp cualquier solución de la ecuación de
la ecuación no homogénea (19) en el mismo intervalo. Entonces F = c1 f 1  c 2 f 2  ...c k f k  c p f p
es también solución de la ecuación no homogénea en el intervalo para constantes cualesquiera
c1 , c 2 , c3 ,..., c k
9.5 Solución de ecuaciones diferenciales homogéneas con coeficientes constantes

dny d n1 y
 a0 ( x) y  0 25
dy
an ( x) n  an1 ( x) x 1
   a1 ( x)
dx dx dx
Para n=2 tenemos que
ay ' 'by 'cy  0
a, b, c  son constantes
y  e mx ; y '  me mx ; y ' '  m 2 e mx  am 2 e mx  bme mx  ce mx  0  e mx am 2  bm  c  0  
am  bm  c  0
2
 b  b 2  4ac Ecuación auxiliar

m
2a
Teniendo en cuenta el discriminante de solución de la ecuación cuadrática tenemos los siguientes casos
y su respectiva solución:
Darwin D. Peña G.
m1 , m2  R ; m1  m2 m1 , m2  R ; m1  m2
 2  2
Caso I: b  4ac  0 Caso II: b  4ac  0
 
 y  c1e  c2 e  y  c1e  c2 xe
m1x m2 x m1x m1x
m1 , m2  C
 2
b  4ac  0
Caso III: m     i

y  ex (c1 cos  x  c2 senx)
Ejemplo 9.6: Resolver y' ' '3 y' '4 y  0

Solución:
m 3  3m 2  4  m  1m  2  0  m1  1; m2  m3  2
2
y  c1e x  c2 e 2 x  c3 xe 2 x
Ejemplo 9.7: Resolver y' '  4 y'13 y  0 y(0)  1, y' (0)  2
Solución:
4  16  52 4  36 4  6i
m 2  4m  13  0  m  m m
2 2 2
  2;  3
m  2  3i
y  e 2 x c1 cos 3 x  c2 sen3 x ;
p ara las condicione s inciales y (0)  -1
 1  e 0 c1 cos(0)  c2 sen(0)   1  c1
y  e 2 x  cos 3 x  c2 sen3 x  
y '  e 2 x 3sen3 x  3c2 cos 3 x   2e 2 x  cos 3 x  c2 sen3 x 
 2  e 0 3sen(0)  3c2 cos(0)   2e 0  cos(0)  c2 sen(0) 
 2  3c2  2  c2  0
soluciòn : Y  e 2 x cos 3 x
Darwin D. Peña G.
9.6 Operador Diferencial
dny
n
 Dn y
dx
dny d n 1 y dy
a n ( x) n  a n 1 ( x) n 1    a1  a0 ( x) y  0
dx dx dx
a n D n y  a n 1 D n 1 y    a1 Dy  a0 y  0
a 
n D  a D    a D  a y  0
n
 n 1
 
n 1
1 0
P( D)
¿Es conmutativo el operador diferencial?
9.7 Operador Anulador:
1. El operador anulador de las funciones polinómica 1, x, x 2  x n1 , se anula con un operador anulador de
mayor grado Dn .
2. El operador anulador de las funciones exponenciales ex , xe 2 x , x 2ex ,, x n1ex , donde n es un grado
mayor del grado de la variable x y  es el término que acompaña a la variable en el exponencial.
D   n
3. El operador anulador de las funciones ex senx, xe 2 x , x 2 ex senx, , x n1ex senx
y 
ex cos x, xe 2 x cos x, x 2 ex cos x, , x n1ex cos x es: D 2  2D   2   2  n
Ejemplo 9.8: Resolver f ( x)  4 x 2 e 5 x


Solución:   5; n  3  D  5 4 x 2 e 5 x  0
3

Ejemplo 9.9: Resolver f ( x)  e 4 x
Solución:  
D  4 e 4 x  0  De4 x  4e 4 x  0  4e 4 x  4e 4 x  0
Ejemplo 9.10: Resolver f ( x)  e x  2 xe 3 x

Solución: D  1e x  0  D  3 2 xe 3 x   0  D  1D  3 e x  2 xe 3 x   0
2 2
Ejemplo 9.11: Resolver f ( x)  e 2 x sen3x

Solución:
 
  2;   3; n  1  D 2  4D  13 e 2 x sen3x  0 
Darwin D. Peña G.
si   0, n  1, El operador se transforma en
senx
D 2
2  
cos x
9.8 Método por coeficiente indeterminado
an ( x) y n  an 1 ( x) y n 1    a1 y ' a0 y  g ( x) 26
Recordemos que: Para la ecuación diferencial lineal homogénea tenemos una solución general o
complementaria y para la ecuación diferencial lineal no homogénea tenemos una solución general +
solución particular. y  yg  y p
El método da resultado, siempre que g x  tenga una de las formas tales que podamos utilizar los
operadores anuladores.
Ejemplo 9.12: Resolver y' '3 y'2 y  4 x 2

Solución:
Resolvemos la ecuación diferencial
y' '3 y'2 y  0  m 2  3m  2  0  m  2m  1  0  m1  2; m2  1
Raíces distinta  y g  c1e 2 x  c 2 e  x
Para solución particular buscamos el operador que anula a g (x).
3
 
g ( x)  4 x 2  D 4 x 2  0  m 3  0  m1  m2  m3  0
y p  c1e  c 2 xe  c3 x 2 e 0
0 0
y p '  c 2  2c3 x
y p ' '  2c 3
2c3  3c 2  2c1  0  6c3  2c 2  0  2c3  4  c3  2
y p  7  6x  2x 2
Y  c1e  2 x  c 2 e  x  7  6 x  2 x 2
.
n 1
Demuestre que la ecuación diferencial a n y  a n 1 y    a1 y'a 0 y  k
n
Ejemplo 9.13:
k
donde k es constante, a 0  0 tiene una solución particular y p  .
a0
Solución:
Dk   0  m  0  y p  c1e 0  c1
Ahora, y p  c1 , y ' p  0,..., y 0n  0 
k
0  0    0 a0 c1  k  y p 
a0
Darwin D. Peña G.
EJERCICIO 9
9.1 Determina la dependencia e independencia lineal de las siguientes funciones.

1. f1 x   cos 2 x, f 2 x   1, f x   cos 2 x
2. f1 x   e x , f 2 x   e  x , f x   cosh x
3. f1 x   5, f 2 x   sen 2 x, f x   cos 2 x
Halle la segunda solución de cada una de las siguientes ecuaciones diferenciales
1. 1  2 x y' '4 xy '4 y  0; y1  e  2 x
2 .x 2 y ' ' xy ' y  0; y1  cos ln x 

Resuelva las siguientes ecuaciones diferenciales.
1. y ' '4 y  4 cos x  3senx  8
2. y 4   5 y ' '4 y  2 cosh x  6
3. y 4   4 y ' '  5 x 2  e 2 x
4 y" y  4 cos x  senx
5 y" y ' y  xsenx
6 y"4 y  cos 2 x
7 y" '8 y"  6 x 2  9 x  2
8 y" ' y" y ' y  xe x  e  x  7
9 y" '3 y"3 y ' y  e x  x  16
10 2 y" '3 y"3 y '2 y  (e x  e  x ) 2
11 y ( 4 )  2 y" ' y"  e x  1
12 y ( 4 )  4 y"  5 x 2  e 2 x
Darwin D. Peña G.
10. VARIACIÓN DE PARAMETROS
Este método proporciona analíticamente una solución particular de la ecuación
a 2 x  y' ' a1 x  p( x) y' a 0 ( x) y  g ( x) 27

Dividiendo por a 2 x  tenemos
y' ' p( x) y'Q( x) y  f ( x) 28

Ahora bien, encontramos la solución para la ecuación homogénea asociada es decir, la solución general
o complementaria de
y ' ' p( x) y'Q( x) y  0 29

Que es y c  c1 y1  c 2 y 2 30 
Este método consiste en proponer una solución particular de la forma y p  u1 y1  u 2 y 2 31
Donde u1 , u 2 pueden ser funciones de x no constantes. Como hemos supuesto que (31) es solución
de (28) tenemos lo siguiente
y´p  u1 ( x) y1  u1 ( x) y1  u 2 ( x) y 2  u 2 ( x) y 2

condicion para u1 , u 2
y1u1  y 2 u 2  0 32 
y´p  u1 ( x) y1  u 2 ( x) y 2
y p  u1 y1  u1 y1  u 2 y 2  u 2 y 2
reemplazamos en (28)
u1 y1  u1 y1  u 2 y 2  u 2 y 2  P ( x)u1 ( x) y1  u 2 ( x) y 2   Q( x)u1 ( x) y1  u 2 ( x) y 2   f ( x)
u1  y1  P ( x) y1  Q( x) y1   u 2  y 2  P ( x) y 2  Q( x) y 2   u1 y1  u 2 y 2  f ( x)
u1 y1  u 2 y 2  f ( x) (33)
formando un sistema de ecuaciones con (32) y (33)
 y1u1  y 2 u 2  0

u1 y1  u 2 y 2  f ( x)
Resolviendo el sistema por Kramer
y1 y2 0 y2 y1 0
  w( y1 , y 2 ); u1    f ( x) y 2 ; u 2   f ( x) y1
y1 y 2 f ( x) y 2 y1 f ( x)
u1  y 2 f ( x) u  y f ( x)
entonces : u1   ; u 2  2  1
 W ( y1 y 2 )  W ( y1 , y 2 )
Darwin D. Peña G.
 y2 f ( x ) y1 f ( x)
Es decir, u1  ; u2 
W ( y1 y2 ) W ( y1 , y2 )
Finalmente al integrar las expresiones anteriores vemos que
  y 2 f ( x)   y f ( x) 
y p     y1    1  y 2
 W ( y1 2 
y )  W ( y1 , y 2 
)
Ejemplo 10.1: Resolver y' ' y  sec x

Solución:
y ' ' y  0  m 2  1  0  m  i    0;   1
cos x senx 
y g  c1 cos x  c2 senx  W (cos x, senx)     cos 2 x  sen 2 x  1
   senx cos x 
y1 y2
 senx.sec x  senx
u '1    u '1   dx  u1  ln cos x 
1 cos x
cos x.sec x cos x
u '2    u '2   dx  u 2  x
1 cos x
sol : y p  ln(cos x). cos x  xsenx
Y  c1 cos x  c 2 senx  cos x ln(cos x)  xsenx
Para una ecuación de orden 3 se realiza el mismo procedimiento. Veámoslo:

y p  u1 y1  u 2 y 2  u 3 y 3
u '1 y1  u ' 2 y 2  u ' 3 y 3  0
u '1 y '1 u ' 2 y ' 2 u ' 3 y 3  0
u '1 y ' '1 u ' 2 y ' ' 2 u '3 y ' ' 3  f ( x)
 y1 y 2 y 3  0 y2 y3 
 
   y '1 y ' 2 y ' 3   W  y1 , y 2 , y 3 ; u '1  0  y'2 y '3   f ( x) y 2 y ' 3  y 3 y ' 2 
 y ' '1 y ' ' 2 y ' ' 3   f ( x) y ' ' 2 y ' ' 3 
 y1 0 y3   y1 y2 0

u ' 2   y '1 0  
y '3    f ( x) y1 y '3  y 3 y '1 ; u ' 3   y '1 y'2 0   f ( x) y1 y ' 2  y 2 y '1 
 y ' '1 f ( x) y ' '3   y ' '1 y ' ' 2 f ( x)
entonces
u1 f ( x ) y 2 y ' 3  y 3 y ' 2  u 2  f ( x) y1 y '3  y 3 y '1 
u '1   ; u'2  
 W ( y1 y 2 , y 3 )  W  y1 , y 2 , y 3 
u 3 f ( x) y1 y ' 2  y 2 y '1 
u '3   ; luego se Integra para encontrar la solucion
 W  y1 , y 2 , y 3 
Darwin D. Peña G.
Ejemplo 10.2: Resolver: 2 y' ' '6 y' '  x 2
Solución:
2 y ' ' '6 y ' '  0
2m 3  6m 2  0
2m 2 (m  3)  0  m1  m2  0; m3  3
y1  1; y 2  x; y 3  e 3 x
y g  c1  c 2 x  c3 e 3 x ; y p  u '1 y1  u ' 2 y 2  u '3 y 3
y p  u '1 y1  u ' 2 y 2  u '3 y 3
1 x e3x 
  x e3x 

W 1, x, e  0
3x
  
1 3e 3 x   9e 3 x ;W x, e 3 x   3x 
 3 xe 3 x  e 3 x
0 0 9e 3 x  1 3e 

1 e3x 
 
W 1, e 3 x   3x 
 3e 3 x
0 3e 
 x 1
W (1, x)    1
0 1

x 2 3 xe 3 x  e 3 x 
3x 3 x 2 3x 3 x2 1 x4 1 x3 x4 x3
u '1 
9e 3 x

9

9
  1  9 9
u '    u '1   
3 4 9 3 12 27

u'2 

 x 2 3e 3 x 

 x2
  u'2  
 x2 1 x3  x3
 
9e 3 x 3 3 x 3 6
 
2 2 2
x (1) x x 1 1
u '3  3 x  3 x   u '3   3 x   x 2 e 3 x dx   x 2 3e 3 x  2 x9e 3 x  54e  x
9e 9e 9e 9 9
Re emplazamos y obtenemos la solución particular

y p  u '1 y1  u ' 2 y ' 2 u '3 y '3
Darwin D. Peña G.
EJERCICIOS 10
1. y   4 y   4 y  ( x  1)e 4 x
ex
2. y   2 y   y 
1 x2
3. 3 y   6 y   6 y  e x sec x
4. y   4 y   4 y  (12 x 2  6 x)e 2 x
5. y   2 y   8 y  2e 2 x  e x
6. y   y   tan x
7. ( D 2  1) y  csc x
8. ( D 2  4D  3) y  (1  e  x ) 1
9 y"3 y '2 y  1 /(1  e x )
10. y"3 y '2 y  e 3 x /(1  e x )
11 y"3 y '2 y  sene x
12 y"2 y ' y  e x arctan x
13 y"2 y ' y  e x /(1  x 2 )
14 y"2 y '2 y  e x sec x
15 y"2 y ' y  e  x Inx
16 y"10 y '25 y  e 10x / x 2
17 y"2 y '10 y  3e x tan 3 x
18 4 y"4 y ' y  e x / 2 1  x 2
19 y" '2 y" y '2 y  e 3 x
20 2 y" '6 y"  x 2
Darwin D. Peña G.
11. APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE SUGUNDO ORDEN
11.1 MOVIMIENTO ARMÓNICO SIMPLE
Una masa m , sujeta a un resorte flexible suspendido de un soporte rígido hace que se alargue el resorte.
La ley de Hooke dice que el resorte mismo ejerce una fuerza F de restitución opuesta a la dirección de
alargamiento. Si s es el alargamiento total, entonces F  ks , donde k es una constante de
proporcionalidad que depende del material con que se hace el resorte.
Usando la segunda ley de Newton, resulta que el desplazamiento xt  del resorte sigue la ecuación
diferencial
d 2x k
 x  0 34
dt 2 m
k k
Como es positiva, las raíces de la ecuación característica son  i y por tanto la solución
m m
general es
c1 c2
xt   c1 cos t 35 . Si hacemos sen 
k k
t  c 2 sen y cos   , vemos que
m m c1  c 2 c1  c 2
 k 
xt  se puede expresar de la forma alternativa xt   c1  c 2 sen t    . Haciendo
 m 
A  c1  c 2 la solución nos queda xt   Asent   

k
 y
m
Darwin D. Peña G.
Un cuerpo que pesa 2 lb estira un resorte 6 plg. Dicho cuerpo se suelta en t  0 desde un punto que está
a 8 plg bajo la posición de equilibrio, con una velocidad dirigida hacia arriba de 4 pie / seg . Determine
3
la función xt  que describe el movimiento libre resultante.
Solución: Puesto que estamos usando el sistema de unidades inglesas gravitacionales, las magnitudes
dadas en pulgadas deben expresarse en pies: 6 p lg  6

 1 pie, 8 p lg  8  2 pie. Además,
12 2 12 3
debemos convertir las unidades de peso a unidades de masa. Tenemos
2 1
m  slug.
32 16
Además, por la ley de Hooke se tiene
1
2  k   lo que implica que k  4 lb .
2 pie
Por consiguiente, la ecuación nos queda
1 d 2x d 2x
  4 x o bien  64 x  0
16 dt 2 dt 2
Con las condiciones iniciales x0  , x 1 0  

2 4
3 3
xt   cos 8t  sen8t
2 1
Resolviendo obtenemos y aplicando la forma alternativa resulta
3 6
xt   sen8t  1.816

17
6
Un cuerpo que pesa 20lb estira un resorte 6 plg. El peso, en reposo, se suelta desde un punto que esta 6
plg debajo de la posición de equilibrio.
a) Encuentre la posición del resorte en t   2 , t   8 , t   6 , t   4 , t  9 32 segundos.
¿Cuál es la velocidad del cuerpo cuando t  3 16 segundos? ¿En que dirección va el cuerpo en dicho
instante?
Darwin D. Peña G.
20
w  20lb  m 
32
1
6in  ft
2
20
F  kx  k   40
1
2
2
d x 20 d 2 x 20 d 2 x
m   kx   40 x   40 x  0
dt 2 32 dt 2 32 dt 2
d 2x
2
 64 x  0  m 2  64  0  m   64  m  8i
dt
x(t )  c1 cos 8t  c2 sen8t
1 1
x ( 0)   c1 
2 2
1
x(t )  cos 8t  c2 sen8t  x´ (t )  4 sen8t  8c2 cos 8t
2
x (0)  0  8c2  0  c2  0
´
1  1 8 1
x(t )  cos 8tt   x  cos x
2 12 2 12 4
 1 8 1
t x cos x
12 2 12 4
 1 8 1
t  x  cos x
8 2 8 2
 1 8 1
t  x  cos x
6 2 6 9
 1 8 1
t  x  cos x
9 2 9 2
9 1 72 2 1
t  x  cos x 
32 2 32 4 2 2
b.
v?
3
t  seg
16
x ´ (t )  4 sen8t
 3 
x ´  4 sen seg   x ´  4 pie / seg
 2 
Darwin D. Peña G.
Una fuerza de 400N estira un resorte 2 m. Una masa de 50 kg se sujeta al extremo del resorte y le suelta
desde la posición de equilibrio con una velocidad dirigida hacia arriba de 10 m/s. Halle la ecuación del
movimiento.
F  40 N
400
w  41kg
9.8
41
F  kx  k 
2
2
d t d 2t 41
m 2  kx  0  50 2  x0
dt dt 2
d 2t 41 41 41 41 41
2
 x 0    x (t )  c1 cos t  c2 sen t0
dt 100 100 10 10 10
x(0)  0
c1  0
41
x(t )  c2 sen t
10
41 41
x´ (t )  c2 cos t  t  0; x´  10
10 10
41 100 100 41
c2 cos 0t  10  c2   x (t )  sen t
10 41 41 10

Otro resorte cuya constante es de 20N/m, se suspende del mismo soporte rígido , paralelamente al
mismo sistema del ejercicio anterior, una masa de 20kg se sujeta al segundo resorte y ambas se sueltan
desde la posición de equilibrio con una velocidad dirigida hacia arriaba de 10m/s.
Cuál de las dos masas tiene mayor amplitud?
Darwin D. Peña G.
Para ejemplo1 Para ejemplo 2
Para el segundo ejemplo la amplitud es dos veces mayor que para el primero pero la frecuencia del
primero es 2 veces mayor que la de la segunda

Una fuerza de 2lb estira un resorte 1 ft, un cuerpo que pesa 3.2lb se sujeta al resorte y luego el sistema
se sumerge en un medio que comunica una fuerza de amortiguación numéricamente igual a 0.4 veces la
velocidad instantánea.
a) Obtenga la ecuación del movimiento si el cuerpo se suelta a partir del reposo desde un punto que está a
un pie sobre la posición de equilibrio.
Primero calculamos k
Darwin D. Peña G.
11.2 MOVIMIENTO AMORTIGUADO
Resulta que el desplazamiento de la masa sujeta a las condiciones de la primera parte se amortigua por
fuerzas internas como fuerza externas. Por lo menos habrá la fuerza de amortiguación debida al medio
dx
que la rodea. Se ha determinado que las fuerzas de amortiguación están por un múltiplo de , o sea
dt
que esta fuerza es proporcional a la velocidad instantánea de la masa. Llamamos  a la constante de
proporcionalidad.
La segunda ley de Newton ahora dice que

d 2x 
m 2  kx  
dx
36 . Dividiendo por m , tenemos que  2 y haciendo  2 
k
dt dt m m
encontramos que la raices de la ecuación auxiliar son m1, 2    2   2 . De 2   2 ,
distinguimos tres casos posibles.
CASO I: 2   2  0 solución: xt   c1e 1  c 2 e 2

mt m t
Puesto que el coeficiente de amortiguación es grande comparado con la constantes k del resorte,
decimos que el sistema está sobreamortiguado.
CASO II: 2   2  0 . En este caso, xt   c1e m1t  c2 te m1t . Decimos que el sistema está
críticamente amortiguado ya que una pequeña disminución de la fuerza de amortiguación provocaría un
movimiento oscilatorio.
CASO III: 2   2  0 . En este caso, tenemos que:
 
xt   e t c1 cos 2   2 t  c 2 sen 2   2 t . Decimos que el movimiento está subamortiguado
porque  es pequeña comparada con k . El movimiento es oscilatorio, aunque con amplitud cada vez
menor, ya que xt    cuando t  

Darwin D. Peña G.
Un cuerpo que pesa 16 lb sujeta a un resorte de 5 pie de largo. En estado de equilibrio, el resorte mide
8.2 pie. Si el peso se empuja hacia arriba y se suelta, a partir del reposo, desde un punto que esta 2 pie
sobre la posición de equilibrio, determinar los desplazamientos xt  sabiendo además que el medio
ofrece una resistencia numéricamente igual a la velocidad instantánea.
Solución: El alargamiento experimentado por el resorte después que se le sujeta el peso es
8.5  2  3.2 pie , luego por la ley de Hooke se obtiene que 16  k 3.2, k  5 lb . Además
pie
16 1
m  slug y la ecuación diferencial nos queda
32 2
1 d 2x d 2x
 2  10 x  0 sujeta a x0  2, x 1 0  0 . Resolviendo esta
dx dx
2
 5 x  o bien 2
2 dt dt dt dt
 
ecuación bajo estas condiciones nos resulta xt   e t   2 cos 3t  sen3t  . La forma alternativa para
2
 3 
 
esta solución es xt   Ae t sen  2  2 t   . De esta forma tenemos que la solución en la forma
alternativa es: xt   10e t sen3t  4.391

2
3
Una fuerza de 2lb estira un resorte 1 pie. Un cuerpo que pesa 3.2 lb se sujeta al resorte y luego el sistema
se sumerge en un medio que comunica una fuerza de amortiguación numéricamente igual a 0.4 veces la
velocidad instantánea. Obtenga la ecuación del movimiento si el cuerpo se suelta a partir del reposo,
desde un punto que esta 1 pie sobre la posición de equilibrio. Exprese la ecuación del movimiento.
2
F  2lb; x1  1 pie; w  3.2lb;   0.4 
5
F 2 3.2
F  kx  k  k  ; x2   x2  1.6
x 1 2
d 2x dx 1 d 2 x 2 dx
m    kx  0    2x  0
dt 2 dt 10 dt 2 5 dt
d 2x dx  4  16  4(1)(20)
2
4  20 x  0  m 2  4m  20  0  m 
dt dt 2
 4   64  4  8i
m m  m1  2  4i  m2  2  4i
2 2
  2;  4
Darwin D. Peña G.
x(t )  e  2t c1 cos 4t  c 2 sen4t   x(0)  e  2 ( 0 ) c1 cos 4(0)  c 2 sen4(0)  1

c1  1
x(t )  e 2t  cos 4t  c 2 sen4t 
x ´ (t )  e  2t 4 sen4t  4c 2 cos 4t   2 cos 4t  c 2 sen4t e  2t
x ´ (t )  e  2 ( 0 ) 4 sen4(0)  4c 2 cos 4(0)  2 cos 4(0)  c 2 sen 4(0)e  2 ( 0 )
2 1
x ´ (0)  4c 2  2  0  c 2   c2 
4 2
 1 
x(t )  e  2t  cos 4t  sen 4t 
 2 
Después que un cuerpo que pesa 10 lb se sujeta a un resorte de 5 pie de largo, el resorte mide 7 pie. Se
quita el cuerpo de 10 lb y se le reemplaza por uno de 8 lb; el sistema completo se coloca en un medio
que ofrece una resistencia numéricamente igual a la velocidad instantanea.
Obtenga la ecuación del movimiento si el peso se suelta desde un punto que esta ½ pie debajo de la
posición de equilibrio con una velocidad dirigida hacia debajo de 1 pie/s.
Exprese la ecuación el movimiento.
w 8 1
w1  10lb; w2  8lb; m2   m
g 32 4
F 10
F  kx  k  
x1 7
87 28
x2  
10 5
2
d x dx 1 d 2 x dx 10
m    kx  0    x0
dt 2 dt 4 dt 2 dt 7
d 2x dx 40 40
2
4  x  0  m 2  4m  0
dt dt 7 7
16
 4  16 
m 7  m  2  1.3i
1
2
x (t )  e  2t (c1 cos 1.3t  c 2 sen1.3t )
1 1
x (0)   c1 
2 2
Darwin D. Peña G.
 1.3  1 
x´ (0)  e 2t  sen1.3t  1.3c2 cos 1.3t   ( 2)e 2t  cos 1.3t  c2 sen1.3t 
 2  2 
2 20
x´ (0)  1.3c2  1  1  c2   c2 
1.3 13
1 20 
x(t )  e 2t  cos 1.3t  sen1.3t 
2 13 
11.3 MOVIMIENTO FORZADO
Si aplicamos una fuerza externa f t  al sistema, entonces la ecuación diferencial se escribe

d 2x f t 
  2 x  F t  37  donde F t  
dx
2
 2
dt dt m
Para resolver esta ecuación no homogénea usamos los métodos de coeficientes indeterminados o
variación de parámetros. Y si el movimiento carece de amortiguación obtenemos.
d 2x f t 
2
  2 x  F t  38 donde F t  
dt m
En el caso de las oscilaciones subamortiguadas, demuestre que la solución general de la ecuación
diferencial: Muestre que la solución es

Darwin D. Peña G.
Solución:
d 2x dx
 2   2 x  F0 sent
dt 2 dt
d 2x dx
 2   2 x  m 2  2m   2  0
dt 2 dt
 2  42  4 2
m  m    2   2
2
   ;   
2   2  xc  e 2  c1 cos 2   2 t  c2 sen 2   2 t 
( D 2   2 ) F0 sent  0
D 2   2  0  m 2   2  m  i
x p (t )  c3 cos t  c4 sent
x´p (t )  c3 sent  c4 cos t
x´´p (t )   2 c3 cos t   2 c4 sent
  2 c3 cos t   2 c4 sent  2 c3 sent  2 c4 cos t
  2 c3 cos t   2 c4 sent  F0 sent
(  2 c3  2 c4   2 c3 ) cos t  0
(  2 c4  2 c3   2 c4 ) sent  F0 sent
c3 ( 2   2 )  2 c4  0
c4 ( 2   2 )  2 c3  F0
2 c4  2 
 c3    c4 ( 2   2 )  2   2 
  F0
(   )
2 2
 ( 2
  ) 
Darwin D. Peña G.
 ( 2   2 )  42 2 
c4 
( 2   2 )   F0
 
F0 (   )
2 2
c4 
(   2 )  42 2
2
 2
 F0 ( 2   2 ) 
c3   2 2 2 
( 2   2 )
 (   )  4  
2
2 F0
c3  
(   2 )  42 2
2
2 F0 F0 ( 2   2 )
x(t )   cos t  sent
( 2   2 )  42 2 ( 2   2 )  42 2
x(t )  Ae t sen(t   )
F0
x p (t )  sen(t   )
(   2 )  42 2
2
A c32  c42
42 2 F0 ( 2   2 ) 2 F0
2 2
A 
  
42 2  ( 2   2 ) 2 42 2  ( 2   2 ) 2 
F 4   (   ) 
2 2 2 2 2 2
A 0
4   (   ) 
2 2 2 2 2 2
F0
A
4 2 2
 ( 2   2 ) 2 
x(t )  Ae t sen  2   2t     4 
2 2
F0

 ( 2   2 ) 2
sen(t   )
OTROS PROBLEMAS RESUELTOS
1.) Variación de los parámetros. La posición y la aceleración, en función del tiempo, de una masa
puntual moviéndose unidimensionalmente, vienen relacionadas por la ecuación diferencial
x(t )  a(t )  tg t  0 (unidades mks)
Determinar la ecuación del movimiento (posición en función del tiempo) de la partícula si la misma
parte del origen con una velocidad de 3 m/s.
Solución
Expresando la aceleración como la derivada segunda de la posición y reordenando la ecuación tenemos:
x(t )  x(t )  tg t
Darwin D. Peña G.
Hallemos primero la solución de la ecuación homogénea asociada. Es ésta:
x(t )  x(t )  0
La ecuación característica es:
m 2  1  0    m  i  xc  C1 cos t  C2 sent
Ahora debemos hallar una solución particular del problema no homogéneo. Vemos que, por el tipo de
función excitación, deberemos usar En vista de la base de soluciones del problema homogéneo halladas,
será:
x*  v1 cos t  v2 sen t
Hallemos ahora v1 y v2:
 multiplicando arriba por sent
 1v  cos t  v 
2 sen t  0
y abajo por cost
  v1 cos t sen t  v 2 sen 2 t  0
 v  sen t  v  cos t  tg t  sen t  
 v1 sen t cos t  v 2 cos t  sen t
2
 1 2
cos t

Si ahora sumamos las dos ecuaciones de este último sistema tendremos.
Re emplazando en la
primera ecuación

sen 2 t
v 2  sen t  v 2   cos t  v1 cos t sen t  sen 3 t  0  v1   
cos t
tablas
1  cos t2
1 
  cos t   v1  sen t  log sec t  tg t
cos t cos t
Con las funciones v1 y v2 así obtenidas podemos escribir:
x*  v1 cos t  v2 sen t  sen t  log sec t  tg t cos t  cos t sen t  log sec t  tg t cos t
con lo cual la solución general del problema no homogéneo es:
x  xc  x p  C1 cos t  C2 sent  log sec t  tgt cos t
Las condiciones iniciales indican que cuando t es 0 , x = 0 y v = 3. De esta manera, podemos escribir:
Darwin D. Peña G.
x(0)  C1 cos 0 C 2 sen 0  log sec 0  tg 0 cos 0  C1  0

v(0)  x (0)  C1 sen 0 C 2 cos 0  1  log sec 0  tg 0 sen 0 C 2 1  3 C 2  2
Con estos valores de la constante tenemos, finalmente:
x  2 sen t  log sec t  tg t cos t
Que es la solución al problema de valores iniciales planteado.
2.) Coeficientes indeterminados. Resolver el problema de valores iniciales:
x (t )  2 x (t )  2 x(t )  e t cos t

x(0)  x (0)  0
Solución:
Primero resolvemos la ecuación homogénea asociada. La ecuación característica será:
 2  2 2  4 1  2
m  2m  2  0  m 
2
 1  i  xc  C1e t cos t  C 2 e t sent
2
Ahora nos toca hallar la solución particular del problema no homogéneo. Vemos que la función
excitación es un producto de una exponencial por un coseno, por lo cual podemos intentar hallar la
solución particular por el método de los coeficientes indeterminados. La función excitación es e-tcost, de
modo que normalmente propondríamos como función de prueba una combinación lineal de senos y
cosenos, Ae-tcost + Be-tcost. Sin embargo, en este caso particular ésta sería una solución del problema
homogéneo, por lo cual debemos multiplicarla por la variable t. Obtenemos así:
x(t )  Ate  t cos t  Bte t sen t

x (t )  Ae t cos t  Ate t cos t  Ate t sen t  Be t sen t  Bte t sen t  Bte t cos t
x (t )   Ae t cos t  Ae t sen t  Ae t cos t  Ate t cos t  Ate t sen t  Ae t sen t  Ate t sen t 
 Ate t cos t  Be t sen t  Be t cos t  Be t sen t  Bte t sen t  Bte t cos t  Be t cos t 
 Bte t cos t  Bte t sen t 
 2 Ae t cos t  2 Ae t sen t  2 Ate t sen t  2 Be t sen t  2 Be t cos t  2 Bte t cos t
Con estas expresiones podemos reemplazar en la EDO del enunciado y obtener:
Darwin D. Peña G.
x   2 x   2 x  2 Ae t cos t  2 Ae  t sen t  2 Ate  t sen t  2 Be t sen t  2 Be t cos t  2 Bte t cos t 

 2 Ae t cos t  2 Ate t cos t  2 Ate t sen t  2 Be t sen t  2 Bte t sen t  2 Bte t cos t  2 Ate t cos t 
igualando con la
función excitación

 2 Bte t sen t  2 Ae t sen t  2 Be t cos t  e t cos t  A  0, B   12
De modo que una solución particular al problema no homogéneo será:
x p (t )   12 e t tsent  x(t )  xc (t )  x p (t ) C 1e t cos t  C2 e t sent  12 e t tsent
Queda como ejercicio para el lector deducir a partir de las condiciones iniciales que las constantes deben
ser nulas y obtener que x(t )   12 e t t sen t .
Como se ve, este método puede requerir derivaciones algo complicadas.
3.) Sistema masa-resorte. Un sistema masa-resorte está caracterizado por los siguientes valores: masa m
= 0,2; constante del resorte k = 80; constante de amortiguamiento h = 2. Si se aplica una excitación F(t)
= 2cos30t, obtener el estado estacionario de la respuesta, expresándolo en la forma x(t) = Acos(t - ).
Solución:
El cuerpo se mueve bajo la acción de 3 fuerzas: F, la del resorte y la del amortiguador. La del resorte es
proporcional y de signo contrario al desplazamiento, mientras que la de amortiguación es proporcional y
de signo contrario a la velocidad. En otras palabras:
F
i
i  F  Fresorte  Famortiguador  F  (kx)  (hv)  F  kx  hv
Por otro lado, y según la segunda ley de Newton, esta sumatoria debe ser igual a la masa por la
aceleración. Si además tenemos en cuenta que la velocidad y la aceleración son la derivada primera y
segunda de la posición, tendremos:
datos

F  kx  hv  ma  ma  hv  kx  F  mx  hx  kx  F  0,2 x  2 x  80 x  2 cos 30t
Esta última es una EDO de 2º orden a coeficientes constantes, no homogénea. Para resolverla, primero
encaramos el problema homogéneo. Planteando la ecuación característica y resolviendo tenemos:
 2  2 2  4  0,2  80
m  5  375i  xc  C1e 5t cos 375t  C 2 e 5t sen 375t
0,4
Ahora buscaremos una solución particular del problema no homogéneo. Vemos que la función
excitación es cosenoidal, y por ende usaremos como función de prueba una combinación lineal de senos
y cosenos.
Darwin D. Peña G.
x p  A cos 30t  Bsen30t

x p  30 Asen30t  30 B cos 30t
x p  900 A cos 30t  900 Bsen30t
0,2 x p  2 x p  80 x p  180 A cos 30t  180 Bsen30t  60 Asen30t  60 B cos 30t  80 A cos 30t 
 80 Bsen30t  (100 A  60 B) cos 30t  (60 A  100 B) sen30t  F (t )  2 cos 30t 
 100 A  60 B  2 A   340 5
   x p   340
5
cos 30t  340
3
sen30t
  60 A  100 B  0 B  3
340
De modo que la solución general del problema no homogéneo será:
x  xc  x p  C1e 5t cos 375t  C2 e 5t sen 375t  340

5
cos 30t  340
3
sen30t
En el estado estacionario, para tiempos muy grandes, los dos primeros términos tienden a 0 y podemos
escribir:
xest   340
5
cos 30t  340
3
sen 30t
La manera normal de expresar un movimiento oscilatorio de este tipo es de la forma x(t) = Acos(t - ).
Si aplicamos la identidad trigonométrica del coseno de una suma, tendremos:
en nuestro caso

A cos(t   )  A cos t cos   A sen t sen    340
5
cos 30t  340
3
sen 30t 
dividiendo m.a.m.
  30 la 3a. con la 2a.
  340
5
 A cos    340
5
 tg       tg    2,6  A 
3
5
1 3
5  0,01716
cos 2,6
A sen   340
3
De modo que finalmente podemos escribir:
xest = 0,01716cos(30t - 2,6)
Nótese que A, la amplitud del movimiento, es siempre un valor positivo. Por eso elegimos un  tal que
su coseno fuera negativo, de modo que al despejar A nos diera un número positivo.
Darwin D. Peña G.
4.) Circuito eléctrico. Hallar la corriente que circula en un circuito serie RLC, sabiendo que R = 120 ,
L = 10 H, C = 10-3 F si la fuerza electromotriz viene dada por E(t) = 17sen2t V, y si la intensidad cumple
las condiciones i(0) = 0, i (0)   201 . Hallar asimismo la corriente en estado estacionario.
Solución:
La intensidad es la derivada de la carga eléctrica Q. Genéricamente se puede escribir para un circuito

RLC:
en nuestro problema
 Q
2
Q 1 
 2Q Q
L R  Q  E (t )  10  120  1000Q  17 sen 2t
t 2
t C t 2
t
Ésta es una típica ecuación de 2º orden a coeficientes constantes. Hallamos primero la solución general
del problema homogéneo y luego le sumamos una solución particular del problema no homogéneo. Para
lo primero planteamos:
 120  120 2  4  10  1000

m  6  8i  Qc  C1e 6t cos 8t  C 2 e 6t sen8t
20
Para la solución particular del problema no homogéneo, podemos usar una función de prueba:
Q p  Asen 2t  B cos 2t
Q p  2 A cos 2t  2 Bsen 2t
Q p  4 Asen 2t  4 B cos 2t
10Q p  120Q p  1000Q p  40 Asen 2t  40 B cos 2t  240 A cos 2t  240 Bsen 2t 
en este problema

 1000 Asen 2t  1000 B cos 2t  (960 A  240 B) sen2t  (240 A  960 B) cos 2t 
960 A  240 B  17 A  601
 17 sen2t     Qp  1
sen2t  240
1
cos 2t
240 A  960 B  0 B   240
1 60
De modo que la solución general del problema no homogéneo es:
Q  Qc  Q p  C1e 6t cos 8t  C2 e 6t sen8t  601 sen2t  240

1
cos 2t
Introduciendo las condiciones iniciales es:
Q  C1 (6e 6t cos 8t  8e 6t sen 8t )  C 2 (6e 6t sen 8t  8e 6t cos 8t )  301 cos 2t  120
1
sen 2t 
 Q(0)  6C1  8C 2  301  i(0)  0  6C1  8C 2  1
30
Similarmente,
Darwin D. Peña G.
Q   C1 (36e 6t cos 8t  48e 6t sen 8t  48e 6t sen 8t  64e 6t cos 8t )  C 2 (36e 6t sen 8t 
48e 6t cos 8t  48e 6t cos 8t  64e 6t sen 8t )  151 sen 2t  601 cos 2t 
 Q (0)  18C1  96C 2  601  i (0)   201  18C1  96C 2  151
Reuniendo los dos últimos resultados, podemos escribir:
 6C1  8C 2  301 C 2   3600

1
   
18C1  96C 2  15 C1  1350
1 7
De modo que la corriente que circula por el circuito vendrá dada por:
Q  Qc  Q p  1350
7
e 6t cos 8t  3600
1
e 6t sen8t  601 sen2t  240
1
cos 2t
La corriente en estado estacionario vendrá dada por los dos últimos términos de esta expresión.
Darwin D. Peña G.
EJERCICIOS 11
11.1 El periodo de las oscilaciones libres no amortiguadas de una masa sujeta a un resorte es
 segundos. Si la constantes del resorte es 16 lb/pie, ¿cuál es el valor numérico del peso?.
4
11.2 Una objeto de 1kg se sujeta a un resorte cuya constante es de 16N/m y el sistema completo se
sumerge en un líquido que comunica una fuerza de amortiguación numéricamente igual a 10 veces la
velocidad instantánea. Determine las ecuaciones del movimiento si
a) El peso se suelta, a partir del reposo, desde un punto que esta a 1m debajo de la posición de equilibrio.
b) El peso se suelta desde un punto que esta a 1m debajo de la posición de equilibrio con un velocidad
dirigida hacia arriba de 12m/s.
11.3 En las partes (a) y (b) del problema anterior, determine si el cuerpo pasa por la posición de
equilibrio. En cada caso, encuentre el instante en el cual el cuerpo alcanza su desplazamiento externo
desde la posición de equilibrio. ¿Cuál es la posición del citado cuerpo en dicho instante.?
11.4 Un cuerpo que pesa 16lb estira un resorte 8 pie . Inicialmente, el peso parte del reposo desde un
3
punto que está 2pie debajo de la posición de equilibrio y el movimiento posterior se realiza en un medio
que opone una fuerza de amortiguación numéricamente igual a 1 de la velocidad instantanea.
2
Encuentre la ecuación del movimiento si el peso es impulsado por una fuerza exterior igual a
f t   10 cos 3t
11.5 Un cuerpo que pasa 24 lb estira un resorte 4 pie. El movimiento subsiguiente se realiza en un
medio que ofrece una resistencia numéricamente igual a veces ( >0) la velocidad instantánea. Si el
cuerpo parte de la posición de equilibrio con una velocidad dirigida hacia arriba de 2 pie/s y si > ,
demuestre que la ecuación del movimiento es
11.6 Demuestre que la solución del problema de valor inicial

2
F
  2 x  F0 cos t , x0  0, x 1 0  0 es xt   2 0 2 cos t  cos t 
d x
dt 2
 
F0
Calcule el lim
     2
2
cos t  cos t 
Darwin D. Peña G.
12. ECUACIONES LINEALES CON COEFICIENTES VARIABLES
La solución de las ecuaciones con coeficientes variables o ecuación de Cauchy-Euler o también llamada
equidimensional es análoga a las ecuaciones con coeficientes constantes, veámoslo.
dny d n 1 y
an x n n
 a n 1 x n 1 n 1
   a1 x
dy
 a0 y  g ( x) 39
dx dx dx
para n  2
d2y
a2 x 2 2
 a1 x
dy
 a0 y  0 40
dx dx
La solución es de la forma
y  xm
y   mx m 1 ;
y   m( m  1) x m  2
Re emplazando en 40 obtenemos

ax 2 y   bxy   cy  0
     
ax 2 m(m  1) x m  2  bx mx m 1  c x m  0
am  m(b  a )  c  0
2
 (b  a )  (b  a ) 2  4ac
m
2a
m1 , m2  R ; m1  m2 m1 , m2  R ; m1  m2
 
Caso I: (b  a) 2  4ac  0 Caso II: (b  a) 2  4ac  0
 
 y  c1 x  c 2 x  y  c1 x  c 2 x ln x
m1 x m2 m1 m1
m1 , m2  C

(b  a)  4ac  0
2
Caso III: 
m     i
 y  x  c cos ln x   c senln x 
 1 2
Darwin D. Peña G.
Ejemplo 12.1: Resolver xy ' ' y'  x
Solución:
xy ' '  y '  0
y  xm
y '  mx m1
y ' '  m(m  1) x m2
 
x xm  1x m2  x m2  mx m1  0  x m1 mm  1  m  0
m 2  m  m  0  m 2  2m  0  m(m  2)  0
m1  0
m2  2
Lo cual obtenemos la solución general y las funciones y1 , y 2

y g  c1  c 2 x 2
y1  1
y2  x 2
y p  u1 y1  u 2 y 2
1 x 2 
W (1, x 2 )    2x
0 2 x 
 y 2 f ( x )  x 2 .x x2 1
u '1      x3
W 2x 2 6
y f ( x) x 1
u'2  1   x
W 2x 2
1 1 1
y p   x3  x3  x3
6 2 3
1
Y  c1  c 2 x 2  x 3
3
EJERCICIOS 12
1. x 2 y   2 xy   4 y  0
2. 4 x y   8xy   y  0
3. x 3 y   5x 2 y   7 xy   8 y  0
4. x 2 y   2 xy   2 y  x 4 e x
5. xy   4 y   x 4
6. 4 x 2 y   10 xy   2 y  2 x 2  2 x
Darwin D. Peña G.
Capítulo III
Darwin D. Peña G.
13. SOLUCIONES EN SERIES DE POTENCIAS
13.1 Solución en Series de Potencias alrededor de puntos ordinarios, en particular alrededor de cero.
Para este caso vamos suponer una solución en la forma de una serie de potencias de x ,

y   Cn x n 40
n 0
dy
Ejemplo 13.1: Resolver  y0
dx

Solución: y   C n x n
n 0
  
y'   nCn x n1
n 0
Al sustituir en la ecuación se tiene  nCn x n1   Cn x n  0
n 0 n 0
 
0.C0 x 1   nCn x n 1   C n x n  0
n 1 n 0
k  n 1 k n
n  k 1
 
0.C0 x 1   ( k  1)Ck 1 x k   C x x k
k 0 k 0
 k  1C  C k x k  0  k  1C k 1  C k  0  C k 1 

Ck
k 1
k 0 k 1
k  0  c1  c 0
c1 c 0 c 0 c c c c c c
k  1  c2    2; k  2  c 3  2  0  0 ; k  3  c 4  3  0  0
2 2 2! 3 2.3 3! 4 2.3.4 4!
c c0 c
k  4  c5  4   0
5 2.3.4.5 5!

c0 x 2 c0 x 3 c0 x 4 c0 x 5
y   C n x  c 0  c1 x  c 2 x  c 3 x  c 4 x  c 5 x    c 0 c 0 x 
n 2 3 4 5
   
n 0 2! 3! 4! 5!
 x2 x3 x4 x5  
xn
 c 0 1  x        c 0 
 2! 3! 4! 5!  n  0 n!

xn
y  c0   ex
n o n!
Darwin D. Peña G.
Ejemplo 13.2: Resolver ( x  1) y' ' y'  0

Solución:
  
y   C n x n ; y '   nC n x n 1 ; y ' '   n(n  1)C n x n  2
n 0 n 0 n 0
 
( x  1) n(n  1)C n x n  2   nC n x n 1  0
n 0 n 0
  
 n(n  1)C n x n1   n(n  1)C n x n2   nC n x n1  0

n 0 n 0 n 0
 
 n(n  1)  nC
n 0
n x n 1   n(n  1)C n x n  2  0
n 0
 
nC n x n 1   n(n  1)C n x n  2  0

n 0
2
n 0
  
 
k  n 1 k n2
n  k 1 nk  2
 
0.c0 x 1  0.(1)c0 x  2  1(0)c1 x 1   n 2 C n x n 1   n(n  1)C n x n  2  0
     n 1 n2
c0  0 c0  0 c1  0
 
 (k  1) 2 c k 1 x k   (k  2)(k  1)c k  2 x k  0
k 0 k 0
 (k  1) 

2
C k 1  ( k  2)(k  1) x k  0
k 0
 ( k  1) 2 c k 1  ( k  2)(k  1)c k  2  0
( k  1) 2 c k 1 ( k  1)c k 1 ( k  1)c k 1
Ck 2   
( k  2)(k  1) k2 k2
c
k  0  c2  1
2
2c 2 c
k  1  c3   1
3 3
3c3 c1
k  2  c4  
4 4
4c 4 c
k  3  c5   1
5 5
5c5 c
k  4  c6   1
6 6
6c 6 c
k  5  c7   1
7 7
Luego
Darwin D. Peña G.

y   C n x n  c 0  c1 x  c 2 x 2  c 3 x 3  c 4 x 4  c 5 x 5  c 6 x 6  c 7 x 7  
n 0
c1 2 c1 3 c1 4 c1 5 c1 6 c1 7
 c 0  c1 x  x  x  x  x  x  x 
2 3 4 5 6 7
 x 2
x 3
x 4
x 5
x 6
x 7

 c 0  c1  x        
 2 3 4 5 7 7 

xn
y  c 0  c1 
n 1 n
Ejemplo 13.3: Resolver: y” - xy’ - y = 0 , y(0) = 1 , y’(0) = 0
Solución:

y   cn x n
n 0
  
y    nc n x n 1  xy    nc n x n   nc n x n
n 1 n 1 n 0
 
y    n(n  1)c n x n  2   (n  2)(n  1)c n  2 x n
n2 n 0
Reemplazando en la Ecuación Diferencial y tenemos:
  
y   xy   y   (n  2)(n  1)c n  2 x n   nc n x n   c n x n 
n 0 n 0 n 0

  (n  2)(n  1)c n  2  (n  1)c n x n  0  (n  2)(n  1)c n  2  (n  1)c n  c n  2 
cn
n 0 n2
Por otro lado:
y (0)  1  c0  1
y (0)  0  c1  0
c0 1 c c 1 1 1
c2    c4  2  0   c6  ···  ··· c 2 n  
2 2 4 24 24 246 2  4  6     2n
1 1 1
 c2n   n  n
(2  1)  (2  2)  (2  3)     2n 2 (1  2  3     n) 2 n!
Los coeficientes impares, por su parte, serán nulos debido a que c1 es nulo y la fórmula de recurrencia
remite a ese primer coeficiente.
Darwin D. Peña G.
De esa forma, la función que buscamos queda:

x 2n
y n
n  0 2 n!
Ejemplo 13.4 Hallar los cuatro primeros términos no nulos de la solución en series de potencias de la
Ecuación Diferenial:
y’ = x2 + y2 , y(0) = 1
Solución:
Si planteamos (en cada caso los puntos suspensivos indican términos de orden superior):
y  a 0  a1 x  a 2 x 2  a3 x 3  ···; y(0)  1  a 0  1
y   a1  2a 2 x  3a3 x 2  ···(*)
Por otro lado, de la Ecuación Diferencia l :

y   x 2  y 2  x 2  a 0  a1 x  a 2 x 2  a3 x 3  ···  2
 x 2  a 02  2a 0 a1 x  a1 x 2  2a 0 a 2 x 2  ··· a 02  2a 0 a1 x  1  a1  2a 0 a 2 x 2 
 1  2a1 x  1  a1  2a 2 x 2 (**)
Igualando (*) y (**) :
a1  1
2a 2  2a1  2  a 2  1
3a3  1  a1  2a 2   4  a3  4
3
Con lo cual :
y  a 0  a1 x  a 2 x 2  a3 x 3  ··· 1  x  x 2  43 x 3  ···
Darwin D. Peña G.
Ejemplo 13.5 Solucionar en serie de potencias

Solución:
K=0
K=1
K=2
K=3
K=4
K=5
Darwin D. Peña G.
EJERCICIOS 13
Resolver:
1. y   xy  0
2. y   2 xy   y  0
3. y   x 2 y   xy  0
4. y   ( x  1) y   y  0
5. ( x 2  1) y   6e  0
6. ( x 2  2) y   3xy   y  0
7. ( x  1) y   y   0
8. 2 x 2 y   xy   ( x 2  1) y  0
9. 3xy   2 y   x 2 y  0
10. y   2 xy   8 y  0 , y(0)  3; y (0)  0

Darwin D. Peña G.
14. LA TRANSFORMADA DE LAPLACE
Definción 14.1: Sea f t  una función definida para t  0 ; entonces la integral

 b
e
 st
f (t )dt  lim  e st f (t )dt
b
41
0 0
Se llama transformada de Laplace de f siempre que el límite exista.
Simbólicamente  f (t )  f (s)
Ejemplo 14.1: Resolver f (t )  1
Solución:
 b
 e (1)dt  lim  e dt
 st  st
b 
0 0
b
  e  st    e  st 1  1
lim    lim   
b  b  s
 s 0  s s
1 
1
s0
s
Ejemplo 14.2: Resolver f(t) = senkt
Solución:
 b
e senkt.dt  lim  e  st .senkt.dt

 st
b 
0 0
u  senkt
du  k . cos kt.dt
dv 1
  e  st
v s
  e  st .senk (b)  k b 
 lim   0    e  st cos kt.dt 
b 
 s  s 0 
  e  sb senk (b)  k b 
 lim   0    e  st cos kt.dt 
b 
 s  s 0 
k b 
 lim   e  st cos kt.dt 
b  s
 0 
Darwin D. Peña G.
u  cos kt
du  k .senkt.dt
dv  e  st
1
v   e  st
s
 k   e  st cos kt b k b    k   e  sb 1  k b  st 
lim     e senkt.dt    lim  
 st
    e senkt.dt 
    b  s  s
b  s s s0 s s 0
0 
b
k k2
 2  2 lim  e  st senkt.dt  lim  e  st senkt.dt
s s b  0 b 
k
k  k 2  s2  b
 senkt  2
2 k k
 2    lim  e  st senkt.dt  2 s 2  2
s  s
2
 b  0 k s k s 2
k  s2
s2
Condiciones suficientes para la transformada de laplace
1.continua parte por parte para t  0

2.orden Exponencia l
14.1 FORMULAS DE TRANSFORMADAS

1 
1
 
 tn 
n!
S s n 1
 
 e at 
1
S a
senkt 
k
s  k2
2
cos kt  2 senhkt  2

s k1
S  k2 S  k2
 f (t )  g (t )
cosh kt  2
s
S  k2
Ejemplo 14.3: Resolver f (t )  3  t 2  sen2t

Solución:

 f (t )   3  t 2  sen2t 
  sen2t
 3(1)   (t ) 2
3 2 2
  3  2
s s s 4
Darwin D. Peña G.

Ejemplo 14.4: Resolver:  cos 2 t 
Solución:
 cos 2t  1  cos 2t   1 
     
 2   2  2
 cos 2t  1
1 1
2 2
1 s 11 s 1
   
2 s  4 2 s 2 s  4 2s
2

2

14.2 INVERSAS
1   s 
 1    1  1  2 2 
 cos kt
S  S k 
 n!   k 
 1  n 1   t n  1  2 2 
 senhkt
S  S k 
 1   s 
 1  e
at
 1  2 2 
 cosh kt
S  a S  k 
 k 
 1  2 2 
 senkt
S  k 
1
Ejemplo 14.5: Resolver f ( s )  
( s  1)( s  2)( s  4)
Solución:
1 A B C
f ( s)    
( s  1)( s  2)( s  4) s  1 s  2 s  4
1 1 1
15  6  10
s 1 s  2 s  4
 1  1   1  1   1  1 
          
15  s  1   6  s  2  10  s  4 
1  1  1  1  1  1 
       
15  s  1 6  s  2  10  s  4 
1 t 1  2t 1  4t
e  e  e
15 6 10
Darwin D. Peña G.
EJERCICIOS 14.1
1. Aplicando la definición de transformada de Laplace, evalúe

1. f (t )  te 4t
2. f (t )  e t sent
3. f (t )  t 2 e  2t
4. f (t )  t cos t
2. halle la transformada inversa

 S 
1.  1  
 ( S  2)( S  4) 
2
 S 
2.  1  2 
 S  6 S  11
 3S  2 
3.  1  3 2 
 S ( S  4) 
 6S  3 
4.  1  4 
 S  5S  4 
2
14.3 Teorema de Translación

 e at f (t )    f (t )  ss a  F (s) s s a  F (s  a)
Ejemplo 14.6: Resolver:  e at f t   

Solución:

 e 5t t 3 

 e e 
b
 st 5t 3
t dt  lim  e t 5 s t 3 dt
b 
0 0
3t 2 e  s 5 t  6te  s 5 t  6e  s 5 t
b b b b
t 3 e  ( s 5)t  6
   
5  s 0 s  52 0
s  5 3
0
s  5 4
0
s  54
Otra forma, por el primer teorema
 
 e at f t    f t s  s  a
e t  t 
5t 3 3
s  s 5
a5
3 6

s4 s  s 5 s  54
Darwin D. Peña G.
Ejemplo 14.7: 
 e 2t cos 4t 
Solución:

 e 2t cos 4t 
a  2
f t   cos 4t
s2
cos 4ts s  2 
s

s  16 s  s  2
2
s  22  16
Forma inversa del primer teorema de traslación:
1  f s  a ¨ 1  f s s s a  e at f (t )
Ejemplo 14.8: Resolver  f s  a ¨ e at f (t )
 s 
Solución:  1  2 
 s  4s  5 
 s 
 1  2 
 
 s  4 s  4  1
 s  2  2   s2   2 
 1   
1
 
1

 s  2  1  s  2  1  s  2  1
2 2 2
 s   1 
 1  2   2 1  2 
 s  1 s  s  2  s  1 s  s  2
a  2
e  2t cos t  2e  2t sent
 e  2t cos t  2sent 
Darwin D. Peña G.
EJERCICIOS 14.2
 1 
1.  1  4
 ( s  1) 
 2s  5 
2.  1  2 
 s  6s  34 
 2s  5 
3.  1  2
 ( s  3) 
 2s  1 
4.  1  2 3
 s ( s  1) 
 1 
5.  1  2 
 s  6s  10 
 1 
6.  1  2 
 s  2s  5 
7. e 2t (t  1) 2 
 
8.  e t cos 2 3t
9. t (e  e ) 
t 2t 2
10. t e 
3  2t
11. (1  e  3e t  4t
) cos 5t 
12. t e 
10  7 t
14.4 Transformada De Una Derivada
 

 f (t ) =  e  st f (t )dt  e  st f (t ) 0  s  e  st f (t )dt   f (0)  s  f (t )

0 0
Para el caso de la primera derivada

 f (t )  sF (s)  f (0)
Para el caso de la segunda derivada
 
 f (t )   e  st
f (t )dt  e  st

f (t ) 0  s  e  st f (t )dt   f (0)  s  f (t )

0 0
 ssF ( s )  f (0)  f (0)

 f (t )  s 2 F ( s )  sf (0)  f (0)
Para el caso general
f n (t ) s n F (s)  s n1 f (0)  s n22 f (0)  .....  f n 1
(0)
Donde F (s)   f (t )
Darwin D. Peña G.
dy
Ejemplo 14.9: Resolver  y  1 ; y(0)  0
dx
 dy 
   y  1
 dx 
1
sY ( s )  y (0)  Y ( s ) 
s
Y ( s )s  1  0 
1
s
1
Y (s) 
s ( s  1)
hallamos la transformada inversa
 1 
 1 Y ( s )   1  
 s ( s  1) 
por fracciones parciales tenemos
1 1 1
 
s ( s  1) s s 1
entonces :
 1  1  1 
 1 Y ( s )   1    
 s  1 s
y (t )  e  1
t
EJERCICIOS 14.3
1. y   5 y   4 y  0 y 0  1; y 0  0
2. y   4 y   6e 3t  3e t y 0  1; y 0  1
3. y   3 y  13sen2t y 0  6
4. y   y  2 sen 2t y 0  10; y 0   0
5. y   2 y   y   2 y  sen3t ; y 0  0; y 0  0; y 0   1
Darwin D. Peña G.
15. SISTEMAS DE ECUACIONES LINEALES
Recordemos que una ecuación diferencial es una igualdad matemática en la cual se vincula la o las
variables independientes con la variable dependiente o función y sus sucesivas derivadas. En cursos
básicos de Cálculo generalmente se estudian ecuaciones diferenciales tales que, en una sola ecuación, se
expresa la variación de una variable dependiente, en función de una sola variable independiente. Este
análisis es adecuado para algunos modelos dinámicos, pero hay otros que necesitan analizar el
comportamiento de dos o más variables dependientes. Para ello es necesario plantear varias ecuaciones
simultáneas, que forman un sistema de ecuaciones diferenciales.
Además siempre existe la posibilidad de resolver cualquier ecuación diferencial ordinaria de orden n
transformándola en un sistema de n ecuaciones diferenciales de primer orden.
Abordaremos el problema de encontrar solución a estas ecuaciones diferenciales simultáneas, así como
en Álgebra se resuelven ecuaciones algebraicas simultáneas. Existen distintos métodos de resolución y
la mayor o menor dificultad de los mismos depende de la estructura lineal o no lineal de las ecuaciones.
El desarrollo del estudio de sistemas de ecuaciones diferenciales lineales, ha sido muy beneficiado con
el aporte del Álgebra Lineal. El empleo de la notación matricial y/o vectorial facilita el trabajo y,
especialmente, facilita su generalización para plantear y resolver sistemas n x n. Dado un sistema
dx1
 a11 (t ) x1  a12 (t ) x 2   a1n (t ) x n
dt
dx 2
 a 21 (t ) x1  a 22 (t ) x 2   a 2 n (t ) x n
dt
 
dx n
 a n1 (t ) x1  a n 2 (t ) x 2   a nn (t ) x n
dt
 x1   a11 (t ) a12 (t )  a1n (t )   x1 

    
d  x 2   a 21 (t ) a 22 (t )  a 2 n (t )   x 2 
 42
dt          
    
 x   a (t ) a (t )  a (t )   x 
 n n1
 n2
nn   n
x A( t ) x
podemos expresar el sistema de ecuaciones diferenciales anterior en forma matricial:

dX
 A(t ) X donde: A = [aij ] es la matriz de coeficientes del sistema, X = [ xi(t) ] es el vector
dt
 dx 
columna que expresa las n funciones incógnitas, dX/dt =  i  es el vector columna de las derivadas
 dt 
de primer orden es una ecuación lineal.
La solución de esta ecuación lineal seria X  Ce At , ¿Esto tiene sentido?, siendo

Darwin D. Peña G.
 a11 (t ) a12 (t )  a1n (t )   x1 

   
 a 21 (t ) a 22 (t )  a 2 n (t )   x2 
A ; X   
  
   
 a (t ) a (t )  a (t )  x 
 n1 n2 nn   n
La idea es encontrar un operador transformador que nos solucione este problema es decir,
¿Cuáles vectores v son tales que cuando sufren una transformación A se convierten en un nuevo vector
paralelo al vector original, teniendo magnitud  veces la del vector original?.
Por supuesto, al responder la pregunta anterior estamos buscando vectores no cero o no nulos, a los
cuales se los llama vectores propios, y valores  que se los denomina valores propios asociados a la
matriz A.
Veamos como se soluciona este sistema lineal homogéneo.
Teorema: Sean 1 , 2 ,..., n n valores propios tales que de la matriz A de los coeficientes del sistema
homogéneo (42), y sean K1 , K 2 ,..., K n los correspondientes vectores propios. Entonces la solución de
(42) en el intervalo    t   esta dada por
X  c1 K1e 1t  c2 K 2 e 2t  ...  cn K n e nt 43
Ejemplo 15.1 Resolver (valores propios reales)

dx
 x  3y
dt
Esta ecuación la podemos escribir de forma matricial
dy
 5x  3 y
dx
d  x  1 3   x 
    
dt  y   5 3   y 
1 3 
A   
5 3 
x
X   
 y
Una solución general a este sistema homogéneo es de la forma
X  c1K1e1t  c2 K 2e2t
Donde K1, K2 son los vectores propios asociados a la matriz A y 1 ,  2 son los valores propios asociados
a los vectores propios de la matriz A. como hallar los vectores y los valores propios, veámoslo (repaso
de algebra lineal.
Para hallar los valores propios recurrimos a la ecuación característica
Darwin D. Peña G.
det  A  I   0
1 3 1 0 
det        0
 5 3 0 1
1   3 
det    0
 5 3   
 2  4  12  0
(  2)(  6)  0
Los valores propios son   2,   6
Para hallar los vectores propios asociados a estos valores propios hacemos lo siguiente
Para   2, recurrimos al sistema  A  I K  0 es decir
1 3  1 0    k 1   0 
   2       
  5 3   0 1   k 2   0 
 3 3   k1   0 
      
 5 5  k 2   0 
Es decir
3k1  3k 2  0
5k1  5k 2  0
Luego k1 = -k2 k1 =1, k2 = -1 entonces el vector propio asociado al valor propio   2 es
1 
K 1   
  1
Ahora para   6 tenemos que  A  I K  0 es
1 3  1 0    k 1   0 
   6       
 5 3   0 1  k 2   0 
  5 3   k1   0 
      
 5  3   k2   0
Es decir
 5k1  3k 2  0
3
5k1  3k 2  0  k1  k2
5
3
De aquí tenemos si k2 = 1, k1  entonces el vector propio asociado al valor propio   6 es
5
3
 
K2   5 
1 
 
Entonces la solución del sistema es
X  c1 K1e 1t  c2 K 2 e 2t
Es decir
Darwin D. Peña G.
3
x 1   2 t  
   c1  e  c 2  5 e 6t
 y   1 1 
 
3
x(t )  c1e  2t  c 2 e 6t
Luego 5
 2t
y (t )  c1 e c 2 e 6t
Teorema: Sea A la matriz de los coeficientes del sistema homogéneo (42), y supongamos que A tiene
elementos reales; sea K el vector propio correspondiente al valor propio complejo 1    i ,  y
 reales. Entonces
X 1  K1e 1t y X 2  K 2 e 1t
Son soluciones de (42)
Teorema: Sea 1    i un valor propio complejo de la matriz A de los coeficientes del sistema
1
 i

homogéneo (42), y sean B1  K1  K1 y B2  K1  K1 , entonces
2 2
 
X 1  B1 cos t  B2 sent et y X 2  B1 cos t  B2 sent et son linealmente independientes de
(42)..
Ejemplo 15.2 Resolver (valores propios complejos)

dx
 6x  y
dt
dy
 5x  4 y
dx
det  A  I   0
 6   - 1
det    0
  5 4   
2  10  29  0
Los valores propios son 1  5  2i, 2  5  2i
Para hallar los vectores propios asociados a estos valores propios hacemos lo siguiente
Para   5  2i, recurrimos al sistema  A  I K  0 es decir
1  2i k1  k 2  0
5k1  1  2i k 2  0
Darwin D. Peña G.
1 
Como k 2  1  2i k 2 y después de elegir k1  1 tenemos el primer vector propio K 1  
1  2i  y de
 
1 
la misma manera para   5  2i, tenemos el segundo vector propio K 2  
1  2i  . Por lo tanto la
 
solución es:
 1  5 2i t  1  52i t
X  c1  e  c 2  e
1  2i  1  2i 
Reescriba la anterior solución de la siguiente manera
X  B1 cos t  B2 sent   B2 cos t  B1 sent 
Donde B1 
1
2

K1  K 2  y B2 
i
2

K1  K 2 
Teorema: Si m es un entero positivo y   1 m es un factor de la ecuación característica, pero
  1 m 1
no es un factor, entonces se dice que 1 es un valor propio de multiplicidad m.
Distinguimos dos posibilidades:
a) Puede ser que para algunas matrices A de n x n sea posible encontrar m vectores propios
K1 , K 2 ,..., K m , linealmente independientes, que correspondan a un valor propio 1 de multiplicidad
m  n. En este caso, la solución general del sistema contiene la combinación lineal
c1 K1e  c2 K 2 e 1t  ...  cm K m e 1t .
1t
b) Si al valor propio 1 de multiplicidad le corresponde solamente un vector propio, entonces siempre

es posible encontrar m soluciones linealmente independientes de la forma
X 1  K 11e 1t
X 2  K 21te 1t  K 22e 1t

t m 1 t m2
X m  K m1 e 1t  K m 2 e 1t    K mm e 1t
m  1! m  2!
En donde los K ij son vectores columna.
Ejemplo 15.3 Resolver (valores propios repetidos)

dx
 3 x  18 y
dt
dy
 2x  9 y
dx
Darwin D. Peña G.
det  A  I   0
 3   - 18 
det    0
 2  9   
  32  0
 3   3t
Los valores propios son 1  2  3 luego la solución es X 1   e .
1 
Una segunda solución es de la forma X 2  Kte 1t  Pe 1t 44 donde K y P son vectores columnas.
dX
Para ver que (44) es solución los sustituimos en el sistema  A(t ) X y simplificando tenemos
dt
 AK  1 K te  t   AP  1 P  K e  t
1
0 1
Como esta última igualdad debe cumplirse para todos los valores de t, se debe tener
 A  1 I K  0
 A  1 I P  K 45
 3   3t  3 p 
En el ejemplo 15.2 tenemos que X 1  1 e identificamos K  1  y P   1  se tiene de
     p2 
(45) que  A  3I P  K o bien
 6  18  p1   3  6 p1  18 p 2  3
      es decir
 2  6  p 2   1  2 p1  6 p 2  1
1
   3 1
haciendo p1 
1
tenemos p 2  0 , P  2 luego X 2   te 3t   2 e 3t por lo tanto la
0   
2   1 0
 3   3t  3  1 
solución general es c 
X  1 
1  e   3t 
+ c 2  te  2 e   3t 
   1    
 0 
Ejemplo 15.4. El tanque A contiene 50 galones de agua en los cuales se disuelven 25 libras de sal, un
segundo tanque B contiene 50 galones de agua pura. Se bombea liquido desde y hacia los tanques en las
proporciones indicadas en la figura deduzca las ecuaciones diferenciales para x(t), y(t).
Agua pura mezcla

gal gal
3 1
min min
Darwin D. Peña G.
Mezcla mezcla
gal gal
1 3
min min
dA
 (Rapidez con que entra la sal) – (Rapidez con que sale la sale)
dt
dx1  gal   gal  x2 lb   gal  x1 lb  x2 2

 3 0  1     4     x1
dt  min   min  50 gal   min  50 gal  50 25
dx2 gal x1 lb  gal x2 lb gal x2  2 2

4   1  3    x1  x2
dt min 50 gal  min 50 gal min 50  25 25
dx1 2 1 dx2 2 2
  x1  x2  x1  x2
dt 25 50 dt 25 25
 2 1  2 1
   0   
A   25 50  I    dt  A  I   25 50 0
 2
  
2 0  2

2

 25 25  25 25
 2   2    2  1   2  2  1
                0          0
 25   25    25  50   25  25  625
4 4 1 4 3
   2   0  2    0
625 25 625 25 625
4 16 12 4 4 4 2 1
       1  
 25 625 625  25 625   25 25  25
2 2 2 3
2  
25
Darwin D. Peña G.
 1 1 
 
Con   
1
 25 50  k1    0 
 2 1    
25
   k2   0 
 25 25 
1 1
 k1  k2  0 k2  2k1 Con k1  1 y k2  2
25 50
2 1 1
k1  k2  0 k1   
25 25 0
 1 1 
3   k   0 
Con     25 50  1    
25  2 1  k2   0 
 
 25 25 
1 1
k1  k2  0 k2  2k1 Con k1  1 y k2  2
25 50
2 1  1 
k1  k2  0 k2   
25 25   2
 x1  1   1  
t 3t
   C1  e 25  C2  e 25
 x2   2   2
t 3t
 
x1 t   C1e 25
 C2e 25
t 3t
 
x2 t   2C1e 25
 2C2e 25
25  C1  C2 2 50  2C1  2C2

0  2C1  2C2 1 0  2C1  2C2
25 25
50  4C1 C1  C1  C2 
2 2
t 3t t 3t
25  25 25  25  
x1 t   e  e ; x2 t   25e 25
 25e 25
21 2
Ejemplo 15.4. Resolver el siguiente circuito.

Darwin D. Peña G.
Tomando a E como αX (t) y a I como y(t) procedemos a obtener el sistema de ecuaciones

Malla exterior
Malla interior
Sistema de ecuaciones
Ceros de la función:
Calculando K
Darwin D. Peña G.
Para λ=0 y k2=1
Para λ y k2=1
Para los valores

α = 100v
l1 = 10mH
l2 = 12,5 mH
R = 5Ω
Darwin D. Peña G.
15. EJERCICIOS
Resolver
 dx
 dx  dt  4 x  y  z
 dt  2 x  3 y 
 dy
1.  2.   2  5 y  z
 dy  2 x  y  dt
 dt  dz
 dt  y  3z

 dx
 dt  x  y  z  dx
  dt  x  2 y
 dy
3.   2 y 4. 
 dt  dy  4 x  3 y
 dz  dt
 dt  y  z

 dx  dx
 dt  6 x  y  dt   x  3 y
5.  6. 
 dy  5 x  2 y  dy  3 x  5 y
 dt  dt
Darwin D. Peña G.
BIBLIOGRAFIA
Dennis G. Zill, Ecuaciones diferenciales con aplicaciones. Grupo Editorial Iberoamericana. 1986
Frank Ayres, Jr, Ecuaciones diferenciales. Mc Graw Hill. 1991.
George F. Simmos, Jhon S. Robertson, Ecuaciones diferenciales con aplicaciones y notas históricas. Mc
Graw Hill. 1993.

Módulo Ecuaciones

Cargado por

Copyright:

Formatos disponibles

Módulo Ecuaciones

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Módulo Ecuaciones

Cargado por

Copyright:

Formatos disponibles

Introducción a las Ecuaciones

Compilado por el docente

UNIVERSIDAD DEL MAGDALENA

Las Ecuaciones Diferenciales se encuentran inmersas en una rama de las matemáticas

1. ORIGEN DE LAS ECUACIONES DIFERENCIALES

1.1 Problemas geométricos

1.2 problemas físicos

En la figura (2) vemos que la componente tangencial de la fuerza debida al peso W es mg

Ejemplo 1.4 halle la variación de la atmósfera con la altura en la atmósfera terrestre

para hallar la densidad sabemos que pV  Rnt ,luego:

p1  y1 Son condiciones iniciales, por tanto p1  p o y y1  0

Ejemplo 1.7 Obtener la ecuación diferencial asociada con la primitiva

Solución: Derivando tenemos

1.1 ¿Cuál es la ecuación diferencial de la velocidad v de un cuerpo de masa m que cae

1.3 Halle la ecuación diferencial de la familia de circunferencias de radio r cuyos centros

1.5 Obtenga la ecuación diferencial de la familia de curvas. x 2  y 2  1

1.6 Obtenga la ecuación diferencial de la familia de curvas. y  e x cos x  e x sin x

1.7 Hallar la ecuación diferencial de todas las circunferencias del plano.

2.1 Clasificación de las ecuaciones diferenciales según su tipo

2.2 Clasificación de las ecuaciones diferenciales según el orden

2.3 Clasificación de las ecuaciones diferenciales según la linealidad o no linealidad.

Se dice que una ecuación diferencial es lineal si tiene la forma

Solución: Derivando 2 veces tenemos

Ejemplo 2.8 Verificar si x = 4t +1 (2), y = t2 – 2 (3) es una solución para la ecuación

Reemplazando en (5) obtenemos

3. ECUACIONES DIFERENCIALES DE PRIMER ORDEN

El problema de las ecuaciones diferenciales consiste esencialmente en encontrar la

Definición 3.1 Se dice que una ecuación diferencial de la forma

Ejemplo 3.2 Resolver (1 + x )dy – ydx = 0

Ahora al sustituir x  0, y  0 se tiene

250  300  270e T

El sistema se puede escribir en la forma:

La ecuación diferencial es de variables separables y su solución general viene dada por:

Definición 4.1 Una función f (x, y) se llama homogénea de grado n si

Ejemplo 4.1 f ( x, y)  x 4  x 3 y es homogénea de grado 4, ya que

Ejemplo 4.2 f ( x, y)  x 2  sin x cos y no es homogénea, ya que

Ahora por la homogeneidad de M y N se es posible escribir

x n M (1, u )dx  x n N (1, u )udx  xdu   0

4.1 Forma alternativa

Ejemplo 4.4 Utilizando la forma alternativa en el ejemplo anterior tenemos

4.2 Resuelva las siguientes ecuaciones diferenciales

Una ecuación diferencial de la forma

5. ECUACIONES DIFERENCIALES EXACTAS Y REDUCIBLES A EXACTAS

La condición necesaria y suficiente de que

M ( x, y)dx  N ( x, y)dy  0 (6)

5.1 Método de solución

Dada la ecuación M ( x, y)dx  N ( x, y)dy  0

Ejemplo 5.1 Resolver 2 xydx  ( x 2  1)dy  0

5.2 Factores integrantes

Es un factor integrante. Luego

5.1 Resolver las siguientes ecuaciones exactas

5.2 Encuentre el factor integrante para las siguientes ecuaciones y resolverlas

La forma general de la ecuación lineal de orden n es :

Veamos el caso n = 1 es decir

Ejemplo 6.2 Hallar una solución continua que satisfaga

Solución: Resolvemos el problema en dos partes

Como es continua, por lo visto en el calculo diferencial

Problemas de crecimiento y decrecimiento

15. y  (tan x) y  cos 2 x, y(0)  1