Ecuaciones diferenciales

ordinarias
EDO

Carlos Yanes Pérez

 Índice general

I Tema 1: Estudio preliminar

1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

2 Estudio cualitativo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

II Tema 2: EDOs que sabemos resolver

3 EDOs de variables separadas (VS) . . . . . . . . . . . . . . . . . . . . . . 12
3.1 Ley del enfriamiento de Newton 13

4 EDOs homogéneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

5 EDOs lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
5.1 Lineales homogéneas 15
5.2 EDOs lineales completas (o afines) 15
5.2.1 Interpretación de las EDOs lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

6 EDOs de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
6.1 Bernoulli generalizado 20

7 EDOs de Ricatti . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
 3

8 EDOs exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

9 EDOs reducibles a exactas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

III Tema 3: Existencia y unicidad de EDOs

10 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28

11 Iterantes de Picard . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

Anexo 1
12 Ejercicios examen temas 1 y 2 . . . . . . . . . . . . . . . . . . . . . . . . . . 33
I
Tema 1: Estudio preliminar

1 Introducción . . . . . . . . . . . . . . . . . . . . . . . 5

2 Estudio cualitativo . . . . . . . . . . . . . . . . 9
 1. Introducción

Definición 1.0.1 Recibe el nombre de Ecuación Diferencial Ordinaria (EDO) toda

