Ecuaciones Diferenciales Tomo I
Ecuaciones Diferenciales Tomo I
Ecuaciones Diferenciales Tomo I
y sistemas dinámicos
Jorge Augusto Pérez Alcázar
Tomo I
Tomo I
Ediciones
Edición
Gerencia de Investigaciones
Gerente de Investigaciones
Carolina Mejía Corredor
Coordinadora de Publicaciones
Laura Cediel Fresneda
Revisor de estilo
Juan Carlos Velásquez
Diseño, diagramación
Cesar Rubiano
ISBNe: 9789587566055
© Universidad EAN, El Nogal: Cl. 79 No. 11 - 45. Bogotá D.C., Colombia, Suramérica, 2019
Prohibida la reproducción parcial o total de esta obra sin autorización de la Universidad
EAN©
© UNIVERSIDAD EAN: SNIES 2812 | Personería Jurídica Res. nº. 2898 del Minjusticia -
16/05/69| Vigilada Mineducación. CON ACREDITACIÓN INSTITUCIONAL DE ALTA CALI-
DAD, Res. N° 29499 del Mineducación 29/12/17, vigencia 28/12/21
Producido en Colombia.
Contenido
JUSTIFICACIÓN 11
CAPÍTULO 1
INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES
CAPÍTULO 2
CONCEPTOS BÁSICOS (IDEAS CUALITATIVAS)
CAPÍTULO 3
ECUACIONES DIFERENCIALES AUTÓNOMAS
6
Contenido
CAPÍTULO 4
SOLUCIONES NUMÉRICAS Y CUANTITATIVAS DE UNA
ECUACIÓN DIFERENCIAL
CAPÍTULO 5
MÁS SOBRE MÉTODOS CUANTITATIVOS
7
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
CAPÍTULO 6
ECUACIONES DIFERENCIALES LINEALES DE PRIMER
ORDEN Y MODELOS
CAPÍTULO 7
APLICACIONES DE LAS ECUACIONES DE PRIMER ORDEN
8
Contenido
CAPÍTULO 8
ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN
Y DE GRADO SUPERIOR
9
JUSTIFICACIÓN
JUSTIFICACIÓN
11
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
12
JUSTIFICACIÓN
13
CAPÍTULO 1.
INTRODUCCIÓN
A LAS
ECUACIONES
DIFERENCIALES
Contenido
Competencias
1. Identificar el concepto de cambio
y sistemas dinámicos.
Introducción
16
Introducción a las ecuaciones diferenciales
17
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Hace muchos siglos, las personas creían que la energía solar era
causada por un dragón tratando de tragarse el sol y que los terremotos
eran causados por dioses descontentos en el inframundo. Estas
teorías fueron construcciones mentales que les ayudaron a explicar
el mundo que les rodeaba. La historia humana es un largo recuento
de esas visiones interesantes de la realidad, que ilustran una vieja
creencia de que los fenómenos pueden entenderse en comparación
con el sistema más sencillo basado en la versión limitada del objeto
real. Dragones chinos míticos y dioses mayas del inframundo eran
modelos religiosos. Por otra parte, los antiguos griegos, después de
observar el cielo desde hace muchos siglos construyeron un modelo
del universo físico que consiste en una gran esfera hueca con la tierra
en su centro y los cuerpos celestes que se mueven a lo largo de varios
caminos en su superficie. Si bien este modelo mostró razonablemente
buen acuerdo con los datos disponibles en ese momento, desde
entonces ha sido sustituido por una descripción más precisa y com-
pleta del universo físico.
18
Introducción a las ecuaciones diferenciales
19
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
20
Introducción a las ecuaciones diferenciales
Aunque la física moderna sugiere que todos los sistemas son en última
instancia discretos y los cambios ocurren en saltos en puntos distintos
en el tiempo, la escala de tales cambios los hace más comprensibles
al ser tratados como que cambian de forma continua con el tiempo.
Incluso una población de moscas de la fruta es tratada por el biólogo
como un cambio continuo. Más comúnmente estos sistemas de tiempo
continuo son modelados por ecuaciones diferenciales, ejemplos de
estos hechos se verán a lo largo de este libro. Cuando los intervalos
de tiempo son muy grandes, los cambios pueden ser tratados como
si fueran separados a saltos; variaciones diarias, semanales o anuales
en las variables económicas tales como precios de las acciones o
de los ingresos fiscales, o variables ecológicas como las poblaciones
anuales de animales o plantas. Para este tiempo discreto o datos
muestreados, los sistemas de un modelo matemático útil son la
ecuación en diferencia.
21
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
22
Introducción a las ecuaciones diferenciales
A. Modelos de mecánica
Figura 2. Galileo Galilei (1564-1642).
23
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
En mecánica
Isaac Newton fue el responsable de un gran número de descubri-
mientos en física y matemáticas, pero quizás los tres más importan-
tes son los siguientes:
F
m
d dv
(mv)=m =ma (1)
dt dt
(3)
m
m
M
M
1 La primera edición del Principia fue muy pequeña, no superó 200 ediciones.
25
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
x
Fuente. Elaborada por el autor.
- mg = ma = m =m (5)
(6)
26
Introducción a las ecuaciones diferenciales
v=
27
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
B. Modelo poblacional
Consideremos una población P (t) que va variando en el tiempo(2).
Matemáticamente se puede decir que la razón por la cual la población
está cambiando respecto al tiempo está dada por la derivada:
(7)
2 Por el momento, la población puede ser cualquier cosa: humanos, paramecios, mariposas, etc.
28
Introducción a las ecuaciones diferenciales
ó (8)
29
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
D. Modelos geométricos
Algunos problemas de naturaleza geométrica requieren de la solución
de una ecuación diferencial para su respuesta. Recordemos que para
una curva plana con ecuación y = f(x) el ángulo para el cual la
tangente a la curva en el punto (x, y) forma con el eje x, está dado por
la relación:
tan =
30
Introducción a las ecuaciones diferenciales
P (x,y)
1
y– 2
y
M
x
T x
(x, f (x))
tan0)
y=f(x0 ) + m(x–x =
(x0 , f (x0 ))
y = f (x)
Y la ecuación diferencial del problema toma la forma:
31
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(9)
R=k (10)
32
Introducción a las ecuaciones diferenciales
(11)
F(x) = (12)
33
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
20
+
15
10
5
x
-20 -15 -10 -5 5 10 15 20
-5
-10
-3 -2 -1 1 2 (13) 3
-1
Y queremos encontrar la antiderivada más general, F(x) de la
-2
F(x) = ln +C
O como
F(x) =
34
y 25 Introducción a las ecuaciones diferenciales
20
+
15
10
translación
-20 vertical.-10Lo mismo
-15 -5
es verdadero
5
para x10 < 0. 15 20
-5
-10
Figura 10. Antiderivadas de la función f(x) = 1/x.
-3 -2 -1 1 2 3
-1
-2
-3
T - Tm o, (14)
35
1.4 Ejercicios
La frase «y es proporcional a x» implica que y está relacionado
con x mediante una ecuación de la forma y = k x, donde K es una
constante. De manera similar «y es proporcional al cuadrado de x
implica que: y = kx2, y es proporcional al producto de x y z» implica
y = k x z, y «y es inversamente proporcional al cubo de x» implica
y = k/x3. Por ejemplo, cuando Newton propuso que la fuerza de
atracción de un cuerpo a otro es proporcional, el producto de las
masas e inversamente proporcional al cuadrado de la distancia entre
ellos, podemos escribir inmediatamente:
F=
1.5 La derivada
Figura 11. Pierre de Fermat (1601-1665).
39
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
f(x)
C 0
X 1
nxn-1
cos (x) - sen (x)
Sen (x) Cos (x)
ex ex
1/x x≠0
A. La razón de cambio
En cálculo aprendemos que una función tiene una razón de cambio
instantánea, y esta razón es igual a la derivada. Por ejemplo, si tenemos
40
Introducción a las ecuaciones diferenciales
la función distancia x (t) medida desde un punto fijo sobre una línea,
entonces la razón a la cual cambia x (t) con respecto al tiempo es la
velocidad v. Entonces se tiene que:
v (t) =
a (t) = .
Sin embargo, hay casos donde las cosas pueden cambiar de-
pendiendo de otros parámetros, como se verá más adelante. En este
texto veremos estos ejemplos y muchos otros. El punto es que todas
las cantidades cambien, la razón por la cual cambia es la derivada de
dicha cantidad. Es un hecho importante hacer el proceso de modela-
miento y realizar el estudio de la ecuación diferencial, por esta razón
41
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y=f (15)
42
Introducción a las ecuaciones diferenciales
se van acercando.
L(x) = (16)
Este límite nos dice que R(x) es cada vez más pequeño a medida
que . En efecto, este hecho nos dice que R(x) =0 es mucho más
pequeño que , ya que dicho cociente tiende a cero. Esto también
nos dice que la función L(x) definida anteriormente, es la única función
lineal que tiene esta propiedad. Esto es lo que significa cuando decimos
que L es la mejor función lineal que aproxima a la función no lineal f.
También se puede observar que la gráfica de la recta tangente en la
figura 12 es la gráfica de L(x). La figura 12 también nos muestra que
L(x) es una buena aproximación para f(x), siendo x cercano a x0.
43
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
P (x,y)
1
D. El límite de cocientes de diferencias
2
y– y
M
x
Consideremos el cociente de diferencias:
T (18) x
Este cociente es igual a la pendiente de la recta que pasa por los
puntos (x0, f(x0)) y (x, f(x)), tal como se muestra en la siguiente figura.
(x, f (x))
y=f(x0 ) + m(x–x0)
(x0 , f (x0 ))
y = f (x)
(19)
44
Introducción a las ecuaciones diferenciales
E. La tabla de derivadas
Mediante la memorización de una tabla de derivadas y el manejo de
unas pocas reglas, especialmente la regla de la cadena, podemos
lograr las habilidades de la diferenciación. No es difícil ver que se
puede calcular la derivada de cualquier función. Esta habilidad es muy
importante. Muchas veces, es claro que esta definición formulista de
derivada es diferente de las dadas previamente.
45
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
46
1.6 Ejercicios
Recordemos las siguientes reglas de derivación de su curso de cálculo.
Sean f y g funciones diferenciables, entonces:
d.
7. f(x) = 8. F(x) =
=0
Para cada una de las siguientes funciones, trace la gráfica
de f y su Linealización L(x) = f(x0) + f/(x0) (x – x0).
16. f(x) = , si x0 = 1
18. f(x) = , en x0 = 1
20. f(x) = , en x0 = 0
49
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
f(x) f(x)
0 C Cos(x) Sen(x) + C
1 X+C Sen(x) -cos(x) + C
Ln |x| + C
50
Introducción a las ecuaciones diferenciales
B. La antiderivada
f / = g sí y solo sí (20)
(22)
C. La tabla de fórmulas
51
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
52
Introducción a las ecuaciones diferenciales
y/ = f(t, y) (23)
y/ = f(t)
y (t) = (24)
Ejemplo 1:
y/ = cos t
Solución:
De acuerdo a (24), la solución es:
y (t) = (25)
53
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
54
Introducción a las ecuaciones diferenciales
Ejemplo 2:
Solución:
Este es un ejemplo de un problema con condición inicial. Este
tipo de problema siempre requiere encontrar una solución particular
que satisface la condición inicial y (0) = 2. De acuerdo a (21) la solución
general para la ecuación diferencial está dada por:
y (t) = (26)
(27)
55
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
2 = y (0) = e0 (0 – 1) + C = - 1 + C
y (t ) = et (t – 1) + 3 (29)
y (t) – y(0) =
Entonces:
ueu du
=2+
=
56
Introducción a las ecuaciones diferenciales
Ejemplo 3:
y/ = , con y(1) = 3
Solución:
Aquí, tomamos a x como la variable independiente. Integrando,
encontramos que:
y (x) = +C
3 = y (1) = ln (1) + C = C.
Entonces, C = 3.
57
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Podemos resolver esta ecuación usando el método dado en esta
sección. Primero introducimos la velocidad, para reducir la ecuación de
segundo grado a un sistema de dos ecuaciones de primer orden, así:
(30)
v (t) = - g t + C1
58
Introducción a las ecuaciones diferenciales
59
1.9 Ejercicios
En el siguiente ejercicio, encuentre la solución general de la ecuación
diferencial. En cada caso chequee al menos seis de las ecuaciones de
las curvas solución.
1. y/ = 2 t + 2 2. y/ = 3 t2 + 2 t + 3.
5. y/ = 6. y/ =
7. y/ = 8. y/ = t cos(3t)
13. r/ = 14. y/ =
Observación: para resolver los ejercicios 13 y 14 se debe usar
fracciones parciales.
29-1 29-2
29-3 29-4
29-5 29-6
29-7 29-8
29-9 29-10
29-11 29-12
29-13 29-14
29-15 29-16
29-17 29-18
29-19 29-20
29-21 29-22
29-23 29-24
29-25 29-26
29-27 29-28
29-29 29-30
29-31 29-32
29-33 29-34
29-35 29-36
29-37 29-38
29-39 29-40
29-41 29-42
29-43 29-44
29-45 29-46
29-47 29-48
29-49 29-50
29-51 29-52
29-53 29-54
29-55 29-56
29-57 29-58
29-59 29-60
29-61 29-62
29-63 29-64
29-65 29-66
29-67 29-68
29-69 29-70
29-71 29-72
29-73 29-74
29-75 29-76
29-77 29-78
29-79 29-80
29-81 29-82
29-83 29-84
29-85 29-86
29-87 29-88
29-89 29-90
29-91 29-92
29-93 29-94
29-95 29-96
29-97 .
29-98
29-99 29-100
29-101 29-102
29-103 29—104
29-105 29-106
29-107 29-108
29-109 29-110
20-111
29-112
29-113
CAPÍTULO 2.
CONCEPTOS
BÁSICOS
(IDEAS CUALITATIVAS)
Contenido
Competencias
1. Deducir el concepto de
ecuación diferencial.
68
Conceptos Básicos (Ideas Cualitativas)
a. . b. c. = y.
69
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Introducción
La idea de esta sección es formalizar el concepto de lo que es una
ecuación diferencial y desarrollar métodos para analizar ecuaciones
diferenciales ordinarias, aprovechando todas las ideas dadas en la
sección anterior, en donde se vio la potencia del cálculo diferencial
e integral. Lograremos este objetivo abordando los diferentes
temas desde los puntos de vista geométrico (gráfico), numérico y
analítico. Debido a que cada una de estas perspectivas en ocasiones
proporciona información incompleta, siempre necesitamos comparar
la consistencia de los resultados.
70
Conceptos Básicos (Ideas Cualitativas)
Ejemplo 1:
(1)
y(x) = (2)
(3)
Donde g(x, y) es una función de las variables x, y, y.
72
Conceptos Básicos (Ideas Cualitativas)
(4)
Observaciones:
Debido a que una solución explícita tiene una derivada en el
intervalo a < x < b, debe ser continua en ese intervalo, ya que tiene
derivada. Una solución explícita nunca posee una tangente vertical (ya
que tiene derivada en dicho intervalo).
73
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
74
Conceptos Básicos (Ideas Cualitativas)
75
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y (x) = (5)
(6)
y (0) = 0 (7)
(8)
76
Conceptos Básicos (Ideas Cualitativas)
Ejemplo 3:
dy / dx = y/ = , (9)
Solución:
Al tener la derivada 1/x de la función que estamos buscando, y
recordando nuestro curso de cálculo podemos analizar:4
y = ln + C, o y (x) =
4 Este método cualitativo no es más que recordar la aplicación de la derivada del curso de cálculo diferencial.
77
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
78
Conceptos Básicos (Ideas Cualitativas)
dy / dx = y/ =
79
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y// =
Solución:
• Monotonía: la derivada de y que es para todo valor
de x representa un numero positivo, por lo tanto la función y
es siempre creciente.
• Concavidad: si derivamos de nuevo a dy/dx se tiene:
81
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
82
Conceptos Básicos (Ideas Cualitativas)
y/ = f (x, y)
5 Estos conceptos son clave en el análisis cualitativo de una ecuación diferencial ordinaria. Estos son los mis-
mos elementos que usamos para la gráfica de funciones en el Cálculo. Una completa referencia de estos métodos
y pruebas son dadas por J. H. Hubbard & B. H. West Differential Equations: A Dynamical Systems Approch. Springer
Verlang, 1989.
83
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 5:
(10)
Solución:
Esta ecuación describe la función cuya razón de cambio es
siempre 1. Con los conocimientos de cálculo, tracemos algunas curvas
solución de (10). Debido a que el lado de la derecha es positivo, se
tiene que las soluciones son crecientes en todo punto de su dominio.
También se tiene que para todos los puntos (x, y), la solución de esta
ecuación diferencial tiene una línea tangente cuya pendiente es 1.
Para transferir esta información a una gráfica podemos seleccionar
varias coordenadas (x, y) y trazar segmentos cortos de recta con
pendiente 1, tal como se muestra en la siguiente figura.
m=1 P (x, y)
84
Conceptos Básicos (Ideas Cualitativas)
85
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 6:
(11)
Solución:
Esta ecuación diferencial modela la situación en donde la razón
de cambio de la función desconocida es igual al valor de la variable
independiente.
87
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
88
Conceptos Básicos (Ideas Cualitativas)
Ejemplo 7:
(12)
89
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
De acuerdo al análisis previo, sabemos que la gráfica de
cualquier curva solución, siempre tiene concavidad hacia abajo, forma
decreciente para x < 0 y con forma creciente para x > 0. También
sabemos que la familia de soluciones es simétrica respecto al eje y.
90
Conceptos Básicos (Ideas Cualitativas)
91
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Es decir:
92
Conceptos Básicos (Ideas Cualitativas)
En resumen:
Definición
Una isóclina correspondiente a la pendiente m de la ecuación
diferencial y/ = g(x, y), es la curva caracterizada por la ecuación
g(x, y) = m.
Procedimiento
Escriba:
1. g (x, y) = m ( 13 )
2. Donde m es una constante.
3. Dé valores diferentes para m. Para cada valor de m resuelva
dicha ecuación para y en términos de x y de m, o para x en
términos de y, y, m. Si no se puede resolver, intente identificar
las curvas definidas implícitamente por la ecuación (13).
Esto da la isoclina correspondiente a la pendiente de m.
4. Para cada valor de m, grafique la isóclina correspondiente a
la pendiente m.
Si se construyen campos de pendientes a mano, dibuje seg-
mentos de recta cortos con Pendiente m, que crucen a la
isóclina apropiada.
93
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La isóclina correspondiente a la pendiente m está dada por:
(14)
Figura 14. Campo de pendientes para la función de error con algunas isóclinas.
94
2.5 Ejercicios
Resuelva las siguientes ecuaciones diferenciales de primer orden.
Esboce la solución explícita para tres diferentes valores de la constante
arbitraria C. Luego encuentre el valor específico de C y la fórmula y(x)
para la solución particular que pasa por el punto dado.
1. P (1,1)
2. P (1,1)
3. P (0,0)
4. P (0,2)
5. P (0,1)
6. P (1,1)
7. P (-1,1)
8. P (1, )
9. P (2,1)
10. P (1,1)
11. P (0,1)
12. P (0,1)
G(x) =
a) P (1, 1)
b) P (1, 1)
c) P (0, 0)
d) P (0, 2)
e) P (0, 1)
f) P (1, 1)
g) P (-1,1)
h) P(1, )
i) P (2, 1)
j) P (1,1)
k) P (0, 1)
l) P (0, 1)
es errónea.
19. Para la ecuación diferencial:
y/ =
y/ = , y(x0) =0
y/ = G(x) =
a) P (1, 1)
b) P (1, 1)
c) P (0, 0)
d) P (0, 2)
e) P (0, 1)
f) P (1, 1)
g) P ( -1,1)
h) P(1, )
i) P (2, 1)
j) P (1,1)
k) P (0, 1)
l) P (0, 1)
y/ = x + 1, y/ = x – 1, y/ = ln , y/ = x2 – 1.
Donde g = 9.8 y k = 0.2 seg-1. Aquí y(x) es la distancia
recorrida en caída libre al tiempo x; de modo que y(x) = 0. Si
este objeto cae desde 5000 metros de altura, calcule cuántos
segundos recorre antes de que impacte con la superficie de la
tierra por medio de:
a) Campos de pendientes, monotonía y concavidad.
b) Encontrando la solución explícita.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Definición
Una ecuación diferencial de primer orden de la forma:
Ejemplo 9:
(15)
104
Conceptos Básicos (Ideas Cualitativas)
Solución:
Veamos qué información cuantitativa podemos obtener de las
soluciones de (15), aplicando los conceptos del cálculo.
y// = - (16)
y// = (17)
105
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
106
Conceptos Básicos (Ideas Cualitativas)
Figura 19. Campo de pendientes de la ecuación autónoma y´= con algunas curvas
soluciones.
y (x) =
y y/ = 1
107
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y (x) = (18)
108
Conceptos Básicos (Ideas Cualitativas)
Figura 20. Campo de pendientes de la ecuación autónoma y´= con algunas cur-
vas soluciones particulares.
Ejemplo 10:
(19)
Solución:
Observe que el lado derecho no depende de la variable indepen-
diente x.
109
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y// = (20)
y// = (21)
= 16.y. (y - 1) (y – ½)
110
Conceptos Básicos (Ideas Cualitativas)
111
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
112
Conceptos Básicos (Ideas Cualitativas)
6 También llamado Modelo de Malthus, del crecimiento poblacional. Este artículo es incluido por J. R. N Newman,
en El mundo de la Matemática, volumen 2, 1996.
113
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Nuestra hipótesis no dice que la razón de cambio de la población
solo depende del tamaño de la misma. En este modelo de crecimiento
no se tiene en cuenta las limitaciones de espacio o recursos naturales.
Por tal razón, a este modelo se le puede considerar muy incipiente.
Ahora, si llamamos:
(22)
114
Conceptos Básicos (Ideas Cualitativas)
115
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
ln (23)
P (t) = (24)
116
Conceptos Básicos (Ideas Cualitativas)
117
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Definición
Si la ecuación diferencial y/ = g(x, y) tiene una solución de la
forma y(x) = constante, la solución recibe el nombre de solu-
ción de equilibrio.
Definición
Dada la ecuación diferencial y/ = g(x, y), se trata de una ecua-
ción que tenga una solución de equilibrio, es decir, una solu-
ción de la forma y(x) = constante. Una solución de equilibrio
se denomina estable si todas las soluciones de la ecuación
118
Conceptos Básicos (Ideas Cualitativas)
Ejemplo 12:
Solución:
Si tomamos k = 1 entonces nuestra ecuación autónoma se
convierte en:
119
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Para este modelo tenemos de nuevo:
7 El modelo original de Malthus es usado solo para situaciones muy simples. El modelo logístico nos permite
representar situaciones más reales.
120
Conceptos Básicos (Ideas Cualitativas)
(algo) =
(25)
Este modelo es llamado el modelo logístico de la población con
velocidad de crecimiento k y capacidad de soporte M. Es una ecuación
diferencial de primer orden. Se dice que no es lineal porque su lado
derecho no es una función lineal de P (t), como lo era el modelo de
crecimiento exponencial visto en el ejemplo 11.
121
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
, de donde,
122
Conceptos Básicos (Ideas Cualitativas)
O M
P (t)
- - + +
+ + + -
(1 - )
- + -
123
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Figura 27. Gráfica del comportamiento de las curvas solución para la función
logística.
124
Conceptos Básicos (Ideas Cualitativas)
(26)
Solución:
El campo de pendientes de (26) se muestra en la siguiente
figura.
126
Conceptos Básicos (Ideas Cualitativas)
(27)
127
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y (x) =
128
Conceptos Básicos (Ideas Cualitativas)
129
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 15:
(28)
Solución:
Aquí g(x, y) = y2, de modo que . Por lo tanto g(x, y)
y, son continuas para todo x y toda y. El teorema establece que
existe una solución única que satisface la condición inicial y(x0) =y0.,
en otras palabras, existe una curva única que pasa a través del punto
y(x0) = y0. Esto significa que las curvas solución no se intersecan,
debido a que de otra manera habría dos soluciones distintas con la
misma condición inicial en el punto donde se intersectarían. Ahora la
función y(x) = 0 es una solución de (28) de modo que ninguna otra
solución puede intersectarla.
130
Conceptos Básicos (Ideas Cualitativas)
Ejemplo 16:
Solución:
La figura 31 ilustra el campo de pendientes de (28), y la figura
32, una curva solución a través del punto (1,1).
131
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
132
Conceptos Básicos (Ideas Cualitativas)
y (x) = 1
133
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 17:
y (x0) = y0 (30)
y// =
134
Conceptos Básicos (Ideas Cualitativas)
= (31)
y (x) = (32)
135
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Figura 35. Gráfica del campo de pendientes con dos curvas solución para (32).
136
Conceptos Básicos (Ideas Cualitativas)
Para ilustrar lo que ocurre con otro punto inicial, notemos que las
soluciones explícitas para la condición inicial y (2) = 1 es , válida
para x < 4. Esta solución también se muestra en la figura (32). Este
ejemplo nos muestra que se debe tener cuidado cuando se resuelve
y(x). Es un error común introducir funciones que no constituyen
soluciones.
Definición
Una solución particular de y/ = g(x, y), que no se obtiene de
la familia de soluciones (que contenga la constante arbitraria
C) seleccionando un valor finito para C, se denomina solución
singular.
Ejemplo 18:
y/ = (33)
137
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Tenemos que, tanto como su derivada - y , son
continuas en la vecindad de (0, 0), así que no se tiene ningún problema
con la unicidad de la solución.
Arcsen y = x + C, (34)
Arcsen y = x (35)
De donde se obtiene:
y = sen(x) (36)
138
Conceptos Básicos (Ideas Cualitativas)
y = sen(x) donde .
139
2.9 Ejercicios
Análisis cualitativo
1. Para cada una de las siguientes ecuaciones diferenciales:
a. y/ = - y2 b. y/ = -1 - y2
c. y/ = y2 d. y/ = y1/5
e. y/ = y3/2 f. y/ =
Figura 38. Campo de pendientes para alguna de las funciones del ejercicio 4.
c. y/ = - d. y/ = -
Figura 39. Campo de pendientes para alguna de las funciones del ejercicio 5.
Figura 40. Campo de pendientes para alguna de las funciones del ejercicio 5.
Modelos simples
Existencia y unicidad
23. ¿Tiene y/ = 3y – 2 algunas soluciones de equilibrio? ¿Para
qué condición(es) inicial(es) sería esta solución de equilibrio la
única solución de un problema con condición inicial?
ECUACIONES
DIFERENCIALES
AUTÓNOMAS
Contenido
3.1 Comportamiento cualitativo de
soluciones utilizando puntos de equilibrio
y líneas de fase 156
3.2 Ejercicios 185
3.3 Diagramas de bifurcación 194
3.4 Ejercicios 208
Competencias
1. Analizar una ecuación diferencial
autónoma de manera cualitativa.
Introducción
Existe otra manera de analizar las ecuaciones diferenciales autóno-
mas, aplicando lo que se conoce como análisis de la línea de fase.
Este método no ofrece una solución explícita, pero proporciona el
comportamiento cualitativo de las soluciones con muy pocos cálculos.
La clave de esta técnica es el teorema de existencia y unicidad, el cual
garantiza que las soluciones no pueden intersecarse. Veamos unos
ejemplos.
155
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
3.1 Comportamiento
cualitativo de soluciones
utilizando puntos de
equilibrio y líneas de fase
y/ = y (2- y) (1)
Solución:
Primero hallamos la solución de equilibrio de (1), las cuales son y(x)
= 0, y y(x) = 2. Por el teorema de existencia y unicidad sabemos que
ninguna otra solución puede cruzarse con cualquiera de estas dos
soluciones, por esta razón nos concentramos en las tres regiones
formadas por estas soluciones de equilibrio, a saber y < 0, luego 0 < y
< 2 y finalmente y > 2.
156
Ecuaciones Diferenciales Autónomas.
y(x) = 0 y(x) = 2
157
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y
y=2
y=0
x
En resumen:
• Una manera simple de construir la línea de fase para y/ = g
(y), consiste en trazar la gráfica de g(y), para el ejemplo g(y)
= y(2 – y), tomando y como el eje horizontal.
158
Ecuaciones Diferenciales Autónomas.
g(y)
1
–1 0 1 2 3 y
–1
–2
159
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
160
Ecuaciones Diferenciales Autónomas.
y/ = y2 (3)
161
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Observación
Cuando se analice la línea de fase, hay que darse cuenta si el
teorema de existencia y unicidad se satisface, ya que ninguna solución
interseca las soluciones de equilibrio.
• Dibuje la línea.
• Encuentre los puntos de equilibrio (los números tales que f
(y) = 0, y márquelos sobre la línea.
• Encuentre los intervalos de valores y, para los cuales f (y)
> 0, y dibuje flechas que señalen hacia arriba (o hacia la
derecha) sobre esos intervalos.
• Encuentre los intervalos de valores y, para los que f (y)
< 0, y dibuje flechas que señalen hacia abajo (o hacia la
izquierda) en esos intervalos.
162
Ecuaciones Diferenciales Autónomas.
y/ = f(y) = (4)
163
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
f (y) < 0 en ( - ,
f( y ) > 0 en (2, ).
164
Ecuaciones Diferenciales Autónomas.
(5)
Solución:
La línea de fase para esta ecuación diferencial está dada en la
figura 10. Considérese que los puntos de equilibrio
•y=2
•y=0
•y=- ,
w
w=2
w=0 t
–2 2 4
166
Ecuaciones Diferenciales Autónomas.
w
w=2
jo corresponden a las soluciones decrecientes. Las gráficas de las
soluciones no se cortan de acuerdo con el teorema de existencia y
unicidad. En particular, no pueden cortar las gráficas de las solucio-
nes de equilibrio. Además, las soluciones deben continuar creciendo
w=0
o decreciendo hasta que se acercan a una solución de equilibrio.tPor
–2 2 4
consiguiente, podemos esbozar muchas soluciones con condiciones
iniciales diferentes fácilmente. La única información que no se tiene
es qué tan rápido crecen o decrecen.
w
w = 2π
w=π
w=2
w=0 t
–2 2
w = –π
Fuente. Elaborada por el autor con el programa Derive.
167
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Entonces:
(6)
Solución:
Esta ecuación tiene su punto de equilibrio en y = - 1 y dy / dt >
0 para todo otro valor de y.
168
Ecuaciones Diferenciales Autónomas.
t
–4 –2 2 4
y = –1
Si calculamos de forma
y explícita la solución, determinamos que
esas soluciones son de la forma:
y (t) = - 1 +
para alguna
y = 1 constante C. Como se supone que y (0) > -1,
debemos tener:
t
y (0) = -1-
169
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 7:
(7)
t
Solución: –4 –2 2 4
y = –1
Si y > 1, dy/dt es negativa; y si y < 1, dy/dt es positiva. Si y = 1,
dy /dt no existe. La línea de fase tiene un agujero. En y = 1. No hay
una manera estándar para denotar tal punto sobre la línea de fase,
pero usaremos un pequeño círculo vacío para marcarlos (figura 14).
y=1
t
170
Ecuaciones Diferenciales Autónomas.
(8)
Solución:
f (y)
y=c
y y=b
a b c
y=a
171
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y=c
y=b
t
y=a
172
Ecuaciones Diferenciales Autónomas.
Solución:
Esta ecuación diferencial (9) es autónoma y satisface las condi-
ciones del teorema de existencia y unicidad. De manera que tenemos
soluciones únicas que nos permite analizar su línea de fase. Primero
tracemos la gráfica de f (y) = (ver figura 16), pero,
¿qué dominio debería considerarse? Después de observar f(y) vemos
que f (0) = 0, lo que nos da una solución de equilibrio, y (t) = 0. ¿Exis-
ten otra? Si pensamos en f (t) para valores grandes de y, vemos que
el término –5e–y sería despreciable cuando se compara con los otros
dos términos de la expresión general, así, para y muy grande, f (y)
, lo cual es cero cuando y = 10, de manera que podemos
anticiparnos y decir que otra solución de equilibrio aproximadamente
es f (t) = 10. Haciendo uso de un método. Para encontrar raíces y re-
solver f (y) = 0 cerca de 10, encontramos que y ; así que
la segunda solución de equilibrio se halla en f (t) 9.999546.
Además
(10)
173
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y=0 y = 9.999546
174
Ecuaciones Diferenciales Autónomas.
Figura 19. Campo de pendientes con una solución numérica de la ecuación difer-
encial autónoma.
175
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
176
Ecuaciones Diferenciales Autónomas.
El teorema de Linealización
Suponga que y0 es un punto de equilibrio de la ecua-
ción diferencial dy / dt = f (y) donde f es una función
diferenciable continuamente. Entonces:
Si f/(y0) < 0, entonces y0 es un punto estable (o sumidero).
Si f/(y0) > 0, entonces y0 es un punto inestable (o fuente).
Si f/(y0) = 0, o si f/(y0) no existe, entonces y0 puede ser esta-
ble, inestable o semiestable (o nodo).
f (y)
177
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
f (y)
Ejemplo 10:
Solución:
Tenemos que f (y) = y (a – y) = ay – y2 y las soluciones de equilibrio
son y (t) = 0 y y (t)= a. Como f/(y) = a – 2y; encontramos f/ (0) = a, y
f/ (a) = - a. Así, si a < 0 la solución de equilibrio y (t) = 0, mientras que
la solución de equilibrio f (t) = a es inestable. Si a > 0, la solución de
equilibrio f (t) = 0 es inestable, mientras que la solución de equilibrio
f (t) = a es estable. Si a = o, entonces f/ (0) = 0, y la prueba no nos
proporciona conclusión alguna (de hecho, en este caso la solución de
equilibrio es semiestable).
178
Ecuaciones Diferenciales Autónomas.
Teorema
Las soluciones de no equilibrio de y/ = f (y) son monótonas si
f (y) es continua.
179
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
El modelo
180
Ecuaciones Diferenciales Autónomas.
181
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
«algo» =
182
Ecuaciones Diferenciales Autónomas.
Figura 23. Gráfica del comportamiento de las soluciones alrededor de los puntos
de equilibrio.
183
3.2 Ejercicios
En los siguientes ejercicios esboce las líneas de fase para la ecuación
diferencial dada. Identifique los puntos de equilibrio como estables
(sumidero), inestables (fuente), semiestables (nodo).
1. 2.
3. 4.
5. 6.
7. 8.
9. 10.
11. 12.
Figura 24. Gráficas de las funciones f(y) para las ecuaciones diferenciales autónomas.
27 28
29 30
31 32 33 34
, y (0) = ½
y/ = sen con y (0) = y0, con y (0) = y0, donde 0 < y0 <
2? ¿Qué pasa si y0 > 2?
y/ = (y – a)(y – b)(y – c)
a. y/ = a – y2 b. y/ = y(a – y2)
a. y/ = ey b. y/ = y
c. y/ = 1 – y2 d. y/ = y (1 – y2)
e. y/ = 1 – y4 f. y/ = y ey
g. y/ = 1 + y2 h. y/ = (1 – y)(2 – y)(3 – y)
i. y/ = cosh y j. y/ = y3
k. y/ = -y lny l. y/ = sen y
47. 48.
49. 50.
51. 52.
53. 54.
55. 56.
En los ejercicios 57 y 58, determine la estabilidad usando la
línea de fase unidimensional.
57.
58.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
3.3 Diagramas de
bifurcación
Figura 26. Stephen Smale (1930).
194
Ecuaciones Diferenciales Autónomas.
195
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
fc (y) = y2 - 2y + c
f2 (y) = y2 - 2y + 2
196
Ecuaciones Diferenciales Autónomas.
y/ = fa(y) = y(a – y)
Figura 27. Gráficas del comportamiento de curvas solución para fc a medida que
varía c.
y y y
y=2
y=0
x x x
y = -2
a = -2 a=0 a=2
y=0
x x x
y = -2
a = -2 a=0 a=2
197
y
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Se observa que a medida que crece a, de negativo a positivo
pasando por cero, la solución de equilibrio y (x) = a se eleva
gradualmente. Nótese que si a < 0, la solución de equilibrio y (x) =
0 es estable (sumidero), mientras que la solución de equilibrio y (x)
= a es inestable (fuente). Si a = 0, existe solamente una solución
de equilibrio, a saber y (x) = 0, y es semiestable (nodo). Si a > 0 la
solución de equilibrio y (x) = 0 es inestable (fuente), mientras que la
solución de equilibrio
y y (x) = a es estable
y (sumidero).y De este modo, el
comportamiento de la solución cambia para a = 0. y=2
Si en la figura 27
x agregamos flechas a lo largo de xlos ejes
y=0
x
verticales
y = -2 para producir una línea de fase para cada uno de estos tres
casos, obtenemos la figura 28.
a = -2 a=0 a=2
Figura 28. Líneas de fase para los tres casos anteriores.
y y y
y=2
y=0
x x x
y = -2
a = -2 a=0 a=2
198 estable
y=0
x x x
Ecuaciones Diferenciales Autónomas.
y = -2
a = -2 a=0 a=2
y
y2 = a
estable y=0
inestable a
estable
199
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y y
4 y=a y=a
estable
–2
inestable
–4
200
μ < 1 μ = 1 μ > 1
y=0
x x x
Ecuaciones Diferenciales Autónomas.
y = -2
a = -2 a=0 a=2
y
y2 = a
estable y=0
inestable a
estable
Solución:
Como no conocemos aún los valores más interesantes de c, se
escoge para empezar valores enteros, digamos = - 4, = -2, = 0,
=2y = 4. Esto nos da cinco ecuaciones diferenciales autónomas
diferentes con cinco líneas de fase diferentes. Una de estas es:
201
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
estable
Los puntos
2 de equilibrio son y = 1 - y, y =1 + . Entre ambos
valores la función f-2 es negativa, por arriba y debajo de ellos f-2 es
positiva. Por lo tanto y =y =10 - es un punto
estableestable (sumidero)y =y,0 y =
1
–4 + –2es un punto
0 inestable
2 4(fuente).
a Con esta información inestable
podemosa
graficar la línea de fase. Para los otros valores de procedemos de
–2
manera similar y dibujamos las líneasinestable
de fase como se muestran en
la figura 32.–4
= -4 = -2 =0 =2 =4
Cada una de las líneas de fase es algo diferente de las otras. Sin
embargo, la descripción básica para = -4, = -2 y = 0 es la misma:
hay exactamente dos puntos de equilibrio; el menor es un punto
estable (sumidero) y el mayor es un punto inestable (fuente). Aunque
la posición exacta de estos
μ < 1 puntos
μ =cambia
1 cuando
μ > 1 t se incrementa, su
posición relativa y tipo no cambia. Las soluciones de esas ecuaciones
con grandes valores iniciales explotan en un tiempo finito, conforme
t aumenta y tienden a un punto de equilibrio cuando t decrece. Las
202
Ecuaciones Diferenciales Autónomas.
203
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y2 – 2y + =0
204
–2 Ecuaciones Diferenciales Autónomas.
inestable
–4
μ < 1 μ = 1 μ > 1
205
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 13.
Solución:
En esta ecuación a es el parámetro. Si hacemos fa (y) = 0, tenemos
. Por lo tanto, se tienen tres puntos de equilibrio si a > 0,
y = 0, ; pero hay un solo punto de equilibrio si a 0, es y = 0. Por
lo tanto, ocurre una bifurcación cuando a = 0. Para entender mejor lo
que sucede, trazamos un diagrama.
206
Ecuaciones Diferenciales Autónomas.
207
3.4 Ejercicios
c. y/ = a +y + y2 d. y/ = a – y -y2
e. y/ = a + y2 f. y/ = y2 + 3y + a
g. y/ = y2 – ay + 1 h. y/ = cos y + a
a. y/ = 1 + cx + x2 b. y/ = x – cx (1 – x)
c. y/ = x2 – 2x + c d. y/ = x + c x3
8. En los ejercicios a al d localice los valores de bifurcación
para la familia paramétrica y dibuje las líneas de fase para
los valores del parámetro ligeramente menores que, o
ligeramente mayores que , y el valor de bifurcación.
a. α b. α
c. d. α
CAPÍTULO 4.
SOLUCIONES
NUMÉRICAS Y
CUANTITATIVAS
DE UNA ECUACIÓN
DIFERENCIAL
Contenido
4.1 Soluciones numéricas: método de Euler 215
4.2 Pasos a lo largo del campo de pendientes 216
4.3 Método de Euler 218
4.4 Ejercicios 228
4.5 Un método cuantitativo (analítico)
Separación de variables 232
4.6 Problemas de aplicación 243
4.7 Ejercicios 258
Competencias
1. Aplica el método de Euler para
resolver problemas con valor inicial.
Introducción
Ya se ha dicho que hay básicamente tres formas de resolver una
ecuación diferencial, como son: la forma numérica, la forma cualitativa
y la forma cuantitativa; en esta sección queremos mostrar cómo se
resuelve numéricamente un problema con valor inicial mediante el
método numérico de Euler; para luego trabajar el método cuantitativo
más importante para resolver ecuaciones diferenciales lineales,
llamado el de variables separables.
214
Soluciones numéricas ycuantitativas de una ecuación diferencial
215
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
216
Soluciones numéricas ycuantitativas de una ecuación diferencial
217
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y k+1 = tk +
(tk+1, yk+1)
(tk, yk)
218
Soluciones numéricas ycuantitativas de una ecuación diferencial
O: yk+1 = yk +
yk+1
yk
Δt Δt
tk tk+1 tk+2 t
219
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
En resumen:
Ejemplo 1:
, y (0) = 1
Solución:
En este ejemplo, f (t, y) = 2y - 1, por lo que el método de Euler
está dado por:
Yk+1 = y k + (2yk- 1)
1
Si se desea consultar una fuente muy amena de métodos de aproximación para PVI, véase el libro del renomb-
rado analista numérico estadounidense Shampine. L. F. (1994). Numerical Solution of Ordinary Differential Equa-
tions. Nueva York: Charman & Hall.
220
Soluciones numéricas ycuantitativas de una ecuación diferencial
221
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
K Tk yk f(t k, y k)
0 0 1 1
1 0.1 1.1 1.2
2 0.2 1.22 1.44
3 0.3 1.364 1.73
4 0.4 1.537 2.07
5 0.5 1.744 2.49
6 0.6 1.993 2.98
7 0.7 2.292 3.58
8 0.8 2.650 4.3
9 0.9 3.08 5.16
10 1.0 3.596
222
Soluciones numéricas ycuantitativas de una ecuación diferencial
223
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ahora se tiene el punto (t2, y2) = (1.1, 1.105). Este tipo de cálculo
se vuelve tedioso muy rápidamente, pero afortunadamente este tipo
de operaciones son perfectas para un computador o una calculadora.
Para = 0.05, los resultados del método de Euler se muestran en la
tabla 2.
Tabla 2. Tabla para la solución del método de Euler con un paso más pequeño.
K T k
y k
f(t k, y k)
0 0 1 1
1 0.05 1.050 1.100
2 0.1 1.105 1.210
3 0.15 1.166 1.331
: : : :
: : : :
19 0.95 3.558 6.116
20 1.00 3.864
, y (0) = 1
Solución:
El valor de la solución en t = 2 es y (2) = 1/5 si se calcula la
integral y se evalúa en este valor. Nuevamente, es interesante ver qué
tanto nos acercamos a este valor con varias selecciones de . La
fórmula para el método de Euler es:
224
Soluciones numéricas ycuantitativas de una ecuación diferencial
K Tk yk
0 0 1
1 0.1 1.0000
2 0.2 0.9800
3 0.3 0.9416
: : :
19 1.9 0.2101
20 2.0 0.1933
225
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 3:
, y (0) = 1.
Solución:
Se tiene que f (t, y) = t y + y, tenemos también que t0 = 0, y0 =
y (0) = 1 de la ecuación de Euler, se tiene que yk+1 = yk + f (tk, yk) =
yk + (tk +tk yk) .
226
Soluciones numéricas ycuantitativas de una ecuación diferencial
227
4.4 Ejercicios
En los siguientes ejercicios del 1 al 8, use el método de Euler
con el tamaño de paso para aproximar la solución al problema de
valor inicial, en el intervalo de tiempo especificado. Su respuesta debe
incluir una tabla de los valores aproximados de la variable dependiente.
Trate de incluir un croquis de la gráfica de la solución aproximada.
1. , y (0) = 3,
2. , y (0) = 1,
3. , y (0) = 2,
4. y (0) = 1,
5. , y (0) = 4,
6. , y (0) = 0,
7. , y (0) = 2,
8. , y (1) = 2,
y (0) = 1
Usando el método de Euler, calcule tres soluciones
diferentes aproximadas correspondientes a
sobre el intervalo . Grafique
las tres soluciones. ¿Cuáles son sus predicciones acerca
de la solución real al problema de valor inicial?
a. y/ = 1 + y2 y (0) = 0
b. y/ = e-y y (0) = 0
c. y = 2t – 3y
/
y (1) = 1
, y (0) = 1.
Fuente. Mûlges y Uhlig (1996).
Carl Runge pasó sus primeros años en La Habana, donde su padre
Julius Runge ejercía como cónsul danés. La familia se trasladó más
adelante a Bremen, donde Julius murió prematuramente en 1864.
231
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Introducción
Un núcleo inestable es radioactivo. En cualquier instante puede emitir
partículas, transformándolo en un núcleo diferente en el proceso.
Por ejemplo 238U es una alfa emisor en decadencia exportando
de acuerdo al esquema 238
U 234
Th + 4
He, donde 4
He es la alfa
partícula. En una muestra de 238
U, un cierto porcentaje del núcleo
decae durante un periodo de observación . Si en un tiempo t la
muestra contiene N (t) núcleo radioactivo, entonces esperamos que la
cantidad del núcleo radioactivo en un tiempo es aproximadamente
proporcional a N (t) y . En forma simbólica:
) (1)
232
Soluciones numéricas ycuantitativas de una ecuación diferencial
N / (t) = (2)
(3)
o´
Ln (4)
Como e C
y son ambos positivos, entonces tenemos dos
casos:
N (t) =
2 Si usted no está conforme con todos estos pasos, puede verificar la solución por diferenciación y sustitución.
Se puede dar una justificación formal mediante la regla de la cadena.
233
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
N (t) = A (5)
Definición
Una ecuación diferenciable de la forma:
dy/dt = f (t, y)
se dice que es de variables separable si la función f (t,
y) puede escribirse como el producto de dos funcio-
nes: una que dependa solo de t y otra que dependa solo
de y. Es decir, una ecuación diferencial de variables se-
parables sí se puede escribir en la forma:
3 También es cierto que si es de la forma dy/dt=f(t)/h(t) se puede ver como variable separable, escribiéndola
como dy/dt =f(t). (1/h(t)).
234
Soluciones numéricas ycuantitativas de una ecuación diferencial
Ejemplo 4:
dy/dt = y.t
Solución:
Es claramente separable ya que podemos escribir g (t) = t, y h
(y) = y. por lo tanto:
Ejemplo 5:
dy/dt = y + t
Solución:
No es separable, ya que no hay forma de expresar y + t como
una función g (t) y otra h (y) tal que = g (t).h (y).
Ejemplo 6:
Solución:
En esta ecuación tenemos que trabajar un poco más para poder
evidenciar que sí es separable, así:
235
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
g(t) = y, h (y) =
Observación:
Dos importantes tipos de ecuaciones diferenciales se presentan
si t o y están ausentes en el lado derecho de la ecuación.
Ejemplo 7:
236
Soluciones numéricas ycuantitativas de una ecuación diferencial
Solución:
El lado derecho solamente depende de la variable dependiente
P (t), por lo que esta ecuación es otro ejemplo de una ecuación que
es autónoma y por lo tanto, se puede ver como una ecuación de
variables separables.
Ejemplo 8:
Solución:
Efectuando algo de algebra informal y reescribiendo esta
ecuación en la forma:
237
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
-2
-3
-4
238
Soluciones numéricas ycuantitativas de una ecuación diferencial
Ejemplo 9:
(7)
Solución:
Esta ecuación es separable, luego haciendo la separación de
variables se tiene:
Integrando se tiene:
Resolviendo para y:
(8)
y=
Ejemplo 10:
(9)
239
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
De nuevo esta ecuación es separable, luego haciendo la
separación de variables se tiene:
Integrando, se tiene:
-1 = A (y – 1) + B y (10)
Integrando:
O equivalentemente:
(11)
. Esto implica
240
Soluciones numéricas ycuantitativas de una ecuación diferencial
y = C1 .(y – 1).et
Despejando y, se tiene:
(12)
Como y (y – 1) = 0 si y = 0, o si y = 1, se sigue que y = 0 y, y = 1
son soluciones de (9). La solución y = 0 corresponde a k =0. Además,
y = 1 no está incluido (11) salvo si tomamos C1 = . Luego la solución
de (9) es (11) y, y = 1. La solución general involucra una constante C1,
la cual puede ser determinado de alguna condición inicial.
Ejemplo 11:
(13)
Solución:
De nuevo separando variables, colocando todas las expresiones
que involucran la función desconocida en el lado izquierdo y las que
involucran la variable dependiente en el lado derecho, resulta:
(14)
4 Como uno puede sospechar el método de variables separables se remonta casi a los comienzos del cálculo.
Gottfried Leibniz (1646-1716) usó esto implícitamente, cuando resolvía el problema inverso de la tangentes en
1691. Esto fue precisado por Johann Bernoulli en 1594 en una carta dirigida a Leibniz.
241
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Figura 10. Familia solución para la ecuación diferencial del ejercicio ejemplo 11.
4
-4 -3 -2 -1 1 2 3 4
-1
-2
-3
-4
242
Soluciones numéricas ycuantitativas de una ecuación diferencial
Ejemplo 12:
Solución:
a. Si A (t) es la cantidad de dinero en el tiempo t, entonces
dA/dt denotará el incremento del capital en cada instante
t, es claro que en incremento del capital es directamente
proporcional al capital en cada momento, luego tenemos la
ecuación diferencial:
243
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Dólares.
(16)
A (10) = 8244
244
Soluciones numéricas ycuantitativas de una ecuación diferencial
Es decir:
Entonces en t = 10:
Esto quiere decir que estamos retirando a una tasa que excede
a la tasa a la que ganamos interés. Como dA/dt en t = 10 es negativa,
A disminuirá y 0.05A – 1000 permanece negativa para t > 10. Si
0.05A– 1000 es negativa, entonces:
O:
245
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Vemos que:
A (11) 7641
A (12) 7008
246
Soluciones numéricas ycuantitativas de una ecuación diferencial
Solución:
Para elaborar el modelo, tenemos en cuenta que t es el tiempo medido
(la variable independiente). Para la variable dependiente, realmente
se tendrían dos opciones; una escoger la cantidad total de sal S (t)
en el tiempo t medida en cucharadas, o bien C (t), la concentración
de sal en el tanque en el tiempo t medida en cucharadas por galón.
Desarrollaremos el modelo para S (t), dejando el modelo para C (t)
como un ejercicio. Usando la sal total S (t) en el tanque como variable
dependiente, la razón de cambio de S (t) es la diferencia entre la
247
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
cantidad de sal que se añade y la cantidad de sal que sale del tanque.
La sal que entra al tanque llega por los tubos A y B y se puede calcular
fácilmente, multiplicando el número de galones por minuto de la
mezcla salmuera que entra al tanque por la cantidad de sal por galón.
La cantidad de sal que sale del tanque por el tubo C depende de la
concentración de sal en el tanque en ese momento. La concentración
está dada por S (t) / 100, por lo que la sal que sale del tanque es del
número de galones que salen por minuto (3 galones por minuto) y la
concentración (S/100). El modelo es:
Es decir:
248
Soluciones numéricas ycuantitativas de una ecuación diferencial
Encontramos que:
2000 – 3S = C4e-0.03t
S (t) =
249
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Lo primero que notamos es que se está agregando más líquido
por minuto que el que se está extrayendo, de modo que el nivel del
líquido asciende y el número de galones en el contenedor aumenta. De
hecho, podemos utilizar un argumento tipo conservación para advertir
que la velocidad de cambio del volumen V de líquido en el contenedor
(en galones por minuto) es igual a la velocidad de agregación (3
galones por minuto) menos la velocidad de drenado (2 galones por
minuto). De esta manera:
V (t) = t + C galones
250
Soluciones numéricas ycuantitativas de una ecuación diferencial
(18)
251
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Figura 12. Gráfica del campo de pendientes con una solución particular del prob-
lema 15.
252
Soluciones numéricas ycuantitativas de una ecuación diferencial
Integrando:
O:
ln y = - 2ln (t + 200) + C
y (t) =
y (t) = (19)
253
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 15:
(20)
Integrando:
O:
(21)
254
Soluciones numéricas ycuantitativas de una ecuación diferencial
(22)
(23)
255
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
=0
(24)
256
Soluciones numéricas ycuantitativas de una ecuación diferencial
257
4.7 Ejercicios
a. x2y/ - y2 = 0 b. y/ - sen(x/y) = 0
c. y y/ - exy = 0 d. y y/ -ex+y = 0
e. y/ - x + y = 0 f. y/ - y sen x = 0
g. y/ + 4 + y2 = 0 h. y/ - ln x + ln y = 0
i. y/ - ty = 0 j. y/ - t4 y2 = 0
k. y/ - 2y – 1 = 0 l. y/ + y -2 = 0
m. y/ -t2y – 1 - y - t2 = 0 n. y/ - 1/(ty + t +y + 1)
2. 3.
4. 5.
6. 7. .
8. 9.
10. 11.
16. 17.
18.
19. 20.
21. 22.
23. 24.
25.
a. x y y/ - y = 0 y (e) = 1
b. x y/ + y2conx = 0 y (0) =
c. x y y/ - (4 – y2)1/2 = 0 y (8) = 1
d. (1 – x2)1/2 y y/ + y3 = 0 y (1) = 1
h. 3ye y y/ + 1 = 0 y (1)= 2
j. y/ - 2yx/(x2-1) = 0, y (2) = -6
k. y/ - 2yx/(x2-1) = 0 y (0) = 2
l. y/ - 2yx/(x2-1) = 0 y (-2) = -6
28. 29.
30. 31.
32.
33. Una cubeta de 5 galones está llena de agua pura.
Suponga que empezamos a añadir sal a la cubeta a razón
de ¼ libra por minuto. Además, abrimos el grifo de manera
que salga ½ galón por minuto de la cubeta, y agregamos
agua pura para mantener llena la cubeta. Si la solución
de agua salada está siempre bien mezclada, ¿cuál es la
cantidad de sal en la cubeta después de?:
37. Una curva pasa por el punto (1, 1). En cada punto (x,
y) de esta curva, la recta tangente pasa por el punto (2x,
-y). Pruebe que la curva satisface la ecuación diferencial:
Y encuentre la ecuación de la curva.
38. Un esquiador localizado en el punto W, se sujeta a
mediante una cuerda de longitud L a un bote de motor
(M), como se muestra en la figura 15.
M x
a.
b.
c.
Donde k es una constante. Si T (t) > Tm, entonces el
objeto se está enfriando y, en consecuencia, k debe ser
negativa. La misma conclusión se cumple para k si T(t) <
Tm, de manera que en cualquier caso se hace K negativa.
Utilizando variables separables demuestre que:
T (t) = Tm +C e kt
MÁS SOBRE
MÉTODOS
CUANTITATIVOS
Contenido
Competencias
1. Identifica ecuaciones que
se pueden llevar a variables
separables.
268
Más sobre métodos cuantitativos
Introducción
Volvemos en este capítulo a trabajar de manera más amplia el método
cuantitativo o analítico, para resolver ecuaciones diferenciales. El
método de separación de variables se aplica a ciertas ecuaciones que
cumplen como se vio anteriormente, la condición de poder separar las
variables en y, y en x. Pero este método es realmente muy importante
ya que una gran cantidad de ecuaciones mediante cambios de
variables se pueden llevar a variables separables como veremos en
este capítulo.
269
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
1 Polyanin, A. y Zaitsev, V. (1995). Exact Solutions for Ordinary Differential Equations. Boca Ratón, Florida: CRC
Press.
2 Este uso del término Euler – homogénea es diferente del término lineal – homogénea que se usará más adelante.
270
Más sobre métodos cuantitativos
Ejemplo 1:
(1)
Solución:
Empleando la sustitución u = x + y, donde derivando respecto a x se
tiene , es decir que , luego la ecuación (1) se nos
transforma en:
Integrando:
Se tiene: Ln (u + 1) = x + C
Ejemplo 2:
271
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Lo primero que se hace es separar la fracción del lado derecho, es
decir:
(2)
Cancelado, u se tiene:
Separando variables:
Integrando, tenemos:
u = ln x + C
Sustituyendo u:
272
Más sobre métodos cuantitativos
u = x + y (3) u = x – y (4)
u = 2x + 3y – 1. (5)
F (u) =
F (u) = - u2
F (u) = ln u.
) (6)
F (u) =u Ln u
273
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
O bien:
(8)
d u = a dx + b dy
O bien:
x+C= (9)
274
Más sobre métodos cuantitativos
Solución:
Para la ecuación (3) se tiene que F (u) = , por tanto, utilizando
(8), tenemos:
x+C=
d y = x du + u dx
Y finalmente:
(9)
275
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Teorema 1
Si una ecuación tiene la forma:
Ejemplo 4:
Solución:
Ln x + C =
Ln x + C =
z - 1 = C x, luego
ln u – 1 = C x, es decir ln u =C x + 1
Ejemplo 5:
Resolver la ecuación .
Solución:
Haciendo u = , entonces se tiene:
Ln x + C =
276
Más sobre métodos cuantitativos
Ln x + C =
De donde:
Ln x +C = ln (u +
Por lo tanto:
(11)
O:
1 + u2 = (C1x –u)2, entonces 1 + u2 =
Luego:
1= ,ó1=
Por lo tanto:
y=
F(x, y) dy = G(x, y) dx
277
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
F (x, y) dy = G (x, y) dx
H (u) dy = K (u) dx
Ln x + C = (12)
278
Más sobre métodos cuantitativos
Ejemplo 6:
Resolver la ecuación:
Solución:
Observe que F (x, y) = x y – 2x2, y g (x, y) = 5x2 + 2y2 son polinomios
de grado dos. Lugo dividiendo por x2, y haciendo la sustitución u =
y/x, se tiene:
(u -2) dy = (5 + 2u2) dx
O bien en la forma:
Ln x + C =
Ln x + C =
Ejemplo 7:
279
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
En este caso no podemos usar la sustitución u = y/x, ya que
tenemos una constante en el denominador. Sin embargo, si se hace la
sustitución t = x+1, y, u = y - 1, entonces la ecuación toma la forma:
De donde:
+C =
v = ln t + c
O bien:
u = t ln t + C t
Reemplazando u, y v tenemos:
y – 1 = (x + 1) ln (x+1) + C(x+1).
Ejemplo 8:
280
Más sobre métodos cuantitativos
Solución:
Con el fin de encontrar las relaciones entre las variables x, y y
con las nuevas variables t y u, se hace:
x = t + A, y y = u + B
A+B+1=0
2A–B–4=0
x = t + 1, y, y = u – 2
281
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Se convierte en:
Donde v = u/t.
En resumen
El cambio de variable:
x = t + A, y = u + B
Se puede emplear para cualquier ecuación diferencial
de la forma:
(13)
Siempre y cuando se pueda resolver para A, y B que
satisfacen
a A + b B + c = 0
dA+ eB + f = 0.
Ejemplo 9:
Solución:
Al formar el sistema:
A + 2B – 1 = 0
2A + 4B – 3 = 0
282
Más sobre métodos cuantitativos
Este sistema no tiene solución, luego las dos rectas son paralelas.
En este caso la solución se puede encontrar haciendo v = x + 2y – 1,
ya que transforma la ecuación diferencial en:
Ejemplo 10:
Solución:
Esta ecuación a pesar de que f(x, y) no es lineal se puede
resolver por separación de variables haciendo el cambio lineal u = x
+ y, luego du = dx + dy, esto es dy = du – dx, entonces la ecuación se
transforma en:
283
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Es decir:
,o
Ejemplo 11:
Resolver la ecuación:
Solución:
La sustitución racional u = y/x, en este caso nos permite resolver
esta ecuación, ya que y = u x, luego dy = u dx + x du, reemplazando
en la ecuación tenemos:
,ó
, integrando:
Arctan (u) = ,ó
y = x tan
284
Más sobre métodos cuantitativos
En resumen:
Se hace la sustitución racional, siempre y cuando el miembro de
la derecha de la ecuación diferencial se pueda escribir en la forma:
u + G(x).H (u)
(14)
285
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 12:
Resolver la ecuación:
Solución:
En este caso, el término constante 2 sugiere la sustitución lineal.
Si se sustituye u = - 2x + y, entonces –a/b = 2. Entonces du = -2dx +
dy, de donde dy = du +2dx. La ecuación se puede escribir en la forma:
, o separando variables
, integrando tenemos:
Sustituyendo u = -2x-+y
286
5.2 Ejercicios
Resuelva los siguientes ejercicios utilizando una sustitución lineal.
1. 2.
3. 4.
5. 6.
7. 8.
9. 10.
11. 12.
13. 14.
17. 18.
19. 20.
21. 22.
23.
24.
25.
26.
27. (x + y) dx + (x – y) dy = 0.
29.
30.
31.
32.
33.
34.
35.
36. y/ =
37.
292
Más sobre métodos cuantitativos
Ejemplo 13:
Resolver la ecuación:
Solución:
Haciendo z = x2, se tiene que d z = 2x dx. Dividiendo ambos
miembros de la ecuación por 2x, se puede reemplazar x por z como
variable independiente, es decir:
293
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
De donde:
Ln u + z2 + z+C
u = e (z +z+C), o u =
2
Sustituyendo se tiene:
z + 2y + 1 =
Ejemplo 14:
Resolver la ecuación
Solución:
En este caso dividimos por el factor y, y obtenemos:
d u = y-1 dy
, separando variables:
,o
, integrando:
294
Más sobre métodos cuantitativos
Es decir que:
Ln (z + 1) = x + C
O:
z+1=
z= – 1,
x+u=
Y u por su igual:
x + ln y = , o ln y =C1 ex–1–z
y=
Ejemplo 15:
Solución:
Dividiendo por sen x, tenemos:
Haciendo t = - cos x
se tiene d t = sen x dx que produce la ecuación:
295
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
, integrando
y – ln y = t + C sustituyendo se tiene:
y – ln y = -con x + C
La sustitución u = xm yn
Otro método que resulta eficaz con frecuencia para convertir
ecuaciones no separables en separables, involucra una sustitución de
la forma:
u = xm yn
296
Más sobre métodos cuantitativos
Teorema 2
La sustitución u = xm yn transforma una ecuación de la forma:
Donde F (u) puede ser una función cualquiera que depende
de u, en una ecuación separable en términos de x y u.
Ejemplo 16:
Resolver la ecuación:
Solución:
Si factorizamos y/x en la ecuación, se tiene:
d y = x-2du – 2ux-3dx,
297
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Integrando:
Ln u = m ln x + n ln y, derivando tenemos:
,o
, separando variables:
, integrando:
298
Más sobre métodos cuantitativos
Problema 17:
Resolver la ecuación:
Solución:
Si factorizamos y/x en la expresión, tenemos:
Ln x =
Ln x =
Es decir:
Ln x + C = , sustituyendo se tiene:
Ln x + C =
o, Ln x + C =
299
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Problema 18:
Solución:
Reescribiendo la ecuación en la forma:
300
5.4 Ejercicios
En los siguientes ejercicios encuentre una sustitución que convierta
cada ecuación en una que se pueda resolver mediante una sustitución
lineal o racional. Reemplace x o y por las nuevas variables en la
ecuación, y en cada caso indique la relación entre la variable x y la
variable u:
1. 2.
3. 4.
5. 6.
7. 8.
9. 10.
11. 12.
13. 14.
15. 16.
17. 18.
19. 20.
21. 22.
23.
Más sobre métodos cuantitativos
303
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Consideremos la ecuación:
y dx + x dy = 0 (15)
du=0
De lo anterior:
u=C
xy=C
Ejemplo 19:
(y +2) dx + x dy = 0 (16)
(y+2) dx + (x + y2) dy = 0 (17)
(yex+5x4lny)dx + (x5y-1 + ex)dy = 0 (18)
304
Más sobre métodos cuantitativos
Solución:
La única diferencia de la ecuación (15) y (16) es el término extra 2.dx,
por lo tanto, no resulta difícil ver cómo se debe modificar la situación
que se empleó en (15) para resolver (16): si se hace u = x y + 2x,
entonces el miembro izquierdo de (16) es igual du, en consecuencia,
se tiene de nuevo:
d u = 0, o u = C
x y +2x = C
u = x y + 2x+
x y + 2x+ =C
3x.y + 6x+ = 3C
u = y ex +x5ln x
y ex +x5ln x = C
305
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
M(x, y) dx + N(x, y) dy = 0
F(x, y) = C
w = M(x, y) dx + N(x, y) dy
306
Más sobre métodos cuantitativos
Ejemplo 20:
W = x dx + y dy = 0 o,
x2 + y2 = C
Solución:
Para verificar esto basta derivar implícitamente esta ecuación
respecto a x, así:
2x + 2y , o, x dx + y dy = 0
307
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
3
2
-8 -6 -4 -2 2 4 6 8
-1
-2
-3
-4
Definición
Si las soluciones para la ecuación diferencial M(x,y) dx +
N(x,y) dy = 0 estan dadas Implícitamente por la ecuación:
F(x,y) = C
F(x, y) = C (21)
308
Más sobre métodos cuantitativos
(23)
(24)
(25)
(26)
w = M(x, y) dx + N(x, y) dy =
309
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Definición
La diferencial de una función continuamente diferenciable F
es la forma diferencial:
dF =
Una forma diferencial se dice que es exacta si ella es la dife-
rencial de una función continuamente diferenciable.
dF = = M(x, y) dx + N(x, y) dy
(27)
310
Más sobre métodos cuantitativos
Ejemplo 21:
Solución:
Usando las ecuaciones (27), podemos encontrar la función F(x,
y) que satisfacen:
F(x, y) = x2 + y4
x2 + y4 = C
-4 -3 -2 -1 1 2 3 4
-1
-2
-3
-4
311
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
M(x) dx + N (y) dy = 0
F(x, y) =
Teorema
Sea w = M(x, y) dx + N(x, y) dy una forma diferencial donde
M y N son continuamente diferenciables.
a. Si w es exacta, entonces:
(28)
b. Si la ecuación (28) es verdadera en una región rectangular
R, entonces w es exacta en R.
Demostración
Para probar la parte a, supongamos que w = d F, entonces:
(29)
312
Más sobre métodos cuantitativos
F(x, y) = (30)
= N(x, y) - = N(x, y) -
313
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Resolver , integrando:
F(x, y) = (31)
314
Más sobre métodos cuantitativos
Ejemplo 22:
Solución:
En este caso tenemos que M(x, y) = sen (x + y), y, N(x, y) =
2y + sen (x + y), y:
y, (32)
F(x, y) =
Luego:
F(x, y) = -cos (x + y) +
= 2y + sen (x +y)
315
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Entonces:
-8 -6 -4 -2 2 4 6 8
-2
-4
-6
-8
F(x, y) =
316
Más sobre métodos cuantitativos
Ejemplo 23:
(3y + dx + (3x+cosy) d y = 0
Solución:
Como M(x, y) = 3y +ex y N(x, y) = 3x + cos y:
y,
F(x, y) = 3yx + ex +
Luego:
Es decir que:
317
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
De donde:
-8 -6 -4 -2 2 4 6 8
-2
-4
-6
-8
318
Más sobre métodos cuantitativos
Definición
Un factor integrante para la ecuación diferencial
, es una función tal que la forma
es exacta.
Ejemplo 24:
Solución:
Como:
- =0
F(x, y) =
, tenemos:
319
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
F(x, y) =
x + y2 + C x2 = 0
-4 -3 -2 -1 1 2 3 4
-1
-2
-3
-4
320
Más sobre métodos cuantitativos
(33)
M(x, y) dx + N(x, y) dy = 0
321
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(34)
(35)
(36)
Ejemplo 25:
Resolver la ecuación:
Solución:
En este caso tenemos que ,y , es decir que no
es exacta. Pero de acuerdo a la ecuación (34), tenemos que:
322
Más sobre métodos cuantitativos
=0
Es decir que:
F(x, y) =
Luego: x–y=x+ ,o
F(x, y) =
g (y) =
323
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 26:
2xy3 dx + (3 x2 y2 + x2 y3 + 1) d y = 0 (37)
Y resuelva la ecuación.
Solución:
Tenemos que, M(x, y) = 2xy3 y, N(x,y) = 3x2y2 + x2 y3 + 1, de donde se tiene que:
Luego:
h=
324
Más sobre métodos cuantitativos
g=
2xy3ey dx + (3 x2 y2 + x2 y3 + 1) ey d y = 0
F(x, y) =
,
de donde se tiene que:
, y por tanto
F(x, y) =
325
5.6 Ejercicios
Calcule la diferencial d F para la función dada F:
3. F(x, y) = 4. F(x, y) =
9. (2x + y) dx + (x – 6y) dy = 0
11.
12.
13.
14.
16. =0
17.
18.
19.
21. (2r + ln y) dr + r y dy = 0.
22. (y2 – x y) dx + x2 dy = 0,
24. 3(y + 1) dx – 2x dy = 0,
25. (x2 + y2 – x) dx – y dy = 0,
ECUACIONES
DIFERENCIALES
LINEALES DE PRIMER
ORDEN Y MODELOS
Contenido
6.1 Resolución de ecuaciones
diferenciales lineales 332
6.2 Solución general de la
ecuación diferencial lineal 334
6.3 Solución de la ecuación diferencial
homogénea de primer orden 337
6.4 Factores integrantes para la ecuación
diferencial lineal de primer orden 341
6.5 Ejercicios 353
6.6 Algunas ecuaciones especiales
reducibles a lineales 359
6.7 Ejercicios 375
Competencias
1. Identificar la forma estándar de
una ecuación diferencial lineal de
primer orden.
2. Reconocer que en toda ecuación
diferencial lineal de primer orden,
su solución general es la suma de
una solución particular y la solución
de la homogénea.
3. Deducir la solución general de una
ecuación diferencial lineal.
4. Aplicar la solución general de la
ecuación diferencial lineal.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
330
Ecuaciones diferenciales lineales de primer orden y modelos
Introducción
La apariencia sencilla de las ecuaciones diferenciales de primer
orden puede confundir a algunas personas y hacerles creer que son
fáciles de resolver. A veces lo son, tal como se vio en el capítulo 4 y
5, pero con frecuencia, resolverlas no es un desafío menor. No existe
un método general para resolver todas las ecuaciones diferenciales
de primer orden. Los métodos existentes de resolución son aplicables
a ciertos tipos de ecuaciones diferenciales y, por tanto, es necesario
clasificarlas y estudiarlas en grupos separados. En este capítulo
aprenderemos a reconocer las ecuaciones diferenciales lineales de
primer orden, ya que estas siempre se pueden resolver usando un
enfoque sistemático, es decir, mediante una fórmula general, tal como
veremos.
331
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(1)
(2)
(3)
332
Ecuaciones diferenciales lineales de primer orden y modelos
333
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(4)
(5)
334
Ecuaciones diferenciales lineales de primer orden y modelos
(9)
335
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
336
Ecuaciones diferenciales lineales de primer orden y modelos
(10)
337
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Es decir:
ln
(11)
y=
y(x) = C (12)
y1(x) = (13)
(14)
338
Ecuaciones diferenciales lineales de primer orden y modelos
Ejemplo 1:
Solución:
Como la ecuación es homogénea, por separación de variables tenemos:
Integrando:
y (x) = C ,
Ejemplo 2:
Solución:
Separando variables, obtenemos:
339
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Integrando:
ln +
340
Ecuaciones diferenciales lineales de primer orden y modelos
(15)
2 Los dos grandes matemáticos del siglo XVIII, Leonard Euler (1707-1783) y Joseph Louis
Lagrange (1736-1818) hicieron la mejor contribución para la solución , Euler
resolvió la ecuación homogénea, introduciendo la expresión c para p una constante, y Lagrange
resolvió la ecuación no homogénea usando el método de variación de parámetros.
3 El método de variación de parámetros es muy útil ya que esa misma idea se utiliza para
ecuaciones lineales de grado superior.
341
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Asociando, tenemos:
(18)0,2
O en su forma separable:
(19)
(20)
(21)
342
Ecuaciones diferenciales lineales de primer orden y modelos
u (x) = (22)
y(x) = ,ó (24)
y(x) =
Yh(x) = , y
(25)
yp(x) =
343
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 3:
Para x > 0.
Solución:
Para empezar, debemos llevar la ecuación diferencial a su forma
estándar, es decir debemos dividir por x la ecuación dada:
(26)
u(x) =
(27)
344
Ecuaciones diferenciales lineales de primer orden y modelos
Despejando y, tenemos:
y(x) =
Ejemplo 4:
Para x > 0.
Solución:
Para empezar, debemos llevar la ecuación diferencial a su forma
estándar, es decir debemos dividir por x la ecuación dada:
(28)
u(x) =
(29)
345
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Despejando y, tenemos:
y (x) =
Ejemplo 5:
Para x > 0.
Solución:
Para empezar, debemos llevar la ecuación diferencial a su forma
estándar, es decir debemos dividir por x la ecuación dada:
(30)
u(x) =
(31)
346
Ecuaciones diferenciales lineales de primer orden y modelos
Despejando y, tenemos:
y (x) = (32)
Figura 2. Gráfica de las curvas solución de 32.
+ y 4
x
-8 -6 -4 -2 2 4 6 8
-1
-2
-3
-4
347
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La gráfica de la función Q(x) es la siguiente:
348
Ecuaciones diferenciales lineales de primer orden y modelos
y
1
0 1 x
U(x) =
,o
Integrando:
Despejando y, se tiene:
y=1+
349
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
O:
y(x) =
Y es continua en
4 Recuerde que los símbolos (x0+) y (respectivamente y (x0-) denotan el límite de y(x) cuando
x se aproxima a x0 con valores de x mayores o menores que x0.
350
Ecuaciones diferenciales lineales de primer orden y modelos
-1.5
Solución:
La ecuación ya se encuentra en la forma estándar, por lo tanto, se
tiene que:
, de donde se obtiene y =
351
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y=y= o, y(x) =
352
6.5 Ejercicios
En los ejercicios 1-10, encuentre la solución general mediante
la separación de variables. Problemas donde se necesita una
integral definida se indican con asterisco.
1. 2.
3. 4.
5. 6. 2x
7. 8.
9. 10.
17. 18.
19. 20.
21. 22.
23. 24.
25. 26.
27. 28.
29. 30.
31. 32.
33. 34.
35. 36.
37. 38.
41. 42.
43. 44.
45. 46.
47. 48.
49. 50.
51. 52.
53. 54.
55. 56.
57. 58.
59. 60.
61.
62.
63.
64.
65.
66.
67.
68. 4
72. 73.
74. 75.
76. 77.
78.
79. P(x) =
82. P(x) =
83. P(x) =
84. P(x) =
Ecuaciones diferenciales lineales de primer orden y modelos
Donde n es un entero
Se llama una ecuación de Bernoulli5.
5 Los Bernoulli fueron una gran familia de famosos matemáticos suizos durante cuatro
generaciones. Esta ecuación es nombrada por los dos hermanos Jakob (1654-1705) y Johann (1667-
1748) quienes trabajaron desde 1695 hasta 1697 sobre la solución de esta ecuación en competencia
con Leibniz.
359
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Esta es una ecuación de Bernoulli con P(x) = -1, Q(x) =x, y n = 2. Al
dividir por y2, la ecuación se obtiene:
360
Ecuaciones diferenciales lineales de primer orden y modelos
(33)
Solución:
Esta ecuación se puede escribir en la forma:
(34)
,ó
Luego si hacemos
se tiene:
z=
361
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
De donde se tiene:
(35)
362
Ecuaciones diferenciales lineales de primer orden y modelos
Ejemplo 10:
Resolver la ecuación:
(36)
Solución:
La ecuación se puede escribir, después de dividir por en la forma:
(37)
363
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
u=
O:
Luego:
De donde:
u=
Reemplazando u, tenemos:
y=
-8 -6 -4 -2 2 4 6 8
-2
-4
-6
-8
(38)
(39)
365
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
a. Para expresar el peso del pez en función del tiempo, entonces
tenemos de la ecuación (39), despejando L que:
L (t) =
366
Ecuaciones diferenciales lineales de primer orden y modelos
367
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(41)
-3 +0.6u = 1, ó
Donde:
(43)
368
Ecuaciones diferenciales lineales de primer orden y modelos
369
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
370
Ecuaciones diferenciales lineales de primer orden y modelos
u=
y=
371
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(38)
Solución:
Se hace el cambio de variable u = , luego
z = -2
-2
372
Ecuaciones diferenciales lineales de primer orden y modelos
y= y=
,
Ejemplo 13:
(39)
Solución:
Como es una ecuación de Riccati, podemos hacer la sustitución:
Simplificando obtenemos:
Es decir:
z = -x + C
Sustituyendo z, se tiene:
373
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y=x+
374
6.7 Ejercicios
Las siguientes son ecuaciones de Bernoulli, resolverlas haciendo
el cambio de variables necesario.
1. 2.
3. 4.
5. 6.
7. 8.
9. 10.
11. 12.
13. 14.
15. 16.
17. 18.
19.
20. , dado que y = 2x es una solución.
22. , es una
solución.
Y = 2 arctan u
Convierte la ecuación:
Sen y= ,y
Cos y =
APLICACIONES
DE LAS
ECUACIONES DE
PRIMER ORDEN
Contenido
Competencias
1. Identifica problemas de la vida real
que se resuelven mediante ecua-
ciones lineales de primer grado.
2. Aplica la ecuación diferencial lineal
a modelos de crecimiento y decre-
cimiento.
3. Aplica la ecuación diferencial a
modelos de mezclas.
4. Aplica la ecuación diferencial a
modelos de la física.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Introducción
Algún matemático, mencionaba que de alguna manera la matemática
es el arte de resolver problemas de la vida real, es por esto que en
este capítulo se plantean una serie de modelos básicos que se pue-
den resolver mediante ecuaciones diferenciales lineales de primer
orden. Se han tratado de mostrar los modelos más importantes que
cualquier ingeniero, físico, químico, biólogo, economistas o matemá-
tico debe conocer.
380
Aplicaciones de las ecuaciones de primer orden
Solución:
De acuerdo a nuestra hipótesis del modelo, tenemos que si P (t) es
el tamaño de la población en el tiempo t, entonces tenemos la ley del
crecimiento población es aplicable, por lo tanto:
7 Para un químico esta ecuación es llamada una ley de tasas de primer orden, en tanto que
las leyes de tasas como dP/dt=-kP2 son de segundo orden.
381
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Luego:
N (t) = 125000
Así:
140000 = 125000
Por consiguiente:
P (t) = 125000
= 125000
=125000
P = P (40) = 125000
382
Aplicaciones de las ecuaciones de primer orden
Solución:
El modelo del decaimiento radioactivo, es muy similar al del
decrecimiento exponencial, solo que en este caso se sabe la razón
a la que un elemento radioactivo decae en un tiempo cualquiera en
forma proporcional a la cantidad presente en ese tiempo.
0.6 ,o
0.6 =
-50k = ln (0.6)
K= = 0.01022
8 La científica polaca Marie Curie (1867-1934) recibió el premio Nobel de física en 1903 y
en 1911 por sus experimentos precursores con el radio y otras sustancias radioactivas.
383
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
N0/2=N0
tm =
384
Aplicaciones de las ecuaciones de primer orden
7.2 Problemas
de la ecuación logística
Es claro que, bajo un crecimiento exponencial, la población llegaría a
ser infinitamente grande con el transcurso del tiempo. Sin embargo,
en la realidad cuando una población llega a ser suficiente existen
factores ambientales que hacen más lenta la razón de crecimiento.
Ejemplos de estos factores podrían ser: la disponibilidad de alimentos,
los depredadores, la falta de espacio suficiente, etc. Estos factores
ocasionan que en algún momento comience a decrecer y en
este caso, se puede suponer que el tamaño de la población esté
limitado a cierto número máximo M, donde o < P < M y cuando P M,
0, y el tamaño de la población tiende a ser estable.
385
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 3:
Tabla 1. Tabla de valores del crecimiento de una planta de Girasol, tomada en días.
386
Aplicaciones de las ecuaciones de primer orden
Solución:
Primero tratemos de hacer una gráfica con dichos datos:
Solución cualitativa
Una manera de describir este crecimiento mediante ecuaciones
diferenciales es agregar un término al miembro de la derecha de la
ecuación diferencial el cual disminuirá la tasa de crecimiento
para valores grandes de y. Esto debería hacerse de tal manera que
la tasa de cambio de altura por altura unitaria (la tasa de crecimiento
relativo) disminuye a medida que la altura aumenta.
La ecuación:
387
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(40)
388
Aplicaciones de las ecuaciones de primer orden
Las isóclinas
Las isóclinas de (40) están dadas por:
Ay (b - y) = m (41)
Si m=m0, vemos que los únicos lugares donde las curvas solución
tienen tangentes horizontales se encuentran a lo largo de las líneas
horizontales y = 0, y = b. Para considerar otros valores de m, volvemos
a escribir (41) como una ecuación cuadrática en y, así:
(42)
389
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Monotonía
El campo de pendientes también sugiere que las soluciones de (40)
son, ya sea creciente o decreciente, dependiendo del valor inicial de
y. esto puede verificarse a partir de (40); es decir, y/= a. y. (b – y),
en la que la derivada es positiva si 0 < y < b y negativa si y > b.
Luego, si comenzamos con un valor y0 para x = 0, donde y0 > b, la
solución disminuirá hasta el valor límite b. De manera semejante, si
comenzamos un valor inicial de y0, donde 0 < y0 < b, la solución se
incrementará hasta su valor límite b. En contraste con el modelo de
crecimiento de poblaciones, las soluciones de la ecuación diferencial
logística para y0 > 0 se aproximan al valor de b para valores grandes
del tiempo. Este valor límite de b se denomina a menudo capacidad
de transporte de una población específica para situaciones en que y
denota una población.
Soluciones de equilibrio
Existen dos condiciones iniciales específicas para las cuales el análisis
anterior se demora, a saber, y0 = 0, y y0 = b. Si y0 = b, entonces de la
ecuación diferencial vemos que su tasa de cambio es siempre cero,
(recordemos que y = b era una isóclina para pendiente cero). Advierta
390
Aplicaciones de las ecuaciones de primer orden
Puntos de inflexión
La figura 3 muestra una ampliación del campo de pendientes entre las
dos soluciones de equilibrio, donde vemos que hay puntos de inflexión.
Para descubrir la ubicación exacta de los puntos de inflexión, derivamos
la ecuación diferencial original (40) con respecto a t para obtener:
391
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Concavidad
Debido a que y es monótonamente creciente para 0 < y < b, acabamos
de descubrir que las soluciones de la ecuación logística para y > 0
tendrán un punto de inflexión solamente cuando o < y0 < b/2. En las
otras soluciones las curvas solución son crecientes cóncavas hacia
abajo (para b/2 < y0 <b) o bien decrecientes y su concavidad hacia
arriba (para y0 > b). Algunas curvas solución trazadas a mano para la
ecuación logística se ilustran en la figura 4. Observe que se han dado
muchas características esenciales de las soluciones de la ecuación
logística sin tener la solución explícita.
(43)
Y (t) = 0 y, y (t) = b
392
Aplicaciones de las ecuaciones de primer orden
Integrando, se tiene:
ln (44)
(45)
(46)
393
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
C=
y (t) = (47)
394
Aplicaciones de las ecuaciones de primer orden
Solución:
Observe que el crecimiento del número de socios del club, tiene un
tope de 800. Si llamamos N (t) el número de socios en un tiempo t.
Entonces la cantidad de cupos para nuevos socios será 800 – N (t)
y así en número de no socios que se quieren afiliar por unidad de
tiempo serán: (800 – N(t))/800, y por lo tanto la razón de cambio de
los socios será directamente proporcional entre los que son socios y
los que se quieren afiliar, es decir:
395
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Tomando exponencial:
Despejando N (t),
N (t) =
50 = , o
1+b =
b = 15
396
Aplicaciones de las ecuaciones de primer orden
Así,
N (t) = (48)
200 =
1+15
N (t) =
N (t) =
397
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Sea N (t) el número de estudiantes que conocen el rumor en t
semanas después de que comienza la investigación. Entonces, como
el crecimiento es de tipo logístico podemos usar la ecuación, que
obtuvimos en el ejemplo 5, es decir:
N (t) =
300 = ,o
1+b=
b = 149
Por lo tanto:
N (t) =
900 =
1 + 149 =
N (t) =
Cuanto t = 4,
398
Aplicaciones de las ecuaciones de primer orden
P/ (t) =
P (t) = C ert
P (t) = P0 ert
399
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 6:
Solución:
Para este caso tenemos que P0= 1000, y r = 0.05. Luego:
P (40) = 1000
400
Aplicaciones de las ecuaciones de primer orden
Una vez más, el interés ganado está dado por (50). Por otro
lado, si W es la cantidad retirada por unidad de tiempo, la cantidad
retirada en el tiempo es:
Luego
P (t+ - P (t)
= rP – W (51)
= rP + D (52)
Ejemplo 7:
401
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Aplicando nuestro modelo (52) y tomando P (t) en miles de dólares,
tenemos:
Luego:
O,
0 = P (0) = - 40 + C
Luego:
C = 40
P (t) = 40 (53)
402
Aplicaciones de las ecuaciones de primer orden
Ejemplo 8:
Solución:
Sea P (t) el saldo de su cuenta de jubilación en el tiempo t después de
la jubilación. Sea P0 el saldo de su cuenta de jubilación cuando se retire.
Entonces P (0) = P0. El problema es encontrar P0 con P (30) . De
nuevo usamos en miles de dólares como nuestra unidad.
(t) =
403
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
De donde se tiene:
Integrando, se tiene:
o,
Por lo tanto:
C = P0 - 1000
P (t) = 1000 +
Puesto que usted quiere retirar $50.000 cada año hasta que muera
por 30 años después de retirarse, tendrá que ser P (30) . Si usted
tiene su cuenta en cero al cabo de los 30 años, entonces se tendrá que:
0 = P (30) = 1000 +
404
Aplicaciones de las ecuaciones de primer orden
Ejemplo 9:
Solución:
En primer lugar, se da cuenta de que su sueldo actual de $35.000
por año no se quedará para siempre en ese nivel, con suerte por no
más de un año. Necesitamos un modelo de cómo su salario crecerá
con el tiempo. Vamos a suponer que su salario se incrementará en
4 % por año. Eso es solo un poco más que la tasa de inflación. Este
pensamiento conduce a la ecuación diferencial S/ = 0.04 S (t) donde
S (t) es el salario anual en miles de dólares. Estamos familiarizados
con la ecuación exponencial, por lo que es fácil resolver esta ecuación
para obtener:
405
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
P (t + ) – P (t)
Por lo tanto:
Integrando, se tiene:
O:
P (t) = -3500
406
Aplicaciones de las ecuaciones de primer orden
0 = P (0) = - 3500 + C
P (t) = 3500
Ejemplo 10:
407
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
El modelo en este caso es similar al del ejemplo 9, pero tiene que cam-
biar para adaptarse a la variabilidad de su tasa de ahorro, es decir:
0.05P (t)
También:
O:
408
Aplicaciones de las ecuaciones de primer orden
Integrando:
Ln(T-Tm) = kt + C1, o
T – Tm=
409
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 11:
Solución:
Sea T (t) la temperatura del objeto en el tiempo t. Escogemos por
conveniencia el origen de la escala del tiempo t0= 0 a las 11: 05, de
manera que T0= 60 0F. Debemos determinar el tiempo t tal T (t) =
72 0F. Al sustituir en nuestra ecuación (55), T0= 60 0F y Tm= -20 0F,
tenemos que:
T (t) = -20 + C
T (t) = -20 + 80
T (2) = 50 = -20 + 80
Despejando k:
410
Aplicaciones de las ecuaciones de primer orden
K=-
T (t) = -20 + 80
T (t) = 72 = -20 + 80
Despejando t:
Ejemplo 12:
411
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Sea T (t) la temperatura en grados Fahrenheit en el tiempo t, Tm= 70 0F,
también se tiene que T (10) = 100 0F y tomamos T (0) = 120 0F. Por la
fórmula de Newton se tiene que:
T (t) = 70 + C (56)
Despejando:
412
Aplicaciones de las ecuaciones de primer orden
T (t) = 70 +
100 = T (20) = 70 + ,
O:
30 =
T (0) = 70 + F
0
413
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
414
Aplicaciones de las ecuaciones de primer orden
415
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
V0 V1
416
Aplicaciones de las ecuaciones de primer orden
v = f (i, t)
417
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
q = C v, o (57)
v=
i= (58)
I=
Inductor
Un inductor almacena energía en campo magnético. La tensión-
corriente actual de un inductor lineal (llamada ley de Faraday y Lenz)
es:
418
Aplicaciones de las ecuaciones de primer orden
419
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
=0
420
Aplicaciones de las ecuaciones de primer orden
(59)
I= (60)
O:
421
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
E L E C
(61)
Ejemplo 13:
Solución:
Por la ecuación (61), tenemos que:
, I (0) = 0
422
Aplicaciones de las ecuaciones de primer orden
I (t) = 2
(62)
Ejemplo 14:
Solución:
Con los datos dados de la ecuación (62), se tiene:
423
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Integrando, se tiene:
I (q) =
I (t) =
424
Aplicaciones de las ecuaciones de primer orden
Ejemplo 15:
Solución:
La corriente i es la misma en cada rama. En este caso, por la ley
del voltaje, directamente obtenemos la ecuación diferencial para la
corriente:
I (t) = C e-3t
I (t) = 6 e-3t
425
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
426
Aplicaciones de las ecuaciones de primer orden
7.6 Mezclas
Vamos a ver otros ejemplos sobre este tema de mezclas, conviene
recordar el ejemplo 13 del capítulo 2. Profundizaremos más sobre
este tópico por considerarlo de mucha aplicación.
(63)
427
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 16:
428
Aplicaciones de las ecuaciones de primer orden
Solución:
Con el fin de ilustrar los principios generales, y dado que este es
nuestro primer problema de la mezcla, vamos a incluir algunos pasos
más que son necesarios para resolver este problema en particular.
429
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
=0.75 -
430
Aplicaciones de las ecuaciones de primer orden
(64)
Q (t) = 15 + 45
431
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
432
Aplicaciones de las ecuaciones de primer orden
433
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
V (t) = V (0) +
434
Aplicaciones de las ecuaciones de primer orden
Ejemplo 17:
Solución:
Este problema se gráfica con la figura 19. El volumen de agua con sal
está creciendo y en tiempo t este volumen está dado por:
V (t) = 50 + 2t
435
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(66)
Integrando:
C (t) =
436
Aplicaciones de las ecuaciones de primer orden
437
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La clave de este tipo de problemas de mezclas, es la relación que se
ha visto en los problemas 15 y 16, esto es:
Como las dimensiones del local son: (12.5). (10). (3.2) = 400
m , la expresión para la cantidad de monóxido en este espacio en el
3
438
Aplicaciones de las ecuaciones de primer orden
Integrando, obtenemos:
C (t) =
Donde
0.004 = 0.12
1/30 = 1 -
439
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
-0.0001t = ln(29/30)
t=
Solución:
Sea C (t) la concentración del medicamento en el instante t (en mg/
cm3). De manera que las condiciones de este tipo de problema nos
conducen a la relación:
440
Aplicaciones de las ecuaciones de primer orden
O:
(68)
O:
441
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Es decir que:
C (t) = 20(1 -
442
Aplicaciones de las ecuaciones de primer orden
Posición = y (t)
Velocidad= v (t) = ,
Aceleración =
.= F (y, , t)
443
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 20:
Solución:
A continuación, se sabe que F = -mg, donde g es la aceleración debida
a la gravedad. El signo menos se introduce porque la gravedad actúa
hacia abajo desde la superficie de la tierra. Aquí estamos tomando
el sistema de coordenadas de modo que y aumenta hacia el cielo. La
magnitud de la fuerza debida a la gravedad, mg, se llama el peso del
cuerpo. Cerca de la superficie del planeta tierra g es aproximadamen-
te g = 9,8 m / seg2 en el sistema mks, utilizado por la mayoría del
mundo (g = 32 pies / seg2 cuando se utilizan los pies como la unidad
de longitud en lugar de metros). Si suponemos que estamos intere-
sados en una masa que se encuentra lo suficientemente cerca de la
superficie de la tierra, entonces g se puede aproximar por la constan-
te. Si la única fuerza es la gravedad, entonces nuestra ecuación se
convierte en:
O de forma equivalente:
(69)
v (t) = - gt +
444
Aplicaciones de las ecuaciones de primer orden
y=-
y=-
Supongamos que se lanza una pelota hacia arriba desde el nivel del
suelo con velocidad v0 y la única fuerza que actúa sobre ella es la
gravedad. ¿Qué tan alto la bola alcanza antes de caer hacia el suelo?
Solución:
La ecuación diferencial que nos modela esta situación es como antes
la dada por (69), es decir:
445
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Y por lo tanto:
Y= (70)
Solución:
Sea y el recorrido a medir una vez que se aplican los frenos. La
ecuación diferencial es:
Integrando obtenemos:
(71)
446
Aplicaciones de las ecuaciones de primer orden
(72)
76 = -4(3)3 + C1 = - 32 + C1
0 = - 24 + C1 t + C2 = - 16 + 216 + C2
447
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
F=ma
m. v/= F (t, v)
448
Aplicaciones de las ecuaciones de primer orden
Solución:
La fuerza que actúa sobre el objeto es:
F = - mg + F1 (69)
F = -mg – kv (70)
F = ma = mv/
m v/= -mg – kv
O:
449
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(v. )/= -g
Integrando:
v. =
v=- (71)
V0 = -
Entonces:
C = v0+
v=
450
Aplicaciones de las ecuaciones de primer orden
Solución:
Como mg = 640 y g = 32, se tiene entonces que m = 640/30 = 20. La
resistencia de la atmósfera es -3.v lb, si v está expresada en pies por
segundo. Luego:
20 v/= -640 – 3v
(72)
451
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Vh=C
v = vh+ vp= -
v=-
Ejemplo 25:
Solución:
Como el objeto está cayendo, la dirección de la resistencia del aire es
hacia arriba (positiva). De donde:
453
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Integrando se tiene:
Luego:
C=
454
Aplicaciones de las ecuaciones de primer orden
C. Velocidad de escape:
Debido a los trabajos de los astrónomos y las leyes empíricas de Tycho
Brache (1546-1601) y de Johannes Kepler (1571-1630) relativos a las
órbitas de la luna y los planetas, Newton centró su atención en una
sola fuerza, tal como se vio en las secciones anteriores: la gravedad.
Su ley de gravitación universal explica el efecto gravitacional de un
cuerpo sobre otro.
455
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 65:
456
Aplicaciones de las ecuaciones de primer orden
F=
457
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
O,
(75)
Si:
458
Aplicaciones de las ecuaciones de primer orden
459
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
460
Aplicaciones de las ecuaciones de primer orden
461
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Resuelva para :
462
Aplicaciones de las ecuaciones de primer orden
Ejemplo 27:
Solución:
De acuerdo al paso 1, primero escribimos la familia de curvas en la
forma t2 + 2x2 = c. Luego calculamos la derivada con respecto a t:
Integrando:
463
1
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
x
-8 -6 -4 -2 2 4 6 8
-1
-2
X
K=2
c=8 K=1/2
c=4 K=1/8
c=1
K=0 Y
K= -1/8
K= -1/2
K= -2
464
Aplicaciones de las ecuaciones de primer orden
Solución:
Al derivar con respecto a t, se obtiene:
3t2 + 6x
Luego:
Ejemplo 29:
Solución:
Esta familia de curvas son círculos con centros fuera del origen, véase
la figura 28, se puede ver que esta familia de figuras son simétricas al
eje x y, al origen (basta ver que la familia dada se puede escribir como
). Vamos a ver que las trayectorias ortogonales van
a tener la misma simetría. Derivamos la ecuación diferencial dada, de
manera implícita para obtener:
465
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y
10
-10 -5 0 5 10 x
-5
-10
y
10
Despejando c de la ecuación
5 original, se tiene:
-10 -5 0 5 10 x
-5
Reemplazando en la ecuación anterior, tenemos:
-10
Despejando , se tiene:
466
Aplicaciones de las ecuaciones de primer orden
5
O bien:
-10 -5 0 5 10 x
-5
Regresando a nuestras variables originales y comple-
tando el cuadrado perfecto
-10 se obtiene la familia ortogonal
, que es la familia ortogonal a la
dada (figura 29).
y
10
-10 -5 0 5 10 x
-5
-10
467
7.9 Ejercicios
7.9.1 Problemas de crecimiento y
decrecimiento
V L
La corriente inicial es cero. Encuentre i(t) para
y grafique el resultado.
i
Figura 31. Gráfica de un circuito RC.
R
+
-
i
Fuente. Elaborada por el autor.
7.9.6 Mezclas
79. x = t + k 80. X = k t
81. x t = k 82. x2 = k (t + x)
85. x = t2 + k 86. x3 (t + 1) = k
87. t = (x – k)2 88. X = e k t
93. x = k t2 94.
SOLUCIÓN DE
LAS ECUACIONES
DIFERENCIALES DE
SEGUNDO ORDEN Y DE
GRADO SUPERIOR
Contenido
8.1 Solución general de una ecuación lineal
de segundo orden 494
8.2 Ejercicios 500
8.3 Problemas de valor inicial
(para ecuaciones homogéneas) 502
8.4 Ejercicios 510
8.5 Reducción del orden 514
8.6 Ejercicios 523
8.7 Ecuaciones diferenciales lineales de
segundo orden homogéneas con
Competencias coeficientes constantes 525
8.8 Ejercicios 540
1. Identificar ecuaciones 8.9 Ecuaciones diferenciales lineales
diferenciales lineales homogéneas de orden n con
homogéneas y no coeficientes constantes 544
homogéneas. 8.10 Ejercicios 547
2. Deducir cómo se resuelven 8.11 Ecuaciones diferenciales lineales
las ecuaciones diferenciales no homogéneas 550
lineales homogéneas y no 8.12 Ejercicios 571
homogéneas con coeficientes 8.13 Método de variación de parámetros
constantes. (ecuaciones de segundo orden) 580
8.14 Ejercicios 593
3. Aplicar las ecuaciones 8.15 Método de Variación de
diferenciales lineales parámetros (de orden n) 595
homogéneas y no 8.16 Ejercicios 598
homogéneas a problemas de 8.17 Ecuación de Cauchy-Euler 600
la vida cotidiana. 8.18 Ejercicios 615
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
492
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Introducción
Introducimos un modelo central en física e ingeniería, el oscilador
armónico, la ecuación diferencial lineal de segundo orden prototípico
de estas situaciones. Se derivan problemas de valor inicial que
describen las vibraciones mecánicas y el comportamiento de los
circuitos eléctricos, amortiguados y no amortiguados, que ilustran
trayectorias en el plano de fase.
493
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y (x) = y p( x) + c y h (x)
494
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(1)
La solución general es de la forma:
y(t) = yp( x) + c1y1(x) + c2y2(x) (2)
Donde:
· yp(x) es una solución particular de dicha ecuación.
· y1 y, y2 son soluciones para la ecuación diferencial homogé-
nea asociada:
(3)
495
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(4)
Ejemplo 1:
496
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
a. La verificación de que cos x es una solución de la ecuación
homogénea es sencillo, pero ofrece una buena
práctica. Y = cos x, su primera derivada es y/ = -senx, su segunda
derivada es y//= -cos x, y . Similarmente,
se puede mostrar que y = sen x es otra solución. También
tenemos que es una ecuación de segundo orden,
entonces por (5) {sen x, cos x} es un conjunto fundamental
de soluciones y entonces y = c1 sen x + c2 cos x da todas las
soluciones de .
b. Es fácil verificar que y p(x) = e-3x satisface la ecuación
Observe que y, , luego ,
1 = y/ (0) = - 3 + c1
497
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(12)
Es de la forma:
y (x) = y p + c1 y1(x) + c2 y2(x)
(13)
(14)
498
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Introduciendo y h
= y – y p, la diferencia entre y y, yh, se tiene
que y = yp + yh. Para calcular la ecuación diferencial que satisface yh,
calculamos como sigue:
c1y1(x) + c2y2(x)
(15)
499
8.2 Ejercicios
En los ejercicios 1 al 8, verifique que el conjunto dado es un conjunto
de soluciones asociados a la ecuación homogénea, y verifique que yp
es una solución particular. Luego dar la solución general de la ecuación
diferencial, y resolver el problema de valor inicial.
1.
2.
3.
4.
5.
6.
7.
8.
Y:
y2(x) = d1 (1 + x) + d2 (1 – x) +x2 + 3x +1
(16)
(16)
y( x0 ) = u0
y/(x0) = v0
502
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Asumamos que:
Y: (17)
v 0 = c1
503
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(18)
504
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
El Wronskiano
Definimos el Wronskiano1 W de dos funciones y1, y2 como:
(19)
505
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(20)
(21)
Solución fundamental
Hemos mostrado en la sección anterior que c1y1 +´c2y2 resuelve
la ecuación diferencial lineal homogénea. Sin embargo, es posible
que no se pueda resolver el problema con valor inicial para esta
solución. Para poder resolver este problema con valor inicial, lo
primero que se tiene que dar es que el Wronskiano de y1, y2 tiene
que ser diferente de cero. Entonces dos soluciones y1, y2 de una
ecuación diferencial lineal homogénea de segundo orden forman el
conjunto fundamental y se puede mostrar que ellas son linealmente
independientes, sí y solo sí W(x0) . Entonces el criterio para saber
si un conjunto de soluciones es fundamental basta con ver que su
506
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
y h = c1y1 + c2y2
Solución:
a. Primero debemos verificar que y1 = sen x, y y2 = cos x satisfa-
cen la ecuación (22). Se omite este paso, ya que este hecho
fue verificado en el ejemplo 1 para mostrar que {senx, cos x}
es un conjunto fundamental de soluciones, debemos ver que
W [y1, y2] ≠ 0 en x = 0. Observe que .
Entonces:
507
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y / = c1cos x – c2 sen x
y (0) = 6 = c2
y/ (0) = 3 = c1
(24)
508
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
509
8.4 Ejercicios
y // + x y/ + y = 0, y(0) = 1, y/(0) = 2
y, y2 la solución de:
Y3 = c1 y1 + c2 y2
y // + p y/ + q y = 0, y(x0) = 1, y/(x0) = 0
y, y2 la solución de:
y // + p y/ + q y = 0, y(x0) = 0, y/(x0) = 1
y // + p y/ + q y = 0, y(x0) = α, y/(x0) = β
y // +y = 0
16. Verifique que {sen 2x, cos 2x} forman un conjunto fun-
damental de soluciones para y// +4y = 0.
Verifique y3(x) = sen (2x + ) es otra solución de y// +
4y = 0. Encuentre c1 y c2 tal que:
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(25)
514
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
y (27)
haciendo:
w = u/ (28)
515
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(29)
Solución:
Como y1 = x, entonces se tiene que:
Luego:
516
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
517
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y 20
15
10
+ x
-20 -15 -10 -5 5 10 15 20
-5
-10
-15
-20
y // -x y/ + y = 0
Solución:
Como y1 = x, entonces se tiene que:
Luego se tiene: -x + x = 0.
518
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
x. u// + (2 – x2) u/ = 0
Finalmente tenemos:
519
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Ejemplo 5:
Solución:
Sea y = erx y sustituimos en:
y // - 2y / + y = 0
Para obtener:
r 2 e rx — 2re rx +e rx =0
O:
r 2 — 2r+1
Luego:
520
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
u // e x = 0, luego u // = 0
u/ = c1
u = c 1 x + cc
Entonces:
Y = u y 1 = (c 1 x + c 2 ) e x = c 1 xe x +c 2 e x
521
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
522
8.6 Ejercicios
En los ejercicios 1 al 8, verifique que la función dada y1 es una solución
de la ecuación diferencial homogénea dada. Luego use el método
de reducción del orden para encontrar el conjunto fundamental de
soluciones y la solución general para la ecuación diferencial. Los
ejercicios 7 y 8 requieren integrales definidas.
4. x y// - (1 + x) y/ + y = 0, y1 = ex-
7. y// +x y/ - y = 0, y1 = x
8. x2y// +x y/ - 2(1+x) y = 0, y1 = x2
9. x2y// - 3xy/ + 4y = 0
19. y// - 4y = 0
ay// + by/ + c y = 0
y = c1y1 + c2y2
525
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
a (e rx)// + b ( e rx)/+ c e ex = 0
Y derivando, se tiene:
ar2erx + b r e rx + c r rx = 0
a r2 + b r + c = 0
a y// + b y/ + c y = 0
a y// + b y/ + c y = 0
526
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Sustituyendo y = e rx, en la ecuación diferencial, obtenemos la ecuación
característica r2 – r – 20 = 0. Factorizando esta expresión obtenemos:
r2 – r – 20 = (r — 5)(r + 4) = 0
3 Las soluciones exponenciales inicialmente las quiso resolver Leonhard Euler (1707-1783),
uno de los más grandes matemáticos de todos los tiempos, además de los más prolíferos a todas las
ramas de la matemática.
527
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y = c1e5x + c2e-4x,
y// + 4 y/ = 0
Solución:
Por sustitución y = erx encontramos que la ecuación característica
es r2+4r = r(r + 4) = 0. Luego hay dos soluciones reales distintas r =
0, r = - 4. Un conjunto fundamental de soluciones de esta ecuación
será {1, e-4x} ya que e0.x = 1. La solución general es:
528
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
y // - 6 y/ + 9y = 0
Solución:
Sustituyendo y = erx en y// - 6 y/ + 9y = 0, encontramos que la ecuación
característica es r2 – 6r + 9 = (r – 3)2 = 0. Luego r = 3 es una raíz
529
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
O, derivando, se tiene:
u//e3x = 0 o, u// = 0
u = c1x + c2
530
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Sustituyendo por y = r rx, encontramos que la ecuación característica
es: r2 + 2r + 1 = (r + 1)2 = 0. Luego r = - 1 es una raíz de multiplicidad
dos. Una solución de la ecuación homogénea es entonces y1 = e-x. Como
la raíz es repetida, la segunda solución de acuerdo al resultado anterior
es x e-x. Luego la solución general de la ecuación homogénea es:
y = c1e-x + c2xe-x
Donde
531
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La ecuación característica es:
y = c1 + c2x + c3ex + c4 x ex + c5 x2 ex
532
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
ar2 + b r + c = 0
r1 = α+ iß (i2 = -1)
r2 = α — iß
y= (32)
Sin embargo esta no es una solución muy útil para algunos pro-
e r1x, como e r2x involucran números com-
blemas de la física, ya que
plejos. Podemos reemplazar {e r1x, e r2x} por un conjunto fundamental
de soluciones diferentes. Las constantes en (32) son denotadas por
y para distinguirlas de c1 y c2, usadas en las otras soluciones.
533
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
O:
(33)
Formula de EULER
La relación entre la exponencial compleja y seno y coseno se sigue de
la fórmula de Euler4:
(34)
(35)
4 Recordemos que la fórmula de Euler se derivó de los desarrollos de las series de Taylor de
cada una de las funciones que intervienen.
534
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(36)
(37)
Donde:
535
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
y // + y/ + y = 0
536
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Sustituyendo y = erx, vemos que la ecuación característica es r2 + r +
1 = 0. Por la fórmula cuadrática, las raíces son:
Como , tenemos:
537
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Sustituyendo a y= erx, vemos que la ecuación característica es r2 + 4
= 0, luego las raíces son r = α ± i ß = ± 2 i. Luego α = 0, ß = 2, y el
conjunto fundamental de soluciones es:
y = c1cos2x + c 2sen2x.
538
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
, es una solución general si r = -3 y, r = -2 son las raíces
de la ecuación cuadrática. En este caso la ecuación característica
tiene que ser de la forma:
y // + 5y/ + 6y = 0
539
8.8 Ejercicios
En los ejercicios 1 al 33, resuelva la ecuación diferencial. Determine la
solución general si no se dan las condiciones iniciales.
1. y//+ y/ - 6y = 0
2. y// - y = 0
3. y// + y = 0
4. y// + 4y/ + 4y = 0
5. y// + 4y/ + 5y = 0
6. y// - 2y/ + y = 0
7. y// - 3y/ + y = 0
8. 2y// + 2y/ + y = 0
9. y// - y/ = 0
11. 3y// = 0
17. 2y// + 4y = 0
21. y// + 10 y/ + 25 y = 0
26. y//+ y/ + y = 0
27. y//- 12 y = 0
30. y//+ 10 y = 0.
31. y//+ 8 y = 0
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(40)
544
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
La ecuación característica es , que se puede
encontrar por la sustitución . Luego son raíces complejas
conjugadas de multiplicidad 2. En la notación de la forma se
tiene:
545
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Como se tienen 7 raíces, la solución general es una combinación lineal
de 7 soluciones homogéneas usando siete constantes arbitrarias:
546
8.10 Ejercicios
En los ejercicios 1 al 23, encuentre la solución general. Aquí y(n) denota
la derivada n – éssima de y.
2. y///+ 5y// - 4y = 0
3. y/// - 5y// + 4y = 0
5. y///+ y// - 2y = 0
6. Y(4)- 2y/// = 0
9. y/- 3y = 0
10. y/+ 4y = 0
19. y(4) - y = 0
22. 3y/+ 4y = 0
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
a y// + b y/ + c y = f(x)
550
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
1. F(x) = polinomio.
2. F(x) = exponencial = eax.
3. F(x) = sinusoidal = sen ßx, o cos ßx.
4. F(x) = productos de funciones de tipo 1, 2, 3.
5. F(x) = combinación lineal de 1, 2, 3, 4.
551
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
de la ecuación característica.
552
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Primer paso: la ecuación característica es r2 – 3r + 2 = (r - 2)
(r – 1) = 0, tiene dos raíces r =1, y r = 2) Entonces yh = c1ex + c2 y2x.
yp = A0 + A1x + A2x2
O:
-6A2 + 2A1 = o
2A2 = 2
A0 = 4, A1 = 3, A2 = 1
y = yp + yh = 4 + 3x + x2 + c1ex + c2 y2x
553
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Primer paso: la ecuación característica es r2 – 3r = (r – 3) r = 0,
tiene una raíz r = 0 y r = 3. Entonces la solución de la ecuación homo-
génea es: yh = c1 + c2 e3x.
554
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
y// + y =3 e—x
Solución:
Primer paso: por sustitución de y = erx en la ecuación diferencial
homogénea o por inspección, vemos que la ecuación característica es
r2 + 1 = 0, la cual tiene como raíces r = i, y r = - i.
555
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
yp = A e—x
y// - y = 3 e—x
Solución:
La ecuación característica es r2 – 1 = 0, la cual tiene las raíces r = 1
y r = -1. Aquí f(x) = 3 e—x es una función exponencial. Por (44), y p
=
Ax e k —x
donde k es el número de veces en que α = –1, es una raíz
de la ecuación característica r2 – 1. Como r = -1 es una raíz una vez,
tenemos que k = 1.
y p = A x e—x
556
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
El principio de superposición
557
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La ecuación característica es r2 – r = r(r – 1) = 0, la cual tiene los ceros
r = 0, r = 1, luego yh = c1 + c2ex. La función forzada es:
y// - 2y/ + y = 7x ex
558
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
La ecuación característica es r2 – 2r + 1 = (r – 1)2 = 0 tiene la raíz r = 1
de multiplicidad 2. Entonces yh = c1 ex + c2 x ex. La función forzada es
de la forma p(x) eαx, donde p(x) = 7x es un polinomio de primer grado y
r = α = 1. Como r = 1 es una raíz de multiplicidad 2, por (43) con k = 2,
m = 1, se tiene que yp tiene la forma:
Solución:
Primero se requiere encontrar la solución de la ecuación homogénea
asociada y/// - 6y// + 12 y/ - 8y= 0. La ecuación característica es r3 –
6r2 + 12 r – 8 = (r – 2)3 = 0, luego r = 2 es una raíz de multiplicidad 3.
Un conjunto fundamental de soluciones para la ecuación homogénea
asociada es {e2x,xe2x,x2e2x}, luego yh = c1 e2x + c2xe2x + c3x2e2x. El término
correspondiente a la función forzada es f(x) = xe2x. Por (43) con α
= 2, m = 1, xp tiene que ser de la forma xk[A0 + A1x]e2x, donde k es la
multiplicidad de 2 que es la raíz de la ecuación característica. Entonces
k = 3 y yp es de la forma:
559
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Dar la forma de y p si
Solución:
La ecuación característica es r2 + 2r + 2 = 0. Sus raíces son r = - 1 + i
y r = - 1 – i. Entonces yh = C1 e–x + C2 e–x. Ahora vamos a determinar yp.
Observe que f es la suma de dos términos, 3e–x y, 4 cos x. Como – 1
no es una raíz de la ecuación característica. El caso (43) nos dice que
yp incluye un término de la forma A1 e–x. Ahora consideremos a 4 cos
x. Como i no es una raíz de la ecuación característica, por el caso 4
con ß = 1 nos dice que yp incluye términos de la forma A2 cos x + B2
sen x. Entonces:
560
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
La ecuación característica es r2 + 4 = 0, cuyas raíces son r = ßi ±2i.
Entonces y h = C1 cox 2x + C2 sen 2x. El término de la función forzada
sen (2x) es sen (ßi, donde ß=2. Como r = ßi es una raíz de la ecuación
característica de multiplicidad 1, tenemos que k = 1 y por el caso 4 la
forma de la solución particular es:
y p = x [A cos(2x) + B sen(2x)]
561
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
La ecuación característica es r2 + 4 = 0 tiene como raíces a ±2i y
entonces y h = C1 cos (2x) + C2 sen (2x). El término de forzamiento es
de la forma:
562
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Primero resolvemos la ecuación homogénea asociada:
y (4) + 8y // + 16 y = 0
(45)
(46)
563
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Si ,
en la ecuación diferencial
, entonces yp tiene la forma:
564
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
La ecuación característica es r2 + 2r + 2 = 0, la cual tiene como raíces
r = –1±i , luego se tiene:
Solución:
Esta ecuación tiene la misma ecuación característica del ejercicio 22,
luego las raíces son r = –1±i . La función forzada en esta ecuación es
la suma de tres tipos de términos:
565
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Si ,
donde P(x) es un polinomio
de grado m y Q(x) es un
Este caso es el caso más general. Este incluye todas los casos
anteriores para apropiadas escogencias de α, ß, m y n.
y // + 3 y/+ 2y = x e-xcos2x
566
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
La ecuación característica es r2+ 3r + 2 = (r + 1) (r + 2) = 0, luego las
raíces son r = - 1 y r = - 2. Entonces, dos soluciones independientes
de la homogénea son e-x y e-2x. El término forzado es de la forma P(x)
eαxcosßx, donde P(x) = x, es de primer grado y α + iß= - 1 + 2i. Como – 1
+ 2i no es una raíz de la ecuación característica, por el caso 7 (con k
= 0), la forma de yp es:
Solución:
La ecuación característica es r2+ 2r + 5 = 0, y sus raíces son r = - 1 + 2i.
Entonces, las dos soluciones linealmente independientes homogéneas
son . Los términos de la función forzada son
de la forma , donde es un polinomio de tercer
grado y . Como r = - 1 + 2i es una raíz de la ecuación
característica de multiplicidad 1, por el caso 7 (con k = 1), la forma de
yp es:
567
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Resolver:
Solución:
Primero encontramos la solución general de . En el ejemplo
14 encontramos que , y que la ecuación característica es
r - 1 = 0, tiene como raíces r = 1, - 1. Entonces la solución general es:
2
568
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
El método es el siguiente:
(1)
569
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
570
8.12 Ejercicios
En los ejercicios 1 al 12, mencione si el método de coeficiente
indeterminado se puede aplicar a las ecuaciones diferenciales dadas.
Si no se puede, explicar la razón.
1. y//+ y = x tan x
4. y//+ y = ex +1
5. y//+ y =
6. y//+ y/ + y = cosh x
7. y//+ x y = 3e2x
8. y//+ y = x senh 2x
9. y//+ y = x –1ex
14. y// - 4y =
50. 2 y/+ 4y = x
y p= Q(x)
581
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(52)
(53)
(54)
(55)
582
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(56)
(57)
(58)
Y:
583
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Tenemos que y1 = x, y2 = x3. Dividiendo por x2 la ecuación diferencial
se tiene:
O:
Y resolviendo para :
Entonces:
584
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Primero se obtiene la solución de la ecuación diferencial homogénea
asociada que es: y//+ y = 0. Sustituyendo por y = erx, obtenemos la
585
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(59)
(60)
(61)
(62)
Como , tenemos:
586
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(63)
Función de influencia
A pesar de que la solución general en (63) es correcta y satisfactoria
para la ecuación diferencial dada, algunas manipulaciones algebraicas
nos permiten dar un importante e interesante resultado.
587
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
(64)
(65)
588
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
589
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Observación:
1. Si constantes arbitrarias son introducidas en el paso c, cuando
encontramos v1 y v2, entonces y = v1 y1 + v2 y2, será la solución
general.
2. Los coeficientes de en (66) son las entradas de la
matriz:
(67)
(68)
590
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Observe que la función x-1ex no es del tipo de función forzada a la cual
se le pueda aplicar el método de los coeficientes indeterminados. Sin
embargo, la ecuación diferencial tiene coeficientes constantes, luego
conocemos cómo resolver la ecuación homogénea asociada. Podemos
resolver esta ecuación por el método de variación de parámetros.
Escribimos la ecuación dada como:
Luego:
591
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Paso 3:
Paso 4 y 5:
592
8.14 Ejercicios
En los ejercicios 1 al 14, encuentre la solución de la ecuación
diferencial dada por el método de variación de parámetros. Decida
dónde el método de coeficientes indeterminados puede ser usado.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15. Muestre que una solución particular de y// + P(x)
y/+ Q(x) y = f(x), dada por variación de parámetros
puede ser escrita por:
17.
18.
19.
20.
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(69)
(70)
595
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
En general:
Resolver:
Solución:
La ecuación diferencial dada tiene coeficientes constantes y haciendo
y = erx, obtenemos la ecuación característica r3 – r2 = r2(r – 1) = 0 con
raíces r = 0, 0, 1, entonces {1, x, ex} forma un conjunto fundamental
de soluciones de la ecuación diferencial homogénea y/// - y// = 0. Sea
y1 = 1, y2 = x, y3 = ex. Entonces a3(x) = 1,n = 3, y f (x) = ex. El sistema de
ecuaciones dado por (69) para esta ecuación diferencial es:
596
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
(71)
Entonces:
Y:
597
8.16 Ejercicios
En los ejercicios 1 al 4, resuelva la ecuación diferencial por el método
de variación de parámetros y dar la solución general.
1.
2.
3.
4.
5. Verifique que
6.
7.
8.
9.
10.
11.
12.
14.
15.
16.
17.
18.
19.
20.
21.
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
donde a3, a2, a1, a0 son constantes con a3 ≠ 0. Hay siempre una ecuación
de Euler de primer orden,
De donde:
601
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Es bueno para el estudiante que obtenga la ecuación de índices
las primeras veces. Basta hacer y = xm en la ecuación diferencial dada
y luego dividir por xm para obtener la ecuación de índices:
O:
602
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
O de forma equivalente:
Solución:
Haciendo y = xm se obtiene la ecuación de índices 9.m. (m – 1) + 15 m
+ 5 = 9.m2 + 6m + 5 = (3m + 1)2 + 4 = 0, cuyas raíces son:
603
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
604
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Sea y = xm, esto nos lleva a la ecuación de índices m (m – 1) – 7m
+ 16 = m2 – 8m + 16 = (m – 4)2 = 0, con lo cual m = 4 es una raíz
de multiplicidad dos. Como tenemos que una solución es y = x4, la
segunda solución para el conjunto fundamental de soluciones puede
ser encontrada por reducción del orden, sea y = v.x4.
Separando variables:
605
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Integrando, se tiene:
De donde:
Solución:
Haciendo y = xm tenemos la ecuación de índices
la cual tiene raíz
repetida a 1. Entonces:
606
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
está dada como sigue. Por sustitución y = xm, obtenemos la
ecuación de índices am (m – 1) + b. m + c = 0. Encontrar las
raíces m1, m2:
2. Si , entonces:
3. Si m1 = m2, entonces:
607
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Haciendo y = xm y calculando las tres primeras derivadas, se tiene:
608
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
Solución:
Como se trata de una ecuación diferencial no homogénea,
primero resolvemos la ecuación diferencial homogénea asociada a
la ecuación diferencial lineal dada. Resolveremos primero la ecuación
diferencial x2y// - 3xy/ + 3y = 0. De la ecuación de índice se tiene que
(m -1) (m- 3) = 0, luego yh = C1 x + C2 x3. Ahora recordemos que
para hallar yp debemos encontrar v1 y v2 tales que yp= v1y1 + v2y2 y
que , donde W1, W2 y W son los determinantes que se
definieron cuando definimos el método de variación de parámetros.
Por lo tanto, dividiendo la ecuación dada por x2 tenemos que:
609
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
610
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
En resumen
Dada la ecuación diferencial de Cauchy-Euler de la forma:
Donde a, b y c son constantes, entonces:
1. Si x > 0, transforme la variable independiente de la x
a r, donde: x = e r
De modo que:
Y:
(Si x < 0, utilice la transformación x = - er de modo que la
transformación nos lleva a r = ln (-x) y las otras ecuaciones
quedan iguales).
2. Sustituir (72) y (73) en la ecuación de Cauchy-Euler
dada, y simplificar los términos semejantes, así se tiene:
La cual es una ecuación diferencial lineal de segundo grado
con coeficientes constantes.
3. Resuelva (74) para y = y(r).
4. Utilice el cambio de variables dado inicialmente y la solu-
ción y = y(r) obtenida en el paso 3 para encontrar la solu-
ción final y = y(r).
611
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
Solución:
Cambiamos la variable independiente de x a r usando el cambio de
variable x = er. Sustituyendo (72) y (73) en la ecuación diferencial
dada, obtenemos:
Solución:
El primer paso es resolver la ecuación diferencial homogénea asociada:
612
Aplicaciones de las ecuaciones diferenciales lineales de segundo orden y de grado superior
tiene:
613
Ecuaciones diferenciales y sistemas dinámicos Jorge Augusto Pérez Alcázar
614
8.18 Ejercicios
En los ejercicios 1 al 12, resuelva la ecuación diferencial de segundo
orden de Cauchy-Euler para x > 0. Si no hay condición inicial, dé la
solución general.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
Se convierte en:
Es decir, la ecuación de Euler permanece inalterada
mediante este cambio de escala de la variable indepen-
diente.
b. Esto contrasta con lo que sucede cuando se realiza
el mismo cambio de escala en la ecuación de coeficien-
te constante a y// + b y/ + c y = 0.
En:
d. Verifique que si k > 0 y m1 es un número real constan-
te, entonces x = k s cambia:
En:
e. Verifique que, si K > 0 y α, ß son constantes reales,
ß≠0, entonces x = ks cambia:
En:
En la ecuación de segundo orden con coeficientes cons-
tantes:
17.
18.
19.
20.