Ejercicios resueltos

Ecuaciones Diferenciales Ordinarias

Hernán González ([email protected])

Due on: 07/14. Last Update: 03/21

Ecuaciones Diferenciales Ordinarias [email protected]

Índice
Ecuaciones diferenciales de primer orden 5
Variables separables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problem 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problem 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Problem 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Modelos y aplicaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Problem 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Problem 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Problem 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Problem 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Problem 8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Factor integrante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
Problem 9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Problem 10 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Ejercicio 11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Problem 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Problem 13 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
Problem 14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
Problem 15 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Ricatti . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Problem 16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Problem 17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Ecuaciones exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Problem 18 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
Problem 19 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
Problem 20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Problem 21 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
Problem 22 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
Problem 23 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Problem 24 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
Cambios de variable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Problem 25 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Problem 26 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Homogéneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Problem 27 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Problem 28 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Problem 29 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Problem 30 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Problem 31 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Existencia y unicidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Ejercicio 32 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Problem 33 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Problem 34 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Mixtos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Problem 35 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Problem 36 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

Page 1 of 116
Ecuaciones Diferenciales Ordinarias [email protected]

Problem 37 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
Problem 38 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
Problem 39 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
Problem 40 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
Problem 41 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Problem 42 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Problem 43 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Problem 44 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Problem 45 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Problem 46 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Problem 47 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51

Ecuaciones de segundo orden 52

Reducción de orden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
Problem 48 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
Polinomio caracterı́stico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Problem 49 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Cauchy-Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Problem 50 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Variación de parámetros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Problem 51 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Problem 52 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
Problem 53 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
Coeficientes indeterminados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Problem 54 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Aniquiladores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
Problem 55 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
Modelos y aplicaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
Problem 56 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
Problem 57 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
Cambio de variable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
Problem 58 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
Problem 59 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
Problem 60 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
Existencia y unicidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
Ejercicio 61 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
Mixtos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
Problem 62 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
Problem 63 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
Problem 64 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
Ejercicio 65 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
Ejercicio 66 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
Ejercicio 67 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
Ejercicio 68 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
Ejercicio 69 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
Ejercicio 70 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
Ejercicio 71 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
Ejercicio 72 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72

Page 2 of 116
Ecuaciones Diferenciales Ordinarias [email protected]

Ejercicio 73 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Ejercicio 74 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Ejercicio 75 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
Problem 76 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
Problem 77 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
Ejercicio 78 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
Problem 79 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
Problem 80 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
Problem 81 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77

Sistema de ecuaciones 78
Cambio de variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
Problem 82 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
Problem 83 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
Diagonalización . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Problem 84 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Problem 85 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
Problem 86 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
Problem 87 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
Problem 88 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
Problem 89 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
Matriz fundamental . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
Problem 90 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
Principio de Duhamel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
Problem 91 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
Problem 92 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
Coeficientes indeterminados . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
Problem 93 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
Modelos y aplicaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
Problem 94 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
Existencia y unicidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
Mixtos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
Problem 95 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
Problem 96 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
Problem 97 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

Series 98
Series de potencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
Problem 98 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
Series de Frobenius . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
Problem 99 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
Problem 100 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
Problem 101 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100

Page 3 of 116
Ecuaciones Diferenciales Ordinarias [email protected]

Laplace 101
Inversa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
Problem 102 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
Problem 103 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
Delta de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
Problem 104 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
Heaviside . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
Problem 105 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
Convolución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
Problem 106 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
Ecuaciones integro-diferenciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
Problem 107 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
Mixtos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
Problem 108 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
Problem 109 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
Problem 110 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
Problem 111 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
Problem 112 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
Problem 113 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111

Sistemas dinámicos 112

Hartman-Grobman . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
Problem 114 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
Sistemas Hamiltonianos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
Mixtos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
Ejercicio 115 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
Ejercicio 116 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
Problem 117 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
Problem 118 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Problem 119 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Problem 120 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Problem 121 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Problem 122 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Problem 123 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116

Page 4 of 116
Ecuaciones Diferenciales Ordinarias [email protected]

Ecuaciones diferenciales de primer orden

Variables separables
Problem 1
Resolver
dy
ln (y x ) = 3x2 y
dx

Operando, separando variables e integrando ambos miembros de la ecuació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
p
3x2 +c
y = e

Problem 2
ex
Resolver en terminos de una integral definida, y 0 ex + yex = con y(2) = 1
x+1

ex ex
y 0 ex + yex = ) (yex )0 =
x+1 x+1
d(yex ) ex ex
= ) d(yex ) = dx
dx x+1 x+1
Si x = 2, yex = e2 Integrando:
Z yex Z x
ex
d(yex ) = dx
e2 2 x+1
Z x x Z x x
e e
yex e2 = dx ) yex = dx + e2
2 x + 1 2 x + 1
Z x x
e
y(x) = e x dx + e2 x
2 x+1

Page 5 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 3

Problem 3

Resolver las siguientes ecuaciones

dy y dy 1
=1 , =1+ 2
dt 3 dt y
Solución:

dy y
La ecuació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 función exponencial a ambos lados de la igualdad para ası́ obtener:¡br /¿

t t
|3 y| = e 3 +C2 =e 3 · e C2

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, llamémosleK, con lo cual
(despejando y) la solución general de la EDO es:
t
y=3 Ke 3

Problem 3 continued on next page. . . Page 6 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 3 (continued)

dy 1 dy y2 + 1
=1+ 2 , =
dt y dt y2

Z
y2
dy = dt/
y2 + 1
Z Z
y2
dy = dt
y2 + 1

Z
y2
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 7 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 3

Modelos y aplicaciones
Problem 4
Un pueblo llamado Melosova, que conlinda con Ponmelov tenı́a una población de 5000 habitantes el año 2000
y una población de 7000 almas el año 2007. Suponiendo que su población siga creciendo exponencialmente
con su ı́ndice constante, ¿Cuál será la población de Melosova el año 2050?

Solució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 añ o 2050 (es decir 50 años después) habran y = 5000 · e0,04806·50 ⇡ 55,281 habitantes.

Problem 5
Una curva de ecuació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.

Solución:

Se puede decir que el á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 8 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 6

Problem 6
En la conservación de alimentos, el azúcar de caña sufre un proceso de inversión y se transforma en glucosa
y fructosa. En soluciones diluidas el ritmo de inversión es proporcional a la concentración del azucar. Si la
concentración al inicio es de 1/50 gr/L y 1/200 gr/L tras tres horas, halle la concentración de azúcar después
de 6 horas.

Solución:

dy
El modelo es = ky
dt
Z Z
dy
= kdt ! ln |y| = kt + c ) y = ekt+c
y
Que con la info, se transforma en 2 ecuaciones con 2 incógnitas

1
= ec
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

Problem 7

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 mó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 será necesario esperar?

dT
= k (T (t) Tm )
dt
Para hallar la función T(t), se resuelve la ecuación diferencial. Primero separamos las variables:
Z Z
dT dT
= kdt ) = kdt ) ln (T (t) Tm ) = kt + C
T (t) Tm T (t) Tm

Problem 7 continued on next page. . . Page 9 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 7 (continued)

Aplicamos exponencial y acomodamos un poco:

kt+C kt
T (t) Tm = e ) T (t) = e · e C + Tm

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

kt
T (t) = e · 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 función ahora nos queda:

kt
T (t) = 80e + 20

Otra condición planteada en el problema es que a los dos minutos la temperatura de la empanada es de
80o C, es decir:

✓ ◆
2k 3 2k 1 3
T (2) = 80 , 80 = 80e + 20 , =e ,k= ln
4 2 4

Luego:

t
✓ ◆ 2t
t
ln( 34 ) ln( 34 ) 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 función en 5:

✓ ◆ 52 p p !
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 10 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 8

Problem 8

Un tanque contiene, en un principio, 40[L] de salmuera (sal muy concentrada en agua). Entra al tanque,
1
a una razón de 8[L/min] , salmuera conteniendo + ↵ de sal por litro. La mezcla, bien agitada, sale del
8
23
tanque a razó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 concentración de sal en el tanque sea de [Kg] de sal por litro.
25
Solució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
= v e ⇤ re v s ⇤ rs
dt

dA 4A 1
+ = 8 ⇤ ( + ↵)
dt 40 + 4t 8
resolviendo

1 1
A= (1 + 8↵)(t + 10) + C(t + 10)
2
tenemos dos condiciones

1 1 23
A(10) = (1 + 8↵)(20) + C(20) =
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 11 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 8

Factor integrante
Como se verá, 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 algún factor, que depende de x , µ(x) siempre mayor que cero.

dy
µ(x) + µ(x)p(x)y = µ(x)q(x)(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 ecuación (2) cumpla con la condició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, usando variables separables)
R
p(x)dx
µ(x) = e
Donde µ(x), es el famoso y popular factor integrante.

Page 12 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 9

Problem 9
Resuelva las siguientes 3 ecuaciones.
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 2 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 integración por partes (Con el gentil auspicio de
Cálculo I, destruyendo cerebros desde que Newton se aburrio un dı́a. Lo veran hasta en la sopa a futuro)

Z Z Z
x dex y x x
e y= dx = e (x + 1)dx = e (x + 1) ex dx = ex (x + 1) ex + C = xex + C.
dx
x
Despejando y de la ecuación, queda finalmente y(x) = e (xex + C) = x + Ce x

Page 13 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 10

Problem 10
Resolver
ds
+ s · tan(t) = 2t + t2 tan t
dt

Z
ln|cos t|
µ(x) = exp tg tdt = e = sec t

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

Page 14 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 10

Bernoulli
Ejercicio 11
Resuelva la ecuación de Bernoulli:
y0 = y + y3

Solución:

Se acomoda la ecuación quedando:

y0 1
=1
y3 y2

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 solución de la EDO lineal es:

R
✓ Z R
◆
2dx dx
v(x) = e c+ 1e

2x
v(x) = ce 1

Pero
1
v=
y2
por lo tanto, la solución queda en:

1
y(x) = ± p
ce 2x 1

Page 15 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 12

