Apuntesode

Mat1640: EDO
Due on Segundo semestre 2014
07/14
Hernan Gonzalez A. , higonzal@uc.cl
1
Hernan Gonzalez A. , higonzal@uc.cl Mat1640 (07/14 ): EDO
Índice
Ejercicio 1 7
Ejercicio 2 8
Ejercicio 3 10
Ejercicio 4 11
Ejercicio 5 13
Ejercicio 6 14
Ejercicio 7 15
Ejercicio 8 17
Ejercicio 9 18
Ejercicio 10 19
Ejercicio 11 20
Ejercicio 12 21
Ejercicio 13 23
Ejercicio 14 24
Ejercicio 15 25
Ejercicio 16 26
Ejercicio 17 27
Ejercicio 18 29
Ejercicio 19 30
Ejercicio 20 32
Ejercicio 21 34
Ejercicio 22 35
Ejercicio 23 37
Ejercicio 24 38
Ejercicio 25 39
Ejercicio 26 40
Ejercicio 27 41
Page 2 of 117
Ejercicio 28 41
Ejercicio 29 42
Ejercicio 30 43
Ejercicio 31 44
Ejercicio 32 47
Ejercicio 33 48
Ejercicio 34 50
Ejercicio 35 52
Ejercicio 36 55
Ejercicio 37 56
Ejercicio 38 57
Ejercicio 39 58
Ejercicio 40 59
Ejercicio 41 60
Ejercicio 42 61
Ejercicio 43 62
Ejercicio 44 62
Ejercicio 45 63
Ejercicio 46 64
Ejercicio 47 64
Ejercicio 48 65
Ejercicio 49 66
Ejercicio 50 67
Ejercicio 51 68
Ejercicio 52 69
Ejercicio 53 70
Ejercicio 54 72
Ejercicio 55 73
Page 3 of 117
Ejercicio 56 74
Ejercicio 57 75
Ejercicio 58 76
Ejercicio 59 77
Ejercicio 60 77
Ejercicio 61 78
Ejercicio 62 79
Ejercicio 63 79
Ejercicio 64 80
Ejercicio 65 81
Ejercicio 66 82
Ejercicio 67 84
Ejercicio 68 85
Ejercicio 69 87
Ejercicio 70 90
Ejercicio 71 91
Ejercicio 72 92
Ejercicio 73 93
Ejercicio 74 93
Ejercicio 75 94
Ejercicio 76 95
Ejercicio 77 96
Ejercicio 78 96
Ejercicio 79 97
Ejercicio 80 97
Ejercicio 81 98
Ejercicio 82 98
Ejercicio 83 100
Page 4 of 117
Ejercicio 84 100
Ejercicio 85 101
Ejercicio 86 101
Ejercicio 87 102
Ejercicio 88 102
Ejercicio 89 103
Ejercicio 90 104
Ejercicio 91 104
Ejercicio 92 105
Ejercicio 93 105
Ejercicio 94 106
Ejercicio 95 106
Ejercicio 96 106
Ejercicio 97 107
Ejercicio 98 107
Ejercicio 99 108
Ejercicio 100 108
Ejercicio 101 109
Ejercicio 102 110
Ejercicio 103 111
Ejercicio 104 112
Ejercicio 105 112
Ejercicio 106 113
Ejercicio 107 114
Ejercicio 108 115
Ejercicio 109 116
Ejercicio 110 116
Ejercicio 111 117
Page 5 of 117
Ejercicio 112 117
Page 6 of 117
Ejercicio 1
Resolver √
du 1 − x2 4 x arctan(x) √
x + u= √ u
dx 1 + x2 1 + x2
Solución:
√
Dividiendo por x u, suponiendo u 6= 0, x 6= 0, se obtiene:
1 du 1 − x2 1 4 arctan(x)
√ + u2 = p
u dx x(1 + x2 ) x(1 + x2 )
√ du
Sea v = u, donde v 0 = 1
√
2 u
, ası́:
dx
1 − x2 2 arctan(x)
v0 + v=p
2x(1 + x2 ) x(1 + x2 )
Calculando el factor integrante:
1 − x2
Z Z
1 x
µ(x) = exp dx = exp − dx
2x(1 + x2 ) 2x 1 + x2
( 12 )!
1 1 x
= exp ln(x) − ln(1 + x2 ) = exp ln
2 2 1 + x2
r
x
=
1 + x2
Multiplicando la ecuación diferencial por µ(x)

r
d x 2 arctan(x)
v =
dx 1 + x2 1 + x2
Integrando se tiene
r
x
v = arctan2 (x) + C
1 + x2
De donde sustituyendo v,
√
r
x
u = arctan2 (x) + C
1 + x2
Se obtiene finalmente:
1 + x2 2
u(x) = arctan2 (x) + C x 6= 0, C ∈ R
x
Page 7 of 117
Hernan Gonzalez A. , higonzal@uc.cl Mat1640 (07/14 ): EDO Ejercicio 1
Ejercicio 2
En un dı́a de Agosto con 20o C temperatura ambiente, el Inombrable compra en el QuickDeli empanadas
napolitanas antes de la ayudantı́a de EDO que iba a impartir al módulo 6, a una temperatura de 100o C y
se quema (se lo merece). Luego de 2 minutos quiso comerlas pero todavı́a estaban muy calientes (80o C),y se
quema nuevamente. ¿Cual es la temperatura luego de 5 minutos? Si se supone que la temperatura adecuada
para comerlas es 40o C ¿Cuanto tiempo será necesario esperar?
Solución:
dT
= −k (T (t) − Tm )
dt
Para hallar la función T(t), se resuelve la ecuación diferencial. Primero separamos las variables:
Z Z
dT dT
= −kdt ⇒ =− kdt ⇒ ln (T (t) − Tm ) = −kt + C
T (t) − Tm T (t) − Tm
Aplicamos exponencial y acomodamos un poco:
T (t) − Tm = e−kt+C ⇒ T (t) = e−kt · eC + Tm
Como la temperatura del ambiente es 20 o C, entonces la función queda definida por:
T (t) = e−kt · eC + 20
Ahora vemos las condiciones iniciales. En primer lugar, la temperatura inicial de la empanada es de
100o C, esto es:
T (0) = 100 ⇔ 100 = e0 · eC + 20 ⇔ 80 = eC
Por lo que la función ahora nos queda:
T (t) = 80e−kt + 20
Ejercicio 2 continued on next page. . . Page 8 of 117

Hernan Gonzalez A. , higonzal@uc.cl Mat1640 (07/14 ): EDO Ejercicio 2 (continued)
Otra condición planteada en el problema es que a los dos minutos la temperatura de la empanada es de
80o C, es decir:

−2k 3 1 3
T (2) = 80 ⇔ 80 = 80e + 20 ⇔ = e−2k ⇔ k = − ln
4 2 4
Luego:
t
2t
t
ln( 34 ) ln( 43 ) 2 3
T (t) = 80e 2 + 20 ⇔ T (t) = 80e + 20 ⇔ T (t) = 80 + 20
4
Finalmente:
2t
3
T (t) = 80 + 20
4
Resolviendo el problema, la temperatura a los 5 minutos queda definida al evaluar la función en 5:
52 √ √ !
3 9 3 45 3
T (5) = 80 + 20 = 80 · · + 20 = + 20 [◦ C] ≈ 58, 25 [◦ C]
4 16 2 2
Para que la empanada se encuentre a 40o C, se deben esperar:
2t 2t
3 1 3 1 t 3
40 = 80 + 20 ⇒ 4 = ⇒ ln = ln ⇒
4 4 4 2 4
ln (1) − ln (4) 4 ln (2)

t=2· ⇒t= [mı́n] ≈ 9, 63 [mı́n]
ln (3) − ln (4) 2 ln (2) − ln (3)
Page 9 of 117
Ejercicio 3
Un tanque contiene, en un principio, 40[L] de salmuera (sal muy concentrada en agua). Entra al tanque,
1
a una razón de 8[L/min] , salmuera conteniendo + α de sal por litro. La mezcla, bien agitada, sale del
8
23
tanque a razón de 4[L/min] . A los 10[min] la cantidad de sal en el tanque es de [Kg] .
2
6
Calcule α de modo que a los 20[min] la concentración de sal en el tanque sea de [Kg] de sal por litro.
25
Solución:
Sea A la cantidad de sal, y V el volumen, ve el volumen de entrada, y vs volumen de salida
V0 = 40
1
ve = 8, re = +α
8
A
vs = 4, rs =
V
V = 40 + 8t − 4t = 40 + 4t
la ecuacion diferencial queda:
dA
= ve ∗ re − vs ∗ rs
dt
dA 4A 1
+ = 8 ∗ ( + α)
dt 40 + 4t 8
resolviendo
1
A= (1 + 8α)(t + 10) + C(t + 10)−1
2
tenemos dos condiciones
1 23
A(10) = (1 + 8α)(20) + C(20)−1 =
2 2
A(20) 1 (1 + 8α)(30) + C(30)−1 6
C(20) = = =
V (20) 2 120 25
resolviendo el sistema obtenemos que
24
α=
125
Page 10 of 117
Ejercicio 4
Resolver la siguiente ecuación diferencial
6xy dx + (4y + 9x2 ) dy = 0
Solución:
Se elige M = 6xy y N = 4y + 9x2 .Notemos que no es Difrencial exacta,

∂M ∂N
6=
∂y ∂x
Se calula
∂M ∂N
− = 6x − 18x = −12x
∂y ∂x
Luego
1 ∂M ∂N 12x
( − )=−
N ∂y ∂x 4y + 9x2
Notar que no nos ayuda mucho, asi que se calcula entonces
1 ∂N ∂M 1 2
( − )= (12x) =
M ∂x ∂y 6xy y
R 2dy
.Por lo tanto el factor adecuado es µ(y) = e = e2ln(y) = y 2 .Multiplico la EDO original por y 2
y
quedando:
6xy 3 dx + (4y 3 + 9y 2 x2 )dy = 0
. Luego esta EDO es exacta y se procede a lo metodico: Existe una Función F tal que
∂F ∂F
= 6xy 3 ∧ = (4y 3 + 9y 2 x2 )
∂x ∂y
Buscamos F por el lado de las x:
Z x
F = 6xy 3 = 3x2 y 3 + C(y)
con C(y), çonstante”que depende de y Buscamos C(y):
∂F
= 4y 3 + 9y 2 x2
∂y
9y 2 x2 + C 0 (y) = 4y 3 + 9y 2 x2 ), C 0 (y) = 4y 3
Por lo tanto, C(y) = y 4 , por lo que finalmente F = 3x2 y 3 + y 4 + k con k constante

Para transformar la ecuación en una diferencial exacta
M (x, y)dx + N (x, y)dy = 0
Utilizamos un factor integrante tal que depende de x, una función µ(x) > 0 (queda como ejercicio
encontrar un factor integrante que dependa de y, igual al ejercicio anterior)
µ(x)M (x, y)dx + µ(x)N (x, y)dy = 0
Para que sea exacta debe cumplirse que
∂(µ(x)M (x, y)) ∂(µ(x)N (x, y))

=
∂y ∂x
Notemos que si solo depende de x, la diferencial parcial respecto a y queda
∂µ(x)M (x, y) ∂M (x, y)

= µ(x)
∂y ∂y
Y por regla de la cadena, la derivada parcial respecto a x
∂µ(x)N (x, y ∂N (x, y) ∂µ(x)

= µ(x) + N (x, y)
∂x ∂x ∂x
Reemplazando en la condición para que sea exacta
∂M (x, y) ∂N (x, y) ∂µ(x)

µ(x) = µ(x) + N (x, y)
∂y ∂x ∂x
Reordenando

∂M (x, y) ∂N (x, y) ∂µ(x)
µ(x) − = N (x, y)
∂y ∂x ∂x
dµ(x) ∂M (x, y) ∂N (x, y)

−
dx = ∂y ∂x
µ(x) N (x, y)
dµ(x)
d ln µ(x)
Notamos que es una EDO de variables separables, ya que = dx
dx µ(x)
∂M (x, y) ∂N (x, y)
Z −
∂y ∂x
ln µ(x) = dx
N (x, y)
Finalmente, aplicando exponencial, el factor integrante es de la forma:
∂M (x,y) ∂N (x,y)
!
Z
∂y − ∂x
µ(x) = exp dx
N (x, y)

Ejercicio 5
y 3 dx + 2(x2 − xy 2 ) dy = 0
Solución:
La idea es mostrar que pueden hacerla por diferenciales exactas, o por separables.
La ecuación es equivalente a
y3
y0 = − .
2x2 − 2xy 2
Reemplazar y 2 = xz, y la EDO queda

√ 0 √
√

xz 1 z z
√ =− √ + z ,
z 2 x 1−z
que es de variables separables.

Para que sean este tipo de ecuaciones sean de variables separables debe cumplirse si y solo si:
1)f (tx, ty) = tm f (x, y), g(tx, ty) = tn g(x, y), y n = m. cuando las ecuaciones diferenciales de este tipo
cumplen eso, se dice que son homogéneas de grado n = m., y su sustitución es y = vx.
2) y 0 = G(ax + by), donde la sustitución es z = ax + by
3)f (x, y) = a1 x + b1 y + c1 y g(x, y) = a2 x + b2 y + c2 , en este caso se dice que la EDO, es de coeficientes
lineales. Por otro lado, si a1 b2 = a2 b1 , entonces se tiene el caso anterior, y si c1 = c2 = 0, se tiene el caso
1.
Page 13 of 117
Ejercicio 6

y 0 − y = e2x y 3
Solución:
Que corresponde a un bernoulli, en este tipo, divides todo por y 3 , quedando
y −3 y 0 − y −2 = e2x
Realizando la sustitución u = −y −2 , u0 = 2y −3 y 0 y normalizando:
u0 + 2u = 2e2x
Esto es:
1
u(x) = c1 e−2x + e2x
2
Volviendo a la variable original:
1
−y −2 = c1 e−2x + e2x
2
Page 14 of 117
Ejercicio 7
Resolver las siguientes ecuaciones
dy y dy 1
=1− , =1+ 2
dt 3 dt y
Solución:
dy y
La ecuación = 1 − es de variables separables, por lo cual se puede resolver ası́:
dt 3
dy y
= 1−
dt 3
dy 3−y
=
dt 3
dy 1
= dt
3−y 3
Integramos a ambos lados de la igualdad:
Z Z
dy 1
= dt
3−y 3
1
− ln |3 − y| = t + C1
3
1
ln |3 − y| = − t + C2 (dondeC2 = −C1 )
3
Aplicamos función exponencial a ambos lados de la igualdad para ası́ obtener:¡br /¿

t t
|3 − y| = e− 3 +C2 = e− 3 · eC2
eC2 es constante, ası́ que podemos hacer C = eC2 y luego:

t
|3 − y| = Ce− 3
t
Al quitar las barras de valor absoluto, se tienen ”dos” soluciones: 3 − y = Ce− 3 y
t
3 − y = −Ce− 3 .Entonces podemos representar las constantes por una sola, llamémosleK, con lo cual
(despejando y) la solución general de la EDO es:
t
y = 3 − Ke− 3

dy 1 dy y2 + 1
=1+ 2 ⇔ =
dt y dt y2
y2
Z
dy = dt/
y2 + 1
y2
Z Z
dy = dt
y2 + 1
y2
Z
dy = t + k
y2
+1
Z 2
y +1 1
2
dy − 2 dy = t + k
y +1 y +1
Z Z
1
dy − dy = t + k
y2 + 1
| {z }
A
donde A es una integral bien conocida, pero se saca haciendo la sustitucion

y = tan θ
entonces :
y − arctan (y) + K = t + k
Page 16 of 117
Ejercicio 8
Una curva de ecuación cartesiana y = f (x)pasa por el origen.Por un punto arbitrario de la curva se trazan
rectas paralelas a los ejes coordenados que forman un rectangulo con ellos.La curva divide tal rectangulo en
dos regiones A y B una de la cuales tiene un area n veces la de la otra.Hallar la funcion f
Solución:
Se puede decir que el área A debajo de la curva es n veces menor a la otra:

Z x Z x
nA = n f (r)dr = xf (x) − f (r)dr
0 0
Derivando la ultima igualdad:
nf (x) = f (x) + xf 0 (x) − f (x)
Finalmente queda una EDO de variables separables:
f 0 (x) n
=
f (x) x
Integrando y aplicando exponencial
f (x) = kxn
Page 17 of 117
Ejercicio 9
Resolver el problema de Cauchy
y 0 + 2y = 1
5
y(0) =
2
Solución:
Se puede resolver con método de variables separables, pero veamos 2 metodos diferentes:
La E.D.O. homogénea asociada es y 0 + 2y = 0, que tiene por solución general yh (t) = ke−2t . Se tiene ası́,
que la solución general de la E.D.O. no homogénea viene dada por y(t) = yh (t) + yp (t), donde yp (t) es
una solución particular de la ecuación. No es difı́cil notar (el ”método”que se uso para hallar soluciones
particulares se llama çoeficientes indeterminados”) que yp (t) = 21 es solución particular de la ecuación
y 0 + 2y = 1. Se tiene ası́, que y(t) = ke−2t + 12 es solución general de la E.D.O. no homogénea. Imponiendo
la condición inicial se obtiene el valor de k para que dicha función satisfaga el problema de Cauchy.
Pueden ocupar el otro método
R
que funciona para las EDOs lineales de primer orden también, que
consiste en multiplicar por e p(x)dx donde p(x) es el término que acompaña a y en la ecuación (siempre y
cuando a y 0 lo acompaña un 1).
La EDO queda
(ye2x )0 = e2x y ahı́ ya es cosa de integrar y despejar.
Page 18 of 117
Ejercicio 10
Determinar la solucion de la ecuacion diferencial
(x + 2y)dx + (2x − 3y)dy = 0
Solución:
La puedes escribir como

x + 2y
y 0 (x) = = f (x, y)
3y − 2x
Y se cumple que
f (λx, λy) = f (x, y)
Por lo tanto con el cambio de variable