identidad que involucre a la variable independiente x perteneciente a un intervalo (a, b)
con −∞ ≤ a < b ≤ ∞ y a la función incógnita y junto con sus sucesivas derivadas
y′ , y′′ , . . . , y(n aunque no necesariamente deben estar todas presentes.
 
F x; y, y′ , y′′ , . . . , y(n = 0, x ∈ (a, b) (∗)

Veamos a continuación algunos ejemplos de EDO donde y = y(x)

y′ = 0 es claro que la solución para esta EDO es y = C ∈ R

y′ = f (x) donde f : I −→ R cuya solución es y(x) = f (x) +C; C ∈ R
R

y′ = y tiene como solución y = kex

2
y′ = 2xy tiene por solución y = Cex
y′′ = y tiene una solución compuesta por la suma de exponenciales y = C1 ex +C2 e−x
y′′ = −y la solución será de funciones trigonométricas y = C1 cos(x) +C2 sin(x)
y′′ + 1 = 0 no tiene solución
(y′ )2 + y2 = 0 tiene como única solución a y = 0

Nota 1.1. Nótese que xy′′ − y = 0 no es, en general, una EDO, pues para x = 0 no
comparece ninguna derivada de y en la ecuación.

En este curso trabajaremos con la ecuación despejando la derivada de mayor orden

obteniendo así la forma explícita o normal de la EDO:
 
y(n = G x; y′ , y′′ , . . . , y(n−1
6 Capítulo 1. Introducción
Definición 1.0.2 Se denomina orden de una EDO al mayor orden de derivación que
comparece en ella.

Definición 1.0.3 Se denomina solución (particular) de la EDO (*) al par (g, J) donde
g : J −→ R; J = (c, d) tal que
 
′ ′′ (n
F x; g(x), g (x), g (x), . . . , g (x) = 0, x ∈ (c, d) ⊆ (a, b)

Para comprobar si una función g es solución de una EDO nos basta con sustituir.
Definición 1.0.4 Al conjunto de todas las soluciones de una EDO se le denomina
solución general. Normalmente, esta es una familia n-paramétrica si la EDO es de orden
n.

Nota 1.2. En algunas ocasiones podrá existir alguna solución que no pueda ser recogida
por la solución general, estas serán las soluciones singulares.

Las EDOs surgen para facilitar el modelaje de fenómenos físicos y su resolución es

interesante para intentar predecir cómo se comportará este en el futuro o cómo se comportó
en un pasado, por ello hay casos en los que conviene partir de un dato experimental, a esto
lo llamaremos observación o valor inicial, para que las predicciones sean más certeras. De
esta forma dada una EDO F(t; x, x′ ) = 0, ∀t ∈ (a, b) mediremos en un instante que nos
interese (t0 ) el valor de la variable objeto de estudio x(t0 ) = x0 este será el valor inicial.
Definición 1.0.5 Al sistema conformado por una EDO F(t; x, x′ ) = 0, ∀t ∈ (a, b) y
un dato inicial x(t0 ) = x0 lo denominaremos problema de valores iniciales (PVI) o
problema de Cauchy (PC). Se tendrá una solución para estos problemas cuando la
solución de la EDO verifique además el dato inicial.

Considérese x′ = f (t, x) donde f : Ω ⊆ R2 → R donde Ω es un dominio (abierto y

conexo) no necesariamente acotado. Entonces el PVI se dará como

x′ = f (t, x)
x(t0 ) = x0

Esto en el caso de primer orden, si no f quedaría en función de las derivadas pertinentes

tras haber despejado la de mayor orden.
Definición 1.0.6 El par (x(t), I) es solución del PVI

x′ = f (t, x)
x(t0 ) = x0

si verifica
1. x es derivable en I
2. (t, x(t)) ∈ Ω; ∀t ∈ I
3. x′ (t) = f (t, x(t)); ∀t ∈ I
4. x(t0 ) = x0
 7

Teorema 1.0.1 — Teorema local de existencia y unicidad.. Sea Ω ⊆ R2 un dominio

y f (t, x) ∈ C (Ω) tal que f admite derivada parcial continua respecto de x en Ω ( fx ∈
C (Ω)). Entonces dado (t0 , x0 ) ∈ Ω existe un intervalo I tal que t0 ∈ I y una única función
definida en dicho intervalo satisfaciendo

x′ = f (t, x)
x(t0 ) = x0

Implícitamente en un PVI (de primer orden) se obliga a que la solución tenga una
pendiente en el valor inicial, además esta se propaga como f (t, x(t)). De esta forma al
darse la unicidad, por el resultado anterior, la solución es de clase C 1 .

Bien es cierto que aunque nuestra EDO sea una función f buena (clase infinito) la
solución puede no ser tan bondadosa (buscar la solución al PVI y′ = y2 , con y(x0 ) ̸= 0).

Hay otros casos, en los que puede existir solución, pero no ser única, véase el siguiente
PVI:
√
x′ = 2 x
x(0) = 0

Claramente x ≡ 0 es solución, pero x = t 2 también es solución. Es más, este PVI tiene

infinitas soluciones ya que x = (t − a)2 , t ≥ a. No hay unicidad, pero hay solución. Por
tanto debemos dar algunos resultados que nos permitan estudiar por separado existencia y
unicidad:
Teorema 1.0.2 — Teorema local de existencia (Cauchy-Peano). Sea Ω ⊆ R2
dominio y f ∈ C (Ω). Entonces para todo (t0 , x0 ) ∈ Ω el problema de valores iniciales
tiene al menos una solución definida en un cierto intervalo I, tal que t0 ∈ I.

Definición 1.0.7 Diremos que una EDO es lineal cuando se tenga linealidad respecto a
la segunda variable

f (t; λ x1 + µx2 ) = λ f (t; x1 ) + µ f (t; x2 ), ∀λ , µ ∈ R

Definición 1.0.8 Diremos que una EDO es autónoma si f (t, x) no depende de t. Es

decir si la EDO queda definida como x′ = f (x).

Teorema 1.0.3 — Lema. x(t) es solución del PVI

x′ = f (t, x(t))
x(t0 ) = x0

si, y solo si
1. x es continua para un cierto intervalo I
2. (t, x(t)) ∈ Ω para t ∈ I
8 Capítulo 1. Introducción
Rt
3. x(t) = x0 + t0 f (s, x(s)) ds

Demostración. Recordemos que por definición el par (x(t), I) es solución del PVI si
verifica
1. x es derivable en I
2. (t, x(t)) ∈ Ω; ∀t ∈ I
3. x′ (t) = f (t, x(t)); ∀t ∈ I
4. x(t0 ) = x0
"⇐"
Los dos primeros puntos del resultado anterior los tenemos por definición falta por ver
el tercero. Es claro que la solución integral verifica que x(t0 ) = x0 falta ver que se verifica
la condición 3. de la definición anterior, pero derivando la ecuación integral se tiene que
derivando x′ (t) = f (t, x(t)) para todo t ∈ I.
"⇒"
En este caso solo tenemos que comprobar que la EDO verifica la ecuación integral,
para ello integramos:
Z t Z t Z t
′
x (s) ds = f (s, x(s)) ds =⇒ x(t) = x0 + f (s, x(s)) ds
t0 t0 t0
■
Dada una familia uniparamétrica de funciones podemos encontrar la EDO que dicha
familia verifica eliminando el parámetro de la ecuación.
■ Ejemplo 1.1 Dada la ecuación x2 + y2 + 2cx = 0 encontrar la EDO asociada. Derivando
implícitamente se tiene que

y2 − x2
x2 + y2 − 2x(x + yy′ ) = 0 =⇒ y′ =
2xy
■
 2. Estudio cualitativo

Aún sin saber resolver una EDO podemos extraer información de la misma, por ejemplo
en un PVI sabemos que se está imponiendo una cierta pendiente para la recta tangente de
la solución en el punto (t0 , x0 ). Pero podemos ir más allá, dada una EDO en su dominio
Ω podemos obtener un campo de direcciones con las pendientes en cada punto; en lugar
de trabajar punto por punto nos preocuparemos de diferenciar regiones en Ω en las que
coincidan las pendientes de la solución para ello trabajaremos con el conjunto de las
isoclinas:

Φm = {(t, x) ∈ Ω : f (t, x) = m}

Donde m indica la pendiente que estamos estudiando, asimismo m ∈ Im(f ). De esta

forma el lugar geométrico obtenido en el estudio de cada isoclina nos indica la pendiente
que tendrá la solución en dicho lugar geométrico.
■ Ejemplo 2.1 Hallar las isoclinas de x′ = x + t 2 = f (t, x). Claramente el dominio de la
función en todo el plano, es más por ser un polinomio en dos variables es de clase C ∞ en
el plano. En este caso la imagen de f es todo R por tanto el valor de m es cualquier valor
real. Calculemos las isoclinas con m = 0 y m = 1:

Φ0 = {(t, x) ∈ Ω : f (t, x) = 0} = {(t, −t 2 ) : t ∈ R} ≡ x = −t 2

Por tanto la única solución de esta EDO (sabemos que es única por el teorema de
existencia unicidad) sobre la parábola x = −t 2 tendrá pendiente cero. Veamos para m = 1

Φ1 = {(t, x) ∈ Ω : f (t, x) = 1} = {(t, −t 2 + 1) : t ∈ R} ≡ x = −t 2 + 1

En este caso la solución tendrá pendiente m = 1 sobre la parábola x = −t 2 + 1. ■

10 Capítulo 2. Estudio cualitativo

Tras dibujar en el plano algunas isoclinas y observar las pendientes obtenida debería-
mos de ser capaces de esbozar el comportamiento de la solución de la EDO (campo de
direcciones en GeoGebra aquí.)
II
Tema 2: EDOs que sabemos
resolver

3 EDOs de variables separadas (VS)

12
3.1 Ley del enfriamiento de Newton

4 EDOs homogéneas . . . . . . . . . . . . . . 14

5 EDOs lineales . . . . . . . . . . . . . . . . . . . . . 15
5.1 Lineales homogéneas
5.2 EDOs lineales completas (o afines)

6 EDOs de Bernoulli . . . . . . . . . . . . . . . . 19
6.1 Bernoulli generalizado

7 EDOs de Ricatti . . . . . . . . . . . . . . . . . . . 21

8 EDOs exactas . . . . . . . . . . . . . . . . . . . . 23

9 EDOs reducibles a exactas . . . . . 25

 3. EDOs de variables separadas (VS)

Son las EDOs de la forma

x′ = p(t)q(x)

Teorema 3.0.1 — Teorema de existencia-unicidad (VS). Sean p(t) ∈ C ((a, b)) y

q(x) ∈ C ((c, d)) con q(x) ̸= 0 se dará que el PVI

x′ = p(t)q(x)x(t0 ) = x0
tiene solución única en R = (a, b) × (c, d) con t en un cierto intervalo I.

Demostración. Veamos primero a lo que tenemos que llegar. Supongamos que tenemos
solución (x) para el PVI. Usando que q(x) ̸= 0:

x′ (t)
= p(t); t ∈ I, t0 ∈ I
q(x(t))
Sabemos que si x(t) es solución entonces es una función continua, por lo que q(x(t))
es continua por ser composición de funciones continuas, x′ (t) será continua (por el Th.
E-U local), ahora, sabemos que el cociente de funciones continuas donde el denominador
no se anula es continuo. Por último, de nuestras hipótesis tenemos que p(t) es continua
(esto en R). Puesto que tenemos funciones continuas podemos integrar (sentido Riemann):
Z t ′ Z t
x (s)
ds = p(s) ds
t0 q(x(s)) t0

Aplicando el cambio de variable η = x(s) se tiene

 3.1 Ley del enfriamiento de Newton 13

Z x(t) Z t
dη
= p(s) ds
x0 q(η) t0
Ya hemos visto a dónde tenemos que llegar, pero hemos asumido que teníamos solución
(usamos la tesis para darnos una idea de lo que hacer, ahora debemos justificar que lo
anterior es correcto). Entonces proponemos la función
Z x Z t
dη
F(t, x) = − p(s) ds ≡ 0
x0 q(η) t0

Sabemos que el conjunto de ceros de F es no nulo, pues F(t0 , x0 ) = 0. Veamos ahora

que se cumplen las hipótesis del Teorema de la función implícita para garantizar que
x = x(t):
F es continua (?)
∂F
∂ x (t0 , x0 ) ̸= 0 (?)
La respuesta a la primera pregunta es clara, F es continua. Para la segunda basta con
derivar

∂F 1
(t, x) = ̸= 0, ∀(t, x) ∈ R2
∂x q(x)
Por lo que el teorema nos permite concluir que efectivamente se puede dar x en función
de t en un cierto entorno: x = x(t), t ∈ I.
■

3.1 Ley del enfriamiento de Newton

Para modelizar el cambio de temperatura de un cierto objeto podemos usar el siguiente
PVI:

T ′ = −k(T − Ta ) (k > 0)
T (0) = T0

Nótese que dependiendo de si el proceso físico a modelizar es de aumento de tempera-

tura o de disminución de la misma deberemos de alterar el planteamiento de la ecuación.

Realmente podríamos enunciar esta ley como

T ′ = k(T − Ta ) (k ∈ R \ {0})
T (0) = T0

Siempre que seamos capaces de inferir el comportamiento de la k para cada caso.

 4. EDOs homogéneas

Son EDOs de la forma

x
x′ = f (t ̸= 0)
t
x
Para la resolución de estas ecuaciones haremos uso del cambio de variable: u = t

x = tu =⇒ x′ = u + tu′
Por lo que se da

f (u) = u + tu′
Ahora tenemos una EDO en u con variable independiente t, despejando la derivada:
1
u′ =
· ( f (u) − u) = p(t)q(u)
t
Tenemos ahora una EDO de variables separadas.

Nota 4.1. Para resolver EDOs homogéneas las transformamos en variables separadas.

Teorema 4.0.1 — Teorema de existencia-unicidad (hom). Sean f ∈ C (R) con

t ̸= 0 se dará que el PVI
x
′
x =f x(t0 ) = x0
t
tiene solución única en un rectángulo, R, con t en un cierto intervalo I que no contiene
al cero.
 5. EDOs lineales

5.1 Lineales homogéneas

Diremos que una EDO es lineal homogénea cuando f sea lineal respecto de la segunda
variable. Además, la EDO tendrá la forma

x′ = a(t)x
Como vemos f (t, x) es una recta que pasa por el origen.
Para que se dé el teorema local sabemos que tendremos solución única debemos exigir
que a(t) ∈ C ((a, b)).
De esta forma se garantizará que f ∈ C 1 ((a, b) × R) siendo nuestro dominio una banda.
Es más, veremos que la solución será maximal en el intervalo de continuidad de a(t).
A ojo de buen cubero podemos afirmar que la solución de estas EDOs será
R
a(s) ds
x(t) = k · e (k ∈ R)
Para el caso de un PVI: x(t0 ) = x0 se tendría que
Rt
t0 a(s) ds
x(t) = x0 · e
Con la solución siempre definida allí donde lo esté a(t).

5.2 EDOs lineales completas (o afines)

Diremos que tenemos una EDO lineal completa si

x′ = a(t)x + b(t)
16 Capítulo 5. EDOs lineales

Donde a(t), b(t) ∈ C ((a, b)), basta con esto para que se verifique el teorema de
existencia-unicidad local. Ahora obtendremos la solución para estas ecuaciones.

Analizando la solución homogénea:

Rt
− t0 a(s) ds
x(t) · e = x0
Ahora, derivando
Rt
− t0 a(s) ds
(x′ − a(t)x) · e =0
Por tanto hemos obtenido la derivada de un producto:

− tt a(s) ds ′
 R 
x(t) · e 0 =0
Como vimos, el factor
Rt
− t0 a(s) ds
e
fue crucial para obtener la solución en el caso homogéneo, veamos si ahora nos es de
utilidad; multipliquemos nuestra EDO por dicho factor:
Rt Rt Rt
− t0 a(s) ds − t0 a(s) ds − t0 a(s) ds
x′ e = a(t)xe + b(t)e
Reordenando:
Rt Rt Rt
− t0 a(s) ds − t0 a(s) ds − t0 a(s) ds
x′ e − a(t)xe = b(t)e
Ahora nos percatamos de que en el primer miembro tenemos la derivada de un producto,
concretamente se da que

− tt a(s) ds ′ − t a(s) ds
 R  R
x(t) · e 0 = b(t) · e t0
Si integramos tenemos
Z t Rs ′ Z t Rs
− t0 a(u) du − t0 a(u) du
x(s) · e ds = b(s) · e ds
t0 t0
entonces
Rt Z t Rs
− t0 a(s) ds − t0 a(u) du
x(t) · e − x0 = b(s) · e ds
t0
Despejando, la solución es
Rt Rt Z t Rs
t0 a(s) ds t0 a(s) ds − t0 a(u) du
x(t) = x0 · e +e b(s) · e ds
t0
Por tanto
Rt Z t Rt
t0 a(s) ds a(u) du
x(t) = x0 · e + b(s) · e s ds
t0

La solución general de la EDO lineal completa que estará definida allí donde lo estén
a(t) y b(t).
 5.2 EDOs lineales completas (o afines) 17

Nota 5.1. Debemos percatarnos que al expresar la EDO lineal como:

− t a(s) ds ′ − t a(s) ds
 R  R
x(t) · e t0 = b(t) · e t0

Podemos llegar intuitivamente a la fórmula de la solución, sin memorizarla.

Teorema 5.2.1 — Teorema de existencia-unicidad (lin, GLOBAL). Dados a(t), b(t) ∈

C ((a, b)) el PVI

x′ = f (t, x) = a(t)x + b(t)

x(t0 ) = x0

tiene solución única en la banda B = (a, b) × R, dicha solución viene dada como
Rt Z t Rt
t0 a(s) ds a(u) du
x(t) = x0 · e + b(s) · e s ds
t0

5.2.1 Interpretación de las EDOs lineales

Puesto que tenemos una función f lineal podemos construir una aplicación lineal:
L: C 1 ((a, b)) −→ C ((a, b))
x −→ L(x) := x′ − a(t)x
Donde el núcleo es

N = ker(L) = (x ∈ C 1 ((a, b)) : L(x) = 0



Es decir, las soluciones homogéneas de la EDO. Además, se da que

Rt
t0 a(s) ds
N =< e >
Que es un subespacio vectorial de dimensión uno.
Por otro lado, resolver la EDO lineal completa equivale a encontrar L(x) = b(t), es
decir, calculamos el conjunto imagen de L. De esta forma se tendrá que una solución
general de la completa es de la forma

xg = xgh + x p
Es decir, la solución general de la completa es la suma de la solución general de
la homogénea más una solución particular de la completa; pensando en la fórmula de
las dimensiones de espacios vectoriales encontramos cierta similitud. . . Pero es más,
conocemos la forma de xgh y de x p , de modo que
Rt Z t Rt
t0 a(s) ds a(u) du
xg = k · e + b(s) · e s ds (k ∈ R)
t0
18 Capítulo 5. EDOs lineales

Ejercicio 5.1 Dadas y1 e y2 soluciones particulares de una EDO lineal, obtener la

solución general para dicha ecuación.

SOLUCIÓN:
Dada la EDO lineal general L[y] = b(x) sabemos que como y1 e y2 son soluciones
particulares se dará que L[y1 ] = b(x) y que L[y2 ] = b(x) por tanto, restando

L[y1 ] − L[y2 ] = 0

Por la linealidad del operador L:

L[y1 − y2 ] = 0

Entonces
k · (y1 − y2 ), k ∈ R
es solución de la EDO lineal homogénea, por tanto la solución general será

yg = k · (y1 − y2 ) + y1

Nota 5.2. También serían soluciones válidas las siguientes yg = k · (y1 − y2 ) + y2 ,

yg = k · (y2 − y1 ) + y1 o yg = k · (y2 − y1 ) + y2

RESUMEN: Para resolver las EDOs lineales debemos seguir los siguientes pasos Si
nuestra ecuación es x′ = a(t) · x + b(t)
1. Despejar x′ si es necesario (la derivada tiene que estar sola).
2. Pasar al mismo miembro de la igualdad lo que dependa de la segunda variable:
x′ − a(t) · x = b(t)
3. Multiplicar todo por una función ρ(t) > 0
ρx′ − ρa(t) · x = ρb(t)
4. Buscamos expresar el lado izquierdo como la derivada de un producto:
(x · ρ)′ = ρb(t) =⇒ x′ ρ + xρ ′ = ρb(t)
5. Comparamos la ecuación del paso 3. Y la del 4. vemos que el término que acompaña
a la derivada coincide y que el lado derecho del igual también, lo único que falta es
forzar que ρ ′ = ρa(t). R
6. Resolver la EDO ρ ′ = ρa(t). Como es lineal se tendrá que ρ = e a(t) dt recordemos
que ρ > 0 entonces si al obtener la función ρ tenemos algo negativo debemos
modificar el resultado de forma oportuna para arreglarlo.
7. RSustituimos en 4. eR integramos:
x
)′ = xx0 ρb(t) dt
R
a(t) dt
x0 (t · e
8. Simplificamos la solución.
 6. EDOs de Bernoulli

Diremos que una EDO es de Bernoulli si es de la forma

x′ = a1 (t)x + an (t)xn (n > 1, n ∈ Z)

donde a1 y a2 son coeficientes continuos en (a, b) ⊆ R.
Para resolver este tipo de ecuaciones debemos de garantizar que x ̸= 0, ya que el
método empleado consiste en dividir toda la EDO por la mayor potencia de x (despejamos
an (t)):

x′
n
= a1 (t)x1−n + an (t)
x
Como vemos aplicando el cambio de variable u = x1−n se logra una EDO lineal, ya
que derivando u se tiene:

x′
u′ = (1 − n)
xn
Entonces

u′
= a1 (t)u + an (t)
1−n
Despejando la derivada tenemos (a = a1 y b = a2 )

u′ = (1 − n)a(t)u + (1 − n)b(t)
tenemos pues una EDO lineal completa.
 20 Capítulo 6. EDOs de Bernoulli

Teorema 6.0.1 — Teorema de existencia-unicidad (ber). Dados a(t), b(t) ∈ C ((a, b)),
con x ̸= 0 entonces se tendrá existencia y unicidad de soluciones en (a, b) × (−∞, 0) o
(a, b) × (0, ∞).

Notemos que si n = 1 se tiene directamente una EDO de variables separadas. Por otro
lado si el dato para un problema de Cauchy es x(t0 ) = 0 entonces la solución de la EDO es
la trivial.

6.1 Bernoulli generalizado

Diremos que una EDO es de Bernoulli (generalizada) si tiene la forma

x′ = a1 (t)x + an (t)xr (r ∈ R \ {1})

para resolver estas EDOs procedemos de forma idéntica al caso anterior: despejando
an (t).

RESUMEN: para resolver EDOs de tipo Bernoulli debemos hacer lo que sigue
1. Despejar la x′ en caso de que sea necesario.
2. Dividir por xr .
3. Tomar u igual a lo que acompaña a la a1 (t).
4. Resolver la EDO lineal obtenida.
 7. EDOs de Ricatti

Diremos que una EDO es de Ricatti si es de la forma

x′ = a0 (t) + a1 (t)x + a2 (t)x2

con ai (t) ∈ C ((a, b)). Además, este tipo de ecuaciones solo se podrá resolver si previa-
mente conocemos una solución particular, x p , pues en tal caso haremos la sustitución

x = u + xp
Así pues llegaremos a una EDO de Bernoulli, que ya sabemos resolver. Tras realizar el
cambio de variables y notando que x′ = u′ + x′p sustituyendo y organizando de forma ade-
cuada los términos (en adelante si se tiene ai entiéndase que es un coeficiente dependiente
de t, las constantes serán indicadas de forma adecuada):

u′ + x′p = a0 + a1 (u + x p ) + a2 (u + x p )2
= a0 + a1 x p + a2 x2p + a1 u + 2x p u + a2 u2

Ahora, nos faltaría desarrollar el lado izquierdo de la igualdad, en este momento usamos
que x p es solución particular y que, por tanto

x′p = a0 + a1 x p + a2 x2p
Por tanto simplificando nos queda una EDO de Bernoulli:

u′ = a1 u + 2x p u + a2 u2 = (a1 + 2x p )u + a2 u2

Nota 7.1. En la década de 1760, Euler propuso el cambio de variable

1
x = xp +
u
22 Capítulo 7. EDOs de Ricatti

con el que la EDO que se obtiene tras las simplificaciones y cálculos pertinentes es
lineal.

RESUMEN: para resolver EDOs de Ricatti seguir los siguientes pasos

1. Despejar x′ .
2. Obtener una solución particular x p (se suele proporcionar la forma que tiene), para
ello se impone la EDO.
3. Hacer el cambio x = u + x p y también derivarlo.
4. Sustituir en la EDO lo obtenido en el paso anterior.
5. Tras operar debe llegarse a una Bernoulli.
 8. EDOs exactas

Diremos que una EDO de la forma siguiente

M(t, x)
x′ = (1)
N(t, x)
y verificando que M, N ∈ C (Ω) con Ω ⊆ R2 dominio y N ̸= 0 es exacta si, y solo si
para cualquier punto (t, x) de Ω existe ω ⊆ Ω entorno de (t, x) y u ∈ C 1 (ω) tal que

ut (t, x) = −M(t, x)
ux (t, x) = N(t, x)

Donde u es denominada función potencial.

Nota 8.1. Si M y N son más regulares se dará que Mx + Nt = 0 o, en otras palabras: el

laplaciano de u se anula (∆u = 0) es decir u es armónica.

Teorema 8.0.1 — Teorema (condición de exactitud). Sean Ω ⊆ R2 dominio y

M, N ∈ C 1 (Ω) (N ̸= 0) tales que Mx + Nt = 0 entonces la EDO (1) es exacta.

Teorema 8.0.2 — Teorema de existencia-unicidad (exactas). Dados la EDO exacta

(1) en el dominio Ω y la condición de Cauchy x(t0 ) = x0 admite una única solución
definida localmente mediante u(t, x) = u(t0 , x0 ).

Demostración. Como la EDO es exacta se dará que ut = −M y que ux = N. Ahora

procederemos a buscar u(t, x) = u(t0 , x0 ) para un cierto entorno del punto.
Queremos ver que a la función auxiliar
24 Capítulo 8. EDOs exactas

H(t, x) = u(t, x) − u(t0 , x0 )

le podemos aplicar el Teorema de la función implícita, para concluir que x es una
función de t:
Vemos que H(t, x) se anula, al menos en (t0 , x0 ), veamos si el jacobiano (gradiente en
este caso) de la función es no nulo en el punto inicial:

Ht = ut (t, x) = −M(t, x)Hx = ux (t, x) = N(t, x)

Por hipótesis N(t, x) ̸= 0 por lo que en el valor inicial el jacobiano es no nulo y, por
tanto se puede dar x en función de t además se puede dar u(t, x(t)) = u(t0 , x0 ) localmente,
cerca del valor inicial.
■

Nota 8.2. Si expresamos la EDO como

−M(t, x)dt + N(t, x)dx = 0

esta será exacta si −Mx (t, x) = Nt (t, x)

 9. EDOs reducibles a exactas

Dada la EDO

P(t, x)dt + Q(t, x)dx = 0

si resulta que Px (t, x) ̸= Qt (t, x), es decir que la EDO no es exacta, quizás podamos
encontrar una función ρ(u) no nula (donde u es una combinación de las variables t y x),
tal que al multiplicar la expresión anterior por él nos quede una EDO exacta, es decir:

ρ(u)P(t, x) dt + ρ(u)Q(t, x) dx = 0
| {z } | {z }
P̃(t,x) Q̃(t,x)
reescribiendo la expresión

P̃(t, x)dt + Q̃(t, x)dx = 0

donde ahora P̃x (t, x) = Q̃t (t, x), esto es, mediando el factor integrante ρ hemos obtenido
una EDO exacta de una que no lo era.

¿Por qué no buscamos factores integrantes que dependan de forma independiente de t

y de x? porque si así lo hiciésemos llegaríamos al imponer la condición de exactitud a una
EDP (ecuación en derivadas parciales) complicada de resolver.

Por lo general se suele probar a buscar factores integrantes que solo dependan de una
variables (o de las dos pero en forma de combinación de estas), por ejemplo busquemos
ρ(u) con u = t. De esta forma se tendrá que

P̃(t, x) = ρ(t)P(t, x)
Q̃(t, x) = ρ(t)Q(t, x)
26 Capítulo 9. EDOs reducibles a exactas

Entonces imponiendo la condición de exactitud:

P̃x (t, x) = Q̃t (t, x) ⇐⇒ Px (t, x)ρ(t) = ρ ′ (t)Q(t, x) + ρ(t)Qt (t, x)

Hemos obtenido una EDO, que reescribiéndola:

ρ ′ (t) Px (t, x) − Qt (t, x)

=
ρ(t) Q(t, x)
Así pues, existirá factor integrante ρ si, y solo si, el lado derecho de la igualdad solo
dependa de t (¿por qué?). Análogamente, si buscáramos factor integrante ρ(u) con u = x
llegaríamos a

ρ ′ (x) Qt (t, x) − Px (t, x)

=
ρ(x) P(t, x)
con la salvedad de que en este caso debemos verificar que, al menos en un entorno
del valor inicial P(t, x) es no nula. Nuevamente, existirá factor integrante ρ si, y solo si,
el lado derecho de la igualdad solo depende de x. En caso de que no podamos encontrar
factor integrante con u = t o u = x se tendrá que hacer u = f (t, x) (función escalar).

Nota 9.1. En caso de haber elegido ρ(u) con u = f (t, x) a la hora de resolver la EDO
resultante para obtener ρ debemos comprobar que el lado de la derecha solo dependa de
u, si no, no tendremos factor integrante y tendríamos que elegir u de otra forma.
III
Tema 3: Existencia y unicidad
de EDOs

10 Introducción . . . . . . . . . . . . . . . . . . . . . . 28

11 Iterantes de Picard . . . . . . . . . . . . . . . 30
 10. Introducción

En este capítulo trabajaremos en R2 , pero todo esto es extrapolable a n variables con

la sutileza de que en lugar de usar el valor absoluto tomaremos la norma que más nos
convenga (en un espacio de Banach todas las normas son equivalentes).
Dado un problema de Cauchy

x′ = f (t, x)
x(t0 ) = x0

estudiaremos diferentes causísticas, como por ejemplo que nuestros datos iniciales es-
tén perturbados (continuidad respecto a los parámetros) o bien f depende de un parámetro
λ ∈ R, esto será estudiado al final de esta sección.

Para tener existencia-unicidad la f no necesariamente debe de ser buena, sino que la f

y el pertinente dominio deben ser complementarios de forma que aunque la f no sea muy
buena si el dominio lo permite se podrá dar la exitencia y la unicidad.

Necesitaremos conceptos como condición de Lipschitz o iterantes de Picard.

Definición 10.0.1 Sea I ⊆ R intervalo, una función f : I → R se dice que el lips-

chitziana en I o que cumple una condición de Lipschitz en I cuando existe L > 0 tal
que
| f (x) − f (y)| ≤ L|x − y|
Donde la constante L es la constante de Lipschitz para f en I. Lo denotaremos como
f ∈ Lip(I).

Desarrollando la definición de Lipschitz se tiene que

f (x) − L|x − y| ≤ f (y) ≤ f (x) + L|x − y|

 29

Suponiendo que x − y < 0 entonces

f (x) − L(y − x) ≤ f (y) ≤ f (x) + L(y − x)

Por tanto la función queda acotada por rectas (considerando x − y > 0) se verá que f
está en una región con forma de pajarita por lo que la función no puede tener asíntotas
(por ejemplo).

Sabemos que
1. Si f ∈ Lip(I) entonces f ∈ UC(I) (uniformemente continua en I)
2. Si f es derivable con derivada acotada, entonces es Lipschitz.
3. Si L < 1 entonces la función es contractiva y, por tanto continua.

Nota 10.1. No toda función Lipschitz es derivable, basta tomar |x|.

¿Cómo afecta la condición de Lipschitz a un problema de Cauchy? Solo nos interesa

en la segunda variable (x) entonces se tiene

Definición 10.0.2 Sea Ω ⊆ R2 dominio y f : Ω → R. Se dice que f es lipschitziana en

Ω respecto de la segunda variable x si existe L > 0 tal que

| f (t, x1 ) − f (t, x2 )| ≤ L|x1 − x2 |, (x1 , x2 ) ∈ Ω

Esto lo denotaremos como f ∈ Lip(x, Ω).

Teorema 10.0.1 Sean Ω ⊆ R2 dominio convexo y f una función que admite parcial
respecto de x acotada, entonces f es Lipschitz en la segunda variable en Ω.

Demostración. Definimos
g(x) = f (t, x)
por tanto
g′ = fx (t, x)
y por hipótesis
|g′ (x)| = | fx (t, x)| ≤ L
Por el teorema del valor medio aplicado en g y teniendo que Ω es convexo (pensamos en
un segmento):
g(x) − g(y) = g′ (θ )(x − y)
entonces
|g(x) − g(y)| ≤ L|x − y|
para cuales quiera (t, x), (t, y) ∈ Ω ■
El recíproco teniendo que la función es derivable también es cierto.
 11. Iterantes de Picard

En capítulos anteriores vimos un resultado que podemos resumir en

Teorema 11.0.1 Sea Ω ⊆ R2 dominio y f ∈ C (Ω). Sea (t0 , x0 ) ∈ Ω e I un intervalo tal

que t0 ∈ I. Una función x : I → R con su gráfica en Ω ((t, x(t)) ∈ Ω), es solución del
problema de Cauchy

x′ = f (t, x)
x(t0 ) = x0

si, y solo si x es una función continua en I verificando que

Z t
x(t) = x0 + f (s, x(s)) ds
t0

Sabemos que el operador

T: C (I) −→ C (I)
x −→ Tx
donde
Tx : I −→ R R
t −→ x0 + tt0 f (s, x(s)) ds
para esto es necesario que la gráfica de la función pertenezca al dominio de la función.

Dado el problema

x′ = x
x(0) = 1
 31

buscaremos un método iterativo que nos aproxime la solución. Para empezar necesita-
mos una función de partida, por ejemplo x0 = 1, ahora x1 = T x0 por tanto
Z t Z t
x1 = x0 + f (s, x0 ) ds = 1 ds = 1 + t
t0 =0 t0

por su lado
t2
Z t
x2 = x0 + f (s, x1 (s)) ds = 1 + t +
0 2
siguiendo el proceso

t2 t3
Z t
x3 = x0 + f (s, x2 (s)) ds = 1 + t + +
0 2 3!
Vemos que en estas iteraciones tenemos vemos que tenemos las sumas parciales del
desarrollo de la exponencial, por inducción debemos ver que para xn se tiene el desarrollo
de la exponencial, que converge absoluta y uniformemente, por tanto la solución del
problema es la exponencial. Pues si los interantes convergen deben hacerlo al punto fi-
jo que es la exponencial, además los iterantes están bien definidos por el dominio de f (t, x).

Veamos el iterante n + 1 con

n
tk
xn (t) = ∑ k!
k=0
entonces Z t Z t
xn+1 = x0 + f (s, xn (s)) ds = 1 + xn (s) ds
0 0
por tanto
Z t n k Z t
s 1
xn+1 = 1 + = 1+∑
∑ k! ds C.U.A. sk ds
0 k=0 k k! 0

Así el resultado es
n n+1 k
t k+1 t
xn+1 (t) = 1 + ∑ =∑
k=0 (k + 1)! k=0 k!

entonces
n
tk
xn (t) = ∑
k k!

Nota 11.1. En general los iterantes no tienen interés práctico, solo teórico.
Anexo 1
 12. Ejercicios examen temas 1 y 2

En esta sección se dará una recopilación de ejercicios como entrenamiento para el

examen de los temas 1 y 2.
Ejercicio 12.1 Identificar y resolver la siguiente ecuación diferencial:

y2 + (xy′ − x3 )y = 0

SOLUCIÓN:

Despejando la derivada:
−y
y′ =
+ x2
x
Vemos claramente que tenemos una EDO lineal completa, con factor integrante
R 1
ρ(x) = e x dx =x
Por lo que
(yρ)′ = x2 ρ
entonces
x3 x4
yx = +C ⇒ y = + xC
3 3

Ejercicio 12.2 Se considera el siguiente problema de Cauchy

(
y′ = y2 − ax2 , a ∈ R
y(x0 ) = y0
34 Capítulo 12. Ejercicios examen temas 1 y 2

1. Determina el dominio en el que el PVI anterior tiene garantizada la existencia y

la unicidad.
2. Determina la regularidad de la solución en función del dato inicial (x0 , y0 ).
3. Describir el campo de direcciones asociado a la EDO anterior y hacer un esbozo
de las soluciones para los distintos valores de a.
HINT: Haz tres estudios distintos, para a < 0, a = 0 y a > 0.
4. Obtén la solución del problema de Cauchy para a = 0. Contrasta el resultado con
el análisis anterior. Por otro lado, para a ̸= 0 identifica el tipo de EDO que se
tiene.
■

SOLUCIÓN:

1. Para ver en qué dominio se tiene existencia y unicidad debemos ver en qué dominio
es
f (x, y) = y2 − ax2
derivable.

Sabemos que por ser un polinomio en dos variables f ∈ C ∞ (R2 ) por tanto tenemos
existencia y unicidad en todo el plano.
2. El dato, en este caso no afectará a la regularidad, ya que nuestra EDO viene dada
por un polinomio.
3. Comencemos el estudio:
a) Para a < 0:
Entonces nuestra EDO queda como y′ = y2 + ax2 ahora con a > 0. Vemos que
Im( f ) = R+ por tanto solo podemos probar con pendientes positivas para las
isoclinas:
Φc = {(x, y) ∈ R2 : y2 + ax2 = c}
Para que el estudio sea más sencillo fijemos a = 1 (después discutiremos que
pasa si este valor crece).
Para c = 0:
Φ0 = {(x, y) ∈ R2 : y2 + x2 = 0} ≡ y2 + x2 = 0
Como vemos, el lugar geométrico donde la pendiente de f es nula es el punto
(0, 0).
Para c = 1:
Φ1 = {(x, y) ∈ R2 : y2 + x2 = 1} ≡ y2 + x2 = 1
El lugar geométrico de los puntos donde la pendiente es 1 es la circunferencia
de radio 1 y centro cero.
Para c = 4:
Φ4 = {(x, y) ∈ R2 : y2 + x2 = 1} ≡ y2 + x2 = 4
El lugar geométrico de los puntos donde la pendiente es 4 es la circunferencia
de centro cero y radio 2.
 35

Figura 12.1: En amarillo los lugares geométricos y en cian el esbozo de la solución.

Como vemos para a = 1 tenemos circunferencias. Para a > 1 se tendrán elipses

cuyo eje mayor estará en el eje Y y para 0 < a < 1 se tendrá una elipse con el
eje mayor en el eje X.

b) Para a = 0:
Entonces nuestra EDO queda como y′ = y2 . La imagen de f es [0, ∞). Veamos
las isoclinas.
Para c = 0:

Φ0 = {(x, y) ∈ R2 : y2 = 0} ≡ y2 = 0 ≡ y = 0

Para c = 1

Φ1 = {(x, y) ∈ R2 : y2 = 1} ≡ y2 = 1 ≡ y = 1 ∨ y = −1

Por tanto para a = 0 tenemos rectas como isoclinas.

Figura 12.2: En amarillo los lugares geométricos y en cian el esbozo de la solución.

c) Para a > 0: Entonces nuestra EDO queda como y′ = y2 − ax2 , además vemos
que y = 0 es solución. La imagen de f vuelve a ser R+ . Veamos las isoclinas
para a = 1 (después discutimos el resto de valores): Para c = 0:

Φ0 = {(x, y) ∈ R2 : y2 − x2 = 0} ≡ y2 = x2 ≡ y = −x ∨ y = x
36 Capítulo 12. Ejercicios examen temas 1 y 2

Para c = 1:

Φ1 = {(x, y) ∈ R2 : y2 − x2 = 1} ≡ y2 − x2 = 1

Para c = 2:

Φ2 = {(x, y) ∈ R2 : y2 − x2 = 2} ≡ y2 − x2 = 2

Tenemos hipérbolas centradas en el origen.

Figura 12.3: En amarillo los lugares geométricos y en cian el esbozo de la solución.

Para valores de 0 < a < 1 se tendrá una mayor excentricidad y para valores de
a > 1 un excentricidad cada vez menor.
4. Para a ̸= 0 se tiene una EDO de Ricatti. Ahora, hallemos la solución para a = 0, en
este caso la EDO era y′ = y2 por ejemplo, mediante variables separadas se obtiene
que
1
y=
c−x
con c ∈ R (además, recordemos que y = 0 es solución) vemos que el estudio previo
con las isoclinas cuadra con el resultado pues estas soluciones variando c quedan
como

Figura 12.4: Cada color es la solución con un c distinto.

 37

Ejercicio 12.3 Es posible resolver la ecuación

y′ = y
p
|x|

para el dato y(0) = 0. ¿Es la solución única? ■

SOLUCIÓN:
p
Sí, pues la función que define la EDO f (x, y) = y |x| que es continua en todo el plano,
por tanto existe solución. En particular para el dato en el (0, 0) vemos que y ≡ 0 es solución.

La unicidad vendrá dada por la diferenciabilidad de f respecto de y:

p
fy = |x|

y vemos que esta parcial es continua en todo el plano. Por tanto hay unicidad de
solución en todo R2 .
Ejercicio 12.4 Dibuja las isoclinas y esboza las soluciones de

y′ =
p
|y|

Luego calcula la solución de la EDO. ■

SOLUCIÓN:

Como tenemos una EDO autónoma los ceros de la f son soluciones, por tanto y ≡ 0 es
solución. Diferenciaremos ahora para y > 0 e y < 0.

Si y > 0:
√
Nuestra EDO queda como y′ = y. Vemos que la imagen de f es R+ , veamos las isoclinas,
para c = 0:
√
Φ0 = {(x, y) ∈ R × R+ : y = 0} ≡ y = 0
Para c = 1:
√
Φ1 = {(x, y) ∈ R × R+ : y = 1} ≡ y = 1
Para c = 4:
√
Φ1 = {(x, y) ∈ R × R+ : y = 4} ≡ y = 2
Por lo que las isoclinas en el semiplano superior son rectas. Para el semiplano inferior
√
y < 0 tenemos la EDO y′ = −y por lo que las isoclinas son Para c = 1:
√
Φ1 = {(x, y) ∈ R × R+ : −y = 1} ≡ y = −1

Para c = 4:
√
Φ1 = {(x, y) ∈ R × R+ : −y = 4} ≡ y = −2
Tenemos rectas en el semiplano inferior. Parecido al caso b) del ejercicio anterior.
38 Capítulo 12. Ejercicios examen temas 1 y 2

Para el semiplano superior resolvemos la EDO

√
y′ = y
y para el semiplano inferior resolvemos la EDO
√
y′ = −y
Ahora, por variables separadas se tendrá para el primer caso
1
y(x) = (2c1 x + c21 + x2 ), c1 ∈ R
4
Para el otro caso será
1
y(x) = − (2c1 x + c21 + x2 ), c1 ∈ R
4

Ejercicio 12.5 Inicialmente un tanque contiene 10L de agua. Se vierte en el tanque

salmuera a razón de 1 litro por minuto. El agua se mezcla bien y se drena a 1 litro por
minuto. En 20 minutos hay 15 gramos de sal en el tanque. ¿Cuál es la concentración de
sal en la salmuera entrante? ■

SOLUCIÓN:

Tenemos un problema de mezclas, cuya ecuación es

S′ (t) = ci − co
S(0) = 0
S(20) = 15
Donde S es la cantidad de sal en el tanque en el tiempo t y ci , co son las concentraciones
de sal que entran y salen del tanque respectivamente. En este caso nos piden ci . La
concentración de salida será
S(t)
co =
V
en este caso el volumen es constante 10L por tanto tenemos la EDO lineal
S
S′ = ci −
10
resolviendo mediante factor integrante llegamos a que
S = ci · 10 + k · e−t/10
por tanto imponiendo que S(0) = 0 = ci · 10 + k ⇒ k = −10ci entonces
S = S = ci · 10 − 10ci · e−t/10 = 10ci (1 − e−t/10 )
como sabemos que
15
S(20) = 15 = 10ci (1 − e−2 ) ⇒ ci = ≈ 1,73 (g/L)
10(1 − e−2 )
 39

Ejercicio 12.6 Bob se hace un café todas las mañanas, pero solo se lo bebe cuando
está a 60oC. La primera vez que Bob toma la temperatura del café esta es de 89oC. Un
minuto más tarde es de 85oC. Sabiendo que la habitación en la que está el café tiene
una temperatura de 22oC ¿Cuándo debería empezarse Bob su café? ■

SOLUCIÓN:

Tenemos que resolver la EDO

T ′ = −k(T − Ta )
T (0) = 89
T (1) = 85
T (t0 ) = 60

Por tanto resolviendo la EDO para Ta = 22 y con k > 0 (notemos que T − 22 > 0) se
da que T = e−kt ·C + 22
Imponiendo T (0) = 89:
T (0) = 89 = C + 22 ⇒ C = 67
Por lo que T = e−kt · 67 + 22 ahora, como T (1) = 85:

63
T (1) = 85 = e−k · 67 + 22 ⇒ e−k =
67
entonces la solución es  t
63
T= · 67 + 22
67
Ahora imponemos que T (t0 ) = 60 por lo que
 t0
63
T (t0 ) = 60 = · 67 + 22
67

Despejando llegamos a que

ln (38/67)
t0 = >0
ln (63/67)
así que Bob deberá tomarse el café una vez transcurridos t0 ≈ 9,21 minutos desde que lo
hizo.

Ejercicio 12.7 Comprobar que la siguiente EDO no es exacta

(y + x3 y + 2x2 )dx + (x + 4xy4 + 8y3 )dy = 0

Pero sabemos que existe un factor integrante dependiente de xy, una vez hallado resuelve
la EDO. ■

SOLUCIÓN: me da pereza. Pero es fácil.

40 Capítulo 12. Ejercicios examen temas 1 y 2

Ejercicio 12.8 Sabemos que si dada la EDO lineal

x′ = a(t)x + b(t)

se da que a(t) es una función periódica de periodo T y b(t) = 0 entonces decir que la
solución de la EDO es periódica equivale a ver que
Z T
a(t) dt = 0
0

Ahora, si tenemos que b(t) ̸= 0 es periódica de periodo t y

Z T
a(t) dt ̸= 0
0

¿Podemos concluir que x(t) es solución periódica? ¿Para este caso se verifica que la
integral de 0 a T de b(t) es nula? Justifica tus respuestas. ■

SOLUCIÓN:

Similar a lo que se hizo en clase, solo

Rt
hay que resolver por factor integrante y percatarse
t0 a(s) ds
de que ρ(t) = ρ(t + T ) con ρ(t) = e .
Ejercicio 12.9 Comprueba que

e2x (x − 1) + x + 1
y(x) =
e2x + 1
es solución del PVI

y′ = (x − y)2
y(0) = 0

SOLUCIÓN:

Lo primero que miramos es que nuestra función sea derivable en un entorno del (0, 0)
(el dato), evidentemente lo es por ser un cociente de funciones donde el denominador no
se anula y el numerador es infinitamente diferenciable.

Ahora veamos que cumple el dato

−1 + 1
y(0) = =0
1

por lo que se cumple el dato.

 41

Nos falta por ver que satisface la EDO:

(2e2x x + e2x − 2e2x + 1)(e2x + 1) − (e2x (x − 1) + x + 1)(2e2x )

y′ =
(e2x + 1)2

Ahora operando
(e2x − 1)2
y′ =
(e2x + 1)2
Sustituyendo en la EDO:
2
(e2x − 1)2 (?) e2x (x − 1) + x + 1 (e2x − 1)2

= x − =
(e2x + 1)2 e2x + 1 (e2x + 1)2
Por lo que verifica la EDO. Por tanto es solución.