Problem 12
Resuelva la ecuación de Bernoulli:
x2 y 0 + 2xy = 5y 3

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
-5
, x u’ - 4x u = - 10x - 6
-4

d x - 4u
, = - 10x - 6
dx Z
, x u = - 10 x - 6 dx
-4

x-5
, x - 4 u = - 10 +c
-5
2
, u= + cx4
x
2 + cx5
, y-2 =
x
x
, y2 =
2 + cx5

r
x
) |y| =
2 + cx5

Problem 13

Resolver la siguiente ecuación diferencial

y0 y = e2x y 3
Solución:

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

3 0 2
y y y = e2x

Realizando la sustitución u = y 2
, u0 = 2y 3 0
y y normalizando:

u0 + 2u = 2e2x

Esto es:

2x 1
u(x) = c1 e + e2x
2
Volviendo a la variable original:

2 2x 1
y = c1 e + e2x
2

Page 16 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 14

Problem 14
Solve
dy 4 p
= y+x y
dx x

Solución:

dy 4 p
= 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 = y2 ! = 2v
dx x 2 dx dx

✓ ◆
dv 4 1
! 2v + v2 = xv / ·
dx x 2v
✓ ◆ R
dv 2 x 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 1
·x 2+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 17 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 15

Problem 15
e3x y 3 2y
Sea la EDO y 0 (x) = (x y). Resuelva la EDO con condición inicial y(0) = 1.
y x
Solución:

En este caso, la EDO ya viene ”normalizada”, por lo tanto, sólo debemos darnos cuenta que es de
Bernoulli. Veamos que para R2 {y = x} la EDO tiene solució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,
u0 = 2y 3 · y 0 (x), este resultado nos dice que debemos multiplicar la EDO de Bernoulli por 2y 3 , se
sigue que:

y 0 (x) 2y = e3x y 3 () 2y 3 · y 0 (x) + 4y 2

= 2e3x () u0 (x) + 4u = 2e3x .

Notemos que esta última

R
EDO es lineal, con p(x) R
= 4 y q(x) = 2e3x , luego, considerando el Factor de
p(x)dx 4dx
Integración µ(x) = e , se tiene que µ(x) = e = 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 solució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
región positiva, ya sea el 1er. o 2do cuadrante, de ahı́ que y > 0, por tanto, la solución de la EDO o del
PVI es y(x) = 1.

Page 18 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 15

Ricatti
Problem 16
1
Resolver conciendo la solución particular es y1 =
x
2
y0 = y2 , y(1) = 15
x2

Solución:

1
Efectuamos el cambio de variable y = v + , entonces:
x
dv 2v dv
= v2 + () 2x 1
v = v2
dx x dx
Se tiene una ecuació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
2 du
u 2u 1 x 1 = u 2
dx
du
+ 2x 1 u = 1
dx
Ahora se tiene una ecuació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 ecuación por este factor y volviendo a nuestras variables originales se tendrá:

x2 u0 + 2xu = x2
(x2 u)0 = x2
x3
x2 u = +c
3
3cx 2 x
u =
3
3
v =
3cx 2 x
3 1
y = +
3cx 2 x x

Problem 16 continued on next page. . . Page 19 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 16 (continued)

Finalmente imponemos la condició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= + =) =
u x dx u2 dx x2
Luego al reemplazar nos queda

1 du 1 1 2 1 2
= + +
u2 dx x2 u2 ux x2 x2
1 du 1 2
= +
u2 dx u2 ux
du 1
=) = 1 2ux
dx
du
) + 2x 1
u = 1
dx
Esta es la conocida ecuación de Ricatti

dy
+ p(x)y + q(x)y 2 = f (x)
dx
Esta ecuación se resuelve si previamente se conoce una solución particular, sea y1 (x) Conocida dicha
solución, encontrada por inspección, se hace el cambio:

y(x) = z(x) + y1 (x)

Page 20 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 17

Problem 17
Resolver, sabiendo que la solución particular es de la forma aekx

y0 = e x 2
y + 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)
x x - 2x
, e u’ - e u=e
d x
, e u = e - 2x
dx Z
, e xu = e - 2x dx

x 1 - 2x
, e 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
2x
ce - 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 21 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 17

Ecuaciones exactas
Para una ecuación exacta, debe cumplirse que

@M @N
M (x, y)dx + N (x, y)dy = 0, = .
@y @x
@f @f
Y significa que la ecuación proviene f, tal que df = M (x, y)dx + N (x, y)dy, que además df = @x dx + @y dy.
(esto se conoce como campo conservativo, donde f es el potencial, será visto en cálculo III.)

Problem 18
Resolver EDO (todos los ejercicios del curso son asi)

3x(xy 2)dx + (x3 + 2y)dy = 0

Solución:

M (x, y) = 3x(xy 2), N (x, y) = x3 + 2y.

Notamos que se cumplı́o la condición de ecuació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 22 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 19

Problem 19
Resuelva la ecuación diferencial sujeta a la condición dada:

(x + y)2 dx + (2xy + x2 1)dy = 0; y(1) = 1

Solució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 función F tale que :

@F @F
=M^ =N
@x @y

con F (x, y) = 0 solución de la EDO

Por lo tanto, si se hace la integración
Z a N , respecto de Y , se obtiene que:
y
N =F

Z y
2xy + x2 1 = xy 2 + yx2 y + C(x)

Con C(x), una ç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 solución (implicita)de la EDO es :

x3
xy 2 + yx2 y+ =k
3
Por último, basta reemplazar x e y dados por la condición inicial y(1) = 1 para asi calcular k:

1 4
1+1 1+ =k!k=
3 3
Por lo tanto, al solución implicita al PVI es
x3 4
xy 2 + yx2 y + =
3 3

Page 23 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 20

Problem 20

[n cos(nx + my) msen(mx + ny)]dx + [m cos(nx + my) nsen(mx + ny)]dy = 0

Solució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 solución de la ecuació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 solución es:

sin(nx + my) + cos(mx + ny) = C

Problem 21

y + senx cos2 xy x
dx + ( 2 + seny)dy = 0
cos2 xy cos xy

y x
M (x, y) = + sin(x) N (x, y) = + sin(y)
cos2 (xy) cos2 (xy)

@M 1 2(xy) tan(xy) @N
= + =
@y cos2 (xy) cos2 (xy) @x
Luego existe F (x, y) = C solución de la EDO.
Z ✓ ◆
x
F (x, y) = + sin(y) dy = tan(xy) cos(y) + h(x)
cos2 (xy)

@F y y
= 2
+ h0 (x) = M (x, y) = 2
+ sin(x)
@x cos (xy) cos (xy)
) h0 (x) = sin(x) ) h(x) = cos(x) + c1
Luego la solución es: F (x, y) = tan(xy) cos(y) + cos(x) = C

Problem 21 continued on next page. . . Page 24 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 21 (continued)

Transformar la ecuación en una diferencial exacta M (x, y)dx + N (x, y)dy = 0.

(Queda como ejercicio encontrar un factor integrante que dependa de y)
Utilizamos un factor integrante tal que depende de x, una función µ(x) > 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))

=
@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 condició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:
Z @M (x,y) @N (x,y)
!
@y @x
µ(x) = exp dx
N (x, y)

Page 25 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 22

Problem 22

Resolver la siguiente ecuación diferencial

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

Solució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 Funció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), ç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

Page 26 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 23

Problem 23
Resolver
(3xy + y 2 )dx + (x2 + xy)dy = 0

Solució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 + 12 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 27 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 24

Problem 24
Para x > 1 encuentre la solución general implı́cita de la ecuación

x2 + y 2 + 1 dx (xy + y) dy = 0
n
Sabiendo que tiene factor integrante de la forma µ(x, y) = (x + 1)

Solución:

M (x, y) = x2 + y 2 + 1 N (x, y) = (xy + y)

utilizando el factor integrante

µM dx µN dy = 0
2
n
(x + 1) (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 n
(x, y) = 2y(x + 1)
@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 n
(n 1)(x + 1) y = 2y(x + 1)
n 1=2
n=3
3
µ(x, y) = (x + 1)

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

Problem 24 continued on next page. . . Page 28 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 24 (continued)

2
3
(x + 1) (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) 3
= (x + 1) (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 ⇣ ⌘
) g(x, y) = + ln(x + 1) + 2x+1
(x+1) 2
2(x + 1)2
Las soluciones estan definidas implicitamente por:

g(x, y) = C C 2 R

Finalmente se obtiene :
1
y(x) = (2(x + 1)2 ln(x + 1) + 2(2x + 1) C) 2

Page 29 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 24

Cambios de variable
Problem 25
Resolver
y(ln(y) ln(x) + 1)
y0 =
x

Reescribiendo la ecuació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

Problem 26

(x + ey )y 0 = xe y
1

Solución:

x
xe y 1 x ey 1 ey 1 1
y0 = y
= · = x ·
x+e x + ey ey ey + 1 ey

x x u u xu0
u= y
) y = ln ) y 0 = ·
e u x u2
]
✓ ◆
u u xu0 u 1 u u xu0 u 1 u 1
· = · ) = ) xu0 = u u2
x u2 u+1 x u2 u+1 u+1
Z Z
1
) du = xdx
u (u + 1)

Z ✓ ◆ ✓ ◆
1 1 u+1 x2 1 x2 1 x
du = ln = + C ) 1 + = Ke 2 ) u = x2
=
u u+1 u 2 u Ke 2 1 ey
⇣ x2
⌘ h ⇣ x2
⌘i
) ey = x Ke 2 1 ) y = ln x Ke 2 1

Problem 26 continued on next page. . . Page 30 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 26 (continued)

Homogéneas
Problem 27

p
xy 0 = y x2 + y 2

Solución:

Notemos
r que reordenando
y ⇣ y ⌘2 y
y0 = 1 , luego sea u =
x x x
p
xu0 + u = u 1 u2
du dx
p =
1 u2 x
1
sen (u) = ln(x) + C
Volviendo a la variable original
⇣y⌘
sen 1 = ln(x) + C
x

Page 31 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 28