y
z=
x
La transformas en una edo de variables separables Quedando como
1 + 2z
z + xz 0 =
3z − 2
Y notemos que es exacta (creo que mas bien, notemos que es homogénea, todos vieron que era exacta)
M (x, y) = x + 2y N (x, y) = 2x − 3y
∂M ∂N
⇒ =
∂y ∂x
Page 19 of 117
Ejercicio 11
y 00 + tan(x)y 0 + cos2 (x)y = sin(x)cos2 (x)
Solución:
dt dt p
Sea x = arcsin(t). Entonces dx = √ , es decir, = 1 − t2 .¡Por otro lado,
1 − t2 dx
dy dy dt p dy
= · = 1 − t2
dx dt dx dt
y
d2 y

d dy d p 2
dy d p 2
dy dt
= = 1−t = 1−t ·
dx2 dx dx dx dt dt dt dx
d2 y p

−t dy p
= √ + 1−t 2 1 − t2
1 − t2 dt dt2
dy d2 y
= −t + (1 − t2 ) 2
dt dt
t √
Como t = sin(x), entonces tan(x) = √ y cos(x) = 1 − t2 . Reemplazando, y usando la notación
1 − t2
y = y(t),
t p
−ty 0 + (1 − t2 )y 00 + √ 1 − t2 y 0 + (1 − t2 )y = t(1 − t2 )
1 − t2
y 00 + y = t
Se resuelve la homogenea y 00 + y = 0, que es de coeficientes constantes. La ecuación caracterı́stica
asociada es k 2 + 1 = 0 ⇔ k = ±i. Entonces, la solución general de la homogenea es¡br /¿
yh (t) = C1 cos(t) + C2 sin(t)
Se encuentra una solución particular de la no homogenea por el método de los coeficientes

indeterminados. Dado que g(t) = t (lado derecho de la ecuación), entonces la forma de la solución
0 00
particular de la no homogenea,ynh (t), es ynh (t) = A + Bt. Ası́, ynh (t) = B, ynh (t) = 0. Reemplazando,
A + Bt = t ⇔ A = 0, B = 1. Entonces, ynh (t) = t.Finalmente,
y(t) = yh (t) + ynh (t)

y(t) = C1 cos(t) + C2 sin(t) + t
y como t = sin(x), la solución general de la no homogenea es
y(x) = C1 cos(sin(x)) + C2 sin(sin(x)) + sin(x)
Page 20 of 117
Ejercicio 12
Resuelva las siguientes 3 ecuaciones. Luego deduzca el método de factor integrante
dy dy dy
= cos x, x + y = cos x, + y = 1 + x.
dx dx dx
Respecto al primero, con variables separables
Z Z
dy
dx = cos xdx
dx
y(x) = sin x + C, donde C ∈ R.
El segundo, notemos que por regla de la cadena queda
dy d(xy)
+y =
dx dx
Volviendo al problema original y realizando variables separables
d(xy)
= cos x
dx
Z Z
dxy
xy = dx = cos xdx = sin x + C
dx
1
y(x) = (sin x + C) .
x
Para la tercera, notemos que al multiplicar por ex ( Siempre > 0, por lo que no altera la igualdad)
dy
ex + ex y = ex (1 + x)
dx
dy d(ex y)
Tenemos nuevamente por regla de la cadena ex + ex y =
dx dx
d(ex y) dy
= ex + ex y = ex (x + 1)
dx dx
Resolvemos la ecuacion de variables separables, mas integración por partes (Con el gentil auspicio de
Cálculo I, destruyendo cerebros desde que Newton se aburrio un dı́a. Lo veran hasta en la sopa a futuro)
dex y
Z Z Z
x x x
e y= dx = e (x + 1)dx = e (x + 1) − ex dx = ex (x + 1) − ex + C = xex + C.
dx
Despejando y de la ecuación, queda finalmente y(x) = e−x (xex + C) = x + Ce−x

Como se ve, la idea es formar una regla de la cadena de alguna forma, para llevarla a variables
separables. Es decir, sea la EDO de la forma
dy
+ p(x)y = q(x), (1)
dx
Multiplicando por algún factor, que depende de x , µ(x) siempre mayor que cero.
dy
µ(x) + µ(x)p(x)y = µ(x)q(x)c(2)
dx
Siguiendo la idea anterior, supongamos que se cumple que: (pensar en condiciones para que se cumpla)
dµ(x)y dy dµ
= µ(x) +y (3)
dx dx dx
Sin embargo, queremos que la parte de la izquierda de la ecuación (2) cumpla con la condición (3)
dy d(µ(x)y)
µ(x) + µ(x)p(x)y = . (4)
dx dx
Reemplazando (3) en (4)
dy dy dµ(x)
µ(x) + µ(x)p(x)y = µ(x) +y
dx dx dx
Simplificando, nos queda una EDO de variables separables
dµ(x)
µ(x)p(x)y = y
dx
Suponiendo que y > 0
dµ
= µ(x)p(x)
dx
Resolviendo (lo dejo como tarea para el estimado lector)
R
p(x)dx
µ(x) = e
Donde µ(x), es el famoso y popular factor integrante.
Page 22 of 117
Ejercicio 13
Resolver EDO (todos los ejercicios del curso son asi)
3x(xy − 2)dx + (x3 + 2y)dy = 0
Solución:
Para una ecuación exacta, debe cumplirse que
∂M ∂N
M (x, y)dx + N (x, y)dy = 0, = .
∂y ∂x
Y significa que la ecuación proviene f, tal que df = M (x, y)dx + N (x, y)dy, que además
∂f
df = + ∂f
∂x dx ∂y dy. (esto se conoce como campo conservativo, donde f es el potencial, será visto en
cálculo III. Aprovecho el corte comercial para el aviso publicitario: Jueves Módulo 4 para que vayan a la
otra ayudatı́a que realizo del cálculo III. Si, tienen premios, juegos, ayudantı́as resueltas y comentarios
rı́diculos como estos en las resoluciones)
M (x, y) = 3x(xy − 2), N (x, y) = x3 + 2y.
Notamos que se cumplo la condición de ecuación exacta
∂M ∂N
= 3x2 = .
∂y ∂x
Es decir, viene de un potencial f tal que
∂f ∂f
df = dx + dy.
∂x ∂y
∂f ∂f
= 3x(xy − 2), = x3 + 2y.
∂x ∂y
Integrando para encontrar el potencial
Z Z Z
∂f
f (x, y) = dx = 3x(xy − 2)dx = 3x2 y − 6xdx = x3 y − 3x2 + h(y)
∂x
Ahora utilizando la derivada parcial respecto a y
∂f ∂x3 y − 3x2 + h(y)

x3 + 2y = = = x3 + h0 (y)
∂y ∂y
Por ende h0 (y) = 2y → h(y) = y 2 + C

Finalmente
f (x, y) = x3 y − 3x2 + h(y) = x3 y − 3x2 + y 2 + C
Page 23 of 117
Ejercicio 14
e3x y 3 − 2y
Sea la EDO y 0 (x) = (x − y). Resuelva la EDO con condición inicial y(0) = 1.
y−x
Solución:
En este caso, la EDO ya viene ”normalizada”, por lo tanto, sólo debemos darnos cuenta que es de
Bernoulli. Veamos que para R2 − {y = x} la EDO tiene solución. Ası́, se tiene y 0 (x) = 2y − e3x y 3 , y es
claro observar la forma y 0 (x) − 2y = −e3x y 3 , vale decir, es una EDO de Bernoulli con n = 3. Hagamos el
cambio definido por u = y 1−n , en este caso, para n = 3, se tiene ası́, u = y −2 . Por otra parte,
u = −2y −3 · y 0 (x), este resultado nos dice que debemos multiplicar la EDO de Bernoulli por −2y −3 , se
0
sigue que:
y 0 (x) − 2y = −e3x y 3 ⇐⇒ −2y 3 · y 0 (x) + 4y −2 = 2e3x ⇐⇒ u0 (x) + 4u = 2e3x .
Notemos que esta última

R
EDO es lineal, con p(x)R
= 4 y q(x) = 2e3x , luego, considerando el Factor de
Integración µ(x) = e p(x)dx , se tiene que µ(x) = e 4dx = e4x . Luego, multiplicando la EDO Lineal por µ,
se tiene:
u0 (x) + 4u = 2e3x
e4x · u0 (x) + 4e4x u = 2e7x
(u · e4x )0 = 2e7x
Z Z
(ue4x )0 dx = 2e7x dx
2 7x
ue4x + C1 = e + C2
7
2 3x
u = e + e−4x C
7
2
Remember that u = y −2 , ası́, y −2 = e3x + e−4x C. Despejando y, se sigue que la solución es
r 7
7
y(x) = ± . Ahora bien, como y(0) = 1, sabemos que el punto (0, 1) pertenece a una
2e3x + 14e−4x C
región positiva, ya sea el 1er. o 2do cuadrante, de ahı́ que y > 0, por tanto, la solución de la EDO o del
PVI es y(x) = 1.
Page 24 of 117
Ejercicio 15
Resolver
(3xy + y 2 )dx + (x2 + xy)dy = 0
Solución:
Notemos que no es exacta, por ende, debemos encontrar un factor integrante.
µ(x)(3xy + y 2 )dx + µ(x)(x2 + xy)dy = 0
Para que sea exacta queremos que:
∂µ(x)(3xy + y 2 ) ∂µ(x)(x2 + xy)

=
∂y ∂x
µ(x)(3x + 2y) = µ0 (x)(x2 + xy) + µ(x)(2x + y)

µ0 (x)(x2 + xy) = µ(x)(3x + 2y) − µ(x)(2x + y)
µ0 (x)(x2 + xy) = µ(x)(x + y)
µ0 (x)(x(x + y)) = µ(x)(x + y)
xµ0 (x) = µ(x)
µ0 (x) 1
=
µ(x) x
d ln µ 1
=
dx x
ln µ(x) = ln x
µ(x) = x
x(3xy + y 2 )dx + x(x2 + xy)dy = 0
∂f ∂f
df = (3x2 y + xy 2 )dx + (x3 + x2 y)dy, df = dx + dy
∂x ∂y
∂f ∂f
= 3x2 y + xy 2 y = x3 + x2 y
∂x ∂y
Z Z
∂f 1
f (x, y) = dx = 3x2 y + xy 2 dx + h(y) = x3 y + x2 y 2 + h(y)
∂x 2
∂f ∂x3 y + 21 x2 y 2 + h(y)
x3 + x2 y = = = x3 + x2 y + h0 (y)
∂y ∂y
h0 (y) = 0, h(y) = k, k = constante
1
f (x, y) = x3 y + x2 y 2 = C
2
Page 25 of 117
Ejercicio 16
Resolver y 00 − xy 0 + y = −x cos x, dada las condiciones iniciales y(0) = 0, y 0 (0) = 2
Todo lo que necesitamos saber son las condiciones Iniciales de la ecuación para el mismo punto, x=0 en
nuestro caso, al menos para y 0 y y 00 y el Teorema Taylor, es decir saber que al rededor de un punto
dado(o sea, en un entorno del punto) x=a una función y=f(x) admite ser expresada como solucion en
serie de potencias:
∞
X
y(x) = an (x − a)n
n=0
ahora, operemos un poco con la EDO dada y la condiciones de la misma:
y 00 − xy 0 + y = −x cos x
y 00 (x) = xy 0 (x) − y(x) − x cos x.....(1)y 0 (0) = 2, y(0) = 0.....(c)

hacemos x=0
y 00 (0) = (0), y 0 (0) − y(0) − (0) cos(0) = 0.....(2)
derivando a ambos lados de (1) respecto a x
y 000 (x) = xy 00 (x) − cos x + x sin x
y (4) (x) = y 00 (x) + xy 000 (x) + 2 sin x + x cos x

y (5) (x) = 2y 000 (x) + xy (4) (x) + 3 cos x − x sin x
y (6) (x) = 3y (4) (x) + xy (5) (x) − 4 sin x − x cos x
y (7) (x) = 4y (5) (x) + xy (6) (x) − 5 cos x + x sin x
y reemplazando las condiciones (c) en (1) y sucesivamente observamos que
y (0) (0) = 0, y (1) (0) = 2, y (2) (0) = 0, y (3) (0) = −1
y (4) (0) = 0, y (5) (0) = 1, y (6) (0) = 0, y (7) (0) = −1

Si suponemos que la solución es de la forma :
∞
X y k (0)
y(x) = xk
k!
k=0
notamos que este patron obtenido obecede a la serie taylor:
y (0) (0)x0 y (1) (0)x0 y (k) (0)xk

+ + ... + + ......
0! 1! k!
Reordenando:
x3 x5 x7 (−1)k x2k+1
x + (x − + − + ... )
3! 5! 7! (2k + 1)!
la parte entre paréntesis corresponde a la expanción en serie de potencias Mclaurin entorno a cero, ası́
finalmente: y(x) = x + sin x
Page 26 of 117
Ejercicio 17
1
Resolver conciendo la solución particular es y1 =
x
2
y0 = y2 − , y(1) = 15
x2
Solución:
1
Efectuamos el cambio de variable y = v + , entonces:
x
dv 2v dv
= v2 + ⇐⇒ − 2x−1 v = v 2
dx x dx
Se tiene una ecuación de Bernoulli, luego efectuando el cambio definido por u = v 1−2 = v −1 y
diferenciando, tendremos:
du 1 dv
= −
dx v 2 dx
dv du
= −u−2
dx dx
du
−u−2 − 2u−1 x−1 = u−2
dx
du
+ 2x−1 u = −1
dx
Ahora se tiene una ecuación diferencial lineal de 1er. orden de la forma y 0 + P (x)y = f (x), donde
P (x) = 2x−1 , luego el factor integrante es:
Z
2
µ(x) = exp 2x−1 = eln x = x2
Amplificando la ecuación por este factor y volviendo a nuestras variables originales se tendrá:
x2 u0 + 2xu = −x2
(x2 u)0= −x2
x3
x2 u = − + c
3
3cx−2 − x
u =
3
3
v = −2
3cx − x
3 1
y = +
3cx−2 − x x

Finalmente imponemos la condición dada y(1) = 15, y tendrı́amos:

3
15 = +1
3c − 1
3
= 14
3c − 1
3
3c − 1 =
14
17
3c =
14
17
∴ c =
42
Otra forma
1 1
y= +
u x
y nos olvidamos de Bernoulli
1 1 dy 1 du 1
y= + =⇒ =− 2 − 2
u x dx u dx x
Luego al reemplazar nos queda
1 du 1 1 2 1 2
− − = + + − 2
u2 dx x2 u2 ux x2 x
1 du 1 2
− 2 = +
u dx u2 ux
du
=⇒ = −1 − 2ux−1
dx
du
∴ + 2x−1 u = −1
dx
Esta es la conocida ecuación de Ricatti
dy
+ p(x)y + q(x)y 2 = f (x)
dx
Esta ecuación se resuelve si previamente se conoce una solución particular, sea y1 (x) Conocida dicha
solución, encontrada por inspección, se hace el cambio:
y(x) = z(x) + y1 (x)
Page 28 of 117
Ejercicio 18
Resolver
dy
ln (y x ) = 3x2 y
dx
Solución:
Operando, separando variables e integrando ambos miembros de la ecuación cuando sea necesario,
tendremos:
dy
ln(y x ) = 3x2 y
dx
dy
ln(y) = 3xy
dx
ln y
dy = 3xdx
y
Z Z
ln y
dy = 3xdx
y
ln2 y 3 2
= x + c1
2 2
ln2 y = 3x2 + 2c1
p
ln y = 3x2 + 2c1
√
3x2 +c
y = e
Page 29 of 117
Ejercicio 19
Para x > −1 encuentre la solución general implı́cita de la ecuación
x2 + y 2 + 1 dx − (xy + y) dy = 0

−n
Sabiendo que tiene factor integrante de la forma µ(x, y) = (x + 1)
Solución:
M (x, y) = x2 + y 2 + 1 N (x, y) = −(xy + y)
utilizando el factor integrante
µM dx − µN dy = 0
2
(x + 1)−n (x2 + y + 1) dx + −(x + 1)−n (xy + y) dy = 0
| {z } | {z }
M∗ N∗
Para que la ecuacion sea exacta, se necesita:
∂M ∗ ∂N ∗
(x, y) = (x, y)
∂y ∂x
Entonces
∂M
(x, y) = 2y(x + 1)−n
∂y
∂N ∗
= n(x + 1)−n−1 (xy + y) − (x + 1)−n y
∂x
∂N ∗ ny(x + 1)
= (x + 1)−n ( − y)
∂x (x + 1)
∂N ∗
= (x + 1)−n y(n − 1)
∂x
Se debe cumplir:
(n − 1)(x + 1)−n y = 2y(x + 1)−n

n−1=2
n=3
µ(x, y) = (x + 1)−3
se obtiene
µ(x, y) = (x + 1)−3

2
(x + 1)−3 (x2 + y + 1) dx + −(x + 1)−3 (xy + y) dy = 0
| {z } | {z }
M∗ N∗
Entonces existe g(x,y) tal que:
∂g(x, y) ∂g(x, y)
= M ∗ (x, y) y = N ∗ (x, y)
∂x ∂y
Ası́:
∂g(x, y)
= −(x + 1)−3 (xy + y) /Integrando
∂y
Z Z
∂g(x, y) xy y
dy = − + dy
∂y (x + 1)3 (x + 1)3
xy 2 y2
g(x, y) = − − + c(x)
2(x + 1)3 2(x + 1)3
y2
g(x, y) = − + c(x)
2(x + 1)2
Por otra parte:
∂g(x, y) 2
= M ∗ (x, y) = (x + 1)−3 (x2 + y + 1)
∂x
Entonces:
∂g(x, y) y2
= + c‘(x)
∂x (x + 1)3
R (x2 + 1)2 − 2x
Igualando e Integrando: c(x) = dx
(x + 1)3
c(x) = ln(x + 1) − 2I
−x 1 2x + 1
I= − =−
2(x + 1)2 2(x + 1)
2(x + 1)
2
2x + 1
c(x) = ln(x + 1) +
(x + 1)2
y2
2x+1

∴ g(x, y) = − + ln(x + 1) + (x+1)2
2(x + 1)2
Las soluciones estan definidas implicitamente por:
g(x, y) = C C ∈ R
Finalmente se obtiene :
1
y(x) = (2(x + 1)2 ln(x + 1) + 2(2x + 1) − C) 2
Page 31 of 117
Ejercicio 20
Calcular y(x) que cumpla:(3x + 5y)dx + (4x + 6y)dy = 0
M (x, y) = 3x + 5y, N (x, y) = 4x + 6y
dy −M (x, y)
=
dx N (x, y)
dy − (3x + 5y)
=
dx (4x + 5y)
−x 3 + 5 xy

dy
=
x 4 + 6 xy

dx
− 3 + 5 xy

dy
=
dx 4 + 6 xy
Hacemos el cambio de variable, y resolvemos por ser función homogénea de grado cero
y
z= ⇒ y = zx
x
dy dz
⇒ =x +z
dx dx
dz − (3 + 5z)
x +z =
dx 4 + 6z
dz − (3 + 5z)
x = −z
dx 4 + 6z
dz −3 − 5z − 4z − 6z 2
x =
dx 4 + 6z
dz −3 − 9z − 6z 2
x =
dx 4 + 6z
4 + 6z dx
dz =
−3 − 9z − 6z 2 x
2 (2 + 3z) dx
2
dz =
−3 (2z + 3z + 1) x
3z + 2 3 dx
dz = −
2z 2 + 3z + 1 2 x
3z + 2 3 dx
dz = −
(2z + 1) (z + 1) 2 x

Hemos llegado a una EDO (1) de Variables Separadas, Ahora Integramos:
1 1 3 dx
dz + dz = −
(2z + 1) (z + 1) 2 x
Z Z Z
1 1 3 dx
dz + dz = − + ln (c)
(2z + 1) (z + 1) 2 x
1 3
ln |2z + 1| + ln |z + 1| = − ln |x| + ln (c) / · 2
2 2
ln |2z + 1| + 2 ln |z + 1| = −3 ln |x| + ln (c)
c
2
ln |2z + 1| + ln (z + 1) = ln 3

x
c
2
ln (2z + 1) (z + 1) = ln 3

x
c
ln 2z 3 + 5z 2 + 4z + 1 = ln 3

x
c
⇒ 2z 3 + 5z 2 + 4z + 1 =
x3
Volvamos a nuestras variables originales:
y3 3 y2 y
⇒2 3
x + 5 2 x3 + 4 x3 + x3 = c
x x x
⇒ 2y 3 + 5y 2 x + 4yx2 + x3 = c
2
⇒ (x + y) (x + 2y) = c
Se deja como propuesto realizarlo por diferenciales exactas.
Page 33 of 117
Ejercicio 21
Resolver
(8y + 10x)dx + (5y + 7x)dy = 0
Solución:
Parecido al anterior
(8y + 10x)dx + (5y + 7x)dy = 0

(5y + 7x)dy = −(8y + 10x)dx
dy 8y + 10x
= −
dx 5y + 7x
8y
+ 10
= − x
5y
+7
x
dy dv
Sea y = vx =⇒ = v + x , luego:
dx dx
dv 8v + 10
v+x = −
dx 5v + 7
dv 8v + 10
−x = +v
dx 5v + 7
dv 8v + 10 + 5v 2 + 7v
−x =
dx 5v + 7
1 dv v 2 + 3v + 2
− x =
5 dx 5v + 7
5v + 7 5
dv = − dx
(v + 1)(v + 2) x
Z Z
5v + 7 1
dv = −5 dx
(v + 1)(v + 2) x
−5
2 ln |v + 1| + 3 ln |v + 2| = ln x + ln |c|
2 3
ln |v + 1| + ln |v + 2| = ln x−5 + ln |c|

ln (v + 1)2 (v + 2)3 = ln cx−5

(v + 1)2 (v + 2)3 = cx−5

y 2 y 3
+1 +2 = cx−5
x x
(y + x)2 (y + 2x)3 = c
Page 34 of 117
Ejercicio 22
Hallar la solución general a la ecuación diferencial:
dy
cos x · sen x + y sen x = 1
dx
Solución:
Ya, esta es ”peluda”
dy 1
cos (x) · sen (x) + y · sen (x) = 1/
dx cos (x) sen (x)
dy
+ y sec (x) = sec (x) cos ec (x)
dx
R
e sec(x) = eln|sec(x)+tan(x)| = sec (x) + tan (x)
Multiplicando por el factor integrante
dy sec (x) + tan (x)

(sec (x) + tan (x)) + y sec (x) (sec (x) + tan (x)) =
dx cos (x) s en (x)
d sec (x) + tan (x)

((sec (x) + tan (x)) y) =
dx cos (x) s en (x)
Z
sec (x) + tan (x)
(sec (x) + tan (x)) y = dx
cos (x) s en (x)
| {z }
I
Z
sec (x) + tan (x)
I= dx
cos (x) s en (x)
Z
1
I= dx
sen (x) − sen2 (x)
Sustitución (invocando los poderes de MAT1610):

x
t = tan
2
2t
⇒ sen (x) =
1 + t2
2
dx = dt
1 + t2

Reemplazando:
Z 2
1+t2
I= 2 dt
2t 2t
1+t2 − 1+t2
Z
dt
I=2 4t
2t − 1+t 2
1 + t2
Z
I=2 2 dt
2t (t − 1)
Z !
1 1
I=2 + dt
2t (t − 1)2
2
I = ln (t) − +C
t−1
x 2
I = ln tan − x
+C
2 tan 2 − 1
Volviendo a la EDO:
x 2
(sec (x) + tan (x)) y = ln tan − +C
tan x2 − 1

2
x 2

ln tan 2 − tan( x
+C
2 )−1
y=
(sec (x) + tan (x))
Page 36 of 117
Ejercicio 23
Resuelva la ecuación de Bernoulli:
x2 y 0 + 2xy = 5y 3
Solución:
x2 y’ + 2xy = 5y3
2 5
⇔ y’ + y = 2 y3 , u = y1 - 3 = y - 2
x x
4 1 R 1
⇔ u’ - u = - 10 2 , F = e−4 x dx = x - 4
x x
⇔ x - 4 u’ - 4x - 5 u = - 10x - 6

d x - 4u
⇔ = - 10x - 6
dx Z
⇔ x - 4 u = - 10 x - 6 dx
x-5
⇔ x - 4 u = - 10 +c
-5
2
⇔ u= + cx4
x
2 + cx5
⇔ y-2 =
x
2 x
⇔ y =
2 + cx5
r
x
∴ |y| =
2 + cx5
Page 37 of 117
Ejercicio 24
Resuelva la ecuación de Bernoulli:
y0 = y + y3
Solución:
Se acomoda la ecuación quedando:
y0 1
3
− 2 =1
y y
se hace el tipico cambio de variables (Bernoulli)
1
v=
y2
Luego:
2 0
v0 = − y
y3
queda entonces, una ecuacion lineal:
2v 0 + v = −1
La solución de la EDO lineal es:
R
Z R

v(x) = e− 2dx
c + −1e dx
v(x) = ce−2x − 1
Pero
1
v=
y2
por lo tanto, la solución queda en:
1
y(x) = ± √
ce−2x −1
Page 38 of 117
Ejercicio 25
Muestre que la ecuación diferencial es homogénea, y resuélvala.
(x2 + y 2 )dx + 2xydy = 0
M (x, y) = x2 + y 2 , N (x, y) = 2xy
Ambas son funciones Homogeneas de Grado n=2, Luego:

dy − x2 + y 2
=
dx 2xy
2
dy −1 − xy
=
dx 2 xy
Hacemos el nunca mal ponderado cambio de variable:

y dy dz
z= ⇒ y = zx ⇒ =x +z
x dx dx
Reemplazamos:
dz −1 − z 2
x +z =
dx 2z
dz −1 − z 2
x = −z
dx 2z
dz −1 − z 2 − 2z 2
x =
dx 2z
dz −1 − 3z 2
x =
dx 2z
2z dx
dz = −
1 + 3z 2 x
Hemos llegado a la EDO (1) de Variales Separadas, integramos (You dont say)
Z Z
2z dx
dz = − + ln (c)
1 + 3z 2 x
1
ln 1 + 3z 2 = − ln |x| + ln (c) / · 3

3 c
ln 1 + 3z 2 = ln 3

x
2 c
⇒ 1 + 3z = 3
x
Volvamos a nuestras variables:
2
y c
⇒1+3 = 3
x2 x
⇒ x3 + 3y 2 x = c
p
−3x (x3 − 3c)
∴ y (x) = ±
3x
Page 39 of 117
Ejercicio 26
Resuelva la ecuación diferencial sujeta a la condición dada:
(x + y)2 dx + (2xy + x2 − 1)dy = 0; y(1) = 1
Solución:
Sea M = (x + y)2 ∧ N = 2xy + x2 − 1. Se vee que esta EDO es diferencial exacta pues
∂M ∂N
= = 2x + 2y
∂y ∂x
Por lo tanto existe una función F tale que :
∂F ∂F
=M∧ =N
∂x ∂y
con F (x, y) = 0 solución de la EDO

Por lo tanto, si se hace la integración a N , respecto de Y , se obtiene que:
Z y
N =F
Z y
2xy + x2 − 1 = xy 2 + yx2 − y + C(x)
Con C(x), una çonstante”que depende de x

∂F
Luego como F = xy 2 + yx2 − y + C(x) ∧ =M
∂x
se obtiene que
y 2 + 2yx + C 0 (x) = x2 + 2xy + y 2
C 0 (x) = x2
x3
→ C(x) = +C
3
x3
Por lo tanto, F = xy 2 + yx2 − y + +C
3
Por lo que la solución (implicita)de la EDO es :
x3
xy 2 + yx2 − y + =k
3
Por último, basta reemplazar x e y dados por la condición inicial y(1) = 1 para asi calcular k:
1 4
1+1−1+ =k→k=
3 3
Por lo tanto, al solución implicita al PVI es
x3 4
xy 2 + yx2 − y + =
3 3
Page 40 of 117
Ejercicio 27
Resolver
ds
+ s · tan(t) = 2t + t2 tan t
dt
Solución:
Hallemos el factor integrante definido por
Z
µ(x) = exp tg tdt = e− ln|cos t| = sec t
Entonces tendremos que
s0 sec t + s sec t tg t = 2t sec t + t2 sec t tg t

(s sec t)0 = (t2 sec t)0
Z Z
(s sec t)0 dt = (t2 sec t)0 dt
s sec t = t2 sec t + c
s = t2 + c cos t
Ejercicio 28
Resuelva
x2 y 00 + xy 0 − y = 0
Solución:
Esta es una ecuacion de euler - cauchy; del tipo:

x2 y 00 + axy 0 + by = 0; enestecaso :
a = 1; b = −1
En este caso, realizando el cambio de variable y = xm nos queda:
m2 + (a − 1)m + b = 0
m2 + (1 − 1)m − 1 = 0
m2 − 1 = 0
(m + 1) (m − 1) = 0
m1 = 1m2 = −1
y1 = x, y2 = x−1
Asi, la base de la solucion quedaria:
y = c1 y1 + c2 y2 ; c1 , c2 ,constantes. → y = c1 x1 + c2 x−1
Page 41 of 117
Ejercicio 29
Resolver, sabiendo que la solución particular es de la forma aekx
y 0 = −e−x y 2 + y + ex ; y (0) = 6
Recordando a Ricatti, Trasformandola en una ecuacion de Bernoulli:

w0 + (p − 2gs)w = gw2
-x x
⇔ w’ + ( - 1 - 2[ - e e ]w = - e - x w2
⇔ w’ + ( - 1 - 2[ - 1]w = - e - x w2
⇔ w’ + w = - e - x w2 , E.D.O. bernoulli: p = 1; g = - e - x ; u = w - 1
⇔ u’ - 1pu = - 1g
R R
⇔ u’ - u = e - x /F ; F = e pdx
= e− dx
= e−x (F: factor integrante)
⇔ e−x u’ - e−x u = e - 2x
d −x
⇔ e u = e - 2x
dx Z
⇔ e−x u = e - 2x dx
1 - 2x
⇔ e−x u = - e + c/ex
2
1 -x
⇔ u=- e + cex , u = w - 1 , nueva constante 2c = c
2
2
⇔ w= x
ce - e - x
2
Luego, por y = s + w (recordar Ricatti) : y = ex +
cex - e - x
ex (cex - e - x ) + 2
y=
cex - e - x
ce2x - 1 + 2
y=
cex - e - x
ce2x + 1
y= x
ce - e - x
por la condicion inicial:
y(0) = 6
ce0 + 1
⇔ =6
ce0 - e0
c+1
⇔ =6
c-1
⇔ c + 1 = 6c - 6
7
⇔ c=
5
7 2x
5e +1
∴ y= 7 x
5e - e-x
Page 42 of 117
Ejercicio 30
Resolver
dy
= sen(x + y)
dx
u=x+y
dy du
= −1
dx dx
Reemplazamos en la EDO (1):
du du
− 1 = sen (u) = dx
dx sen (u) + 1
Integrando:
Z Z
du
= dx + c
sen (u) + 1
| {z }
I
Trabajamos la integral, hacemos un cambio de variable:

2t 2
⇒ sen (u) = 2
du = dt
1+t 1 + t2
Reemplazando en I:
Z 2 Z
1+t2 dt
I= 2t dtI = 2 2
1+t2 + 1 (t + 1)
2
I=−
t+1
2
I=− u

tan 2 +1
Volviendo a la EDO (1):

Z Z
du
= dx + c
sen (u) + 1
2
− u
=x+c
tan 2 + 1
u 2
tan =− −1
2 x+c

u 2
= − arctan +1
2 x+c

2
u = −2 arctan +1
x+c
Volviendo a la variable Original:

2 2
x + y = −2 arctan + 1 ∴ y (x) = − x + 2 arctan +1
x+c x+c
Page 43 of 117
Ejercicio 31
1 −e2x
y 0 + e−2x y 2 − (1 + 4x + 2x2 )y = (1 + x + 2x2 + x3 )
x x
se pide encontrar la solución particular de la forma y1 (x) = e2x (Ax + B)
Solución:
Ricatti de nuevo (será un ejercicio propuesto claro)
S (x) = e2x (Ax + B)
es una solucion particular, en este caso este seria nuestro .a reemplazar, encontremos nuestro y’
2
S 0 (x) = e2x A + e2x 2 (Ax + B)
S 0 (x) = e2x (A + 2Ax + 2B)
Ahora es cosa de reemplazar en nuestra EDO de Ricatti, pero ojo, hay q acomodarla un poco, para q se
vea como una ec de ricatti:
Recordemos que la ecuacion de Ricatti es de la forma:
y 0 = P (x) + Q (x) y + R (x) y 2
ahora nuestra ecuacion quedaria de la forma:
e2x 1
y0 = − 1 + x + 2x2 + x3 + 1 + 4x + 2x2 y − e2x y 2

x x

Luego Reemplazamos nuestra solucion particular en la y e y 0 :
e2x 1
e2x (A + 2Ax + 2B) = − 1 + x + 2x2 + x3 + 1 + 4x + 2x2 e2x (Ax + B)

x x
2
−e−2x e2x (Ax + B)
e2x 1
e2 (A + 2Ax + 2B) = − 1 + x + 2x2 + x3 + 1 + 4x + 2x2 e2x (Ax + B)

x x
−e−2x e4x A2 x2 + 2ABx + B 2

e2x 1
e2x (A + 2Ax + 2B) = − 1 + x + 2x2 + x3 + 1 + 4x + 2x2 e2x (Ax + B)

x x
−e2x A2 x2 + 2ABx + B 2

1 1
e2x (A + 2Ax + 2B) = e2x − 1 + x + 2x2 + x3 + 1 + 4x + 2x2 (Ax + B)

x x
2 2 2

− A x + 2ABx + B
1 B
 
2 2
− − 1 − 2x − x + A + + 4Ax + 4B + 2Ax + 2xB
e2x (A + 2Ax + 2B) = e2x  x x 
2 2 2
−A x − 2ABx − B
1
 
2 2

x −1 + 2A − A + x (−2 + 4A + 2B − 2AB) + (−1 + B)
e2x (A + 2Ax + 2B) =  x 
+ −1 + A + 4B − B 2


Ahora con esto podemos formar 4 ecuaciones igualando los terminos correspondientes del lado izquierdo
con el derecho, para asi poder calcular A y B

1 − 1 + 2A − A2 = 0

2 − 2 + 4A + 2B − 2AB = 2A

3 −1+B =0

4 − 1 + A + 4B − B 2 = A + 2B

⇒
3
B=1
1
−1 + 2A − A2 = 0
A2 − 2A + 1 = 0
2
(A − 1) = 0
A=1
Con esto encontramos los valores de A y B, asi nuestra solucion particular es:
S (x) = e2x (x + 1)
Despues haces un cambio de variable, tipico para las ecuaciones de Ricatti:
1
y = S (x) +
z
1
y = e2x (x + 1) +
z
1 dz
y 0 = e2x (2x + 3) −
z 2 dx
Page 46 of 117
Ejercicio 32
Resolver
y(ln(y) − ln(x) + 1)
y0 =
x
Solución:
Reescribiendo la ecuación
y(ln y − ln x + 1) y y
y0 = ⇐⇒ y 0 = ln + 1
x x x
dy y dv
Es de la forma =F . Efectuamos un cambio definido por y = vx =⇒ y 0 = x + v, ası́ que
dx x dx
dv
x +v = v ln v + v
dx
1 1
dv = dx
Z v ln v x
Z
1 1
dv = dx
v ln v x
ln |ln v| = ln |x| + ln |c|
ln v = cx
y = xecx
Page 47 of 117
Ejercicio 33
Resolver con Transformada de Laplace
y” + 4y = 16x, y(0) = 3, y 0 (0) = −6
Solución:
Conociendo algunos valores de la transformada de Laplace
L [a · f (t) + b · g (t)](s) = a · L [f (t)](s) + b · L [g (t)](s)
n!
L [tn ](s) =
sn+1
s
L [cos (a · t)](s) =
s2 + a2
a
L [sen (a · t)](s) =
s2 + a2
L [f n (t)](s) = sn L [f (t)](s) − sn−1 f (0) − sn−2 f 0 (0) ... − s · f n−2 (0) − f n−1 (0)

Resolviendo finalmente
y 00 + 4y = 16t y (0) = 3 y 0 (0) = −6
L [y 00 ](s) + 4L [y](s) = 16L [t](s)
1
s2 ϕ (s) − s · y (0) − y 0 (0) + 4ϕ (s) = 16 · 2

s
16
s2 ϕ (s) − 3s + 6 + 4ϕ (s) = 2

s
16
ϕ (s) s2 + 4 = 2 + 3s − 6
s
16 3s 6
ϕ (s) = + 2 − 2
s2 2
(s + 4) (s + 4) (s + 4)

1 1 s 2
ϕ (s) = 4 − + 3 − 3
s2 s2 + 4 s2 + 4 s2 + 4

1 2 s 2
L−1 [ϕ (s)](t) = 4L−1 − 2L−1 + 3L−1 − 3L−1
s2 (t) s2 + 4 (t) s2 + 4 (t) s2 + 4 (t)
y (t) = 4t − 2sen (2t) + 3 cos (2t) − 3sen (2t)
y (t) = 4t + 3 cos (2t) − 5sen (2t)
Page 49 of 117
Ejercicio 34
(
x0 = x − y x(0) = 4
Resuelva el PVI:
y 0 = 5x − y y(0) = 0
Solución:
Por Laplace
Sea L(x) = X y L(y) = Y
sX − x(o) = X−Y (1)

sY − y(o) = 5X − Y (2)
x(0) = 4, y(0) = 0
Luego:

4 1 s−1 4

0 s+1 −5 0
X= , Y=

s−1 1 s−1
1

−5 s + 1 −5 s + 1
4s 2·2 10 · 2
X= + 2 , Y=
s2 +2 2 s + 22 s2 + 22
Aplicando la transformada inversa . . .
x(t) = 4 cos 2t + 2 sin 2t, y(t) = 10 sin 2t

Ahora reemplazando una ecuación en la otra (creeı́an que solo con el Wroskiano se podı́a resolver)
(
x0 = x − y x(0) = 4
y 0 = 5x − y y(0) = 0
y = x − x0
0
(x − x0 ) = 5x − (x − x0 )
x0 − x00 = 5x − x + x0
x00 + 4x = 0
r2 + 4 = 0
r = ±2i
x(t) = c1 e2it + c2 e−2it
y(t) = (c1 e2it + c2 e−2it ) − (c1 e2it + c2 e−2it )0 =
= c1 e2it + c2 e−2it − 2ic1 e2it + 2ic2 e−2it =
= c1 (1 − 2i)e2it + c2 (1 + 2i)e−2it
(
4 = c1 + c2
0 = c1 (1 − 2i) + c2 (1 + 2i)
c1 = 2 − i
c2 = 2 + i
x(t) = (2 − i)e2it + (2 + i)e−2it
y(t) = (2 − i)(1 − 2i)e2it + (2 + i)(1 + 2i)e−2it =
= −5ie2it + 5ie−2it
x(t) = (2 − i)(cos 2t + i sin 2t) + (2 + i)(cos 2t − i sin 2t)
y(t) = −5i(cos 2t + i sin 2t) + 5i(cos 2t − i sin 2t)
x(t) = 2 cos 2t + sin 2t + 2 cos 2t + sin 2t = 4 cos 2t + 2 sin 2t
y(t) = 5 sin 2t + 5 sin 2t = 10 sin 2t
Page 51 of 117
Ejercicio 35

0
x = x x(0) = 1


Resuelva el PVI, a traves del metodo de los valores/vectores propios: 0
y =x+y y(0) = 1

z 0 = x + y + z

z(0) = 1
Solución:
Sea X la matriz
 
x
X= y 
z
El sistema se escribe matricialmente como:

   
1 0 0 1
X0 =  1 1 0  X; X(0) =  1 
1 1 1 1
La solución de este sistema tiene la forma:
X = φ ert
donde r es el valor propio y φ es el vector propio asociado a r. Los valores propios de esta matriz
se obtienen fácilmente de acuerdo al polinomio caracterı́stico:
p(r) = (r − 1)3
Aquı́ se obtiene r = 1 de multiplicidad 3.

En este caso, hay un único vector propio asociado a dicho valor propio . . .
Entonces, obtenemos el vector propio para r = 1:
 
0 0 0
 1 0 0  φ = 0v
1 1 0
Luego:
     
φ1 0 0
 φ 2  =  0  = φ3  0 
φ3 φ3 1
Podemos considerar:
 
0
φ= 0 
1
y por ahora:
 
0
X(1) =  0  et
1

Ahora tomamos una solución de la forma:
X = φ t et + η et
derivando,
X0 = φ et + φ t et + η t et
Reemplazando en la ecuación inicial . . .
X0 = φ et + φ t et + η t et = Aφ t et + Aη et
Luego:
(A − I) φ = 0V , (A − I) η = φ
La primera ecuación está resulta . . .

De aquı́ se desprende que :
 
0
η= 1 
1
Llegamos a una solución:

   
0 0
X(2) =  0  t et +  1  et
1 1
por último consideramos una solución de la forma:
t2 t
X=φ e + η t et + ϕ et
2
Satisfaciendo las ecuaciones:
(A − I) φ = 0V , (A − I) η = φ , (A − I) ϕ = η
La primera y segunda ecuación ya han sido resueltas . . .

Luego, se encuentra que:


1
ϕ= 0 
1
Entonces:
     
0 2 0 1
t
X(3) =  0  et +  1  t et +  0  et
2
1 1 1
Por lo tanto la solución general del sistema es:

      
0 0 0
X= C1  0  et + C2  0  t et +  1  et  +
1 1 1

      
0 2 0 1
t t t
C3  0  e +  1  te +  0  et 
2
1 1 1
y finalmente, aplicando la condición inicial:
C1 = −1 , C2 = 1 , C3 = 1 ,
     
1 0 0
=
X(t)  1  et +  1  t et +  0  t2 et
1 2 1/2
Page 54 of 117
Ejercicio 36
Resolver
y 00 + 5y 0 + 4y = 17sin(x) + ex
Solución:
En lo personal, no me gusta el método de los aniquiladores, pero bueno aca va.

Sea D una operador lineal tal que Dy = y 0 . Entonces la ecuación
y 00 + 5y 0 + 4y = 17 sin x + ex
se puede escribir como
(D + 4)(D + 1)y = 17 sin x + ex
.
Vemos aquı́ que la solución de la ecuación homogénea es yh = c1 e−4x + c2 e−x . Ahora nos falta encontrar
una solución particular y estamos listos.
Podemos anular el lado derecho de la ecuación aplicando (D − 1) y (D2 + 1).
Entonces:
(D + 4)(D + 1)(D − 1)(D2 + 1)y = 0
.
La solución particular tiene la forma:
yp = c3 ex + c4 sin x + c5 cos x
yp0 = c3 ex + c4 cos x − c5 sin x
yp00 = c3 ex − c4 sin x − c5 cos x

Como es solución de la ecuación original, podemos reemplazarla para obtener condiciones para c3 ,c4 y c5 .
Quedamos con:
10c3 ex + (3c4 − 5c5 ) sin x + (3c5 + 5c4 ) cos x = 17 sin x + ex

1 51 −85
Resolviendo este simple sistema de ecuaciones, averiguamos que c3 = 10 , c4 = 34 y c5 = 34
Como y = yh + yp , entonces
1 x 51 85
y= e + sin x − cos x + c1 e−4x + c2 e−x
10 34 34
, para dos constantes c1 y c2 cualquieras
Page 55 of 117
Ejercicio 37
Muestre que la ecuación
(2x + xy 2 )dx + (4y + x2 y)dy = 0

es de variables separables, exacta y de Bernoulli, luego escoja uno de los casos anteriores, y resuelva.
Solución:
(2x + xy 2 )dx + (4y + x2 y)dy = 0 (3)
Veamos si existe F (x, y) = c
∂F ∂F
p(x, y) = = 2x + xy 2 q(x, y) = = 4y + x2 y
∂x ∂y
Luego:
∂p ∂q
= 2xy = 2xy
∂y ∂x
Se muestra que es exacta.
A partir de la primera ecuación:
2x 1 x
y0 + + y=0
4 + x2 y 4 + x2
x −2x −1
y0 + y= y
4 + x2 4 + x2
Se muestra que es Bernoulli. Nuevamente tomando la primera ecuación:
x y
dx + dy = 0
4 + x2 2 + y2
Se muestra que es de variables separables. Usando el primer método:
∂F
= 2x + xy 2
∂x
1
⇒ F = x2 + x2 y 2 + g(y)
2
∂F
= xy 2 + g 0 (y) = 4y + xy 2
∂y
⇒ g(y) = 2y 2 + c1
Finalmente la solución general es:
1
x2 + x2 y 2 + 2y 2 + c = 0
2
Page 56 of 117
Ejercicio 38
Pruebe que, si y1 es solución de y10 + p(x)y1 = q1 (x) y si y2 es solucion de y20 (x) + p(x)y2 = q2 (x), entonces
la EDO y 0 + p(x)y = (α)q1 (x) + (β)q2 (x) tiene por solucion a y(x) = (α)y1 (x) + (β)y2 (x)
Solución:
y10 + p(x)y1 = q1 (x) ⇒ αy10 + αp(x)y1 = αq1 (x)
y20 + p(x)y2 = q2 (x) ⇒ βy20 + βp(x)y2 = βq2 (x)
Sumando:
αy10 + αp(x)y1 + βy20 + βp(x)y2 = αq1 (x) + βq2 (x)
(αy10 + βy20 ) + (αp(x)y1 + βp(x)y2 ) = αq1 (x) + βq2 (x)
(αy10 + βy20 ) + (αy1 + βy2 )p(x) = αq1 (x) + βq2 (x)
Como α, β constantes:
(αy1 + βy2 )0 + (αy1 + βy2 )p(x) = αq1 (x) + βq2 (x)
Luego la solución de y 0 + p(x)y = αq1 (x) + βq2 (x) es:
y = αy1 + βy2
Page 57 of 117
Ejercicio 39
Resuelva:
Z t
f (u)
1+t= √ du
0 t−u
Hint:
Γ (α + 1)
L [tα ](s) = ;s > 0 ; α > −1
sα+1
√

1
Γ = π
2
Γ (n + 1) = nΓ (n)
Z ∞
Donde Γ es la función gamma Γ(n) = tn−1 e−t dt
0
Solución:
Notar la semejanza de la función Gamma con la transformada de Laplace. Utilizando la convolución:
Z t
1 1
f (u) · √ du = f∗√ (t)
0 t−u t
Aplicando la transformada de Laplace a ambos lados de la ecuación:

1 1 1
+ 2 = L(f ) · L √
s s t
1 1 Γ(1/2)
+ 2 = L(f ) · 1/2
s s s
1 1 1 1
L(f ) = +
Γ(1/2) s1+(−1/2) Γ(1/2) s1+(1/2)
1 Γ(1 − 1/2) 2 Γ(1 + 1/2)

L(f ) = + 2
Γ2 (1/2) s1+(−1/2) Γ (1/2) s1+(1/2)
1 −1/2 2 1/2
f (t) = t + t
π π
Page 58 of 117
Ejercicio 40
1
Resolver en terminos de una integral definida, y 0 + y = con y(2) = 1
x+1
Solución:
1 ex
y0 + y = ⇒ y 0 ex + yex =
x+1 x+1
ex
⇒ (yex )0 =
x+1
d(yex ) ex ex
= ⇒ d(yex ) = dx
dx x+1 x+1
Si x = 2, yex = e2 Integrando:
yex x
ex
Z Z
d(yex ) = dx
e2 2 x+1
x x
ex ex
Z Z
yex − e2 = dx ⇒ yex = dx + e2
2 x+1 2 x+1
x
ex
Z
y(x) = e−x dx + e2−x
2 x+1
Page 59 of 117
Ejercicio 41
[n cos(nx + my) − msen(mx + ny)]dx + [m cos(nx + my) − nsen(mx + ny)]dy = 0

Solución:
Si m, n son constantes reales.

Sea M (x, y) = n cos(nx + my) − m sin(mx + ny), N (x, y) = m cos(nx + my) − n sin(mx + ny).
Se cumple que
∂M ∂N
= = −mn sin(nx + my) − mn cos(mx + ny)
∂y ∂x
por lo que existe F (x, y) = c solución de la ecuación.
Z Z
F (x, y) = n cos(nx + my)dx − m sin(mx + ny) = sin(nx + my) + cos(mx + ny) + h(y)
∂F
= m cos(nx + my) − n sin(mx + ny) + h0 (y) = N (x, y)
∂y
⇒ h0 (y) = 0 ⇒ h(y) = c0
La solución es:
sin(nx + my) + cos(mx + ny) = C
Page 60 of 117
Ejercicio 42
Usando la transformacion t = sen (x) convierta la ecuación diferencial
y 00 + tan (x) y 0 + cos2 (x) y = 0

en una ecuación con coeficientes constantes. Utilice esta transformacion para resolver la ecuacion diferencial.
Solución:
Si t = sin(x)
dy dy dt dy
= = cos(x)
dx dt dx dt
d2 y

d dy dt d dy dt d dy
2
= = = cos(x)
dx dx dx dx dt dx dx dt dt
d2 y dt d2 y d2 y

dt dy dx dy dy
= 2
cos(t) − sin(t) = 2
cos(t) − sin(t) = 2 cos2 (t) − sin(t)
dx dt dt dt dx dt dt dt dt
Reemplazando en la ecuación inicial:
d2 y

dy dy
cos2 (t) − sin(t) + tan(x) cos(x) + cos2 (t)y = 0
dt2 dt dt
d2 y
cos2 (t) + cos2 (t)y = 0
dt2
Considerando que
(1 + 2n)π
x 6=
2
sólo queda:
d2 y
+y =0
dt2
Aplicando transformada de Laplace L(y) = Y :
B As
Y s2 − sA − B + Y = 0 ⇒ Y (s2 + 1) = B + As ⇒ Y = +
s2 + 1 s2 + 1
Aplicando la transformada inversa:
y(t) = B sin(t) + A cos(t) ⇒ y(x) = B sin(sin(x)) + A cos(sin(x))

displaymath39
Page 61 of 117
Ejercicio 43
Resolver el PVI, donde δ corresponde al delta de dirac
y 0 + y = δ(t − 1), y(0) = 1
Solución:
y 0 + y = δ(t − 1) y(0) = 1
Se sabe que L(δ(t − t0 )) = e−st0

Sea L(y) = Y (s), Aplicando la transformada de Laplace:
L(y 0 ) + L(y) = L(δ(t − 1))
sY (s) − y(0) + Y (s) = e−s

e−s + 1 e−s 1
Y (s) = = +
s+1 s+1 s+1
Aplicando la transformada inversa:
y(t) = e−(t−1) u1 (t) + e−t
y(t) = e1−t u1 (t) + e−t
Ejercicio 44
Resolver con Laplace:
∞
X
y 00 + y = δ2jπ (t), y(0) = y 0 (0) = 0
j=0
e2πs 1
L {y(t)} = = L{x(t)}L{sin t}.
e2πs − 1 s2 + 1
, notamos que es una suma geometrica e−2πs < 1
 
e2πs 1 X
−2πks
X 
= = e = L δ −2kπ (t) .
e2πs − 1 1 − e−2πs  
k≥0 k≥0
Ahora, simplemente utilizamos convolución:

 
X X X
y(t) =  δ−2kπ (t) ∗ sin t = (δ−2kπ (t) ∗ sin t) = sin(t − 2kπ)
k≥0 k≥0 k≥0
Page 62 of 117
Ejercicio 45
d2 y

3 dy 1
4x4 2 + 8x + y = tan
dx dx 2x
1
Indicacion: use el cambio de variable x =
t
Solución:
1
Sea t = , luego:
x
dy dy dt dy −1
= = ·
dx dt dx dt x2
d2 y

d dy dt d dy dt d dy −1
= = = ·
dx2 dx dx dx dt dx dx dt dt x2
2
d2 y 1

dt d y −1 dy 2 dx dy 2
= 2
· 2
+ 3
· = 2
· 4+ ·
dx dt x dt x dt dt x dt x3
Reemplazando en la ecuación inicial:
d2 y 1

dy 2 dy −1 t
4x4 · + · + 8x 3
· + y = tan
dt2 x4 dt x3 dt x2 2

t
4y 00 + y = tan
2
La ecuación homogenea asociada es:
4y 00 + y = 0
que tiene por solución:
yh = c1 cos(t/2) + c2 sin(t/2)
Usando el método de variación de parámetros, se puede hallar la solución particular:
yp = u1 cos(t/2) + u2 sin(t/2)
donde u1 = u1 (t), u2 = u2 (t).

Finalmente se obtiene:
" #

t
Z
t

t t t 1 + sin t
2
u2 (t) = −4 cos u1 (t) = −2 sin tan dt = 4 sin sin − ln

2 2 2 2 2 cos 2t
La solución general de la ecuación es:

" #

t

t

t t t 1 + sin t
2 2
y(t) = c1 cos + c2 sin − 4 cos + 4 sin sin − ln

2 2 2 2 2 cos 2t
1
Y volviendo a t = se obtiene la solución original.
x
Page 63 of 117
Ejercicio 46
Sea y 0 = 1 + sin(xy) y D = {(x, y) ∈ R2 : 0 < x < 1, −∞ < y < ∞} Demuestre que para cualquier punto
(a, y0 ) con 0 < a < 1, el PVI asociado a algun intervalo centrado en a tiene solucion unica
Solución:
Tenemos el PVI:
y 0 = 1 + sin(xy) = f (x, y)
y(a) = y0
f (x, y) es continua en todo R2, en particular, también en D. Por lo tanto existe AL MENOS UNA
solucion definida en un intervalo abierto J, que contiene a x=a
Además:
∂f (x, y)
= −x cos(xy)
∂y
La cual tambien es continua en todo R2, en particular D, entonces la solución también es UNICA en un
intervalo abierto J0 que contiene a x=a.El teorema de existencia y unicidad nos asegura esto.
Ejercicio 47
Resuelva, y diga que sucede en este caso
10xy 0,4
y0 = , y(0) = −1
3
Por variables separables (do it!) y evaluando la condicion inicial queda:
y = (x2 − 1)5/3
El problema radica que hay un problema de continuidad, ya que
10xy 0,4
f (x.y) =
3
tiene diferente rango que la siguiente función, contradiciendo la unicidad, y no la existencia
4xy −0,6
fy =
3
Page 64 of 117
Ejercicio 48
Determine los valores de k ∈ R , para los que y = ekx es solucion de la ecuacion

diferencial y 000 − 6y 00 + 11y 0 − 6y = 0 .Si k1 , k2 , k3 son los valores encontrados ,
verifique que independientemente de los valores de las constantes c1 , c2 , c3 , la
funcion y = c1 ek1 x + c2 ek2 x + c3 ek3 x tambien satisface la ecuacion dada .
Solución:
Es fácil ver que
y = ekx ⇒ y 0 = kekx (4)

00 2 kx
⇒y =k e (5)
000 3 kx
⇒y =k e (6)
Luego, la ecuación a resolver la podemos reescribir como
k 3 ekx − 6k 2 ekx + 11kekx − 6ekx = 0 (7)
Dado que ekx > 0, ∀x ∈ R, podemos simplicar la ecuación (4) y reescribirla como
k 3 − 6k 2 + 11k − 6 = 0 (8)
Pero, al factorizarla obtenemos
(k − 3)(k − 2)(k − 1) = 0 (9)
De lo cual se concluye que k1 = 1, k2 = 2 y k3 = 3.

.
Ahora, para la segunda parte, dado que
y = c1 ek1 x + c2 ek2 x + c3 ek3 x (10)

0 k1 x k2 x k3 x
⇒ y = c1 k1 e + c2 k2 e + c3 k3 e (11)
00
⇒y = c1 k12 ek1 x + c2 k22 ek2 x + c3 k32 ek3 x (12)
000
⇒y = c1 k13 ek1 x + c2 k23 ek2 x + c3 k33 ek3 x (13)
Luego, la ecuación diferencial podemos reescribirla como
c1 ek1 x (k13 − 6k12 + 11k1 − 6) + c2 ek2 x (k23 − 6k22 + 11k2 − 6) + c3 ek3 x (k33 − 6k32 + 11k3 − 6) = 0
Como cada uno de los paréntesis es cero, se concluye que no importa el valor de los coeficientes c1 , c2 , c3 .
Page 65 of 117
Ejercicio 49
Sean y1 (x), y2 (x) dos soluciones de la EDO en un intervalo abierto, donde A, B, C son continuas y A(x)
nunca se anula.
A(x)y 00 + B(x)y 0 + C(x)y = 0,
Considere
y y2
W (x) = 10 .
y1 y20
Demuestre que
dW (x)
A(x) = A(x)y1 y200 − A(x)y2 y100 .
dx
Sustituya Ay200 y Ay100 en la EDO original para demostrar que
dW (x)
A(x) = −B(x)W (x).
dx
Por último, resuelva esto último para deducir
Z
B(x)
W (x) = K exp − dx .
A(x)
Solución:
W (x) = y1 y20 − y2 y10 ⇒ W 0 (x) = y10 y20 + y1 y200 − y20 y10 − y2 y100 = y1 y200 − y2 y100
∴ A(x)W 0 (x) = A(x)y1 y200 − A(x)y2 y100

Despejando se obtiene:
y1 Ay200 = AW 0 + Ay2 y100
y2 Ay100 = Ay1 y200 − AW 0

Multiplicando la EDO inicial por y1 (la que contiene y2 ) e y2 (la que contiene y1 ) y reemplazando el
último resultado en estas nuevas ecuaciones se obtiene:
Ay1 y200 − AW 0 + By10 y2 + Cy1 y2 = 0
AW 0 + Ay2 y100 + By20 y1 + Cy1 y2 = 0

restando se obtiene:
2AW 0 − A(y1 y200 − y2 y100 ) + B(y1 y20 − y2 y10 ) = 0
2AW 0 − AW 0 + BW = 0
AW 0 = −BW
Resolviendo para W (x):
dW
A(x) = −B(x)W (x)
dx
dW B(x)
=− dx
W A(x)
Z
B(x)
Z
B(x) − dx
ln(W ) = − dx + c ⇒ W (x) = Ke A(x)
A(x)
Page 66 of 117
Ejercicio 50
Utilizando el método de Frobenius, resolver la EDO:
4xy 00 + 2y 0 + y = 0
Solución:
∞
X
Supongamos existe una solución de la forma: y(x) = Ak xk+ρ con la condición de que A0 6= 0
k=0
Le daremos a la ecuación la forma: x2 y 00 + xf (x)y 0 + g(x)y = 0
x x x
⇒ 4xy 00 + 2y 0 + y = 0 / ⇒ x2 y 00 + y 0 + y = 0
4 2 4
(todo lo anterior se hizo dado que x=0 es un pto. singular regular de la EDO) de donde se extrae que:
1 1
f (x) = ⇒ f0 = , fk = 0, ∀k ≥ 1
2 2
x 1
g(x) = ⇒ g0 = 0, g1 = , gk = 0, ∀k ≥ 2
4 4
Ahora, la función caracterı́stica para la EDO viene dada por:I(ν) =ν(ν − 1) + νf0 + g0
ν 1
en nuestro caso: I(ν) = ν(ν − 1) + = ν ν −
2 2
1
Como sabemos debe cumplirse que I(ρ) = 0 ⇒ ρ1 = 0 ∧ ρ2 =
2
Para este tipo de ecuaciones, la ley de recurrencia viene dada por:
n−1
X
− [(k + ρ)fn−k + gn−k ]Ak
k=0
An =
I(n + ρ)
−1 A
⇒ An = · n−1
4 1
(n + ρ) n + ρ −
2
Tomando ρ1 = 0:
−An−1 −An−1
An = =
1 2n(2n − 1)
4n n −
2
Ahora, veamos algunos valores:
−A0 −A1 A0 (−1)k A0
A1 = , A2 = = , Ak =
1·2 4·3 1·2·3·4 (2k)!
∞ k ∞ k √
X (−1) A0 k X (−1) √
Ası́, tenemos: y1 (x) = x = A0 ( x)2k = A0 cos( x)
(2k)! (2k)!
k=0 k=0
1
Ahora, para la segunda solución, tomando ρ2 = :
2
−An−1 −An−1
An = =
1 2n(2n + 1)
4n n +
2
entonces:
−A0 −A1 A0 (−1)k A0
A1 = , A2 = = , Ak =
2·3 4·5 2·3·4·5 (2k + 1)!
∞ ∞
X (−1)k A0 k+ 1 X (−1)k √ 2k+1 √
Luego: y2 (x) = x 2 = A0 ( x) = A0 sin( x)
(2k + 1)! (2k + 1)!
k=0 k=0 √ √
Finalmente, la solución de la EDO queda dada por: y(x) = A0 cos( x) + A0 sin( x)
Page 67 of 117
Ejercicio 51
y + senx cos2 xy x
2
dx + ( 2 + seny)dy = 0
cos xy cos xy
Solución:
y x
M (x, y) = + sin(x) N (x, y) = + sin(y)
cos2 (xy) cos2 (xy)
∂M 1 2(xy) tan(xy) ∂N
= 2
+ 2
=
∂y cos (xy) cos (xy) ∂x
Luego existe F (x, y) = C solución de la EDO.
Z
x
F (x, y) = + sin(y) dy = tan(xy) − cos(y) + h(x)
cos2 (xy)
∂F y y
= + h0 (x) = M (x, y) = + sin(x)
∂x cos2 (xy) cos2 (xy)
⇒ h0 (x) = sin(x) ⇒ h(x) = − cos(x) + c1
Luego la solución es:
F (x, y) = tan(xy) − cos(y) + − cos(x) = C
Page 68 of 117
Ejercicio 52
 Z t (
y 0 − 4y + 8

y(u)du = δ(t − 3) + h(t) 0 t<4
Resuelva 0 donde h(t)=
y(0) = 0
 1 4≤t
Solución:
La ecuación se reescribe como sigue:
y 0 − 4y + 8(y ∗ 1)(t) = δ(t − 3) + u4 (t)
Sea Y = L(y). Aplicando la transformada de Laplace a ambos lados:
1 e−4s
sY − y(0) − 4Y + 8Y = e−3s +
s s
8 se−3s + e−4s
Y (s − 4 + ) =
s s
se−3s e−4s
Y = +
(s − 2)2 + 22 (s − 2)2 + 22
Considerando el primer término:
s−2 2
F1 (s) = +
(s − 2)2 + 22 (s − 2)2 + 22
f1 (t) = e2t cos(2t) + e2t sin(2t)
Sabiendo que: L(uc (t)f (t − c)) = e−cs F (s) donde F (s) = L(f (t)),
la transformada inversa queda: u3 (t)f1 (t − 3). Análogo para el segundo término.
Ası́, se obtiene:
e2(t−4) sin(2(t − 4))

y(t) = u3 (t)e2(t−3) (cos(2(t − 3)) + sin(2(t − 3))) + u4 (t)
2
Page 69 of 117
Ejercicio 53
Dos tanques de 100 galones están inicialmente llenos de agua pura.
lb gal
Agua que contiene sal con una concentración de 1 gal entra al estanque A desde una fuente externa a 1 min .
gal gal
El agua fluye del tanque A al tanque B a 3 min . En el tanque B, el agua se evapora a una tasa de 1 min y
gal
es bombeada al tanque A a 2 min .
Encuentre la cantidad de sal en cada estanque en función del tiempo.
Solución:
Sean xA (t), xB (t) las cantidades de sal en el estanque A y B, respectivamente, a través del tiempo.
Quedan claras las condiciones iniciales. xA (0) = xB (0) = 0.
La variación de sal en el tiempo queda definida por la diferencia entre la tasa de entrada de sal y la tasa
de salida de sal.
Esto es:
lb lb gal
= · ≡ caudal · concentracion
min gal min
Para el estanque A:

dxA lb gal gal xB (t) lb gal xA (t) lb
= 1 ·1 +2 · −3 ·
dt gal min min 100 gal min 100 gal
Para el estanque B:

dxB gal xA (t) lb gal xB (t) lb gal lb
= 3 · − 2 · +1 ·0
dt min 100 gal min 100 gal min gal
Estas dos ecuaciones se pueden reescribir:
2xB 3xA
x0A = 1 + −
100 100
3xA 2xB
x0B =
−
100 100
Aplicando la transformada de Laplace (con las C.I. ya definidas) para ambas ecuaciones y definiendo
L(xA ) = X(s), L(xB ) = Y (s):

3 2Y (s) 1
s+ X(s) − =
100 100 s

3X(s) 2
− + s+ Y (s) = 0
100 100

Luego:
3 100s + 2
Y (s) = X(s) =
s2 (100s + 5) s2 (100s + 5)
Lo que puede escribirse como sigue:
1 3 1 1 1 2 1 1
Y (s) = −12 + + 12 5 X(s) = 12 + − 12 5
s 5 s2 s + 100 s 5 s2 s + 100
Aplicando la transformada inversa de Laplace:
2 −t
xA (t) = 12 + t − 12e 20
5
3 −t
xB (t) = −12 + t + 12e 20
5
Page 71 of 117
Ejercicio 54
Resolver
s 2
dy
1+
d2 y dx
=k
dx2 a−x
donde a, k son constantes tales que a ∈ R, k 6= 0, 1 y x < a
Solución:
k p
y 00 = 1 + y 02
a−x
Sea v = y 0 ,
dv k p
v0 = = 1 + v2
dx a−x
Luego:
1 k
√ dv = dx
1+v 2 a−x
Integrando:
p
k
ln v + 1 + v 2 = ln (a − x)
p
1 + v 2 = (a − x)k − v
1 + v 2 = (a − x)2k − 2(a − x)v + v 2
1 1 1
v = y0 = (a − x)k −
2 2 (a − x)k
Como k 6= 0, 1:
−1 1
y(x) = (a − x)k+1 + (a − x)1−k k 6= ±1
2(k + 1) 2(1 − k)
Page 72 of 117
Ejercicio 55
Resolver
x0 = x cos t + y sin t + z sin t

0
y = y cos t + z sin t
0
z = z cos t
Solución:
La primera ecuación es de variables separables
z0
z 0 = z cos(t) ⇔ = cos(t)
z
Integramos respecto a t llegando a:
ln(|z|) = sin(t) + C
⇒ z = c1 esin(t)
Ahora reemplazamos z(t) en la segunda ecuación y llegamos a:
y 0 − y cos(t) = c1 esin(t) sin(t)
La ecuación anterior es lineal con factor integrante e− sin(t) , por lo que multiplicando por el factor
mencionado llegamos a:
(ye− sin(t) )0 = c1 sin(t)
E integrando llegamos a:
yesin(t) = −c1 cos(t) + c2
⇒ y = esin(t) · (c2 − c1 cos(t))

Para x(t) el procedimiento es análogo
x0 − x cos(t) = sin(t)esin(t) · (c2 − c1 cos(t) + c1 )
⇒ (xesin(t) )0 = c2 sin(t) − c1 cos(t) sin(t) + c1 sin(t)

1
⇒ x = esin(t) · ( c1 · (cos(2t) − 4 cos(t)) − c2 cos(t) + c3 )
4
Es claro que las soluciones son periódicas (sólo hay senos o cosenos en las soluciones).

La Matriz Fundamental no es más que la que satisfaga la relación
E 0 (t) = AE(t)
Donde E 0 (t) y E(t) en este caso vendrı́an a ser x0 (t)...z 0 (t) y x(t)..z(t), respectivamente. Ahora,
respondiendo puntualmente a la preguntas
1) Matriz Fundamental
 
Cos(t) Sin(t Sin(t)
E 0 (t) =  0 Cos(t) Sin(t)  E(t)
0 0 Cos(t)
Donde
 0   
x (t) x(t)
E 0 (t) = y 0 (t) E(t) = y(t)
z 0 (t) z(t)
2) Indique si existen o no. En este caso observamos que la Matriz A es Triangular Superior por lo cual
sabemos que los valores propios del sistema serán λ1 , λ2 , λ3 = Cos(t) Esto producirá una solución general
que puede ser factorizada por ec os(t) Aunque eso es otro cuento. Como te piden que indiques si tiene o
no soluciones periodicas. Ademas de mencionar que todos los valores propios son cosenos puedes indicar
la solución general del sistema que serı́a
E(t)g = eCos(t) V1 C1 + eCos(t) tV2 C2 + eCos(t) t2 V3 C3
Esa es una solución que se comporta perı́odicamente. siendo V1 , V2 , V3 y C1 , C2 , C3 los Vectores

Propios(tb conocidos como eigenvectores) y las constantes.
Ejercicio 56
Probar que las siguientes funciones son son linealmente independientes en cualquier intervalo I.
ekx , xekx
Solución:
El ejercicio es tonto, pero DEBEN HACERLO EN LA I2, VALE PUNTOS, EN SERIO!
kx
ekx + kxekx

ke 2kx
W = kx =e + kxe2kx − kxe2kx = e2kx 6= 0 ∀k, x ∈ R
e xekx
Page 74 of 117
Ejercicio 57
Solve
dy 4 √
= y+x y
dx x
Solución:
dy 4 √
= y+x y
dx x
dy 4 1 1 1 1 dy dv
+y − = xy 2 ,ec. de Bernoulli, n = → sea v = y 1− 2 → v = y 2 → = 2v
dx x 2 dx dx

dv 4 1
→ 2v + v2 − = xv / ·
dx x 2v

dv 2 x R 2
+v − = ,ec. lineal → sea ρ (x) = e − x dx = e−2 ln x = x−2
dx x 2

dv 2 x
→ +v − = / · x−2
dx x 2
dv −2 1
· x + v −2 · x−3 =

dx
| {z } 2x
d −2 )
dx (vx
Z
d 1
vx−2 =

→ / ()
dx 2x
Z
1 1
vx−2 = dx
2 x
1
vx−2 = ln x + c
2
v = x2 (ln x + c)
1 2
∴ y 2 = x2 (ln x + c) ∴ y = x4 (ln x + c)
Page 75 of 117
Ejercicio 58
Sea la siguiente ecuacion diferencial, en donde a y b son constantes con a distinto de b:
y 00 + (a + b) y 0 + aby = 0
A) Demuestre que la sustitucion

z = y 0 + ay
Reduce el orden de la ecuacion dada, obtenga la solucion general.
B) Aplique el procedimiento anterior para obtener la solucion general de:
y 00 + 2ky 0 + k 2 y = 0
Solución:
z0
z 0 = −bz ⇒ = −b ⇒ ln(z) = −bx + ln(c) ⇒ z = Ce−bx
z
(Obviando las constantes de integración, que añadiremos al final)
Ası́:
y 0 + ay = e−bx
de la ec. homogénea:
y 0 + ay = 0 ⇒ y = Ke−ax
Luego, suponemos como solución:
y(x) = u(x)e−ax ⇒ y 0 = u0 e−ax − aue−ax
reemplazando en la ecuación, se llega a:
e(a−b)x
u0 = e(a−b)x ⇒ u =
a−b
Finalmente, la solución general de la ecuación diferencial queda:
ke−bx
y(x) = , k = cte.
a−b
Page 76 of 117
Ejercicio 59
Un pueblo llamado Melosova, que conlinda con Ponmelov tenı́a una población de 5000 habitantes el año 2000
y una población de 7000 almas el año 2007. Suponiendo que su población siga creciendo exponencialmente
con su ı́ndice constante, ¿Cuál será la población de Melosova el año 2050?
Solución:
Sea y el numero de habitantes presentes en el tiempo t, y k la constante de proporcionalidad.

dy dy
→ = ky → = k · dt , luego integramos
dt y
→ ln y = k · t + c → y = cekt
Pero tenemos las siguientes condiciones de borde y(0) = 5000 y y(7) = 7000
→ 5000 = cek·0 ∴ c = 5000
7 7
→ 7000 = 5000 · ek·7 → = e7k → ln = 7k →∴ k ≈ 0, 04806
5 5
∴ el numero de habitantes durante el tiempo transcurrido desde el 2000 esta definido por
y = 5000 · e0,04806t
∴ En el añ o 2050 (es decir 50 años después) habran y = 5000 · e0,04806·50 ≈ 55,281 habitantes.
Ejercicio 60
(x + ey )y 0 = xe−y − 1
Solución:
x
xe−y − 1 x − ey 1 −1 1
y0 = = · = ey
x ·
x + ey x + ey ey ey + 1 ey
x x 0 u u − xu0
u= ⇒ y = ln ⇒ y = ·
ey u x u2
]
u u − xu0 u − xu0

u−1 u u−1 u−1
· 2
= · ⇒ 2
= ⇒ xu0 = u − u2
x u u+1 x u u+1 u+1
Z Z
1
⇒− du = xdx
u (u + 1)
x2
Z
1 1 u+1 1 x2 1 x
− − du = ln = + C ⇒ 1 + = Ke 2 ⇒ u = x2
= y
u u+1 u 2 u Ke 2 − 1 e
x2 x2
h i
∴ ey = x Ke 2 − 1 ⇒ y = ln x Ke 2 − 1
Page 77 of 117
Ejercicio 61
Resolver la EDO (todo el curso es asi, no se asombren) utilizando la sustitución y = z n para algún n
apropiado.
2x4 yy 0 + y 4 = 4x6
Solución:
Sea y = z n , entonces y 0 = nz n−1 z 0 , sustituyendo
2nx4 z 2n−1 z 0 + z 4n = 4x6
esta ecuacion es homogenea si el exponente de cada monomio es 6, es decir

2n − 1 = 2 , 4n = 6 ; que es cierto para n= 23
reemplazando n
3x4 z 2 z 0 + z 6 = 4x6
4x6 − z 6
z0 =
3x4 z 2
z 6
4− z
z0 = x
z 2 , haciendo u =
x
3
x
4 − u6
xu0 + u =
3u2
4 − u6 − 3u3
xu0 =
3u2
3u2 dx
du =
4 − u6 − 3u2 x
3u2 du dx
3 3
=
(1 − u )(u + 4) x
−3 2 3 2
 
u u dx
 5 5
+ 3  du =

 3
u −1 u +4 x
1 1
ln(u3 + 4) − ln(u3 − 1) = ln(x) + C
5 5
3 1
u +4 5
= Cx
u3 − 1
retornando a
 51
y2

+4
 x3
 = Cx

 2
y
−1
x3
15
y 2 + 4x3

= Cx
y 2 − x3
Page 78 of 117
Ejercicio 62
y 00 + 3y 0 − 10y = 6e4x
Solución:
La solucion homogénea
y 00 + 3y 0 − 10y = 0
tiene la forma emx donde la ecuacion caracteristica es
m2 + 3m − 10 = 0, m = −5ym = 2
La solucion complementaria es
y = Ae−5x + Be2x
La solucion particular tiene la forma y = Ce4x , reemplazando en la EDO
16Ce4x + 12Ce4x − 10Ce4x = 6e4x

1
de esto se obtiene C = La solucion es
3
1
y = Ae−5x + Be2x + e4x
3
Ejercicio 63
p
xy 0 = y − x2 + y 2
Solución:
Notemos
r que reordenando
y y 2 y
0
y = − 1− , luego sea u =
x x x
p
xu0 + u = u − 1 − u2
du dx
√ =−
1−u 2 x
sen−1 (u) = −ln(x) + C
Volviendo a la variable original
y
sen−1 = −ln(x) + C
x
Page 79 of 117
Ejercicio 64
xy 0 − 2y = x5
Solución:
xy 0 − 2y = x5 ⇒
xdy − x5 + 2y dx = 0 ⇒

dy x5 + 2y
= ⇒
dx x
Tomando y = u4 x ⇒ dy = 4xu3 du + u4 dx ⇒
4xu3 du + u4 dx
= x4 + 2u4 ⇒
dx
du x4 + u4
= ⇒
dx 4xu3
Tomando u = zx ⇒ du = zdx + xdz ⇒
zdx + xdz 1 + z4
= ⇒
dx 4z 3
4z 3 dx
dz = ⇒
1 − 3z 4 x

ln 1 − 3z 4
− = ln (x) + C ⇒
3
y
ln 1 − 3 5 = −3 ln (x) − 3C
x
Otra forma!
dy x5 + 2y y
xy 0 − 2y = x5 → = = x4 + 2
dx x x
y dy dv
Sea v = → = x+v
x dx dx
dv dv v
→ x + v = x4 + 2v → = x3 +
dx dx x
v dv dt
Sea t = → = x+t
x dx dx
dt
→ x + t = x3 + t → dt = x2 dx , luego integrando
dx
x3
→t= +c
3
v x3 y x3
∴ = +c→ 2 = +c
x 3 x 3
1
∴ y = x5 + cx2
3
Page 80 of 117
Ejercicio 65
En la conservación de alimentos, el azúcar de caña sufre un proceso de inversión y se transforma en glucosa
y fructosa. En soluciones diluidas el ritmo de inversión es proporcional a la concentración del azucar. Si la
concentración al inicio es de 1/50 gr/L y 1/200 gr/L tras tres horas, halle la concentración de azúcar después
de 6 horas.
Solución:
El modelo es sencillo
dy
= kdt
y
Z Z
dy
= kdt → ln |y| = kt + c ⇒ y = ekt+c
y
Que con la info, se transforma en 2 ecuaciones con 2 incógnitas

1 c

=e
50
ln 14
k=
3

1 1
y (6) = e2 ln 4 ·
50
1 1
e2 ln 4 · → Con la propiedad ab = eb ln a

50
2
1 1 1 1 1
· = · =
4 50 16 50 800
1
y (6) =
800
Page 81 of 117
Ejercicio 66
Considere la siguiente EDO:

α(t)
y 0 tanh (y) + cosh(y) = β (t) tq α (t) y β (t) son funciones continuas, resolver la EDO.
Obtener la solucion explicita para α (t) = β (t) = 2t.
Solución:
α (t)
y 0 tanh (y) + = β (t)
cosh (y)
dy sinh (y) α (t)
⇔ · + = β (t)
dt cosh (y) cosh (y)
dx dy
Sea x = cosh (y) ⇒ = senh (y)
dt dt
ergo :
dx 1 α (t)
⇔ · + = β (t)
dt x x
dx
⇔ + α (t) = β (t) x
dt
dx
⇔ − β (t) x = −α (t)
dt Z Z Z
dx
⇔ exp − β (t) dt − β (t) exp − β (t) dt x = −α (t) exp − β (t) dt
dt
Z Z
d
⇔ exp − β (t) dt x = −α (t) exp − β (t) dt
dt
Z Z Z
⇒ exp − β (t) dt x = − α (t) exp − β (t) dt dt + C
Z Z Z
⇔ x = exp β (t) dt − α (t) exp − β (t) dt dt + C
Como x = cosh (y)

Z Z Z
⇔ y = arc cos h exp β (t) dt − α (t) exp − β (t) dt dt + C
Si α (t) = β (t) = 2t
Z
y = arc cos h exp t2 −2 t exp −t2 dt + C

Z Z
t exp −t2 dt |{z}= −2 exp (u) du = 2 exp (u) = 2 exp −t2

y
u=−t2
ergo :
⇔ y = arc cos h exp t2 −4 exp −t2 + C

⇔ y = arc cos h C exp t2 − 4


Por último
x = cosh (y)
ey + e−y
x=
2
2xey = e2y + 1
e2y − 2xey + 1 = 0
2
(ey − x) = x2 − 1
p
ey − x = ± x2 − 1
n p o
y = ln x ± x2 − 1
2
luego con x = Cet − 4
q
t2 2
2
y = ln Ce − 4 ± t
Ce − 4 − 1
Page 83 of 117
Ejercicio 67
Z t
f (u)
1+t= √ du
0 t−u
Solución:
Para resolver la ecuación integral, se debe usar la Transformada de Laplace y el producto de convolución,
además de considerar que:
Γ(x + 1)
L [tx ] (s) =
sx+1
En efecto, al aplicar la transformada en ámbos miembros, considerando que la integral es la convolución
de las funciones f (t) y g(t) = √1t se tiene:

1
L [1 + t] (s) = L f (t) ∗ √ (s)
t

1 1 1
+ = L [f (t)] (s)L √ (s)
s s2 t
Sea L [f (t)] (s) = F (s), ası́:
1 1 Γ( 1 )
+ 2 = F (s) 12
s s s2
r
1 1 π
+ = F (s)
s s2 s
1 1
F (s) = √ 1 +√ 3
πs πs 2
2
√ √
1 π 2 π
F (s) = +
π s 12 π 2s 23
1 Γ(− 21 + 1) 2 Γ( 12 + 1)
F (s) = 1 + 3
π s2 π s2
Finalmente, al aplicar la transformada inversa:
1 −1 2 1
f (t) = t 2 + t2
π π
Page 84 of 117
Ejercicio 68
Dado un cable de longitud L y masa m, sostenido de dos puntos a igual altura y a distancia H, encontrar la
ecuación de su deformada, suponiendo que la única fuerza externa que actúa es el propio peso del cable.
Consideremos un cable que cumple con las condiciones descritas. Para mayor comodidad, llamamos
”w”la carga uniformemente distribuı́da en el cable (es decir, por la acción de su propio peso).
Sea ç.el punto más bajo del cable.
Ampliando lo que ocurre en el pedazo de cable entre ç ”b”, y anotando las fuerzas participantes, se
2
obtiene algo como la imagen mostrada en el recuadro bajo la figura principal.

Donde ”T. la tensión, ”V.es la tensión horizontal y W=wu es la carga distribuı́da que afecta al trozo
es
entre ç ”b”. α es el ángulo que forma ”Tçon la horizontal.

2
Realizando el triángulo de fuerzas se tiene, sin complicaciones por Pitágoras que:

p p
T (u) = V 2 + W 2 = V 2 + w2 u2
V
Llamemos a
= t; entonces
w
p p
T (u) = t2 w2 + w2 u2 = w t2 + u2
Ahora, escogiendo un diferencial en dirección del eje X (la horizontal) se cumple que: dx = du cos(α)
V wt du
⇒ dx = du cos(α) = du · = √ du = r
T 2
w t +u 2
u2
1+ 2
t
Eligiendo un origen ”tünidades más abajo del pto. ç”, quedando el pto. c(0,t) y el pto. b(x,y). Ası́
entonces, integrando la relación anterior, se tiene:
r
u2
Z
du u
x= = t ln + +1
t2
r
u2 t
1+
t2
Cálculo de la integral:
Z
du
r
u2
1+
t2
u du
Sea: z = ⇒ dz =
t t
Entonces:
Z Z
du dz
r =t √ +C
u 2 z2 + 1
1+ 2
t

Y la integral es directa √
y se llega al resultado mostrado.
Y se sabe que ln(a + a2 + 1) = sinh−1 (a) entonces:

−1 u x
x = t sinh ⇒ u = t sinh
t t
De modo similar a cuando notamos que:
dx = du cos(α)
podemos notar entonces que:
dy = du sin(α)
.
Pero
dx
du =
cos(α)
, luego:
wu wu u
dy = dx tan(α) = dx = dx = dx
V wt t
Y recordando que como ya obtuvimos:

u x
= sinh
t t
. Por tanto:
Z y Z x
x x
dy = sinh dx ⇒ dy = sinh dx ⇒
t t 0 t
x
x x x
⇒ y − t = t cosh = t cosh − 1 ⇒ y = t cosh
t 0 t t
Esta es la famosa ecuación de la catenaria.
Page 86 of 117
Ejercicio 69
dS dI dR
(Propuesto) = −βSI, = βSI − γI, = γI
dt dt dt
Sea el anterior modelo de epidemias de 3 ODEs, donde S(t) son personas susceptibles a enfermar, I(t)
infectados y R(t) sanos y que no pueden enfermar. Plantee variaciones del modelo para cuando: a) Hay
muertes y Nacimientos, b) Los individuos sanos pueden enfermar de nuevo, c) Si hay personas portadoras
(no infectan, ni muestran signos) d) Al modelo en (c) agrege a los recobrados e) Si hay un nuevo grupo que
son potenciales portadores. f) Extras! 1
Notación:
S - Individuos susceptibles (ver modelo SIR).
I - Individuos infectados (ver modelo SIR).
R - Individuos recobrados (ver modelo SIR).
N - Población total.
β - Tasa de contagios (probabilidad de que una persona enferme al estar en contacto con un infectado).
1/γ - Tiempo promedio de infección (para un solo individuo).
µ - Tasa promedio de defunciones (probabilidad de que un individuo infectado muera debido a la
enfermedad).
f - Tasa promedio de pérdida de inmunidad en individuos recobrados.
E - Tasa promedio de nacimientos.
1/ Tiempo promedio de incubación.
B - Tasa promedio de nacimientos.
M - Infantes con inmunidad pasiva.
δ - Tiempo promedio de inmunidad temporal.
Modelo SIR con nacimientos y muertes Éste es una ampliación del modelo SIR y considera nacimientos a
lo largo del tiempo t, de modo que hay una constante renovación de individuos susceptibles a la
enfermedad. Del mismo modo, algunos individuos mueren o bien por la enfermedad o de manera natural
al cabo de un tiempo de haber sido recobrados. El flujo de transiciones de un grupo a otro sigue siendo:
S → I → R Las ecuaciones diferenciales que representan el modelo son:
dS
= −βSI + µ(N − S)
dt
dI
= βSI − γI − µI
dt
dR
= γI − µR
dt
Modelo SIS con nacimientos y muertes Simplificación del modelo SIR con nacimientos y muertes, en
donde un individuo recobrado nunca desarrolla inmunidad a la enfermedad, de modo que el flujo de un
grupo a otro resulta: S → I → S Sus ecuaciones diferenciales son:
dS
= −βSI + µ(N − S) + γI
dt
dI
dt
1 Si, me fui en volada, y que? Acostumbrense Haha