Problem 28
Muestre que la ecuación diferencial es homogénea, y resué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
y 2
dy 1 x
=
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 z2
x +z =
dx 2z
dz 1 z2
x = z
dx 2z
dz 1 z 2 2z 2
x =
dx 2z
dz 1 3z 2
x =
dx 2z
2z dx
2
dz =
1 + 3z x
Z Z
2z dx
2
dz = + ln (c)
1 + 3z x
1
ln 1 + 3z 2 = ln |x| + ln (c) / · 3
3
c
ln 1 + 3z 2 = ln 3
x
c
) 1 + 3z 2 = 3
x
Volvamos a nuestras variables:
✓ 2◆
y c
)1+3 2
= 3
x x
) x3 + 3y 2 x = c
p
3x (x3 3c)
) y (x) = ±
3x

Page 32 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 29

Problem 29
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)

y
dy x 3+5 x
= y
dx x 4+6 x

dy 3 + 5 xy
=
dx 4 + 6 xy
Hacemos el cambio de variable, y resolvemos por ser función homogé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

Problem 29 continued on next page. . . Page 33 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 29 (continued)

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)

2 c
ln |2z + 1| + ln (z + 1) = ln
x3

2 c
ln (2z + 1) (z + 1) = ln
x3

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

c
) 2z 3 + 5z 2 + 4z + 1 =
x3
Volvamos a nuestras variables originales:

y3 3 y2 3 y
)2 x + 5 x + 4 x3 + x3 = c
x3 x2 x

) 2y 3 + 5y 2 x + 4yx2 + x3 = c

2
) (x + y) (x + 2y) = c

Se deja como propuesto realizarlo por diferenciales exactas.

Page 34 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 30

Problem 30
Resolver
(8y + 10x)dx + (5y + 7x)dy = 0

Solució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 5
ln |v + 1| + ln |v + 2| = ln x + ln |c|
2 3 5
ln (v + 1) (v + 2) = ln cx
2 3 5
(v + 1) (v + 2) = cx
⇣y ⌘2 ⇣ y ⌘3
5
+1 +2 = cx
x x
(y + x)2 (y + 2x)3 = c

Page 35 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 31

Problem 31

xy 0 2y = x5

Solució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
3
4z dx
4
dz = )
1 3z 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 36 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 31

Existencia y unicidad
Ejercicio 32
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 2 R para los cuales el PVI no tiene solucion. Encuentre todos los valores
de a, b 2 R para los cuales el PVI tiene más de una solución y determı́nelas.

Sabemos que f (x, y) es continua para y 1 y x 2 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 2 R. Luego en esta region, por Teorema de Existencia y Unicidad,
existe solucion única.
Para determinar la solución, vemos que es de variables separables:
dy dx p
p = 2 1 + y = arctan(x) + C
1+y 1 + x2
✓ ◆2
arctan(x)
y(x) = +C 1
2
Luego, la solución al PVI es:

✓ ◆2
arctan(x) p
y(x) = +2 1+b arctan(a) 1
2

Notemos que para b = 1 y a 2 R, el teorema de existencia y unicidad no aplica. Verificamos que la

solución singular y(x) = 1 es solución al PVI. Entonces poseemos más de 1 solución:

✓ ◆2
arctan(x)
y(x) = 1 & y(x) = arctan(a) 1
2

Page 37 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 33

Problem 33
Sea y 0 = 1 + sin(xy) y D = {(x, y) 2 R2 : 0 < x < 1, 1 < y < 1} Demuestre que para cualquier punto
(a, y0 ) con 0 < a < 1, el PVI asociado a algun intervalo centrado en a tiene solucion unica

Solució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, también en D. Por lo tanto existe AL MENOS UNA
solucion definida en un intervalo abierto J, que contiene a x=a
Además:

@f (x, y)
= x cos(xy)
@y
La cual tambien es continua en todo R2, en particular D, entonces la solución también es UNICA en un
intervalo abierto J0 que contiene a x=a.El teorema de existencia y unicidad nos asegura esto.

Problem 34
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 función, contradiciendo la unicidad, y no la existencia

0,6
4xy
fy =
3

Page 38 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 34

Mixtos
Problem 35
Resolver la EDO (todo el curso es asi, no se asombren) utilizando la sustitución y = z n para algún n
apropiado.

2x4 yy 0 + y 4 = 4x6
Solución:
Sea y = z n , entonces y 0 = nz n 1 0
z , sustituyendo
1 0
2nx4 z 2n z + 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= 32
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
0 1
3 2 3 2
u u dx
B 5 C
@ 3 + 35 A du =
u 1 u +4 x

1 1
ln(u3 + 4) ln(u3 1) = ln(x) + C
5 5
✓ ◆1
u3 + 4 5
= Cx
u3 1
retornando a
0 1 15
y2
+4
B x3 C
@ 2 A = Cx
y
1
x3
✓ ◆ 15
y 2 + 4x3
= Cx
y 2 x3

Problem 35 continued on next page. . . Page 39 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 35 (continued)

Problem 36
Muestre que la ecuació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.

Solución:

(2x + xy 2 )dx + (4y + x2 y)dy = 0 (1)

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 ecuación:
2x 1 x
y0 + + y=0
4 + x2 y 4 + x2
x 2x
y0 + y= y 1
4 + x2 4 + x2
Se muestra que es Bernoulli. Nuevamente tomando la primera ecuación:
x y
dx + dy = 0
4 + x2 2 + y2
Se muestra que es de variables separables. Usando el primer mé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 solución general es:
1
x2 + x2 y 2 + 2y 2 + c = 0
2

Page 40 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 37

Problem 37
Hallar la solución general a la ecuación diferencial:
dy
cos x · sen x + y sen x = 1
dx

Solució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)

Sustitución (invocando los poderes de MAT1610):

⇣x⌘
t = tan
2

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

Problem 37 continued on next page. . . Page 41 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 37 (continued)

Reemplazando:

Z 2
1+t2
I= ⇣ ⌘2 dt
2t 2t
1+t2 1+t2

Z
dt
I=2 4t
2t 1+t2

Z
1 + t2
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 x +C
2 tan 2 1

x 2
ln tan tan( x
+C
2) 1
2
y=
(sec (x) + tan (x))

Page 42 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 38

Problem 38

Resolver p
du 1 x2 4 x arctan(x) p
x + u= p u
dx 1 + x2 1 + x2
Solución:

p
Dividiendo por x u, suponiendo u 6= 0, x 6= 0, se obtiene:

1 du 1 x2 1 4 arctan(x)
p + u2 = p
u dx x(1 + x2 ) x(1 + x2 )

p 1 du
Sea v = u, donde v 0 = p
2 u
, ası́:
dx
1 x2 2 arctan(x)
v0 + v=p
2x(1 + x2 ) x(1 + x2 )

Calculando el factor integrante:

✓Z ◆ ✓Z ⇢ ◆
1 x2 1 x
µ(x) = exp dx = exp dx
2x(1 + x2 ) 2x 1 + x2
✓ ◆ (✓ ◆ 12 )!
1 1 2 x
= exp ln(x) ln(1 + x ) = exp ln
2 2 1 + x2
r
x
=
1 + x2

Multiplicando la ecuació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
p x
u = arctan2 (x) + C
1 + x2

Se obtiene finalmente:
1 + x2 2
u(x) = arctan2 (x) + C x 6= 0, C 2 R
x

Page 43 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 39

Problem 39
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 44 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 40

Problem 40

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.

Solució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
y t exp t2 dt |{z} = 2 exp (u) du = 2 exp (u) = 2 exp t2
u= t2

ergo :
⇥ ⇤
, y = arc cos h exp t2 4 exp t2 + C
⇥ ⇤
, y = arc cos h C exp t2 4

Problem 40 continued on next page. . . Page 45 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 40 (continued)

Por ú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
2 2
y = ln Cet 4± Cet2 4 1

Problem 41

1 e2x
y0 + e 2x 2
y (1 + 4x + 2x2 )y = (1 + x + 2x2 + x3 )
x x
se pide encontrar la solución particular de la forma y1 (x) = e2x (Ax + B)

Solución:

Ricatti de nuevo (será 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

Problem 41 continued on next page. . . Page 46 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 41 (continued)

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
2x 2
e e2x (Ax + B)

e2x 1
e2 (A + 2Ax + 2B) = 1 + x + 2x2 + x3 + 1 + 4x + 2x2 e2x (Ax + B)
x x
2x
e 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
A2 x2 + 2ABx + B 2

0 1
1 2 B 2
1 2x x + A + + 4Ax + 4B + 2Ax + 2xB
e2x (A + 2Ax + 2B) = e2x @ x x A
2 2 2
A x 2ABx B

0 1
1
x2 1 + 2A A2 + x ( 2 + 4A + 2B 2AB) + ( 1 + B)
e2x (A + 2Ax + 2B) = @ x A
+ 1 + A + 4B B2

Problem 41 continued on next page. . . Page 47 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 41 (continued)

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 48 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 41

Propuestos
Problem 42
Resolver la siguiente ecuación diferencial

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

La idea es mostrar que pueden hacerla por diferenciales exactas, o por separables.
La ecuación es equivalente a

y3
y0 = .
2x2 2xy 2

Reemplazar y 2 = xz, y la EDO queda

p 0 ✓ p ◆
xz 1 z z p
p = p + 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 homogéneas de grado n = m., y su sustitución es y = vx.
2) y 0 = G(ax + by), donde la sustitució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.

Problem 43
Resolver el problema de Cauchy
y 0 + 2y = 1

5
y(0) =
2

Se puede resolver con método de variables separables, pero veamos 2 metodos diferentes:
La E.D.O. homogénea asociada es y 0 + 2y = 0, que tiene por solución general yh (t) = ke 2t . Se tiene ası́,
que la solución general de la E.D.O. no homogénea viene dada por y(t) = yh (t) + yp (t), donde yp (t) es
una solución particular de la ecuación. No es difı́cil notar (el ”método”que se uso para hallar soluciones
particulares se llama çoeficientes indeterminados”) que yp (t) = 12 es solución particular de la ecuación
y 0 + 2y = 1. Se tiene ası́, que y(t) = ke 2t + 12 es solución general de la E.D.O. no homogénea. Imponiendo
la condición inicial se obtiene el valor de k para que dicha función satisfaga el problema de Cauchy.
Pueden ocupar el otro método
R
que funciona para las EDOs lineales de primer orden también, que
consiste en multiplicar por e p(x)dx donde p(x) es el término que acompaña a y en la ecuación (siempre y
cuando a y 0 lo acompaña un 1).
La EDO queda
(ye2x )0 = e2x y ahı́ ya es cosa de integrar y despejar.