Modelo SIRS Extensión del modelo SIR, en donde los individuos recobrados pueden perder la inmunidad
a la enfermedad y volver a formar parte del grupo de susceptibles. El flujo de un grupo a otro está dado
como sigue: S → I → R → S Sus ecuaciones diferenciales son:
dS
= −βSI + µ(N − S) + f R
dt
dI
dt
dR
= γI − µR − f R
dt
Modelo SEIS El modelo SEIS considera una nueva clase de individuos E (del inglés exposed), es decir,
aquellos que portan la enfermedad pero que al hallarse en su periodo de incubación no muestran
sı́ntomas y no están en condición de infectar a otros. Se supone además que un individuo que ha
enfermado nunca obtiene inmunidad, de modo que el flujo de grupos resulta: S → E → I → S y sus
ecuaciones diferenciales están dadas por:
dS
= B − βSI − µS + γI
dt
dE
= βSI − ( + µ)E
dt
dI
= E − (γ + µ)I
dt
Modelo SEIR Derivado del modelo SEIS pero agregando esta vez a la población de recobrados. En este
caso, al igual que los anteriores, cada grupo es mutuamente excluyente y la suma de todos es la
población total, esto es, N = S(t) + E(t) + I(t) + R(t) El flujo de un grupo a otro es: S → E → I → R y
sus ecuaciones diferenciales resultan:
dS
= B − βSI − µS
dt
dE
= βSI − ( + µ)E
dt
dI
= E − (γ + µ)I
dt
dR
= γI − µR
dt

Modelo MSIR Éste modelo considera una nueva clase de individuos M, los llamados infantes con
inmunidad pasiva que tras cierto tiempo la pierden y están en condición de ser portadores de la
enfermedad (para una descripción más detallada véase Inmunidad (medicina)). Los flujos de grupo
resultan: M → S → I → R y sus ecuaciones diferenciales:
dM
= B − δM S − µM
dt
dS
= δM S − βSI − µS
dt
dI
dt
dR
= γI − µR
dt
Modelo MSEIR Derivado del modelo anterior, en donde entre el periodo de susceptibilidad e infección
hay un periodo de latencia, de modo que el flujo entre un grupo y otro es: M → S → E → I → R y sus
ecuaciones diferenciales:
dM
= B − δM S − µM
dt
dS
= δM S − βSI − µS
dt
dE
= βSI − ( + µ)E
dt
dI
= E − (γ + µ)I
dt
dR
= γI − µR
dt
Page 89 of 117
Ejercicio 70
y 00 + 4y 0 + 4y = cosh x.
y 00 + 4y 0 + 4y = 0. Del polimonio caracteristico:

λ2 + 4λ + 4 = (λ + 2)2 = 0 λ = −2. Como debe tener 2 soluciones, haciendo u1 = e−2x , y u2 = c1 (x)e−2x
obtendremos que u2 = xe−2x , ahora:
−2x
xe−2x

e

−2e−2x −e −2x
(2x − 1)
= −e−2x e−2x (2x − 1) + 2xe−2x e−2x = −e−4x (2x − 1) + 2xe−4x = (−2x + 1 + 2x)e−4x = e−4x .
Como el Wroskiano es distinto que cero, las soluciones son L.I.
Ahora tomando A(x) y B(x) tal que A(x)u1 + B(x)u2
Z Z
1 1
A(x) = − u2 (x)b(x) dx, B(x) = u1 (x)b(x) dx
W W
con b(x) = coshx
Z Z
1 1 x
A(x) = − −4x
xe−2x cosh x dx = − xe2x cosh x dx = − e (9(x − 1) + e2x (3x − 1)) + C1
e 18
Z Z
1 1 x
B(x) = −4x
e−2x cosh x dx = e2x cosh x dx = e (3 + e2x ) + C2
e 6
Page 90 of 117
Ejercicio 71
y 00 − 3y 0 + 2y = 4x + 4
Vemos la homogenea primero, y 00 − 3y 0 + 2y = 0, y calculando el polimonio caracteristico
yh = c1 ex + c2 e2x
Calculamos el Wroskiano quedando : W (ex , e2x ) = e3x

Usamos variaciones de parametros
−y2 f −e2x 4(x + 1)

Z Z
u1 = dx = dx = 4e−x (x + 2)
W (ex , e2x ) e3x
ex 4(x + 1)
Z Z
y1 f
u2 = dx = dx = −e2x (2x + 3)
W (ex , e2x ) e3x
yp = y1 u1 + y2 u2 = ex (e−x 4(x + 2)) + e2x (−e−2x (2x + 3)) = 4(x + 2) − (2x + 3) = 2x + 5
y(x) = yh (x) + yp (x) = c1 ex + c2 e−2x + 2x + 5
Cabe destacar que en este caso, era mas fácil usar coeficientes indeterminados, ya que la función 4x + 4
es sencilla.
Page 91 of 117
Ejercicio 72
y 00 + y = csc(x)
We see that a fundamental solution set for the corresponding homogeneous equation are: y1(x) = cos(x)
y y2(x) = sin(x) And so we take as our particular solution:
yp (x) = v1(x)cos(x) + v2(x)sin(x)
Next, we want to solve the system:
cos(x)v 0 1 + sin(x)v 0 2 = 0
−sin(x)v 0 1 + cos(x)v 0 2 = csc(x)