Page 49 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 44

Problem 44
Determinar la solucion de la ecuacion diferencial

(x + 2y)dx + (2x 3y)dy = 0

Solució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 homogénea, todos vieron que era exacta)

M (x, y) = x + 2y N (x, y) = 2x 3y

@M @N
) =
@y @x

Problem 45
Resolver la ODE, donde P es la población de organismos, t el tiempo r define la tasa de crecimiento, K es
la capacidad de persistencia, y Po la población inicial. Indicar que pasa en un tiempo infinito,
✓ ◆
dP P
= rP 1
dt K

Solución:

La solución general es:

KP0 ert
P (t) =
K + P0 (ert 1)

And when t goes to infinity :

lı́m P (t) = K.
t!1

Page 50 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 46

Problem 46

Pruebe que, si y1 es solució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)
Solució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 solución de y 0 + p(x)y = ↵q1 (x) + q2 (x) es:

y = ↵y1 + y2

Problem 47
Consideremos el siguiente problema con valor inicial, y(a) = b:

dy
x = |y|1/2
dx
,
donde (a, b) 2 R2 .
- Demuestre que si ab 6= 0, entonces el PVI tiene una única solución definida sobre un intervalo abierto que
contiene al punto a.

Page 51 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 47

Ecuaciones de segundo orden

Reducción de orden
Problem 48
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
y1 y2
W (x) = .
y10 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 último, resuelva esto último para deducir
⇢ Z
B(x)
W (x) = K exp dx .
A(x)

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
ú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
Z
Z B(x)
dW B(x) B(x) dx
= dx ) ln(W ) = dx + c ) W (x) = Ke A(x)
W A(x) A(x)

Problem 48 continued on next page. . . Page 52 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 48 (continued)

Polinomio caracterı́stico
Problem 49

Determine los valores de k 2 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 .

Solución:

Es fácil ver que

y = ekx ) y 0 = kekx (2)

00 2 kx
)y =k e (3)
000 3 kx
)y =k e (4)

Luego, la ecuación a resolver la podemos reescribir como

k 3 ekx 6k 2 ekx + 11kekx 6ekx = 0 (5)

Dado que ekx > 0, 8x 2 R, podemos simplicar la ecuación (4) y reescribirla como

k3 6k 2 + 11k 6=0 (6)

Pero, al factorizarla obtenemos

(k 3)(k 2)(k 1) = 0 (7)

De lo cual se concluye que k1 = 1, k2 = 2 y k3 = 3.

.
Ahora, para la segunda parte, dado que

y = c 1 e k 1 x + c 2 e k2 x + c 3 e k3 x (8)
0 k1 x k2 x k3 x
) y = c 1 k1 e + c2 k2 e + c 3 k3 e (9)
00
)y = c1 k12 ek1 x + c2 k22 ek2 x + c3 k32 ek3 x (10)
000
)y = c1 k13 ek1 x + c2 k23 ek2 x + c3 k33 ek3 x (11)

Luego, la ecuació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 paréntesis es cero, se concluye que no importa el valor de los coeficientes c1 , c2 , c3 .

Page 53 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 49

Cauchy-Euler
Problem 50
Resuelva
x2 y 00 + xy 0 y=0

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
2
m + (1 1)m 1=0
2
m 1=0
(m + 1) (m 1) = 0
m1 = 1m2 = 1
1
y1 = x, y2 = x
Asi, la base de la solucion quedaria:
y = c1 y1 + c2 y2 ; c1 , c2 ,constantes. ! y = c1 x1 + c2 x 1

Variación de parámetros
Problem 51

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:

e 2x xe 2x
2e 2x e 2x
(2x 1)
2x 2x 2x 2x 4x 4x 4x 4x
= e e (2x 1) + 2xe e = e (2x 1) + 2xe = ( 2x + 1 + 2x)e =e .
Como el Wroskiano es distinto que cero, las soluciones son L.I. Ahora tomando 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 2x 1 x
A(x) = 4x
xe cosh x dx = xe2x cosh x dx = e (9(x 1) + e2x (3x 1)) + C1
e 18
Z Z
1 2x 1 x
B(x) = 4x
e cosh x dx = e2x cosh x dx = e (3 + e2x ) + C2
e 6

Page 54 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 52

Problem 52

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

Z Z
y2 f e2x 4(x + 1) x
u1 = dx = dx = 4e (x + 2)
W (e , e2x )
x e3x

Z Z
y1 f ex 4(x + 1)
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

Page 55 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 53

Problem 53

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:

v0 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:

v0 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 56 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 53

Coeficientes indeterminados
Problem 54

y 00 + 3y 0 10y = 6e4x

Solución:

La solucion homogé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
5x
y = Ae + 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

5x 1
y = Ae + Be2x + e4x
3

Page 57 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 54

Aniquiladores
Problem 55
Resolver
y 00 + 5y 0 + 4y = 17sin(x) + ex

Solución:

En lo personal, no me gusta el método de los aniquiladores, pero bueno aca va.

Sea D una operador lineal tal que Dy = y 0 . Entonces la ecuació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 solución de la ecuación homogénea es yh = c1 e 4x + c2 e x . Ahora nos falta encontrar
una solución particular y estamos listos.
Podemos anular el lado derecho de la ecuación aplicando (D 1) y (D2 + 1).
Entonces:

(D + 4)(D + 1)(D 1)(D2 + 1)y = 0

.
La solució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 solución de la ecuació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 4x x
y= e + sin x cos x + c1 e + c2 e
10 34 34
, para dos constantes c1 y c2 cualquieras

Page 58 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 55

Modelos y aplicaciones
Problem 56
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 fricción, k constante del resorte, F(t), fuerza externa.

-Si b=0, no hay amortiguació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:
p c2
R= c1 + c2 , = tan↵
c1
Donde, R es la amplitud, ↵ el delay (phase shift)

-Si hay amortiguación, pero F(t)= 0, hay 3 casos

b2 4km = 0, Raı́ces reales iguales, amortiguación crı́tica

b2 4km > 0 Raı́ces reales diferentes, sobre amortiguación

b2 4km = 0 Raı́ces imaginarias diferentes, sub amortiguación

-Si no hay amortiguació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.

Page 59 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 57

Problem 57
Una función no negativa f (x) tiene la propiedad que su conjunto de ordenadas en un intervalo cualquiera
tiene un área proporcional a la longitud de arco de la gráfica correspondiente en el intervalo. Encuentre f (x).

Sea f : I ! R de clase C2 , tal que f (x) > 0. Luego, el área bajo la curva en algún intervalo cualquiera
[x0 , x1 ], está 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 aplicación lineal:

Z x1 Z x1 p
f (t)dt = k 1 + (f 0 (t))2 dt
x0 x0

Luego como ésto es válido en cualquier intervalo, los integrandos deben ser iguales, es decir:

p
f (t) = k 1 + (f 0 (t))2 (12)

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 funció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 solución L.I no siempre es mayor que cero y la
comprobación en (1) debe ser por separado, pues el principio de superposición no se cumple, se deduce
que:
x
f (x) = k cosh( ) k>0
k

Page 60 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 57

Cambio de variable
Problem 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

Solución:

z0
z0 = bz ) = b ) ln(z) = bx + ln(c) ) z = Ce bx
z
(Obviando las constantes de integración, que añadiremos al final)
Ası́:

y 0 + ay = e bx

de la ec. homogénea:

y 0 + ay = 0 ) y = Ke ax

Luego, suponemos como solución:

y(x) = u(x)e ax
) y 0 = u0 e ax
aue ax

reemplazando en la ecuación, se llega a:

e(a b)x
u0 = e(a b)x
)u=
a b
Finalmente, la solución general de la ecuación diferencial queda:

ke bx
y(x) = , k = cte.
a b

Page 61 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 59

Problem 59
Resolver
s ✓ ◆2
dy
1+
d2 y dx
=k
dx2 a x

donde a, k son constantes tales que a 2 R, k 6= 0, 1 y x < a

Solució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
p dv = dx
1 + v2 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 62 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 60

Problem 60

Usando la transformacion t = sen (x) convierta la ecuación diferencial

y 00 + tan (x) y 0 + cos2 (x) y = 0

en una ecuación con coeficientes constantes. Utilice esta transformacion para resolver la ecuacion diferencial.
Solució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

✓ ◆
dt d2 y dy dx dt d2 y dy d2 y dy
= cos(t) sin(t) = cos(t) sin(t) = 2 cos2 (t) sin(t)
dx dt2 dt dt dx dt2 dt dt dt

Reemplazando en la ecuació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
só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))

Page 63 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 60

Existencia y unicidad
Ejercicio 61
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 64 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 61

Mixtos
Problem 62

y 00 + tan(x)y 0 + cos2 (x)y = sin(x)cos2 (x)

Solución:

dt dt p
Sea x = arcsin(t). Entonces dx = p , 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 dy dt
= = 1 t = 1 t2 ·
dx2 dx dx dx dt dt dt dx
✓ ◆
t dy p d2 y p
= p + 1 t2 2 1 t2
1 t2 dt dt
dy d2 y
= t + (1 t2 ) 2
dt dt
t p
Como t = sin(x), entonces tan(x) = p y cos(x) = 1 t2 . Reemplazando, y usando la notación
1 t2
y = y(t),
t p
ty 0 + (1 t2 )y 00 + p 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 ecuación caracterı́stica
asociada es k 2 + 1 = 0 , k = ±i. Entonces, la solución general de la homogenea es¡br /¿

yh (t) = C1 cos(t) + C2 sin(t)

Se encuentra una solución particular de la no homogenea por el método de los coeficientes

indeterminados. Dado que g(t) = t (lado derecho de la ecuación), entonces la forma de la solució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 solución general de la no homogenea es

y(x) = C1 cos(sin(x)) + C2 sin(sin(x)) + sin(x)

Page 65 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 63

Problem 63

2
✓ ◆
4d y dy 1
4x + 8x3 + y = tan
dx2 dx 2x
1
Indicacion: use el cambio de variable x =
t
Solución:

1
Sea t = , luego:
x
dy dy dt dy 1
= = · 2
dx dt dx dt x
✓ ◆ ✓ ◆ ✓ ◆
d2 y d dy dt d dy dt d dy 1
= = = ·
dx2 dx dx dx dt dx dx dt dt x2
✓ ◆
dt d2 y 1 dy 2 dx d2 y 1 dy 2
= 2
· 2
+ 3
· = 2
· 4+ ·
dx dt x dt x dt dt x dt x3
Reemplazando en la ecuación inicial:
✓ ◆ ✓ ◆ ✓ ◆
d2 y 1 dy 2 dy 1 t
4x4 2
· 4+ · 3 + 8x 3
· 2 + y = tan
dt x dt x dt x 2
✓ ◆
t
4y 00 + y = tan
2
La ecuación homogenea asociada es:

4y 00 + y = 0

que tiene por solución:

yh = c1 cos(t/2) + c2 sin(t/2)

Usando el método de variación de parámetros, se puede hallar la solución particular:

yp = u1 cos(t/2) + u2 sin(t/2)

donde u1 = u1 (t), u2 = u2 (t).

Finalmente se obtiene:
✓ ◆ Z ✓ ◆ ✓ ◆ ✓ ◆" ✓ ◆ t
#
t t t t t 1 + sin 2
u2 (t) = 4 cos u1 (t) = 2 sin tan dt = 4 sin sin ln
2 2 2 2 2 cos 2t

La solución general de la ecuación es:

✓ ◆ ✓ ◆ ✓ ◆ ✓ ◆" ✓ ◆ t
#
t t t t t 1 + sin
y(t) = c1 cos + c2 sin 4 cos2 + 4 sin sin ln 2
2 2 2 2 2 cos 2t

1
Y volviendo a t = se obtiene la solución original.
x

Page 66 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 64

Problem 64

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 , quedará c1 (x) = ln(x)
Ahora, normalizando la ODE

3 0 4
y 00 y + 2 y = lnx
x x
Usando variacion de pará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

Page 67 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 65

Ejercicio 65

y 00 + y = t cos t.
(propuesto con C.I.)

La solució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 di↵erential 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 68 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 66

Ejercicio 66

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,

y p = 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 69 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 67

Ejercicio 67

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)
00 0
u sin(2x) + 2u 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 68

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 70 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 69

Ejercicio 69

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 70

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

Di↵erentiate: 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 71 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 71

Ejercicio 71

dy
= y + ex .
dx

Vemos que la solución particular es de la forma yp = Axex ., ya que la solución homogenea es L.D. con la
parte no homonegea de la ecuación.
d
(Axex ) = Axex + ex
dx
Axex + Aex = Axex + ex
A = 1.
La solución general es:

y = c1 ex + xex .

Ejercicio 72

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 =
x x x x
Ae +Be +Be +Bxe 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 72 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 72

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 73 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 73

Ejercicio 73

d2 y dy
2 3y = 0, y(0) = 0, y 0 (0) = 1.
dx2 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 74

2y 00 + 6y 0 + 5y = 0, y(0) = 1, y 0 (0) = 1/2

p
6± 62 4 · 2 · 5
m1,2 = .
2·2
3 1
m1,2 = ±i
2 2
3
⇣ ⇣x⌘ ⇣ x ⌘⌘
x
y(x) = e 2 A cos + B sin .
2 2
✓ ⇣ x ⌘◆
3 3x ⇣ ⇣x⌘ ⇣ x ⌘⌘ 3 A ⇣x⌘ B
y 0 (x) = e 2 A cos + B sin + e 2x 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
✓ ◆
3 0 1 B
y 0 (0) = e (1 · cos(0) + B sin(0)) + e0 sin(0) + cos(0) (13)
2 2 2
3 B 1
= + = . (14)
2 2 2
B=2
3
⇣ ⇣x⌘ ⇣ x ⌘⌘
y(x) = e 2x cos + 2 sin .
2 2

Page 74 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 75

Ejercicio 75

d2 y dy
+ 2y = 0, y(0) = 1, y(1) = 0
dx2 dx

y(x) = Aex + Be 2x

y(0) = Be0 = B = 1,
A=0
2x
y(x) = e .

Problem 76

d2 y
+ 25y = 2 sin(5t) + cos(t).
dt2

Resolviendo la homogé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), (15)
yp00 = 10A sin(5t) + 10B cos(5t) 25At cos(5t) 25Bt sin(5t) C cos(t) D sin(t). (16)

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 75 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 77

Problem 77
Determine el orden mı́nimo que debe tener la ODE asociada para que la solución sea:

y(x) = c1 t3 + c2 tet sin(t)

Octavo orden Aplicando los aniquiladores

D4 ((D 1)2 + 1)2 y(x) = 0.

Ejercicio 78
Resolver
x3 y 000 x2 y 00 + 2xy 0 2y = x3

ez = x
dy dy dz z dy
y0 = = =e
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

Page 76 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Ejercicio 78

Propuestos
Problem 79
2y 1
Resolver y 00 = (1 x)2 conociendo que es una solución homogénea
(1 x)2 (1 x)

Problem 80
Encontrar la ecuació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)

Problem 81
Modelar los siguientes sistemas de resortes

1. Dos resortes en serie, horizontal, sujentando una masa, sin fuerza externas ni amortiguaciones.

2. Dos resortes en paralelo, horizontal, sujentando una masa , sin fuerza externas ni amortiguaciones.

3. Dos resortes en serie, horizontal, sujentando masas diferentes, sin fuerza externas ni amortiguaciones.

4. Dos resortes en paralelo, horizontal, sujentando masas diferentes , sin fuerza externas ni amortiguacio-
nes. ¿Cuando hay resonancia?

5. Generalizar todo caso anterior, para n resortes, y m masas.

Page 77 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 81

Sistema de ecuaciones
Cambio de variables
Problem 82
(
x0 = x y x(0) = 4
Resuelva el PVI: 0
y = 5x y y(0) = 0

Solución:

Ahora reemplazando una ecuación en la otra (creeı́an que solo con el Wroskiano se podı́a resolver)

(
x0 = x y x(0) = 4
0
y = 5x y y(0) = 0
y=x x0
0 0
(x x ) = 5x (x x0 )
x0 x00 = 5x x + x0
x00 + 4x = 0
r2 + 4 = 0
r = ±2i
x(t) = c1 e2it + c2 e 2it

2it 0
y(t) = (c1 e2it + c2 e 2it
) (c1 e2it + c2 e ) =
2it 2it 2it 2it
= c1 e + c2 e 2ic1 e + 2ic2 e =
2it 2it
= c1 (1 2i)e + c2 (1 + 2i)e
(
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

=
2it 2it
= 5ie + 5ie
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 78 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 83

Problem 83
Resolver

x0 = x cos t + y sin t + z sin t

0
y = y cos t + z sin t
0
z = z cos t

Solución:

La primera ecuació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 ecuación y llegamos a:

y0 y cos(t) = c1 esin(t) sin(t)

La ecuación anterior es lineal con factor integrante e sin(t) , por lo que multiplicando por el factor
mencionado llegamos a:

sin(t) 0
(ye ) = 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 aná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 periódicas (sólo hay senos o cosenos en las soluciones).

Problem 83 continued on next page. . . Page 79 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 83 (continued)

La Matriz Fundamental no es más que la que satisfaga la relació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
0 1
Cos(t) Sin(t Sin(t)
E 0 (t) = @ 0 Cos(t) Sin(t) A E(t)
0 0 Cos(t)

Donde
0 0 1 0 1
x (t) x(t)
E 0 (t) = @y 0 (t)A E(t) = @y(t)A
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 serán 1 , 2 , 3 = Cos(t) Esto producirá una solució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 solució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 solución que se comporta perı́odicamente. siendo V1 , V2 , V3 y C1 , C2 , C3 los Vectores

Propios(tb conocidos como eigenvectores) y las constantes.

Page 80 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 83

Diagonalización
Problem 84
Encuentre la diagonalización de Jordan, del siguiente sistema (repaso Lineal)
0 1
1 2 3 0 0
B0 1 2 0 0C
B C
B C
A = B0 0 1 2 0C
B C
@0 0 0 2 1A
0 0 0 0 1

Si no es diagonalizable la matriz, se debe buscar un sistema ”semi”diagonalizable, para utilizar la matriz

exponencial.
El polinomio caracterı́stico es:

PA ( ) = (1 )4 (2 )

Sus raı́ces son 1 = 1 y 2 = 2 con multiplicidades algebraicas 4 y 1 respectivamente. 1 = 1

Comencemos con 1 , tenemos que hallar 4 vectores linealmente independientes, pues la multiplicidad de
1 es 4. Pero no valen 4 vectores cualquiera. Hay que hacer lo siguiente:
2 3
Hallar la cadena de ker de (A 1 Id) , (A 1 Id) , (A 1 Id) ,...hasta que la dimensión del último sea
la multiplicidad de la raı́z (4 en este caso).
0 1
0 2 3 0 0
B 0 0 2 0 0C
B C
B C
B := A 1 · Id = B0 0 0 2 0C
B C
@ 0 0 0 1 1A
0 0 0 0 0

Calculando el rango de esta matriz nos da rg(B)=4, luego su la dimensión del ker es n(B) = dim (R5 ) -
rg(B) = 5-4 =1. Resolviendo el sistema BX=0, obtenemos que todas las coordenadas excepto la primera
han de valer cero. Ası́ pues, los vectores del ker de B son:

Ker(A 1 · Id) = {(x, 0, 0, 0, 0) : x 2 R}

. Como la dimensión del Ker(A 1 · Id)) es 1, cualquier base de Ker(A 1 · Id) estará formada por un
único vector de Ker(A 1 · Id), linealmente independiente. Tomamos para formar la base, por ejemplo,
al vector canónico (1, 0, 0, 0, 0)

Ker(A 1 · Id) =< (1, 0, 0, 0, 0) >

0 1
0 0 4 6 0
B0 0 0 4 0C
B C
B C
C := (A 1 · Id)2 = B0 0 0 2 2C
B C
@0 0 0 1 1A
0 0 0 0 0

Problem 84 continued on next page. . . Page 81 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 84 (continued)

Realizando un proceso análogo al anterior obtenemos que el rango de C es 3, luego su la dimension del
Ker es n(C) = 2. Resolviendo el sistema CX=0 se obtiene que todas las coordenadas de los vectores de
Ker(A 1 · Id)2 valer cero, excepto las dos primeras. Como Ker(A Id) ⇢ Ker(A Id)2 , sabemos
que podemos expandir la base de Ker(A 1 · Id) para obtener una base de Ker(A 1 · Id)2 Elegimos
entonces el vector (0, 1, 0, 0, 0) . Ası́:

Ker(A 1 · Id)2 =< (1, 0, 0, 0, 0), (0, 1, 0, 0, 0) >

0 1
0 0 0 14 6
B0 0 0 4 4C
B C
3 B C
D := (A 1 · Id) = B0 0 0 2 2C
B C
@0 0 0 1 1A
0 0 0 0 0
El rango de esta matriz es rg(D)=2. Su dimension del ker es por tanto 3. Resolvemos el sistema DX=0 y
observamos que las dos últimas coordenadas valen 0. Expandemos la base de Ker(A 1 · Id)2 para
obtener la de Ker(A 1 · Id)3 , por ejemplo con el vector (0, 0, 1, 0, 0):

Ker(A 1 · Id)3 =< (1, 0, 0, 0, 0), (0, 1, 0, 0, 0), (0, 0, 1, 0, 0) >

0 1
0 0 0 14 14
B0 0 0 4 4 C
B C
B C
E := (A 1 · Id)4 = B0 0 0 2 2 C
B C
@0 0 0 1 1 A
0 0 0 0 0
En este caso, la dimension del ker de E es n(E)=4, y como la dimensión de Ker(A 1 · Id)4 ) no puede
ser superior a la multiplicidad algebraica del valor propio 1, que es 4, ya hemos llegado a la dimensión
máxima. Resolvemos el sistema EX=0 y concluı́mos que la suma de las últimas dos coordenadas es cero.
Ahora tomamos un vector v4 2 Ker(A 1 · Id)4 pero que no pertenezca a ninguno de los anteriores. Por
ejemplo,v4 = (0, 0, 0, 1, 1). Obtenemos ası́ la base deKer(A 1 · Id)4 :

Ker(A 1 · Id)4 =< (1, 0, 0, 0, 0), (0, 1, 0, 0, 0), (0, 0, 1, 0, 0), (0, 0, 0, 1, 1) >

Ahora hallar v3 , v2 y v1 es trivial: (y es lo que hacen en sus sistemas de ODEs sin darse cuenta)

v3 = (A 1 · Id)v4 = (0, 0, 2, 0, 0)

v2 = (A 1 · Id)2 v4 = (A 1 · Id)v3 = (6, 4, 0, 0, 0)

v1 = (A 1 · Id)3 v4 = (A 1 · Id)v2 = (8, 0, 0, 0, 0)

2 =2
0 1
1 2 3 0 0
B 0 1 2 0 0C
B C
B C
A 2 · Id = B 0 0 1 2 0C
B C
@ 0 0 0 0 1A
0 0 0 0 1

Ker(A 2 · Id) =< (14, 4, 2, 1, 0) > .

Problem 84 continued on next page. . . Page 82 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 84 (continued)

Ya tenemos los 5 vectores de la nueva base.La matriz de cambio de base es

0 1
8 6 0 0 14
B0 4 0 0 4C
B C
B C
P = B0 0 2 0 2C
B C
@0 0 0 1 1A
0 0 0 1 0

Para hallar la matriz de Jordan sólo hay que hacer las imágenes por A, de los vectores de la base de
Jordan, y expresarlos en dicha base:

Av1 = v1 ) (1, 0, 0, 0, 0)

Av2 = v1 + v2 ) (1, 1, 0, 0, 0)
Av3 = v2 + v3 ) (0, 1, 1, 0, 0)
Av4 = v3 + v4 ) (0, 0, 1, 1, 0)
Av5 = 2v5 ) (0, 0, 0, 0, 2)
Finalmente la matriz de Jordan es:
0 1
1 1 0 0 0
B0 1 1 0 0C
B C
B C
J = B0 0 1 1 0C
B C
@0 0 0 1 0A
0 0 0 0 2
1
Se cumple J = P AP

Page 83 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 85

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

✓ 0◆ ✓ ◆✓ ◆
x 3 4 x
= .
y0 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 2 1
= Aet + Be 5t .
y 1 2

Problem 86
*Resolver el sistema ✓ ◆
3 1
0 2 4
x = 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
t
x2 (t) = (v1 t + v2 )e

Problem 86 continued on next page. . . Page 84 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 86 (continued)

Problem 87
0
Resuelva x (t) = Ax(t),si
0 1
1 1 0 0
B0 1 0 0C
A=B
@0
C
0 2 1A
0 0 0 2

Por ser , la solución es directa

0 1
et tet 0 0
B0 et 0 0 C
eJt =B
@0
C
0 e2t te2t A
0 0 0 e2t

0 1
C1 et + C2 tet
B C3 et C
x(t) = B C
@C4 e2t + C5 te2t A
C6 e2t

Page 85 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 88

Problem 88
8
> 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
Solución:
Sea X la matriz
0 1
x
X=@ y A
z

El sistema se escribe matricialmente como:

0 1 0 1
1 0 0 1
X0 = @ 1 1 0 A X; X(0) = @ 1 A
1 1 1 1

La solució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 fá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 único vector propio asociado a dicho valor propio . . .
Entonces, obtenemos el vector propio para r = 1:
0 1
0 0 0
@ 1 0 0 A = 0v
1 1 0

Luego:

0 1 10 0 1
1 0 0
@ 2
A=@ 0 A= 3
@ 0 A
3 3 1

Podemos considerar:
0 1
0
=@ 0 A
1

y por ahora:

0 1
0
X(1) = @ 0 A et
1

Problem 88 continued on next page. . . Page 86 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 88 (continued)

Ahora tomamos una solución de la forma:

X= t e t + ⌘ et

derivando,

X0 = et + t e t + ⌘ t e t

Reemplazando en la ecuación inicial . . .

X0 = et + t et + ⌘ t et = A t et + A⌘ et

Luego:

(A I) = 0V , (A I) ⌘ =

La primera ecuación está resulta . . .

De aquı́ se desprende que :
0 1
0
⌘=@ 1 A
1

Llegamos a una solución:

0 1 0 1
0 0
X(2) = @ 0 A t et + @ 1 A e t
1 1

por último consideramos una solución de la forma:

t2 t
X= e + ⌘ t e t + ' et
2
Satisfaciendo las ecuaciones:

(A I) = 0V , (A I) ⌘ = , (A I) ' = ⌘

La primera y segunda ecuación ya han sido resueltas . . .

Luego, se encuentra que:
10
1
'=@ 0 A
1

Entonces:
0 1 0 1 0 1
0 2 0 1
t
X(3) = @ 0 A e t + @ 1 A t et + @ 0 A e t
2
1 1 1

Por lo tanto la solución general del sistema es:

0 1 20 1 0 1 3
0 0 0
X = C 1 @ 0 A e t + C 2 4 @ 0 A t et + @ 1 A e t 5 +
1 1 1

Problem 88 continued on next page. . . Page 87 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 88 (continued)

20 1 0 1 0 1 3
0 2 0 1
4 @ A t @ A @
C3 0 t
e + 1 t
te + 0 A et 5
2
1 1 1

y finalmente, aplicando la condición inicial:

C1 = 1, C2 = 1 , C3 = 1 ,
0
1 0 1 0 1
1 0 0
=
X(t) @ 1 A et + @ 1 A t et + @ 0 A t2 et
1 2 1/2

Problem 89
0
Resuelva x (t) = Ax(t), con
0 1
5 3 2
A= @ 8 5 4A
4 3 3

Utilizando el ejercicio anterior, de buscar la base canonica de Jordan

2 3
2 3
0 2 0 21 0 03 6 2 1 0
7
6 4 07 6 1 07
A = P JP 1 = 4 1 5 40 1 15 6 2 0 7
3 4 3 5
2 1 0 0 1 2 1
2 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
2 3
2 3 2 3 2 1 0
0 2 0 1 0 0 6 1 7
6 4 07 6 7
(t) = P eJt P 1 = 4 1 5 e t4
0 1 t 5 6 2 0 07
3 4 3 5
2 1 0 0 1 2 1
2 2
2 3
4t + 1 3t 2t
(t) = et 4 8t 1 6t 4t 5
4t 3t 2t + 1
2 3 2 3
x(t) c1 et (4t + 1) 3c2 et t 2c3 et t
Xh = 4y(t)5 = eAt = 4 8c1 et t c2 et (6t 1) 4c3 et t 5
z(t) 4c1 et t + 3c2 et t + c3 et (2t + 1)

Page 88 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 89

Matriz fundamental
Problem 90
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

✓ t ◆
e e t
(t) = t
e 3e t

Recordando la formula, debemos invertir la matriz fundamental

Z
1
XP = (t) (t)F (t)

✓ ◆ ✓ ◆
1 1 3e t e t 1 3e t e t
(t) = =
3e t et et e t et et 2 et et

✓ ◆✓ ◆ ✓ ◆ ✓ ◆
1 1 3e t e t
et 3e t et + e t et 4
(t)F (t) = = =
2 et et et et et et et 2e2t

Integrando
Z ✓ ◆
1 4t
(t)F (t) =
e2t

Finalmente
Z ✓ t ◆✓ ◆ ✓ ◆
1 e e t 4t 4tet et
Xp = (t) (t)F (t = t t 2t =
e 3e e 4tet 3et

Page 89 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 90

Principio de Duhamel
Problem 91

  
6/7 15/14 e2t 4
x0 = Ax(t) + f (t) = x(t) + , x(0) =
5/7 37/14 e t 1

2 3
6 t/2 1 3t 3 t/2 3 3t
67e + 7
e
7
e e
7 7
eAt = 42 2 3t 1 t/2 6 3t 5
et/2 e e + e
7 7 7 7
 t/2
3e + e3t
eAt x(0) = t/2
e 2e3t
2 3
6 (t s)/2 1 3(t s) 3 (t s)/2 3 3(t s)
6 e + e e e 7
eA(t s) = 4 7 7 7
2 (t s)/2 2 3(t s) 1 (t s)/2 6 3(t s) 5
7
e e e + e
7 7 7 7
2 3
6 (t s)/2 1 3(t s) 3 (t s)/2 3 3(t s)
6 e + e e e 7
eA(t s) · f (s) = 4 7 7 7
2 (t s)/2 2 3(t s) 1 (t s)/2 6 3(t s) 5
7
e e e + e
7 7 7 7
 2s
e
·
e s
2 3
Z t 3 2t 5 t 2 t/2 1 3t
6 7 e e e + e 7
eA(t s) · f (s) ds = 4 10 28 7 28
2t 13 t 2 t/2 1 3t 5
0 e e e e
21 42 21 14
Z t  t/2
At A(t s) 3e + e3t
x(t) = e x(0) + e · f (s) ds = t/2
0 e 2e3t
2 3
3 2t 5 t 2 t/2 1 3t
6 7e 28
e
7
e + e 7
28
+ 4 10 13 2 1 3t 5
2t t t/2
e e e e
21 42 21 14
2 3
3 2t 5 t 19 t/2 29 3t
6 7e 28
e + e + e 7
7 28
= 4 10 13 19 29 3t 5
e2t e t + et/2 e
21 42 21 14

Page 90 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 92

Problem 92
0
Resuelva x (t) = Ax(t) + f (t), con
0 1 0 1
5 3 2 sin(t)
A=@ 8 5 4A f (t) = @ 0 A
4 3 3 2

Z t
At
X(t) = e X0 + eA(t s)
F (s) ds
t0
2 3
2 3
0 2 0 21 0
3
0 6
2 1 0
7
6 1 1
A = P JP 1
=4 4 07
5 40 1
6
15 6 0 07
7
3 4 2 5
2 1 0 0 1 3
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
2 3
2 3 2 3 2 1 0
0 2 0 1 0 0 6 7
6 1
(t) = P eJt P 1 = 4 1 4 075 e 40 1 t 5 6
t 6 0 077
3 4 2 5
2 1 0 0 1 3
2 2 1
2
2 3
4t + 1 3t 2t
(t) = et 4 8t 1 6t 4t 5
4t 3t 2t + 1
2 3
4t1 3t 2t
1
(t) = e t 4 8t 6t + 1 4t 5
4t 3t 1 2t
2 t 3
e (4t + (4t 1) sin(t))
1
(t) · F (t) = 4 8e t t(sin(t) + 1) 5
t
e ( 4 sin(t)t 4t + 2)
Z
1
w= (t)F (t)
3 2
x(t)
X(t) = Xh + Xp = 4y(t)5 = eAt X0 + (t)w =
z(t)
[c1 et (4t + 1) 3c2 et t 2c3 et t 4t(2t + t sin(t) + (t + 1) cos(t) + 1) + 12t(2t + t sin(t) + (t + 1) cos(t) + 2)
1
2 (4t + 1)(8(t + 1) + (4t 1) sin(t) + (4t + 3) cos(t))], [8c1 et t c2 et (6t 1) 4c3 et t 8t(2t + t sin(t) + (t +
1) cos(t) + 1) + 4(6t 1)(2t + t sin(t) + (t + 1) cos(t) + 2) 4t(8(t + 1) + (4t 1) sin(t) + (4t +
3) cos(t))], [ 4c1 et t + 3c2 et t + c3 et (2t + 1) + 2(2t + 1)(2t + t sin(t) + (t + 1) cos(t) + 1) 12t(2t + t sin(t) +
(t + 1) cos(t) + 2) + 2t(8(t + 1) + (4t 1) sin(t) + (4t + 3) cos(t))]

Problem 92 continued on next page. . . Page 91 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 92 (continued)

Coeficientes indeterminados
Problem 93

✓ ◆ ✓ ◆ ✓ ◆
d~x 2 1 1 t 0
= ~x + e + t
dt 3 2 0 1

Ok, recordando coeficientes indeterminados para matrices, con la bases en R2, con la homogenea ya
resuelta.

~x(t) = ~atet + ~bet + ~ct + d~ + f~te t

+ ~g e t

✓ ◆
3/2
~a =
3/2
✓ ◆
~b = c1
c1 (1/2)
✓ ◆
0
~c =
2

✓ ◆
0
d~ =
1

✓ ◆
0
f~ =
0
✓ ◆
c2
~g =
3c2

Page 92 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 93

Modelos y aplicaciones
Problem 94
Dos tanques de 100 galones están inicialmente llenos de agua pura.
lb gal
Agua que contiene sal con una concentració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 función del tiempo.
Solución:

Sean xA (t), xB (t) las cantidades de sal en el estanque A y B, respectivamente, a través del tiempo.
Quedan claras las condiciones iniciales. xA (0) = xB (0) = 0.
La variació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

Problem 94 continued on next page. . . Page 93 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 94 (continued)

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 94 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 94

Existencia y unicidad

Page 95 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 94

Mixtos
Problem 95
Resolver
3 1 ax 2
ẋ = 2a + 2eax y, ẏ = e ax 3
a + (a2 + b)y ae y , x(0) = 0, y(T ) = 0
2 2

Reordenando el sistema en solo uno:

3 2
ẍ aẋ bẋ + 2ab = 0
4
Si z = x0
dz 3 2
az bz + 2ab = 0
dt 4
Dejando z como la variable independiente
dt 1
= 3 2
dz 4 az + bz 2ab

Trivialmente notamos que (super trivial)

⇣ ⌘
2 tanh 1 2p3az+2b
p
2
b 6a +b
t+C = p p
b 6a2 + b
Despejando z
⇣p p ⇣ p p ⌘ ⌘
2 b 6a2 + b tanh 12 b 6a2 + b(C + t) + b
z=
3a
Volviendo a la variable original e integrando trivialmente
⇣ ⇣ ⇣ p p ⌘⌘ ⌘
2 2 log cosh 12 b 6a2 + b(C + t) + bt
x= +D
3a

Page 96 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 95

Propuestos
Problem 96

✓ ◆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 p · ·p
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

Problem 97
x0 = y + x(1 x2 y2 )
(Freak one) Resolver
y0 = x + y(1 x2 y2 )

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 p =t+c
r2 1
r
p = Aet
r 2 1
r2
= Be2t
r2 1
1 2t
1 = Be
r2
1
2t
r = (1 Be ) 2

Las soluciones son:

1
2t
r(t) = (1 + ce ) 2 , ✓(t) = (t ↵).

Page 97 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 97

Series
Series de potencias
Problem 98
Resolver y 00 xy 0 + y = x cos x, dada las condiciones iniciales y(0) = 0, y 0 (0) = 2
Necesitamos las condiciones Iniciales de la ecuació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 en un entorno del punto x=a una función y=f(x)
admite ser expresada como solucion en serie de potencias:
1
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
(6) (4) (5)
y (x) = 3y (x) + xy (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 solución es de la forma :
1
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 que corresponde a la expansión en serie de potencias Mclaurin entorno a cero de: y(x) = x + sin x

Problem 98 continued on next page. . . Page 98 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 98 (continued)

Series de Frobenius
Problem 99
Utilizando el método de Frobenius, resolver la EDO:

4xy 00 + 2y 0 + y = 0

Solución:
1
X
Supongamos existe una solución de la forma: y(x) = Ak xk+⇢ con la condición de que A0 6= 0
k=0
Le daremos a la ecuació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, 8k 1
2 2
x 1
g(x) = ) g0 = 0, g1 = , gk = 0, 8k 2
4 4
Ahora, la funció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:
nX1
[(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)!
X1 X1
( 1)k A0 k ( 1)k p 2k p
Ası́, tenemos: y1 (x) = x = A0 ( x) = A0 cos( x)
(2k)! (2k)!
k=0 k=0
1
Ahora, para la segunda solució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)!
X1 X1
( 1)k A0 k+ 1 ( 1)k p 2k+1 p
Luego: y2 (x) = x 2 = A0 ( x) = A0 sin( x)
(2k + 1)! (2k + 1)!
k=0 k=0 p p
Finalmente, la solución de la EDO queda dada por: y(x) = A0 cos( x) + A0 sin( x)

Page 99 of 116
Ecuaciones Diferenciales Ordinarias [email protected] Problem 99

Propuestos
Problem 100
Resolver el sistema con ecuaciones, o con series.

y (4) + 3y 000 sen(t)y 00 + 8y 0 = t2 , y(0) = 1, y 0 (0) = 2, y 00 (0) = 3, y 000 (0) = 4

Problem 101
Resolver con Frobenius
z 2 f 00 zf 0 + (1 z)f = 0

✓ ◆
1 0 1 z 1 0 1 1
f 00 f + 2 f = f 00 f + f =0
z z z z2 z
1
X
f= Ak z k+r
k=0
1
X
f0 = (k + r)Ak z k+r 1

k=0
1
X
f 00 = (k + r)(k + r 1)Ak z k+r 2

k=0
1
X 1 ✓ ◆X
1
1X 1 1
(k + r)(k + r 1)Ak z k+r 2
(k + r)Ak z k+r 1
+ Ak z k+r
z z2 z
k=0 k=0 k=0
1
X 1 1 1
1X 1 X 1 X
= (k + r)(k + r 1)Ak z k+r 2
(k + r)Ak z k+r 1
+ 2 Ak z k+r Ak z k+r
z z z
k=0 k=0 k=0 k=0
1
X 1
X 1
X 1
X
(k + r)(k + r 1)Ak z k+r 2
(k + r)Ak z k+r 2
+ Ak z k+r 2
Ak z k+r 1

k=0 k=0 k=0 k=0

1
X 1
X 1
X 1
X
= (k + r)(k + r 1)Ak z k+r 2
(k + r)Ak z k+r 2
+ Ak z k+r 2
Ak 1z
k 1+r 1

k=0 k=0 k=0 k 1=0

1
X 1
X 1
X X1
= (k + r)(k + r 1)Ak z k+r 2
(k + r)Ak z k+r 2
+ Ak z k+r 2
Ak 1z
k+r 2

k=0 k=0 k=0 k=1

( 1
) 1
X X
= ((k + r)(k + r 1) (k + r) + 1) Ak z k+r 2
Ak 1z
k+r 2

k=0 k=1
( 1
) 1
X X
= (r(r 1) r + 1) A0 z r 2
+ ((k + r)(k + r 1) (k + r) + 1) Ak z k+r 2
Ak 1z
k+r 2

k=1 k=1
( 1 1
)
X X
= (r 1)2 A0 z r 2
+ (k + r 1)2 Ak z k+r 2
Ak 1z
k+r 2

k=1 k=1
1
X
= (r 1)2 A0 z r 2
+ (k + r 1)2 Ak Ak 1 z k+r 2

k=1
Ak 1
Como r=1 contribuye como raı́z. (k + 1 1)2 Ak Ak 1 = k 2 Ak Ak 1 = 0 ! Ak = k2

Page 100 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 101

Laplace
Inversa
Problem 102
Resolver con Transformada de Laplace

y” + 4y = 16x, y(0) = 3, y 0 (0) = 6

Solució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) =
s 2 + a2

a
L [sen (a · t)](s) =
s2 + a2

2 0
L [f n (t)](s) = sn L [f (t)](s) sn 1
f (0) sn f (0) ... s · fn 2
(0) fn 1
(0)

Problem 102 continued on next page. . . Page 101 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 102 (continued)

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 ·
s2

16
s2 ' (s) 3s + 6 + 4' (s) =
s2

16
' (s) s2 + 4 = + 3s 6
s2

16 3s 6
' (s) = + 2
s2 2
(s + 4) (s + 4) (s2 + 4)

  
1 1 s 2
' (s) = 4 + 3 3
s2 2
s +4 2
s +4 s2 +4

   
1 1 1 1 2 1 s 1 2
L [' (s)](t) = 4L 2L + 3L 3L
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 102 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 103

Problem 103
(
x0 = x y x(0) = 4
Resuelva el PVI: 0
y = 5x y y(0) = 0

Solución:

Por Laplace
Sea L(x) = X y L(y) = Y

sX x(o) = X Y (17)
sY y(o) = 5X Y (18)

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=
s 2 + 22 s + 22 s 2 + 22
Aplicando la transformada inversa . . .

x(t) = 4 cos 2t + 2 sin 2t, y(t) = 10 sin 2t

Page 103 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 103

Delta de Dirac
Problem 104

Resolver el PVI, donde corresponde al delta de dirac

y 0 + y = (t 1), y(0) = 1

Solució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))

s
sY (s) y(0) + Y (s) = e
s
e +1 e s 1
Y (s) = = +
s+1 s+1 s+1
Aplicando la transformada inversa:
(t 1) t
y(t) = e u1 (t) + e

y(t) = e1 t u1 (t) + e t

Page 104 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 104

Heaviside
Problem 105
¿Como resolveria este ejercicio?
(
00 t, 0t<3
y + 4y = , y(0) = 0, y 0 (0) = 0
1, 3t<1

Solución:

Hay 2 formas, tomar Laplace por partes

Z ( ! Z 3 Z 1
1
t, 0  t < 3
F (s) = e st dt = te st dt + e st
dt
0 1, 3  t < 1 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)) = +
s2 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 105 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 105

Convolución
Problem 106
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) =
(s2 + 4)2 s2 + 4 s2 + 4
3s
1 3 3
L () = 3 cos (2t), L 1 ( 2 ) = sin (2t)
+4 s2 s +4 2
Aca la dificultad, hay una convolucion entremedio
Z
24 1 t
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
1 24 1 3
L = 24 · sin(2t 2⌧ ) sin 2⌧ d⌧ = (sin 2t 2t cos 2t).
(s2 + 4)2 4 0 2

Page 106 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 106

Ecuaciones integro-diferenciales
Problem 107

Resuelva:
Z t
f (u)
1+t= p du
0 t u
Hint:

(↵ + 1)
L [t↵ ](s) = ;s > 0 ;↵ > 1
s↵+1

✓ ◆
1 p
= ⇡
2

(n + 1) = n (n)
Z 1
Donde es la función gamma (n) = tn 1 e t dt
0
Solución:

Notar la semejanza de la función Gamma con la transformada de Laplace. Utilizando la convolución:

Z t ✓ ◆
1 1
f (u) · p du = f⇤p (t)
0 t u t
Aplicando la transformada de Laplace a ambos lados de la ecuación:

✓ ◆
1 1 1
+ 2 = L(f ) · L p
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 (1/2)
+ 2 (1/2)
s1+( 1/2) s1+(1/2)

1 1/2 2 1/2
f (t) = t + t
⇡ ⇡

Page 107 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 107

Mixtos
Problem 108
Probar que toda función periodica, con T > 1 tiene el Laplaciano de la forma

RT st
e F (t) dt
0
L{F (t)} = sT
dt
1 e

Z 1 1 Z
X (k+1)T 1 Z
X T
st st s(t+kT )
e F (t) dt = e F (t) dt = e F (t + kT ) dt.
0 k=0 kT k=0 0

Z T Z T 1
X
st sT 2sT st nsT
L{F (t)} = F (t)e dt · 1 + e +e + ... = F (t)e dt · e =
0 0 n=0

Z T
1 st
= sT
F (t)e dt
1 e 0

Problem 109
Resolver con Laplace:
1
X
00
y +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
2⇡s
, notamos que es una suma geometrica e <1
8 9
e2⇡s 1 X <X =
2⇡ks
= = e = L 2k⇡ (t) .
e2⇡s 1 1 e 2⇡s : ;
k 0 k 0

Ahora, simplemente utilizamos convolución:

0 1
X X X
y(t) = @ 2k⇡ (t)
A ⇤ sin t = ( 2k⇡ (t) ⇤ sin t) = sin(t 2k⇡)
k 0 k 0 k 0

Page 108 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 110

Problem 110
8 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

Solución:

La ecuación se reescribe como sigue:

y0 4y + 8(y ⇤ 1)(t) = (t 3) + u4 (t)

Sea Y = L(y). Aplicando la transformada de Laplace a ambos lados:

4s
1 3s e
sY y(0) 4Y + 8Y =e +
s s

3s 4s
8 se +e
Y (s 4+ )=
s s

se 3s e 4s
Y = +
(s 2)2 + 22 (s 2)2 + 22

Considerando el primer té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). Análogo para el segundo té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 109 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 111

Problem 111
Encontrar el Laplace inverso de las siguientes expresiones

1 1 1 s 2 + a2
s4 , , 2, , {s log 2 }
(s 1)2 (s2 + 1) (s2 2
+a ) n s a2

✓ ◆ ✓ ◆
1n! 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 1 t 1
3)L = (sin !t !t cos !t) , cos t + sin t
(s2 + ! 2 )2 2! 3 2 2

1 1 dn 1 est 1 dn 1
est
4)L {F (s)} = lı́m n 1
(s ia)n 2 2 n
+ lı́m (s + ia)n
s!ia (n 1)! ds (s + a ) s! ia (n 1)! dsn 1 (s2 + a2 )n

1 sin(ax) sin(ax)
oL 1 2 2 n
= ⇤ ··· ⇤
(s + a ) | a {z a }
n times

s 2 + a2
Bonus)f (x) : F = s log
s 2 a2
Z x
F s 2 + a2
f (t) dt :
= log 2
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 110 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 111

Propuestos
Problem 112
Resolver la ecuación integro-diferencial
Z t
dy 10w
+ 25 y(t w)e 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) = 2
=
s(s + 5) 5s (s + 5)2 5(s + 5)
1
L (y(s))

Problem 113
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)

Page 111 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 113

Sistemas dinámicos
Hartman-Grobman
Problem 114
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 112 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 114

Sistemas Hamiltonianos

Page 113 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 114

Mixtos
Ejercicio 115
Encontrar una funció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 116
Encontrar puntos de bifuración: x0 = x x3

p
x = 0, x = ±

Si = 0 hay bifuración.

Problem 117
Encontrar puntos de bifurcació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 bifurcación. Pasando a polares par encontrar el ciclio lı́mite estable.
dy
= ur r3 , d✓/dt = 1.
dt
p
De la primera ecuació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).

Page 114 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 118

Problem 118
Encontrar puntos de estabilidad del Oscilador de Van der Pol, y encontrar una funcion Lyapunov débil

d2 x dx
+ (x2 1) +x=0
dt2 dt
2
Vea que una función Lyapunov fuerte es: W (x, y) = (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 función dé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 función fuerte.

Page 115 of 116

Ecuaciones Diferenciales Ordinarias [email protected] Problem 118

Propuestos
Problem 119
Resolver el sistema Hamiltoniano
dx
= y,
dt
dy
= sin x.
dt

Problem 120
Calcular las trayectorias del Depredador-Presa
dx
= x + xy,
dt
dy
= y xy.
dt
¿Cual es la diferencia con el Hamiltoniano en este caso? ¿Como puedo hacerlo Hamiltoniana? Hint: x =
ex1 , y = ex2

Problem 121
Encontrar puntos de equilibrio.

x0 = 2xy
0
y = x + y + xy y3

Problem 122
Indique condiciones para las constantes para que el sistema tenga equilibrios definidos.
dx
= ↵x xy
dt
dy
= xy y
dt

Problem 123
Encontrar equilibrios para el sitema hamitoliano x0 = xsin(y) + 2y, y 0 = cos(y)

Page 116 of 116