From the first equation, we find:
v 0 1 = −tan(x)v 0 2
And so substituting into the second equation, we obtain:
sin2(x)cos(x)v 0 2 + cos(x)v 0 2 = csc(x)
Multiply through by cos(x):
sin2(x)v 0 2 + cos2(x)v 0 2 = cot(x)
Apply a Pythagorean identity on the left:
v 0 2 = cot(x)
Hence, this implies:
v 0 1 = −1
Integrating, we obtain:
v1 = −x
v2 = ln|sin(x)|
And so out particular solution is:
yp(x) = −xcos(x) + sin(x)ln|sin(x)|
And then by superposition, we obtain the general solution:
y(x) = yh(x) + yp(x) = Acos(x) + Bsin(x) − xcos(x) + sin(x)ln|sin(x)|
Page 92 of 117
Ejercicio 73
x2 y 00 − 3xy 0 + 4y = x2 lnx
Este es largo. Con y = xm , tendremos que m = 2, pero necesitamos otra

Haciendo u2 (x) = c1 (x)x2 , quedará c1 (x) = ln(x)
Ahora, normalizando la ODE
3 0 4
y 00 − y + 2 y = lnx
x x
Usando variacion de parámetros , calculamos el Wroskiano, y vemos que son L.I.
W = x3 + 2x3 lnx − 2x3 lnx = x3
Ya, lo de siempre.
(lnx)3
u1 = −
3
(lnx)2
u2 = −
2
Finalmente
1
y = c1 x2 + c2 x2 lnx + (lnx)3
6
Ejercicio 74
Pruebe que y1(t) = 1 − t y y2(t) = t3 no pueden ser soluciones de
y 00 + p(t)y + q(t) = 0
, con p(t), q(t) continuas en [−1, 5]
Tomamos el Wroskiano, y notar que en t = 1, t = 0 son L.D.
W (y1, y2) = (1 − t)(3t2 ) − (t3 )(−1) = 3t2 − 2t3
Page 93 of 117
Ejercicio 75
y 00 + y = t cos t.
(propuesto con C.I.)
La solución particular de es la forma Pn eαt cos βt, por coeficientes indeterminados.

with n=1, α = 0, and β = 1.
Since α + iβ = i is a simple root of the characteristic equation
λ2 + 1 = 0
we should try a particular integral of the form
yp = t[F1 (t)eαt cos βt + G1 (t)eαt sin βt]

= t[F1 (t) cos t + G1 (t) sin t]
= t[(A0 t + A1 ) cos t + (B0 t + B1 ) sin t]
= (A0 t2 + A1 t) cos t + (B0 t2 + B1 t) sin t.
Substituting yp into the differential equation, we have the identity
t cos t = yp00 + yp
= [(A0 t2 + A1 t) cos t + (B0 t2 + B1 t) sin t]00
+ [(A0 t2 + A1 t) cos t + (B0 t2 + B1 t) sin t]
= [2A0 cos t + 2(2A0 t + A1 )(− sin t) + (A0 t2 + A1 t)(− cos t)]
+ [2B0 sin t + 2(2B0 t + B1 ) cos t + (B0 t2 + B1 t)(− sin t)]
+ [(A0 t2 + A1 t) cos t + (B0 t2 + B1 t) sin t]
= [4B0 t + (2A0 + 2B1 )] cos t + [−4A0 t + (−2A1 + 2B0 )] sin t.
Comparing both sides, we have
4B0 =1
2A0 +2B1 =0
−4A0 =0
−2A1 +2B0 =0
which has the solution A0 = 0, A1 = 1/4, B0 = 1/4, B1 = 0. We then have a particular integral
1 1
yp = t cos t + t2 sin t.
4 4
Page 94 of 117
Ejercicio 76
dy
= t2 − y
dt
t2 , is a polynomial of degree 2, so we look for a solution using the same form,
yp = At2 + Bt + C
, where
dyp
= 2At + B
dt
Plugging this particular integral with constants A, B, and C into the original equation yields,
2At + B = t2 − (At2 + Bt + C)
, where
t2 − At2 = 0, −Bt = 2At
and
−C = B
Replacing resulting constants,
yp = t2 − 2t + 2
To solve for the general solution,
y = yp + yc
where yc is the homogeneous solution yc = c1 e−t , therefore, the general solution is:
y = t2 − 2t + 2 + c1 e−t
Page 95 of 117
Ejercicio 77
y 00 + 4y = 3sin2x
r2 + 4 = 0, r = ±2i
yc = c1 cos(2x) + c2 sin(2x)
There is a linearly dependent term in yc to Q(x) in the original DE. Multiplying by a polynomial of
order 1 gives:
yp = (Ax + B) sin(2x), u = Ax + B, u0 = A, u00 = 0
yp0 = u0 sin(2x) + u cos(2x)

yp00 = u00 sin(2x) + u0 cos(2x) + u0 cos(2x) − u sin(2x)
u00 sin(2x) + 2u0 cos(2x) − u sin(2x) + 4u sin(2x)
u00 sin(2x) + 2u0 cos(2x) + 3u sin(2x)
(0) sin(2x) + 2A cos(2x) + 3(Ax + B) sin(2x)
2A cos(2x) + 3(Ax + B) sin(2x) = 3 sin(2x)
Ejercicio 78
y 00 + y = 2xsin(x)
yp = (Ax2 + xB)cosx + (Cx2 + xD)sinx
2A + 2D = 0
−2B + 2C = 0
−4A = 2
4C = 0
C = B = 0, A = −D = −1/2
y(x) = c2 sin(x) + c1 cos(x) − 1/2x2 cos(x) + 1/2xsin(x)
Page 96 of 117
Ejercicio 79
2y 00 + 3y 0 + y = t
y = At + B, y 0 = A, y 00 = 0
2y 00 + 3y 0 + y = t =⇒ 2(0) + 3A + (At + B) = t
=⇒ At + (3A + B) = t
A = 1(At = t)yB = −3(3A + B = 0).
Ejercicio 80
y 00 + y = sec t
The auxilliary equation is: m2 + 1 = 0 with roots m = ±i
The homogeneous solution is: y = C1 cos t + C2 sin t

Differentiate: y 0 = -C1 sin t + C2 cos t ⇒ C10 cos t + C20 sin t = 0 [1]
..and again: y 00 = -C1 cos t − C2 sin t ⇒ -C10 sin t + C20 cos t = sec t [2]
Solve the system of equations: C10 = − tan t, C20 = 1
...Hence: C1 = ln | cos t|, C2 = t
Therefore, the particular solution is: y = ln | cos t| · cos t + t sin t
The complete solution is: y = C1 cos t + C2 sin t + ln | cos t| · cos t + t sin t
Page 97 of 117
Ejercicio 81
dy
= y + ex .
dx
Vemos que la solución particular es de la forma yp = Axex ., ya que la solución homogenea es L.D. con la
parte no homonegea de la ecuación.
d
(Axex ) = Axex + ex
dx
Axex + Aex = Axex + ex
A = 1.
La solución general es:
y = c1 ex + xex .
Ejercicio 82
y 00 + 4y = xex + xsin2x
y 00 + 4y = 0
r2 + 4 = 0
r2 = −4
r = +/ − 2i
Y c = C1(cos(2x) + isin(2x)) + C2(cos(2x) − isin(2x))
For the particular solution:
Y p = Aex + Bxex + x(Ccos(2x) + Dxcos(2x) + Esin(2x) + F xsin(2x))
Y p = Aex + Bxex + Cxcos(2x) + Dx2 cos(2x) + Exsin(2x) + F x2 sin(2x)

y 0 = Aex + Bex + Bxex + Ccos(2x) − 2Cxsin(2x) + 2Dxcos(2x) − 2Dx2 sin(2x) + Esin(2x) +
2Excos(2x) + 2F xsin(2x) + 2F x2 cos(2x)
y 00 =
Aex +Bex +Bex +Bxex −2Csin(2x)−2Csin(2x)−4Cxcos(2x)+2Dcos(2x)−4Dxsin(2x)−4Dxsin(2x)−
4Dx2 cos(2x)+2Ecos(2x)+2Ecos(2x)−4Exsin(2x)+2F sin(2x)+4F xcos(2x)+4F xcos(2x)−4F x2 sin(2x)
y 00 + 4y = Aex + Bex + Bex + Bxex − 2Csin(2x) − 2Csin(2x) − 4Cxcos(2x) + 2Dcos(2x) − 4Dxsin(2x) −
4Dxsin(2x)−4Dx2 cos(2x)+2Ecos(2x)+2Ecos(2x)−4Exsin(2x)+2F sin(2x)+4F xcos(2x)+4F xcos(2x)−
4F x2 sin(2x) + 4Aex + 4Bxex + 4Cxcos(2x) + 4Dx2 cos(2x) + 4Exsin(2x) + 4F x2 sin(2x)] = xex + xsin(2x)
Page 98 of 117
Equate the coefficients of ex from both sides
5A + 2B = 0
Equate the coefficients of xex from both sides
5B = 1
B = 1/5
5A + 2(1/5) = 0
A = −2/25
Equate the coefficients of sin(2x) from both sides
−4C + 2F = 0
Equate the coefficients of x sin(2x) from both sides
−8D − 4E + 4E = 1
−8D = 1
D = −1/8(∗)
Equate the coefficients of cos(2x) from both sides
2D + 4E = 0
2(−1/8) + 4E = 0
E = 1/16(∗)
Equate the coefficients of x cos(2x) from both sides
−4C + 8F + 4C = 0
F = 0, C = 0(∗)
Equate the coefficients of
x2 sin(2x)
from both sides = 0
A = −2/25, B = 1/5, C = 0, D = −1/8, E = 1/16, F = 0
y = C1(cos(2x)+isin(2x))+C2(cos(2x)−isin(2x))−(2/25)ex +(1/5)xex −(1/8)x2 cos(2x)+(1/16)xsin(2x)
Page 99 of 117
Ejercicio 83
d2 y dy
2
−2 − 3y = 0, y(0) = 0, y 0 (0) = 1.
dx dx
m2 − 2m − 3 = 0
(m − 3)(m + 1) = 0
y(x) = Ae3x + Be−x
0 = Ae3·0 + Be0 = A + B = 0
y 0 (x) = 3Ae3x − Be−x
1 = 3A − B
1 1
A= , B=− .
4 4
Ejercicio 84
2y 00 + 6y 0 + 5y = 0, y(0) = 1, y 0 (0) = −1/2
√
−6 ± 62 − 4 · 2 · 5
m1,2 = .
2·2
3 1
m1,2 = − ± i
2 2
3
x x
y(x) = e− 2 x A cos + B sin .
2 2
x
0 3 −3x x x
− 3
x A x B
y (x) = − e 2 A cos + B sin +e 2 − sin + cos .
2 2 2 2 2 2 2
y(0) = 1
y(0) = e0(Acos(0) + Bsin(0)) = A = 1,
1
y 0 (0) = −
2

0 3 0 0 1 B
y (0) = − e (1 · cos(0) + B sin(0)) + e − sin(0) + cos(0) (14)
2 2 2
3 B 1
= − + =− . (15)
2 2 2
B=2
3
x x
y(x) = e− 2 x cos + 2 sin .
2 2

Ejercicio 85
d2 y dy
+ − 2y = 0, y(0) = 1, y(∞) = 0
dx2 dx
y(x) = Aex + Be−2x
y(0) = Be0 = B = 1,
A=0
y(x) = e−2x .
Ejercicio 86
d2 y
+ 25y = 2 sin(5t) + cos(t).
dt2
Resolviendo la homogénea
m2 + 25 = 0
yc = A cos(5t) + B sin(5t)
Ahora, para C.I., es complicado, ya que la homogenea es igual forma que el lado izquierdo de la ODE,
por lo que debemos utilizar la familia de soluciones de este tipo.
yp = At cos(5t) + Bt sin(5t) + C cos(t) + D sin(t)
yp0 = A cos(5t) + B sin(5t) − 5At sin(5t) + 5Bt cos(5t) − C sin(t) + D cos(t), (16)
yp00 = −10A sin(5t) + 10B cos(5t) − 25At cos(5t) − 25Bt sin(5t) − C cos(t) − D sin(t). (17)
Despejando en el ODE
−10A sin(5t) + 10B cos(5t) + 24C cos(t) + 24D sin(t) = 2 sin(5t) + cos(t)
1 1
B = 0, D = 0, A = − , C =
5 24
Quedando finalmente
1 t
yp = cos(t) − cos(5t)
24 5
1 t
y = yc + yp = A cos(5t) + B sin(5t) + cos(t) − cos(5t)
24 5
Page 101 of 117

Ejercicio 87
Probar que toda función periodica, con T > 1 tiene el Laplaciano de la forma
RT
e−st F (t) dt
0
L{F (t)} = dt
1 − e−sT
Z ∞ ∞ Z
X (k+1)T ∞ Z
X T
e−st F (t) dt = e−st F (t) dt = e−s(t+kT ) F (t + kT ) dt.
0 k=0 kT k=0 0
Z T Z T ∞
X
−st −sT −2sT
F (t)e−st dt · e−nsT =

L{F (t)} = F (t)e dt · 1 + e +e + ... =
0 0 n=0
Z T
1
= F (t)e−st dt
1 − e−sT 0
Ejercicio 88
(New one) Determine el orden mı́nimo que debe tener la ODE asociada para que la solución sea:
y(x) = c1 t3 + c2 tet sin(t)
Octavo orden Aplicando los aniquiladores
D4 ((D − 1)2 + 1)2 y(x) = 0.
Page 102 of 117

Ejercicio 89
Una función no negativa f (x) tiene la propiedad que su conjunto de ordenadas en un intervalo cualquiera
tiene un área proporcional a la longitud de arco de la gráfica correspondiente en el intervalo. Encuentre f (x).
Sea f : I → R de clase C2 , tal que f (x) > 0. Luego, el área bajo la curva en algún intervalo cualquiera
[x0 , x1 ], está dada por
Z x1
A= f (t)dt
x0
Por otro lado como f es rectificable, porque es de de clase C2 , la longitud de curva en tal intervalo es:
Z x1 p
L= 1 + (f 0 (t))2 dt
x0
De esta forma como A ∼ L:

Z x1 Z x1 p
f (t)dt = k 1 + (f 0 (t))2 dt
x0 x0
Considerando que la integral es una aplicación lineal:

Z x1 Z x1 p
f (t)dt = k 1 + (f 0 (t))2 dt
x0 x0
Luego como ésto es válido en cualquier intervalo, los integrandos deben ser iguales, es decir:
p
f (t) = k 1 + (f 0 (t))2 (18)
Luego sea y = f (t), sustituyendo, elevando al cuadrado, derivando y cancelando se llega a:

y
y 00 − =0
k2
De donde se deduce que la función f debe satisfacer la siguiente forma:
x x
f (x) = c1 cosh( ) + c2 sinh( )
k k
Comprobando solo c1 cosh( xk ), en (1), ya que la otra solución L.I no siempre es mayor que cero y la
comprobación en (1) debe ser por separado, pues el principio de superposición no se cumple, se deduce
que:
x
f (x) = k cosh( ) k>0
k
Page 103 of 117

Ejercicio 90
Resolver
x3 y 000 − x2 y 00 + 2xy 0 − 2y = x3
ez = x
dy dy dz dy
y0 = = = e−z
dx dz dx dz
dy
xy 0 =
dz
d2 y dy
x2 y 00 = −
dz 2 dz
Si llamo al operador lineal
d
D=
dz
dando
x2 y 00 = D(D − 1)y
x3 y 000 = D(D − 1)(D − 2)y

D(D − 1)(D − 2)y − D(D − 1)y + 2Dy − 2y = e3z
Resuelvo la homogenea
yc = c1 zez + c2 ez + c3 e2z
yc = c1 x log x + c2 x + c3 x2
y la particular es del tipo (por coef indeterminados)
yp = Ax3 + Bx2 + Cx + D
Ejercicio 91
Encontrar la ecuación diferencial no homogenea
y = c1 + c2 x + c3 e2x cos x + c4 e2x sin x − xe−x
y 0000 − 4 y 000 + 5 y 00 = f (x)
Page 104 of 117

Ejercicio 92
Considere el siguiente problema de valores iniciales:
p
1 + y(x)
y 0 (x) = , y(a) = b
1 + x2
Encuentre todos los valores de a, b ∈ R para los cuales el PVI no tiene solucion. Encuentre todos los valores
de a, b ∈ R para los cuales el PVI tiene más de una solución y determı́nelas.
Sabemos que f (x, y) es continua para y ≥ −1 y x ∈ R, entonces en esta region existe solucion. Luego:
∂f 1 1
(x, y) = p ·
∂y 2 1 + y(x) 1 + x2
Y esta es continua en y > −1 y x ∈ R. Luego en esta region, por Teorema de Existencia y Unicidad,
existe solucion única.
Para determinar la solución, vemos que es de variables separables:
dy dx p
√ = 2 1 + y = arctan(x) + C
1+y 1 + x2
2
arctan(x)
y(x) = +C −1
2
Luego, la solución al PVI es:
2
√

arctan(x)
y(x) = + 2 1 + b − arctan(a) − 1
2
Notemos que para b = −1 y a ∈ R, el teorema de existencia y unicidad no aplica. Verificamos que la

solución singular y(x) = −1 es solución al PVI. Entonces poseemos más de 1 solución:
2
arctan(x)
y(x) = −1 & y(x) = − arctan(a) − 1
2
Ejercicio 93
Linealizar
dy dz
= −3y + yz = z + y2
dt dt
Hartman–Grobman, aproximar por el Jacobiano evaluado en los puntos crı́ticos por las variables
respectivas mas el valor inicial.
dy dz
= −3y =z
dt dt
Page 105 of 117

Ejercicio 94
Consideremos el siguiente problema con valor inicial, y(a) = b:
dy
x = |y|1/2
dx
,
donde (a, b) ∈ R2 .
- Demuestre que si ab 6= 0, entonces el PVI tiene una única solución definida sobre un intervalo abierto que
contiene al punto a.
Ejercicio 95
Encontrar una función de Lyapunov para
(
x˙1 = x2
x˙2 = −xn1
Ver las condiciones, debe ser positiva definida, y su derivada negativa definida
1 1
V (x) = xn+1
1 + x22
n+1 2
Ejercicio 96
Encontrar puntos de bifuración: x0 = λx − x3
√
x = 0, x = ± λ
Si λ = 0 hay bifuración.
Page 106 of 117

Ejercicio 97
Encontrar puntos de estabilidad del Oscilador de Van der Pol, y encontrar una funcion Lyapunov débil
d2 x dx
+ β(x2 − 1) +x=0
dt2 dt
Vea que una función Lyapunov fuerte es: W (x, y) = (2 + β 2 )x2 − 2βxy + 2y 2
Primero debemos separar en un sistema
dx dy
= y, = −x + β(1 − x2 )y.
dt dt
Luego linealizarlo en el punto crı́tico (0,0) (Hartman–Grobman), o bien evaluar el Jacobiano
directamente.

0 1
A=
−1 β
Quedando los vectores propios

p
β ± (β 2 − 4)
λ=
2
Donde 0 < β < 2 es una espiral inestable, y un nodo inestable si β > 2
Luego, una función débil Lyapunov es trivialmente V (x, y) = (x2 + y 2 )/2 ≥ 0, ya que:
dV
= 2(xx0 + yy 0 ) = xy + y(−x + β(1 − x2 )y) = y 2 β(1 − x2 ) ≤ 0
dt
para β > 0 y x2 < 1.
Se deja propuesto la función fuerte (realizado en ayudantı́a).
Ejercicio 98
Resolver con Laplace, con C.I. distintas de cero y 00 + 2y 0 + y = 0; y 0 (0) = 2 y y(1) = 2.
L(y 00 + 2y 0 + y) = (s2 y(s) − sy(0) − y 0 (0)) + 2(sy(s) − y(0)) + y(s) = 0
Ahora fijamosy(0) = a
(s2 y(s) − sa − 2) + 2(sy(s) − a) + y(s) = 0

as + 2a + 2
y(s) =
(s + 1)2
y(t) = e−t (at + a + 2t)
Evaluando finalmente y(1) = 2
y(t) = e−t (et + t + e − 1)

Ejercicio 99
Encontrar puntos de bifurcación y ciclo lı́mite
dx dy
= ux − y − x(x2 + y 2 ), = x + uy − y(x2 + y 2 ).
dt dt
Veamos si tiene primero, linealizamos por Hartman–Grobman

u −1
A= . (19)
1 u
Y es cosa de analizar los valores propios.
λ=u±i
Por ende, hay bifurcación. Pasando a polares par encontrar el ciclio lı́mite estable.
dy
= ur − r3 , dθ/dt = 1.
dt
√
De la primera ecuación, el ciclo lı́mite es cuando r = u. Hacer notar que el Teorema de Bendixson no
aplica, por que apesar de no ser nulo:
df1 df2
+ = u − 4(x2 + y 2 )
dt dt
No es simplemente conexo para una bola abierta cerca de (0,0).
Ejercicio 100
Resolver la ecuación integro-diferencial
Z t
dy
+ 25 y(t − w)e−10w dw = 4; y(0) = 0
dt 0
L(y 0 + (25y(t) ∗ e−10t ) = 4)

1 4
sy(s) − 0 + 25y(s) =
s + 10 s

25 4
y(s) s + =
s + 10 s
4(s + 10) 8 4 8
y(s) = = − −
s(s + 5)2 5s (s + 5)2 5(s + 5)
L−1 (y(s))
Page 108 of 117

Ejercicio 101
Encontrar el Laplaciano inverso de las siguientes expresiones
1 1 1 s2 + a2
s4 , , , , {s log }
(s − 1)2 (s2 + 1)2 (s2 + a2 )n s2 − a2

1
n! 1 3! 1 6 1
1)L {tn } = , t3 = = =
sn+1 6 6 s3+1 6 s4 s4
n! 1 n!
2)L (tn ) = , L eat = , L tn eat =

sn+1 (s − a) (s − a)n+1

1 1 t 1
3)L−1 = (sin ωt − ωt cos ωt) , − cos t + sin t
(s2 + ω 2 )2 2ω 3 2 2
1 dn−1 est 1 dn−1 est
4)L−1 {F (s)} = lı́m n−1
(s − ia)n 2 2 n
+ lı́m n−1
(s + ia)n 2
s→ia (n − 1)! ds (s + a ) s→−ia (n − 1)! ds (s + a2 )n

1 sin(ax) sin(ax)
oL−1 2 2 n
= ∗ ··· ∗
(s + a ) | a {z a }
n times
s2 + a2
Bonus)f (x) : F = s log
s2 − a2
x
s2 + a2
Z
F
= log 2
f (t) dt :
0 s s − a2
Z x
d F 2s 2s
−x f (t) dt : = 2 2
− 2
0 ds s s + a s − a2
Z x
2s 2s
−x f (t) dt : 2 2
− 2
0 s + a s − a2
Z x
−x f (t) dt = 2 cos ax − 2 cosh ax (∗)
0
Z x
− f (t) dt − xf (x) = −2a sin ax − 2a sinh 2x
0
x2 f (x) = 2(ax sin ax + ax sinh ax + cos ax − cosh ax)

2
f (x) = (ax sin ax + ax sinh ax + cos ax − cosh ax)
x2
Page 109 of 117

Ejercicio 102
Explique que pasa si resuelvo con Laplace (y resuelvalo claro)
y 00 + 4y = 12sin(2t); y(π) = −3, y 0 (π) = −3
24
L(y 00 ) = s2 y(s) − sy(π) − y 0 (π), L(4y) = 4y(s), L(12sin(2t)) =
s2 + 4
24 3s 3
y(s) = − 2 −
(s2 + 4)2 s + 4 s2 + 4
3s 3 3
L−1 ( ) = 3 cos (2t), L−1 ( 2 ) = sin (2t)
s2
+4 s +4 2
Aca la dificultad, hay una convolucion entremedio
1 t
Z
24
L−1 ( 2 ) = sin(2t − 2τ ) sin(2τ ) dτ
(s + 4)2 4 0
Z t Z t
sin(2t − 2τ ) sin 2τ dτ = (sin 2t cos 2τ − cos 2t sin 2τ ) sin 2τ dτ
0 0
Z t Z t
= sin 2t cos 2τ sin 2τ dτ − cos 2t sin2 2τ dτ
0 0
Z t Z t
1 1
= sin 2t sin 4τ dτ − cos 2t (1 − cos 4τ ) dτ
0 2 0 2
t
1 t 1 1
= − sin 2t cos 4τ |0 − cos 2t τ − sin 4τ
8 2 4 0
1 1 t 1
= − sin 2t cos 4t + sin 2t − cos 2t + cos 2t sin 4t
8 8 2 8
1 1 t
= sin(4t − 2t) + sin 2t − cos 2t
8 8 2
1 t
= sin 2t − cos 2t
4 2
1
= (sin 2t − 2t cos 2t).
4
Finalmente
Z t
24 1 3
L−1 = 24 · sin(2t − 2τ ) sin 2τ dτ = (sin 2t − 2t cos 2t).
(s2 + 4)2 4 0 2
Page 110 of 117

Ejercicio 103
Resolver con vectores propios, para x(0) = y(0) = 1.
dx dy
= 3x − 4y, = 4x − 7y.
dt dt
0
x 3 −4 x
0 = .
y 4 −7 y

3 −4 1 0
det −λ .
4 −7 0 1
λ2 + 4λ − 5 = 0.
(λ − 1)(λ + 5) = 0
Para el primer valor propio

3 −4 α α
=1 .
4 −7 β β
3α − 4β = α, α = 2β

2
v̂1 = .
1
Luego con λ = −5

1
v̂2 = .
2

x
= Aeλ1 t v̂1 + Beλ2 t v̂2 .
y

x t 2 −5t 1
= Ae + Be .
y 1 2
Page 111 of 117

Ejercicio 104
*Resolver el sistema
−3 −1

x0 = 2 4
−1 x
1 2
−3 −1

−r
2
−1
4 ) = (r + 1)2
1 2 − r

1
Solo tiene m.a. 2 y m.g. con el vector propio.
−2
Para calcular el otro vector, utilizamos una reduccion de orden”, con matrices. (A − λI)v2 = v1 ,
quedando finalmente
x2 (t) = (v1 t + v2 )eλt
Ejercicio 105
x0 = y + x(1 − x2 − y 2 )
(Freak one) Resolver
y 0 = −x + y(1 − x2 − y 2 )
Pasando a polares (si, a polares)
r0 = r(1 − r2 )
θ0 = −1,
Z Z
dr
= dt
r(1 − r2 )
Z
1 1 1
− − dr = t + c
r 2(r + 1) 2(r − 1)
1 1
ln r − ln(r + 1) − ln(r − 1) = t + c
2 2

r
ln √ =t+c
r2 − 1
r
√ = Aet
r2 − 1
r2
= Be2t
r2 − 1
1
1 − 2 = Be−2t
r
−1
r = (1 − Be−2t ) 2
Las soluciones son:

−1
r(t) = (1 + ce−2t ) 2 , θ(t) = −(t − α).

Ejercicio 106
0
Resuelva x (t) = Ax(t),si
1 1 0 0
 
0 1 0 0
A=
0

0 2 1
0 0 0 2
Por ser , la solución es directa
et tet 0 0
 
0 et 0 0 
eJt = 
0 0 e2t te2t 
0 0 0 e2t
C1 et + C2 tet
 
 C3 et 
x(t) = 
C4 e + C5 te2t 
2t

C6 e2t
Page 113 of 117

Ejercicio 107
0
Resuelva x (t) = Ax(t), con
−3 −2
 
5
A=  8 −5 −4
−4 3 3
Utilizando el ejercicio anterior, de buscar la base canonica de Jordan

 
−2 0 1 0 0  −2 1 0
 
0 
 1

A = P JP −1 =  1
 −4 0  0 1 1 − 2 0 0
3 3
− 2 1 0 0 1
 
2 −2 1
2
A − I 6= 0, (A − I)2 = 0, (A − I)v3 6= 0
v3 = (0, 0, 1) =⇒ (A − I)v3 6= 0
v2 = (A − I)v3 = (−2, −4, 2)

3
(A − I)v1 = 0 =⇒ v1 = 0, 1, −
2
 

0 −2 0
   −2 1 0
1 0 0  1 
φ(t) = P eJt P −1 =  1
 −4 0 et
0 1 t
− 0 0 
 2 3
  
3

− 0 0 1

2 1 −2 1
2 2
4t + 1 −3t −2t
 
φ(t) = et  8t 1 − 6t −4t 
−4t 3t 2t + 1
c1 et (4t + 1) − 3c2 et t − 2c3 et t

   
x(t)
Xh = y(t) = eAt =  8c1 et t − c2 et (6t − 1) − 4c3 et t 
z(t) −4c1 et t + 3c2 et t + c3 et (2t + 1)
Page 114 of 117

Ejercicio 108
Caracterizar el sistema de masa resorte para todos los casos. Libre, amortiguado y sus casos, y con vibraciones
forzadas
d2 x dx
m 2
+b + kx = F (t)
dt dt
Fome, m= masa, b coeficiente de fricción, k constante del resorte, F(t), fuerza externa.
-Si b=0, no hay amortiguación, y Si F(t)= 0, el sistema es libre. y por ende las soluciones son
imaginarias, de tal forma que
u = c1 cos(w0 t) + c2 sen(w0 t)
r
k
Donde w0 = es la frecuencia natural.
m
Incluso, puede dejar expresado de esta forma,c1 = Rcos(α), c2 = Rsen(α)
u = Rcos(α)cos(w0 t) + Rsen(α)sen(w0 t) = Rcos(w0 t − α)
Es claro que:
√ c2
R= c1 + c2 , = tanα
c1
Donde, R es la amplitud, α el delay (phase shift)
-Si hay amortiguación, pero F(t)= 0, hay 3 casos
b2 − 4km = 0, Raı́ces reales iguales, amortiguación crı́tica
b2 − 4km > 0 Raı́ces reales diferentes, sobre amortiguación
b2 − 4km = 0 Raı́ces imaginarias diferentes, sub amortiguación
-Si no hay amortiguación, pero F (t) = cos(wt) (puede ser sen(wt) ) hay 2 casos.
w diferente a w0 , existe oscilaciones no controladas.
w igual a w0 , hay efecto resonancia.
Link de interes: http://tutorial.math.lamar.edu/Classes/DE/Vibrations.aspx
Page 115 of 117

Ejercicio 109
Resolver el sistema no homogeneo
0 2 −1 1
~x = ~x + et
3 −2 −1
Las soluciones de la homogenea son

1 1
λ = ±1, e1 = , e2 =
1 3
Resolviendo la matriz fundamental
e−t
t
e
φ(t) = t
e 3e−t
Recordando la formula, debemos invertir la matriz fundamental

Z
XP = φ(t) φ−1 (t)F (t)
−t
−e−t 3e−t −e−t

1 3e 1
φ−1 (t) = =
3e−t et − et e−t −et et 2 −et et
−t
−e−t et 3e−t et + e−t et

−1 1 3e 4
φ (t)F (t) = = =
2 −et et −et −et et − et et −2e2t
Integrando
Z
−1 4t
φ (t)F (t) =
−e2t
Finalmente
e−t
t
4tet − et
Z
−1 e 4t
Xp = φ(t) φ (t)F (t = t =
e 3e−t −e2t 4tet − 3et
Ejercicio 110
0
x 1 − cos t x
=
y cos t 1 y
Sea u = cos t , y luego evaluando la integral de la matriz en la exponencial

t −u 1 −i i t − iu 0 1 i 1
exp = exp √ · ·√
u t 2 1 1 0 t + iu 2 −i 1

x(t) t cos(sin t) − sin(sin t) x0
=e .
y(t) sin(sin t) cos(sin t) y0

Ejercicio 111
Resolver la ODE, donde P es la población de organismos, t el tiempo r define la tasa de crecimiento, K es
la capacidad de persistencia, y Po la población inicial. Indicar que pasa en un tiempo infinito,

dP P
= rP 1 −
dt K
Solución:
La solución general es:
KP0 ert
P (t) =
K + P0 (ert − 1)
And when t goes to infinity :
lı́m P (t) = K.
t→∞
Ejercicio 112
¿Como resolveria este ejercicio?
(
00 t, 0≤t<3
y + 4y = , y(0) = 0, y 0 (0) = 0
1, 3 ≤ t < ∞
Solución:
Hay 2 formas, tomar Laplace por partes

( !
∞ Z 3 Z ∞
t, 0 ≤ t < 3
Z
−st −st
F (s) = e dt = te dt + e−st dt
0 1, 3 ≤ t < ∞ 0 3
O bien plantear esto como un Heaviside
f (t) = t − (t − 1)u(t − 3)
1 e−3s e−3s (3s + 1)

L(t − (t − 1)u(t − 3)) = 2
+ −
s s s2
L(y 00 (t)) = s2 y(s) − sy(0) − y 0 (0) = s2 y(s)
L(4y 0 (t)) = 4(sy(s) − y(0)) = 4sy(s)
Page 117 of 117

Apuntesode

Cargado por

Copyright:

Formatos disponibles

Apuntesode

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntesode

Cargado por

Copyright:

Formatos disponibles

Mat1640: EDO

Due on Segundo semestre 2014

Hernan Gonzalez A. , higonzal@uc.cl

Ejercicio 100 108

Ejercicio 101 109

Ejercicio 102 110

Ejercicio 103 111

Ejercicio 104 112

Ejercicio 105 112

Ejercicio 106 113

Ejercicio 107 114

Ejercicio 108 115

Ejercicio 109 116

Ejercicio 110 116

Ejercicio 111 117

Ejercicio 112 117

Calculando el factor integrante:

Multiplicando la ecuación diferencial por µ(x)

Aplicamos exponencial y acomodamos un poco:

T (t) − Tm = e−kt+C ⇒ T (t) = e−kt · eC + Tm

Como la temperatura del ambiente es 20 o C, entonces la función queda definida por:

T (0) = 100 ⇔ 100 = e0 · eC + 20 ⇔ 80 = eC

Por lo que la función ahora nos queda:

Ejercicio 2 continued on next page. . . Page 8 of 117

Resolviendo el problema, la temperatura a los 5 minutos queda definida al evaluar la función en 5:

Para que la empanada se encuentre a 40o C, se deben esperar:

ln (1) − ln (4) 4 ln (2)

Sea A la cantidad de sal, y V el volumen, ve el volumen de entrada, y vs volumen de salida

la ecuacion diferencial queda:

Resolver la siguiente ecuación diferencial

6xy dx + (4y + 9x2 ) dy = 0

Se elige M = 6xy y N = 4y + 9x2 .Notemos que no es Difrencial exacta,

6xy 3 dx + (4y 3 + 9y 2 x2 )dy = 0

con C(y), çonstante”que depende de y Buscamos C(y):

Ejercicio 4 continued on next page. . . Page 11 of 117

Para transformar la ecuación en una diferencial exacta

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

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

Para que sea exacta debe cumplirse que

∂(µ(x)M (x, y)) ∂(µ(x)N (x, y))

∂µ(x)M (x, y) ∂M (x, y)

∂µ(x)N (x, y ∂N (x, y) ∂µ(x)

∂M (x, y) ∂N (x, y) ∂µ(x)

dµ(x) ∂M (x, y) ∂N (x, y)

Finalmente, aplicando exponencial, el factor integrante es de la forma:

Ejercicio 4 continued on next page. . . Page 12 of 117

Resolver la siguiente ecuación diferencial

Reemplazar y 2 = xz, y la EDO queda

que es de variables separables.

Resolver la siguiente ecuación diferencial

Que corresponde a un bernoulli, en este tipo, divides todo por y 3 , quedando

Realizando la sustitución u = −y −2 , u0 = 2y −3 y 0 y normalizando:

Resolver las siguientes ecuaciones

Aplicamos función exponencial a ambos lados de la igualdad para ası́ obtener:¡br /¿

eC2 es constante, ası́ que podemos hacer C = eC2 y luego:

Ejercicio 7 continued on next page. . . Page 15 of 117

donde A es una integral bien conocida, pero se saca haciendo la sustitucion

Se puede decir que el área A debajo de la curva es n veces menor a la otra: