Capı́tulo 2

Ecuaciones Diferenciales
ordinarias de primer orden

En este capı́tulo nos centraremos en el estudio de las ecuaciones diferenciales ordinarias

de primer orden. En particular comentaremos la interpretación geométrica de las soluciones
de estas ecuaciones. Se clasificarán las ecuaciones según su tipo para aplicar el método
más adecuado para resolverlas y se verán métodos de resolución para algunas de estas
ecuaciones.

2.1. Ejemplo inicial: Modelo de un solenoide de arranque

de un motor
Se considera un solenoide de arranque de un motor [2], esto es, un dispositivo electro-

Figura 2.1: Solenoide de arranque. Fuente [1]

magnético que consta de los siguientes elementos:

Una bobina, esto es un alambre conductor enrollado que al recibir energı́a genera un
campo magnético.

Un émbolo o núcleo móvil que es una pieza de material ferromagnético que se mueve
en el interior de la bobina.

El funcionamiento consiste en que el émbolo acciona un mecanismo para engranar el

motor de arranque con el volante del motor.
Para modelar este dispositivo es necesario tener en cuenta que la resistencia R es
la del alambre de la bobina mientras que la inductancia L es consecuencia del efecto
electromagnético. Por tanto este dispositivo se puede representar por medio de un circuito
RL como es de la siguiente figura:

17
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

i(t)
R

vs

Figura 2.2: Circuito RL: solenoide

Usando la segunda ley de Kirchhof o ley del voltaje, que asegura que la suma de los
voltajes alrededor de un circuito cerrado debe ser cero debido a la conservación de la
energı́a, se obtiene la siguiente ecuación diferencial1 que permite modelar la intensidad de
corriente i(t)
di
L + Ri = Vs (2.1)
dt
Asumiendo que vs = V es constante y que i(0) = 0 se obtiene un problema de valor
inicial cuya solución, como veremos en este capı́tulo es

V V
i(t) = i(0)e−Rt/L + (1 − e−R/L ) = (1 − e−Rt/L ).
R R
A lo largo de este tema se considera este y otros tipos de EDOs de primer orden y se
analizarán los métodos de resolución para cada una de ellas.

2.2. Primeras definiciones

En este primer apartado se reescriben, en caso de primer orden, las definiciones que se
han estudiado en el capı́tulo anterior.

Definición 2.1 (EDO de primer orden) Una ecuación diferencial de primer orden
es una ecuación que relaciona una variable independiente x, una variable dependiente
o función incógnita y(x) y la primera derivada y ′ (x) de la variable dependiente.

Las ecuaciones diferenciales ordinarias de primer orden se pueden expresar de diferentes

formas. Sea y(x) una función derivable:

Una ecuación diferencial ordinaria de orden 1 en forma general o implı́cita es

una expresión de la forma
F (x, y, y ′ ) = 0 .

Una ecuación diferencial ordinaria de orden 1 en forma normal o explı́cita es una

expresión de la forma
y ′ = f (x, y) .
1
En este caso es una EDO lineal completa de orden 1.

18
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Una ecuación diferencial ordinaria de orden 1 en forma diferencial es una expresión

de la forma
M (x, y)dx + N (x, y)dy = 0 .

Ejemplo 2.1 Escribir la ecuación diferencial de orden 1 (x−y)y ′ = x2 +y 2 en sus formas

general, normal y diferencial.
La forma normal o explı́cita consiste en despejar y ′ , ası́ la forma normal será

x2 + y 2
y′ = .
x−y
La forma general o implı́cita admite varias expresiones

x2 + y 2
y′ − =0 o bien (x − y)y ′ − x2 − y 2 = 0 .
x−y
Para escribir la forma diferencial es útil considerar que y ′ = dy/dx. Ası́, la ecuación
en forma normal se puede escribir

dy x2 + y 2
=
dx x−y
A partir de esta expresión es sencillo obtener la expresión diferencial de esta ecuación

(x2 + y 2 )dx − (x − y)dy = 0 .

Definición 2.2 (Solución de una EDO de orden 1) Una función ϕ: I ⊂ R → R

es una solución de la ecuación diferencial ordinaria de orden 1

F (x, y, y ′ ) = 0,

si se cumple

ϕ ∈ C 1 (I) y,

ϕ verifica la ecuación diferencial

F (x, ϕ(x), ϕ′ (x)) = 0,

en todo punto x ∈ I.

El intervalo I se denomina intervalo de definición o de existencia de ϕ.

Una EDO de orden 1 tiene, normalmente, una familia de soluciones dependien-

te de un parámetro C. Para cada valor particular de este parámetro se obtiene una
solución particular de dicha ecuación.
Si todas las soluciones de una EDO de orden 1 se pueden escribir como soluciones
particulares a partir de la misma familia de soluciones G(x, y, C) = 0, entonces esta
familia se denomina solución general de la ecuación de primer orden.
Si una ecuación diferencial de primer orden admite una familia de soluciones depen-
diente de un parámetro C, y existe una solución que no se puede escribir como solución
particular de esa familia para ningún valor de C, entonces esa solución es una solución
singular y esta EDO no admite solución general.

19
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Ejemplo 2.2 La familia y = Cex , siendo C una constante arbitraria, es una familia
uniparamétrica de soluciones de la ecuación diferencial y ′ = y.
En efecto, derivando y(x) = Cex se obtiene y ′ (x) = Cex , de donde se comprueba
trivialmente que y ′ = y.
La curva y(x) = 3ex es una solución particular de esta familia, en particular la solución
que pasa por el punto (0, 3).

Ejemplo 2.3 Verificar que y = 0 e y = x4 /16 son soluciones del PVI

 ′ √
y = x y,
y(0) = 0 .
En efecto, si y(x) es la función constante igual a cero, trivialmente se cumple la ecua-
ción diferencial y la condición inicial, por tanto y(x) = 0 es una solución del PVI.
Derivando la función y(x) = x4 /16 se obtiene y ′ (x) = x3 /4. Por otra parte
√ p
x y = x x4 /16 = x3 /4 ,

es decir, se verifica la ecuación diferencial y además y(0) = 0, por tanto es solución del
PVI.

En el ejemplo anterior se pone de manifiesto que un PVI puede tener más de una
solución.

2.3. Interpretación geométrica de una solución

Consideremos ahora una EDO de primer orden F (x, y, y‘) = 0 si está escrita en forma
implı́cita e y ′ = f (x, y) si se escribe en forma normal. Una solución de esta ecuación será
una función real de variable real y(x) que se puede representar en el plano OXY.

Definición 2.3 (Curva integral) Se llama curva integral o curva solución de una
ecuación diferencial ordinaria de orden 1 a la gráfica de una solución ϕ de dicha ecua-
ción.

Ejemplo 2.4 La siguiente figura representa algunas curvas integrales o soluciones de la

ecuación y ′ = x.

Sea y(x) solución de y ′ = f (x, y), la pendiente de la recta tangente2 a la curva solución
en el punto (x0 , y0 ) está dada por el valor f (x0 , y0 ). Teniendo en cuenta este hecho, es
posible representar una aproximación de las soluciones mediante los denominados campos
de direcciones.

Definición 2.4 (Campo de direcciones) Una ecuación diferencial y ′ = f (x, y) de-

fine un campo de direcciones en el plano asociando a cada punto (x0 , y0 ) un vector cuya
pendiente es igual a la de la solución de la EDO que pasa por ese punto, por ejemplo
el vector (1, f (x0 , y0 )).

2
La pendiente de la recta tangente a la gráfica de la función y(x) en un punto (x0 , y(x0 )) viene dada
por y ′ (x0 ).

20
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Figura 2.3: Soluciones de y ′ = x. Elaboración propia

Ejemplo 2.5 Campo de direcciones de la ecuación y ′ = x y algunas soluciones represen-

tadas.

Figura 2.4: Campo de direcciones. Elaboración propia

Siguiendo los campos de direcciones se pueden dibujar curvas solución de una ecuación
diferencial, sin necesidad de conocer la expresión de la solución.

Definición 2.5 (Isoclinas) Se llama isoclina al lugar geométrico de los puntos del
plano en donde las rectas tangentes a las soluciones de la ecuación diferencial y ′ =
f (x, y) tienen la mima pendiente, esto es , la isoclina de nivel k es el lugar geométrico
definido por la ecuación f (x, y) = k.
La representación gráfica de varias isoclinas se denomina gráfica de contorno. Estas
gráficas permiten observar los intervalos de crecimiento y decrecimiento de las solucio-
nes de la ecuación diferencial

21
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Ejemplo 2.6 Ejemplo de campo de direcciones, isoclinas y soluciones asociadas a la ecua-

ción y ′ = x.
Para calcular las isoclinas de nivel k se plantea la ecuación

f (x, y) = k,

En este ejemplo particular se obtiene x = k, es decir, las isoclinas son rectas paralelas
al eje de ordenadas. Obsérvese también sobre este ejemplo que sobre una isoclina el campo
de direcciones está formado por vectores paralelos entre sı́.

Figura 2.5: Isoclinas. Elaboración propia.

Las isoclinas de pendiente nula separan las isoclinas f (x, y) = k, k > 0, donde las
soluciones de la ecuación diferencial tienen pendiente positiva y por tanto son funciones
crecientes de las isoclinas f (x, y) = k, k < 0 en las que las soluciones tienen pendiente
negativa y por tanto son funciones decrecientes. Como consecuencia de este análisis se
puede concluir que el principal interés en el cálculo de las isoclinas es el estudio de las
denominadas isoclinas nulas.
El estudio de los campos de direcciones e isoclinas nos da información sobre las solu-
ciones de una EDO de primer orden, sin necesidad de calcular la solución.

Ejemplo 2.7 Consideremos por ejemplo el siguiente campo de direcciones y algunas iso-
clinas asociadas a la EDO y ′ = sin y.
Las isoclinas nulas, que son las de mayor interés, son las curvas de nivel sin y = 0, es
decir, las rectas de la forma y = nπ, n ∈ Z.
Vamos a fijarnos en lo que ocurre en el intervalo −π < y < π. Este patrón se repite
en los demás intervalos de amplitud 2π.

Las soluciones que pasan por puntos (x, y), con −π < y < 0 son funciones decre-
cientes, además se cumple

lı́m y(x) = 0 y lı́m y(x) = −π

x→−∞ x→+∞

es decir, las rectas y = −π e y = 0 son ası́ntotas de estas soluciones.

22
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Las soluciones que pasan por puntos (x, y), con 0 < y < π son funciones crecientes,
además se cumple
lı́m y(x) = 0 y lı́m y(x) = π
x→−∞ x→+∞

es decir, las rectas y = π e y = 0 son ası́ntotas de estas soluciones.

Figura 2.6: Isoclinas. Elaboración propia

2.4. Métodos analı́ticos de integración

En este apartado veremos los métodos más sencillos para obtener la solución de una
EDO de primer orden, tanto de forma analı́tica como de forma numérica. Recordemos que
el proceso para hallar la solución de una ecuación diferencial se denomina integración.
Existen diferentes métodos de integración según el tipo de ecuación diferencial. Es
importante que una misma ecuación diferencial puede corresponderse con distintos tipos,
en cuyo caso se podrá resolver por más de un método.

2.4.1. Ecuaciones de variables separables

Definición 2.6 (Variables separables) Una ecuación diferencial ordinaria de pri-

mer orden se dice ecuación de variables separables si se puede escribir de la forma
dy
= f (x)g(y) .
dx

23
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Ejemplo 2.8
x−5 1
y′ = 2
es una ecuación en variables separables con f (x) = x − 5 y g(y) = 2 .
y y
−1
2ydx + 3xdy = 0 es una ecuación en variables separables siendo f (x) = y
3x
g(y) = 2y.

y ′ = cos 3x + 5 es una ecuación en variables separables con f (x) = cos 3x + 5 y

g(y) = 1.

y ′ = y + sin x no es una ecuación en variables separables.

Método resolución EDO de variables separables

Para resolver una ecuación en variables separables se procede del siguiente modo:

Escribir la ecuación en la forma

dy
= f (x)dx ,
g(y)

si g(y) ̸= 0 .

Integrar ambos miembros de la ecuación anterior:

Z Z
dy
= f (x)dx ,
g(y)

Si G(y) y F (x) son las integrales indefinidas de 1/g(y) y f (x), la expresión

G(y) = F (x) + C

proporciona una solución implı́cita de la ecuación diferencial de variables separa-

bles

Nota 2.1 (Pérdida de solución) Es importante tener en cuenta que si g(a) = 0, en-
tonces la función y(x) = a será una solución de la ecuación diferencial.
x−5
Ejemplo 2.9 Resolver la ecuación y ′ = .
y2
Esta ecuación es una ecuación en variables separables. Reescribimos en la forma

y 2 dy = (x − 5)dx ,

y se integra a ambos lados de la ecuación

y3 x2
Z Z
y 2 dy = (x − 5)dx ⇒ = − 5x + C1 ,
3 2
siendo C1 una constante arbitraria. Se puede despejar y si se quiere la expresión explı́cita
1/3
3x2

y= − 15x + C2
2

24
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

donde C2 = 3C1 . Habitualmente esta constante se representa simplemente por C por

comodidad. Ası́ la solución de esta EDO será
 2 1/3
3x
y= − 15x + C
2
Ejemplo 2.10 Resolver la ecuación y ′ = f (x).
Este tipo de de ecuación diferencial de primer orden es el tipo más sencillo que se
puede resolver. Se reescribe en la forma
dy = f (x)dx ,
y se integra a ambos lados de la ecuación
Z Z Z
dy = f (x)dx ⇒ y = f (x)dx + C,

siendo C una constante arbitraria.

Ejemplo 2.11 Resolver la EDO y ′ = y 2 − 4.
Esta es una EDO en variables separables con f (x) = 1 y g(y) = y 2 − 4. Se puede
reescribir de la siguiente forma
dy
= dx, y ̸= ±2 .
y2 − 4
Tal como hemos comentado en la nota 2.1, las funciones y(x) = −2 e y(x) = 2 serán
también soluciones de esta EDO.
Integramos a ambos lados de la igualdad (es necesario descomponer en fracciones sim-
1 3
ples 2 ) y se obtiene
y −4
y−2
ln = 4x + C, C constante arbitraria.
y+2
Despejamos y para obtener la expresión explı́cita de la solución y se obtiene
1 + C1 e4x
y=2 , C1 = ±eC ̸= 0 ,
1 − C1 e4x
que es una familia uniparamétrica de soluciones de la EDO y ′ = y 2 − 4 .
Recordemos que existen las soluciones y = 2 e y = −2.
1 + C1 e4x
La solución y = 2 se corresponde con la expresión y = 2 cuando C1 = 0, por
1 − C1 e4x
tanto podemos considerar la familia uniparamétrica de soluciones
1 + C1 e4x
y=2 , C1 constante arbitraria.
1 − C1 e4x
Sin embargo es imposible escribir la solución y = −2 como un elemento de la familia
anterior, por tanto y = −2 es una solución singular de la ecuación y ′ = y 2 − 4.
En la práctica no siempre se pueden resolver de forma analı́tica las ecuaciones de
variables separables ya que no siempre se pueden integrar las funciones que aparecen en
términos de funciones elementales.
3
La descomposición en fracción simples es
 
1 1 1 1
= − .
y2 − 4 4 y−2 y+2

25
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

2.4.2. Ecuaciones homogéneas

La ecuaciones homogéneas son un tipo de EDO que se pueden reducir a ecuaciones de
variables separadas mediante el cambio de variable conveniente.

Definición 2.7 (Función homogénea de grado k) Una función M (x, y) se dice

homogénea de grado k, siendo k un número real, si se verifica

M (αx, αy) = αk M (x, y) ,

para todo α ∈ R y cualquiera (x, y) ∈ Dom(M ) .

a
Hay que evitar problemas con el dominio de M garantizando que si (x, y) ∈ Dom(M ), entonces
(αx, αy) ∈ Dom(M ).

Ejemplo 2.12
La función M (x, y) = x2 + xy + y 2 verifica
M (αx, αy) = α2 (x2 + xy + y 2 = α2 M (x, y)
por tanto es una función homogénea de grado 2.
La función M (x, y) = 3x2 y 3 verifica
M (αx, αy) = 3α2 x2 α3 y 3 = α5 M (x, y)
por tanto es una función homogénea de grado 5.
√
La función M (x, y) = 2x + y verifica
p √
M (αx, αy) = 2αx + αy = αM (x, y)
por tanto es una función homogénea de grado 1/2.

Definición 2.8 (Ecuación diferencial homogénea) Una ecuación diferencial ordi-

naria de primer orden se dice ecuación homogénea si se puede escribir de la forma

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

Siendo M y N dos funciones homogéneas del mismo grado.

Método resolución EDOs homogéneas

Dada una ecuación diferencial homogénea M (x, y)dx + N (x, y)dy = 0

1. Realizar el cambio de variable y = ux (o x = vy) para obtener una ecuación

separable.

Si y = ux entonces dy = xdu + udx (Si x = vy, dx = vdy + ydv)

2. Resolver la ecuación separable que se obtiene.

3. Deshacer el cambio de variable para obtener la solución de la ecuación inicial.

26
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

x+y
Ejemplo 2.13 Resolver la ecuación y ′ = .
x−y
Reescribimos la ecuación en forma diferencial

(x + y)dx − (x − y)dy = 0 ,

siendo M (x, y) = x + y y N (x, y) = −(x − y) = y − x dos funciones homogéneas de grado

1, por tanto es una EDO homogénea. Procedemos a resolverla paso a paso.

1. Consideramos el cambio de variable y = ux, por tanto dy = udx + xdu.

2. Sustituimos en la ecuación diferencial

(x + ux)dx − (x − ux)(udx + xdu) ⇒ x(1 + u2 )dx − x2 (1 − u)du = 0

3. Reescribir la ecuación diferencial en variables separables.

1−u 1
du = dx .
1 + u2 x

4. Integrar
1−u
Z Z
1
du = dx
1 + u2 x
Se obtiene
1
arctan u − ln |1 + u2 |) = ln |x| + C, Cconstante arbitraria .
2

5. Deshacer el cambio de variable u = y/x, Se obtiene

y 1 y2
arctan − ln(1 + 2 ) = ln(x) + C, Cconstante arbitraria ,
x 2 x
lo cual se puede escribir
p
y = x tan (ln x2 + y 2 + C) Cconstante arbitraria ,

Ejemplo 2.14 Resolver (y 2 + xy)dx − x2 dy = 0.

Las funciones M (x, y) = y 2 + xy) y N (x, y) = −x2 son funciones homogéneas de grado
2, por tanto la ecuación diferencial es una ecuación homogénea

1. Consideramos el cambio de variable y = ux, por tanto dy = udx + xdu.

2. Sustituimos en la ecuación diferencial

(u2 x2 + ux2 )dx − x2 (udx + xdu) ⇒ u2 x2 dx − x3 du = 0

Si x ̸= 0 se puede escribir
u2 dx = xdu

3. Reescribir la ecuación diferencial en variables separables.

du dx
2
=
u x

27
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

4. Integrar Z Z
du dx
=
u2 x

Se obtiene
−1
= ln |x| + C, Cconstante arbitraria .
u
5. Deshacer el cambio de variable u = y/x, Se obtiene

−x
= ln |x| + C + C, Cconstante arbitraria ,
y

lo cual se puede escribir

−x
y= Cconstante arbitraria ,
ln |x| + C

2.4.3. Ecuaciones exactas y factores integrantes

En este epı́grafe se analizará la resolución tanto de las ecuaciones diferenciales exactas
como de las reducibles a exactas mediante factor integrante, es decir, aquellas que no
son exacta pero que al multiplicarlas por una función apropiada en cada caso pasan a
convertirse en ecuaciones diferenciales exactas.

Definición 2.9 Una ecuación diferencial de primer orden escrita en la forma

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

se dice ecuación diferencial exacta si existe una función G(x, y) tal que

∂G ∂G
= M (x, y) y = N (x, y) .
∂x ∂y
La solución de una ecuación exacta viene dada de forma implı́cita por las curvas de
nivel
G(x, y) = C .

La condición necesaria y suficiente para que una EDO sea exacta en un rectángulo R
del plano es
∂M ∂N
= , ∀(x, y) ∈ R
∂y ∂x
Dada una EDO de primer orden exacta

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

vamos a ver un método para encontrar la función G(x, y) que nos permite escribir la
solución de forma implı́cita como G(x, y) = C, siendo C una constante arbitraria. Esta
función G(x, y) debe satisfacer las condiciones

∂G ∂G
= M (x, y) y = N (x, y) .
∂x ∂y

28
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Método resolución EDOs exactas

Existen dos posibilidades

Si partimos de la relación
∂G
= M (x, y)
∂x
1. Integrar con respecto a x la ecuación anterior para obtener
Z
G(x, y) = M (x, y)dx + g(y) .

2. Calcular g(y), para ello

• Calcular la parcial respecto a la variable y en ambos lados de la ecuación
anterior
• Realizar el cambio
∂G
= N (x, y)
∂y
para calcular g ′ (y).
• Integrar g ′ (y).
3. La solución general será la familia uniparamétrica G(x, y) = C, siendo C
una constante arbitraria.

Si partimos de la relación
∂G
= N (x, y)
∂y
.

1. Integrar con respecto a y la ecuación anterior para obtener

Z
G(x, y) = N (x, y)dy + h(x) .

2. Calcular h(x), para ello

• Calcular la parcial respecto a la variable x en ambos lados de la ecuación
anterior
• Realizar el cambio
∂G
= M (x, y)
∂x
para calcular h′ (x).
• Integrar h′ (x).
3. La solución de la EDO será G(x, y) = C, siendo C una constante arbitraria.

Ejemplo 2.15 Comprueba que la ecuación diferencial 2x − 2yy‘ = 0 es exacta y determi-

nar su solución general.
En primer lugar es necesario reescribir la ecuación en forma diferencial

2xdx − 2ydy = 0 .

29
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Se obtiene ası́ que M (x, y) = 2x y N (x, y) = −2y .

A continuación se comprueba si esta ecuación es exacta

∂M

 = 0,
∂y


 ∂N = 0 .


∂x
Ambas parciales coinciden por lo que la ecuación diferencial es exacta. Procedemos a
continuación a resolverla.
Partimos de la relación
∂G
= M (x, y) = 2x
∂x
1. Integrar con respecto a x la ecuación anterior para obtener
Z
G(x, y) = 2xdx + g(y) = x2 + g(y) .

2. Calcular g(y), para ello

Calcular la parcial respecto a la variable y en ambos lados de la ecuación ante-

rior
∂G
= g ′ (y)
∂y
Realizar el cambio
∂G
= 2y
∂y
para calcular g ′ (y). Esto es
g ′ (y) = 2y
Integrar g ′ (y), obteniendo g(y) = y 2 .

3. Se obtiene ası́ G(x, y) = x2 + y 2 . La solución de la EDO exacta será x2 + y 2 = C,

siendo C una constante arbitraria.

Ejemplo 2.16 Comprobar que la EDO 2xydx + (x2 − 1)dy = 0 es exacta y resolverla.
En este caso la EDO está escrita en forma diferencial, ası́ es fácil de identificar que
M (x, y) = 2xy y que N (x, y) = x2 − 1.
Comprobemos la condición necesaria y suficiente de EDO exacta,

∂M

 = 2x ,
∂y


 ∂N = 2x .


∂x
Puesto que ambas parciales son iguales la ecuación es exacta. Procedemos a resolverla
a partir de la relación
∂G
= M (x, y) = 2xy
∂x
1. Integrar con respecto a x la ecuación anterior para obtener
Z
G(x, y) = 2xydx + g(y) = x2 y + g(y) .

30
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

2. Calcular g(y), para ello

Calcular la parcial respecto a la variable y en ambos lados de la ecuación ante-

rior
∂G
= x2 + g ′ (y)
∂y
Realizar el cambio
∂G
= N (x, y) = x2 − 1
∂y
para calcular g ′ (y). Esto es

x2 − 1 = x2 + g ′ (y)

Se obtiene ası́ que g ′ (y) = −1.

Integrar g ′ (y) = −1, obteniendo g(y) = −y.

3. Se obtiene ası́ G(x, y) = x2 y − y. La solución de la EDO exacta será x2 y − y = C,

siendo C una constante arbitraria.

Definición 2.10 (Factor integrante) Una función µ(x, y) se dice que es un factor
integrante para una ecuación diferencial de primer orden

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

si la ecuación
µ(x, y)M (x, y)dx + µ(x, y)N (x, y)dy = 0
es exacta.

Ejemplo 2.17

La función µ(x, y) = ex es un factor integrante en R2 de la ecuación y ′ + y = 0 .

La ecuación y ′ + y = 0 se escribe en forma diferencial de la forma ydx + dy = 0, lo
cual no es una ecuación exacta. Si se multiplica por ex se obtiene

ex ydx + ex dy = ,0

Se tiene ası́ las funciones M (x, y) = ex y e N (x, y) = ex .

Comprobemos la condición necesaria y suficiente de EDO exacta,

∂M

 = ex ,
∂y


 ∂N = ex .


∂x
Las dos parciales coinciden por lo que la ecuación ex ydx + ex dy = 0 es exacta.

La función µ(x, y) = 1/y es un factor integrante en R × (0, ∞) de la ecuación

y′ + y = 0 .

31
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

1 1
En efecto, considerando la ecuación dx+ dy = se obtiene M (x, y) = 1 y N (x, y) =
y y
y por tanto 
∂M

 = 0,
∂y


 ∂N = 0 .


∂x

Obsérvese que para que una función µ(x, y) sea un factor integrante de la ecuación
M (x, y)dx + N (x, y)dy = 0 se tiene que cumplir
 
∂(µM ) ∂(µN ) ∂µ ∂µ ∂N ∂M
= ⇐⇒ M −N = − µ
∂y ∂x ∂y ∂x ∂x ∂y

Esta ecuación permite establecer algunos criterios que permiten encontrar factores
integrantes. Por ejemplo,

El factor integrante solo depende de x, es decir µ = µ(x), entonces la condición

anterior se escribe
∂M ∂N
−
µ′
 
∂µ ∂M ∂N ∂y ∂x
N = − µ ⇐⇒ =
∂x ∂y ∂x µ N

∂M ∂N
−
∂y ∂x My − Nx
Por tanto la función solo depende de la variable x, esto es =
N N
f (x). En este caso el factor integrante será
R
(f (x)dx
µ(x) = e

Nx − My
El factor integrante solo depende de y, esto es, µ = µ(y) si g(y) = y el
M
factor integrante será

Nx − My
dy
R
R
g(y)dy M
µ(y) = e =e

Método resolución EDOs por factor integrante

Para resolver una ecuación por el método de los factores integrantes se procede del
siguiente modo:

1. Encontrar el factor integrante.

2. Multiplicar la ecuación por el factor integrante.

3. Verificar que la nueva ecuación es exacta.

4. Resolver la nueva ecuación mediante el método de resolución de las ecuaciones

exactas.

32
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Ejemplo 2.18 Resolver la ecuación diferencial y 2 + (xy + 1)y ′ = 0 sabiendo que admite
un factor integrante que es función de y.
Suponemos que µ(y) es un factor integrante de la ecuación escrita en forma diferencial

y 2 dx + (xy + 1)dy = 0,

Esto significa que la ecuación

µ(y)y 2 dx + µ(y)(xy + 1)dy = 0, ,

es exacta, esto es

∂(µ(y)y 2 ) ∂(µ(y)(xy + 1))

= ⇐⇒ µy y 2 + 2yµ(y) = yµ(y) ⇐⇒ µy y = −µ(y) .
∂y ∂x

µy −1
Integrando la ecuación = se obtiene
µ y

−1
ln |µ(y)| = − ln |y| ⇐⇒ µ(y) = .
y

Una vez calculado el nuevo factor integrante escribimos la nueva ecuación

1
−ydx − (x + )dy = 0
y

1
y comprobamos que sea una ecuación exacta, siendo M (x, y) = −y y N (x, y) = −(x + ).
y

∂M

 = −1 ,
∂y


 ∂N = −1 .


∂x

Las dos parciales coinciden por lo que la ecuación es exacta. Se aplica ahora el método
de resolución de las ecuaciones exactas.
Procedemos a resolverla a partir de la relación

∂G
= M (x, y) = −y
∂x

1. Integrar con respecto a x la ecuación anterior para obtener

Z
G(x, y) = −ydx + g(y) = −yx + g(y) .

2. Calcular g(y), para ello

Calcular la parcial respecto a la variable y en ambos lados de la ecuación ante-

rior
∂G
= −x + g ′ (y)
∂y

33
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Realizar el cambio
∂G 1
= N (x, y) = −(x + )
∂y y
para calcular g ′ (y). Esto es
1
−x − = −x + g ′ (y)
y
1
De aquı́ se obtiene g ′ (y) = −
y
Integrar g ′ (y), obteniendo g(y) = − ln |y|.
3. Se obtiene ası́ G(x, y) = −yx − ln |y|. La solución de la EDO exacta será −yx −
ln |y| = C, siendo C una constante arbitraria.

2.4.4. Ecuaciones lineales de primer orden

Definición 2.11 (EDO lineal de primer orden) Una EDO lineal de primer orden
es aquella que se puede escribir en la forma
dy
a1 (x) + a0 (x)y = g(x) .
dx

Definición 2.12 (Forma canónica o estándar) La forma canónica o estándar de

una EDO lineal de primer orden es aquella en la que a1 (x) = 1 esto es

y ′ + p(x)y = r(x) .

Además,

Si r(x) = 0 se dice EDO lineal homogénea de orden 1.

Si r(x) ̸= 0 se dice EDO lineal completa de orden 1

Ejemplo 2.19
La ecuación y ′ + (sin x + cos x)y = 0 es una EDO lineal homogénea de primer orden.
La ecuación y ′ + (sin x + cos x)y = x2 + 1 es una EDO lineal completa de primer
orden.
La ecuación y ′ + sin xy 2 = ex + 1 no es una EDO lineal.

Método de resolución de ecuaciones lineales

Antes de describir los métodos de resolución de las ecuaciones lineales de primer or-
den comentaremos algunos resultados teóricos de interés. Estos resultados se extenderán
posteriormente a las ecuaciones lineales de orden superior.
Teorema 2.1 (Estructura del conjunto de soluciones de la EDO homogénea) El
conjunto de soluciones de una ecuación lineal homogénea de primer orden tiene estructura
de espacio vectorial de dimensión uno (el orden de la ecuación). Una base de dicho espacio
es una solución no trivial de la ecuación homogénea.

34
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Teorema 2.2 (Solución de la ecuación lineal completa) La solución general de una

ecuación lineal completa se puede escribir como la suma de la solución general de la ecua-
ción lineal homogénea asociada y una solución particular de la ecuación lineal completa.

Resolución de una ecuación lineal homogénea de orden 1.

Consideremos una EDO lineal homogénea de orden 1 escriba en forma canónica, esto
es
y ′ + p(x)y = 0 .
La ecuación anterior es una EDO en variables separables ya que se puede escribir
(y ̸= 0)
dy
= −p(x)dx .
y
Integrando esta última ecuación se obtiene
Z Z Z
dy R
= − p(x)dx ⇐⇒ ln y = − p(x)dx ⇐⇒ y(x) = Ke− p(x)dx .
y
Nota 2.2 En el procedimiento anterior se habı́a omitido el caso y = 0 que es trivialmente
una solución de la EDO lineal homogénea. Sin embargo esta solución se corresponde con
la expresión anterior cuando K = 0. Por tanto, la solución general de la EDO lineal
homogénea es R
y(x) = Ke− p(x)dx

Resolución de la EDO lineal completa de primer orden

Método resolución EDOs lineales completas de orden 1

Dada una EDO lineal completa

y ′ + p(x)y = r(x),

los pasos a seguir para obtener su solución general son:

1. Resolver la EDO lineal homogénea asociada. Se obtiene ası́ la solución

R
yh (x) = Ke− p(x)dx
, K constante arbitraria

2. Encontrar una solución particular yp (x) de la ecuación completa.

3. Escribir la solución general de la completa como suma de las dos soluciones an-
teriores.
y(x) = yh (x) + yp (x) .

Llegados a este punto, es necesario considerar algún método que permita obtener una
solución particular de la EDO lineal completa. Veremos dos alternativas
Método de Lagrange o variación de parámetros.
En este método se supone que la solución particular de la EDO lineal completa
tiene la misma forma que la solución general de la EDO homogénea asociada, pero
suponiendo que la constante es una función de la variable independiente. Esto es

35
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

1. La solución general de la EDO homogénea es

R
yh (x) = Ke− p(x)dx
, K constante arbitraria

2. Suponemos que la solución general de la EDO completa es de la forma

R
yp (x) = K(x)e− p(x)dx

3. Determinamos la expresión de K(x) imponiendo que yp (x) sea solución de la

EDO completa. Se obtiene ası́
R Z
− p(x)dx
yp (x) = e r(x)dx

4. La solución general de la ecuación completa será

R
Z R 
y(x) = yh (x) + yp (x) = e− p(x)dx e p(x)dx r(x)dx + K .

Método del factor integrante

En este caso se considera la EDO completa como una ecuación diferencial que admite
un factor integrante de la forma
R
p(x)dx
µ(x) = e

y se resuelve la EDO exacta que se obtiene con este factor integrante.

R
p(x)dx
e (p(x)y − r(x))dx + dy = 0 .

Aplicando el método de resolución de las EDOs exactas se llega a la solución

R
Z R 
− p(x)dx p(x)dx
y(x) = e e r(x)dx + K .

Nota 2.3 Se puede verificar que la solución uniparamétrica

R
Z R 
− p(x)dx p(x)dx
y(x) = e e r(x)dx + K

no posee soluciones singurales y por tanto es la solución general de la EDO lineal.

Ejemplo 2.20 Calcular la solución general de la EDO

1 x
y′ = y+ .
x 1 + x2
Reescribiendo la EDO anterior en la forma
1 x
y′ − y=
x 1 + x2
1
se observa fácilmente que es una EDO lineal completa de primer orden donde p(x) = −
x
x
y r(x) = .
1 + x2
Resolveremos esta ecuación utilizando el método de Lagrange.

36
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

1
1. Resolver la ecuación homogénea asociada y ′ − y = 0. Esta es una ecuación de
x
variables separables y su solución es
Z
1
R dx
yh (x) = Ke − p(x)dx
= Ke x = Keln x = Kx .

2. Suponemos que la solución particular de la EDO completa es de la forma

yp (x) = K(x)x .,

3. Calculamos K(x) imponiendo que la función yp sea solución de la EDO completa.

Calculamos
yp′ (x) = K ′ (x)x + K(x) .

Sustituimos los valores de yp (x) e yp′ (x) en la EDO completa

1 x 1
yp′ (x) − yp (x) = r(x) ⇐⇒ K ′ (x)x + K(x) − K(x) = 2
⇐⇒ K ′ (x) =
x 1+x 1 + x2
Integrando se obtiene Z
1
K(x) = dx = arctan x
1 + x2

4. Escribir la solución particular de la EDO completa

yp (x) = x arctan x .

5. Escribir la solución general de la EDO completa

y(x) = yh (x) + yp (x) = Kx + x arctan x, K constante arbitraria .

Ejemplo 2.21 Resolver la EDO lineal y ′ + y = x .

En este caso la EDP lineal completa está escrita en forma canónica. Se observa que
p(x) = 1 y r(x) = x .
Vamos a calcular la solución general usando el método del factor integrante.
R R
p(x)dx dx
µ(x) = e =e = ex .

La ecuación exacta que hay que resolver es

ex (y − x)dx + ex dy = 0 ,

donde M (x, y) = ex (y − x) y N (x, y) = ex . Comprobemos en primer lugar que es una

ecuación exacta observando que las siguientes parciales coinciden.

∂M

 = ex ,
∂y


 ∂N = ex .


∂x
t

37
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Buscamos las soluciones G(x, y) = K, donde G(x, y) es una función que cumple

∂G ∂G
= M (x, y) = ex (y − x) y = N (x, y) = ex .
∂x ∂y
Partimos de segunda expresión. Integrando respecto a la variable y se obtiene
Z
G(x, y) = ex dy + f (x) = ex y + f (x) .

Derivando respecto a x
∂G
= yex + f ′ (x) ⇐⇒ ex (y − x) = yex + f ′ (x)
∂x
Por tanto
f ′ (x) = −xex .
Integrando Z
f (x) = − xex dx = −xex + ex .

Ası́ se obtiene la familia de soluciones G(x, y) = K que se escribe en este caso

ex y − xex + ex = K ⇐⇒ y = x − 1 + Kex .

2.4.5. Otras EDOs de primer orden.

En los ejemplos previos hemos visto los principales tipos de EDOs de primer orden. Sin
embargo existen otras muchas ecuaciones que a priori no se adaptan a ninguno de estos
tipos. En muchas ocasiones, tras un cambio de variable adecuado, nos encontramos ante
ecuaciones de los tipos ya estudiados. En este apartado comentaremos algunos ejemplos
aunque existen otros muchos ejemplos de ecuaciones reducibles a un caso ya conocido.

Ecuación de Bernoulli

Definición 2.13 (ecuación de Bernoulli) Una ecuación diferencial se dice de Ber-

noulli si puede escribirse en la forma

y ′ + f (x)y = q(x)y α , α ∈ R, α ̸= 0, 1 .

Toda ecuación de Bernoulli se puede escribir como una ecuación lineal mediante el cambio
de variable v = y 1−α .

Ejemplo 2.22 Resolver la EDO y ′ − 2y = 4xy 3 .

La EDO anterior es una ecuación de Bernoulli con f (x) = −2, q(x) = 4x y α = 3 .
Dividimos la ecuación y ′ − 2y = 4xy 3 por y −3 , obteniendo

y −3 y ′ − 2y −2 = 4x .

Aplicamos el cambio de variable v = y −2 . Entonces,

dv dv dy dy
= = −2y −3 = −2y −3 y ′ .
dx dy dx dx

38
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Sustituyendo en la ecuación previa se obtiene la ecuación lineal

v ′ + 4v = −8x .

La solución general de esta última ecuación es

1
v(x) = ke−4x − 2x + .
2
Finalmente, deshaciendo el cambio de variable se obtiene la familia de soluciones
1
y −2 = ke−4x − 2x + .
2

Ecuación reducible a separable

Las ecuaciones diferenciales de la forma

y ′ = f (Ax + By + C), conA ̸= 0, B ̸= 0 .

se transforman en una ecuación en variables separables mediante el cambio

u = Ax + By + C .

Ejemplo 2.23 Resolver la EDO y ′ = (−2x + y)2 − 7 .

La ecuación anterior es de la forma y ′ = f (Ax + By + C) siendo

A = −2, B = 1, C = 0 y f (z) = z 2 − 7 .

Considerando el cambio de variable u = −2x + y se obtiene u′ = −2 + y ′ y sustituyendo

en la ecuación
u′
u′ + 2 = u2 − 7 ⇐⇒ 2 =1
u −9
que es una ecuación en variables separables.
Integrando se obtiene
1 u
arctan( ) = x + C , C ∈ R.
3 3
Finalmente, deshaciendo el cambio de variable se obtiene la solución de la ecuación
diferencial de partida
1 −2x + y
arctan( ) = x + C , C ∈ R.
3 3

2.5. Teoremas de existencia y unicidad de solución

En las secciones previas hemos visto que la solución general de una EDO es una familia
de funciones. En la mayorı́a de las aplicaciones prácticas, lo que nos interesa es una solución
particular de la familia, aquella que es solución de un PVI o de un problema de contorno.
En este punto surgen una pregunta de modo natural: ¿existe siempre solución de un
PVI?. En caso de existir ya hemos visto algún ejemplo ( Ejemplo 2.2) en el cual existen al
menos dos soluciones al mismo PVI, ¿bajo que condiciones se puede asegurar la unicidad
de solución?
Los siguientes teoremas proporcionan condiciones suficientes (no necesarias) para la
existencia y unicidad de solución de un PVI asociado a una EDO de primer orden.

39
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Teorema de existencia y unicidad de solución (EDOs de orden 1)

Teorema 2.3 (Teorema de Peano-Picard) Se considera el problema de valor ini-

cial en forma explı́cita  ′
y = f (x, y)
y(x0 ) = y0
Si f es una función continua en un rectángulo

R = {(x, y) : a ≤ x ≤ y, c ≤ y ≤ d}

que contiene al punto (x0 , y0 ) en su interior, entonces el PVI tiene al menos una solu-
ción en el intervalo [x0 , x0 + h) para algún h > 0 .
∂f
Si además es continua en R, la solución es única en [x0 , x0 + h) para algún h > 0 .
∂y

Ejemplo 2.24 Estudiar la existencia y unicidad de solución de los siguientes PVI.

y ′ = x2 − xy 3


y(1) = 6

x2 y ′ = −1


y(x0 ) = y0

En el primer apartado, el PVI ya está escrito en forma explı́cita, por tanto la función
f (x, y) a analizar es f (x, y) = x2 − xy 3 . Esta función tiene como dominio todo R2 y es
una función continua en todo punto, por tanto está garantizada la existencia de solución
del PVI en un intervalo [1, 1 + h), .h > 0.
Para estudiar la unicidad de solución se deriva la función f respecto la variable de-
pendiente, esto es
∂f
= −3xy 2 .
∂y
Esta función es continua en todo R2 y por tanto en un rectángulo abierto que contenga
al punto (1, 6). Por tanto, la solución del PVI es única en un intervalo [1, 1 + h), siendo
h > 0 suficientemente pequeño.
En el segundo ejemplo, en primer lugar hay que escribir la EDO en forma explı́cita
para obtener la función f (x, y). Puesto que y ′ = −1/x2 se tiene f (x, y) = −1/x2 .
La función f (x, y) = −1/x2 es continua en su dominio, esto es, en

{(x, y) ∈ R2 |x ̸= 0}

Por tanto, se puede garantizar la existencia de solución en [x0 , x0 + h), h > 0 siempre
que x0 ̸= 0.
Para analizar la unicidad de solución se analiza ∂f /∂y = 3/x4 , que es una función
continua en {(x, y) ∈ R2 |x ̸= 0}, por tanto se puede garantizar la unicidad de solución en
[x0 , x0 + h), h > 0 siempre que x0 ̸= 0 solo si x0 ̸= 0.

Nota 2.4 En el caso lineal, el teorema de existencia y unicidad de solución de un PVI se

reduce al siguiente resultado:

40
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Teorema de existencia y unicidad de soluciones EDOS lineales de orden 1

Dado el PVI
y ′ + p(x)y = r(x)


y(x0 ) = y0 .
Si p(x) y r(x) son funciones continuas en un intervalo abierto I, con x0 ∈ I,
entonces el PVI tiene una única solución.

2.6. Métodos numéricos

Consideremos un PVI de la forma
 ′
y = f (x, y), x ∈ (x0 , x0 + T ]
y(x0 ) = y0 .

en el cual las funciones f y ∂f /∂y son funciones continuas y por tanto está garantizada la
existencia y unicidad de solución.
Aunque se puede garantizar la existencia y unicidad de solución de este PVI, no siempre
es posible expresar esta solución en términos de las funciones elementales o de forma
explı́cita. Para evitar este problema se busca una aproximación de la solución en un número
finito de puntos xk del intervalo [x0 , x0 + T ]. Este procedimiento se denomina resolución
mediante métodos de discretización.
A lo largo de este apartado del tema utilizaremos la siguiente notación:

y(xk ) representa el valor exacto de la solución en el punto xk ,

yk es el valor aproximado de la solución en el punto xk . Este valor se obtiene después

de aplicar un método numérico.

La diferencia |y(xk ) − yk | es el error de discretización.

Por sencillez solo consideraremos métodos numéricos en los que se realiza una partición
uniforme del intervalo [x0 , x0 + T ].

Esquema general de un método de k pasos

Los métodos que vamos a estudiar son métodos de k-pasos. Es esquema general de
estos métodos es el siguiente:

1. Se realiza una partición uniforme del intervalo [x0 , x0 + T ] en N partes iguales, es

decir, con N + 1 puntos. El paso de integración es h = T /N y los puntos de la malla
son
xj = x0 + jh, j = 0, 1, 2, . . . , N .

2. Conocemos los valores y0 , y1 , · · · , yk−1 .

3. El resto de los valores aproximados yk , yk+1 , ... se calculan de forma progresiva en

función de los k valores previos.

41
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

2.6.1. Método de Euler explı́cito

El método de Euler explı́cito es el más sencillo desde el punto de vista conceptual pero
el de menor interés desde el punto de vista práctico.
La idea de este método consiste en aproximar la curva solución (desconocida) por la
recta tangente a la solución.
Una vez considerada la partición del intervalo [x0 , X0 +T ], la idea es aproximar y(x1 ) =
y(x0 + h) por el valor y1 correspondiente a la recta tangente en el punto x1 . Teniendo en
cuenta que
y(x0 + h) − y(x0 )
y ′ (x0 ) = f (x0 , y0 ) = lı́m
h→0 h
se considera la siguiente aproximación

y(x0 + h) − y(x0 )
y ′ (x0 ) = f (x0 , y0 ) ≈ ,
h
la cual es equivalente a decir

y(x1 ) = y(x0 + h) ≈ y(x0 ) + hf (x0 , y0 ) = y0 + hf (x0 , y0 ) .

Denotamos por y1 = y0 + hf (x0 , y0 ). Ahora el proceso se repite comenzando en el

punto (x1 , y1 ) y tomando como derivada de la solución en ese punto f (x1 , y1 ). Repitiendo
este proceso se llega al algoritmo del método de Euler, que es un método de 1 paso.

x0 , y0 , h valores conocidos
yk+1 = yk + hf (xk , yk ), k = 0, 1, . . . , N − 1 .

Se puede demostrar que en este caso el error de discretización satisface

|y(xk ) − yk | ≤ Ch, k = 0, 1, . . . , N

siendo C una constante independiente de h. Esto es, el método de Euler es un método de

orden 1, esto es, si el paso de la malla se reduce a la mitad, los errores se reducen a la
mitad.

2.6.2. Método de Euler mejorado

El método de Euler mejorado intenta disminuir los errores que se producen con el
método de Euler cuando, por ejemplo, la solución en el punto x1 crece mucho más deprisa
que en el punto x0 . La idea del método de Euler mejorado consiste en aproximar la
pendiente de la recta tangente en el punto (xk , yk ) por la media aritmética de las pendientes
de las rectas tangentes en los extremos xk y xk+1 del intervalo a considerar.
Por ejemplo, partimos del punto (x0 , y0 ), la recta tangente en este punto tiene pendiente
′
y (x0 ) = f (x0 , y0 ). El siguiente punto a considerar es (x1 , y(x1 )), como no se conoce el valor
exacto y(x1 ) este se aproxima según la fórmula que nos proporciona el método de Euler,
considerando una predicción
Y1 = y0 + hf (x0 , y0 ) .

A continuación se considera la aproximación

h
y1 = y0 + (f (x0 , y0 ) + f (x1 , Y1 )) .
2

42
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Repitiendo este proceso se llega al siguiente esquema correspondiente con un método

de 2 pasos.


 x0 , y0 , h valores conocidos

Y1 = yk + hf (xk , yk ) (predicción)
 h
k+1 = yk + (f (xk , yk ) + f (xk+1 , Y1 )), k = 0, 1, . . . , N − 1 .

 y
2
En este caso el error de discretización cumple

|y(xk ) − yk | ≤ Ch2 , k = 0, 1, . . . , N

siendo C una constante independiente de h. Esto es, el método de Euler mejorado es un

método de orden 2, esto es, si el paso de la malla se reduce a la mitad, los errores se
reducen a la cuarta parte.

2.6.3. Métodos Runge Kutta

Tanto el método de Euler como el método de Euler mejorado son casos particulares
de una familia de métodos más amplia conocida como métodos Runge Kutta. Estos méto-
dos se basan en considerar distintas aproximaciones basadas en métodos de integración
numérica. Ası́ por ejemplo el método de Euler se corresponde con aproximar la integral
Z x1
f (x, y(x))dx
x0

por la fórmula del rectángulo mientras que en el método de Euler mejorado se considera
la fórmula del trapecio.
Los métodos Runge Kutta en general son más eficientes que el método de Euler o Euler
mejorado. De todos los métodos Runge Kutta el más utilizado es el método Runge Kutta
4 que se basa en el siguiente algoritmo


 x0 , y0 , h datos conocidos


P redicciones




Y1 = yk ,







 h

 Y2 = yk + f (xk , Y1 ),
2



h
 Y3 = yk + f (xk + h/2, Y2 )



 2



 Y4 = yk + hf (xk + h/2, Y3 )




Aprox.




 yk+1 = yk + h (f (xk , Y1 ) + 2f (xk + h/2, Y2 ) + 2f (xk + h/2, Y3 ) + f (xk + h, Y4 )) ,



6
En este caso el error de discretización cumple

|y(xk ) − yk | ≤ Ch4 , k = 0, 1, . . . , N

siendo C una constante independiente de h. Esto es, el método Runge Kutta 4 es un

método de orden 4, esto es, si el paso de la malla se reduce a la mitad, los errores se
reducen a la dieciseisava parte.

43
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

En general, un método Runge Kutta de s etapas se corresponde con las siguientes

fórmulas las cuales se pueden escribir mediante la tabla de Butcher que resume los coefi-
cientes de las fórmulas del método.

c1 a11 a12 · · · a1s

c2 a21 a22 · · · a2s
.. .. .. .
. . . · · · ..
cs as 1 as 2 · · · ass
b1 b2 · · · bs



 x0 , y0 , h datos conocidos


P redicciones




Y1 = yk + a11 hf (xk + c1 h, Y1 ) + a12 hf (xk + c2 h, Y2 ) + . . . + a1s hf (xk + cs h, Ys )






Y2 = yk + a21 hf (xk + c1 h, Y1 ) + a22 hf (xk + c2 h, Y2 ) + . . . + a2s hf (xk + cs h, Ys )






Y3 = yk + a31 hf (xk + c1 h, Y1 ) + a32 hf (xk + c2 h, Y2 ) + . . . + a3s hf (xk + cs h, Ys )


 ..


 .
..






 .



 Ys = yk + as1 hf (xk + c1 h, Y1 ) + as2 hf (xk + c2 h, Y2 ) + . . . + ass hf (xk + cs h, Ys )





 Aprox.
yk+1 = yk + b1 hf (xk + c1 h, Y1 ) + b2 hf (xk + c2 h, Y2 ) + . . . + bs hf (xk + cs h, Ys ),


2.7. Aplicaciones
En esta parte del capı́tulo vamos a ver algunas de las aplicaciones de las EDOs de
primer orden. Como veremos, son muchas las situaciones que se modelan con este tipo de
ecuaciones. En [?] puede encontrar otros muchos ejemplos detallados.
Estas notas se centran en los ejemplos más caracterı́sticos en los que aparecen EDOs
de primer orden, en particular:

Trayectorias ortogonales e isogonales.

Modelos de desintegración radiactiva.

Dinámica de poblaciones.

Ley de enfriamiento de Newton.

Circuitos LR y RC en serie.

Problemas de mezclas.

2.7.1. Trayectorias ortogonales

Se dice que dos curvas planas son ortogonales si las rectas tangentes a ambas en
los puntos de corte son rectas perpendiculares. Existen problemas comunes de fı́sica que
requieren conocer familias de curvas ortogonales a una familia dada. Por ejemplo, en
electrostática, las lı́neas de fuerza son ortogonales a las curvas equipotenciales.

44
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

Trayectorias ortogonales a una familia dada

Dada una familia de curvas planas, se puede construir una familia de curvas con tra-
yectorias ortogonales del siguiente modo:

F (x, y, c) = 0 familia de curvas.

Derivando respecto a x y eliminando c se llega a la EDO f (x, y, y ′ ) = 0 asociada

a la familia.

Cambiando el papel de y ′ por −1/y ′ , es decir , escribiendo f (x, y, −1/y ′ ) = 0 se

obtiene la EDO de la familia de curvas ortogonales.

Resolviendo esta última ecuación se obtiene G(c, y, c̃) = 0, la familia de curvas

ortogonales.

Ejemplo 2.25 Obtener las curvas ortogonales a las hipérbolas de ecuación y = C/x,
siendo C una constante arbitraria.

1. Calcular la EDO de primer orden que tiene a esta familia como solución.

Derivamos respecto a x:
y ′ = −C/x2
A partir del sistema de ecuaciones
 −C
 y′ = 2

x
 y=C

x
Despejando C en la segunda ecuación y sustituyendo en la primera se obtiene
la EDO
−y
y′ = (2.2)
x
que tiene como solución la familia de hipérbolas inicial.

2. Buscamos una familia de curvas cuyas pendientes sean ortogonales a las de la prime-
ra familia, por tanto las curvas que se buscan son las soluciones de la ecuación que
se obtiene al sustituir en (2.2) y ′ por −1/y ′ , esto es, las soluciones de la ecuación
−1 −y
′
= ⇐⇒ yy ′ = x
y x

3. Integrar la ecuación previa para obtener las curvas ortogonales.

y2 x2
Z Z
′
yy dy = xdx ⇐⇒ = + K1 .
2 2

Por tanto, la familia de curvas ortogonales a las hipérbolas de ecuación Y = C/x es

la familia de hipérbolas de ecuación y 2 − x2 = 2K1 .

45
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Figura 2.7: Trayectorias ortogonales

2.7.2. Trayectorias isogonales

Se dice que dos familias de curvas son trayectorias isogonales ( oblicuas) si las rectas
tangentes en cada punto de intersección forman un ángulo constante α ̸= π/2 .

Trayectorias isogonales

Dada una familia de curvas planas, se puede construir una familia de curvas con tra-
yectorias que forman un ángulo α siguiendo el siguiente procedimiento

F (x, y, c) = 0 familia de curvas.

Derivando respecto a x y eliminando c se llega a la EDO f (x, y, y ′ ) = 0 asociada

a la familia.

Cambiandoa el papel de y ′ por

y ′ + tan α
1 − y ′ tan α
se obtiene
y ′ + tan α
f (x, y, )=0
1 − y ′ tan α
la EDO de la familia de curvas que forman un ángulo α con la inicial.

Resolviendo esta última ecuación se obtienen las trayectorias buscadas.

a
Para obtener este cambio nos basamos en que el ángulo entre dos rectas tangentes que se cortan
en un punto cumple,
y ′ − f (x, y)
tan α = .
1 + y ′ f (x, y)

Ejemplo 2.26 Determinar la familia de trayectorias que forman un ángulo π/4 con

1
y= ,
x+c

siendo C una constante arbitraria.

46
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

1. Calcular la EDO de primer orden que tiene a

1
y= .
x+c

como familia de soluciones.

Derivamos respecto a x:
y ′ = −1/(x + C)2

A partir del sistema de ecuaciones

−1

′
 y = (x + C)2



 y= 1


x+C

1
Despejando C en la segunda ecuación C = − x y sustituyendo en la primera
y
se obtiene la EDO
y ′ = −y 2 (2.3)

que tiene como solución la familia inicial.

2. Realizamos el cambio de y ′ por

y ′ + tan α y′ + 1
=
1 − y ′ tan α 1 − y′

donde se ha utilizado que tan (π/4) = 1.

Se obtiene ası́ la ecuación diferencial

y′ + 1
= −y 2 .
1 − y′

3. Resolver la ecuación anterior para obtener la familia isogonal.

La ecuación previa puede escribirse del siguiente modo:

y2 − 1
dy = dx
y2 + 1

por lo que es una ecuación en variables separables. Integrando a ambos lados

y2 − 1
Z Z Z Z
2
dy = dx ⇐⇒ (1 − )dy = dx
1 + y2 1 + y2

Obteniendo,
y − 2 arctan y = x + C .

47
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Figura 2.8: Trayectorias isogonales

2.7.3. Desintegración radiactiva

Una sustancia radioactiva tiende a descomponerse con el tiempo formando nuevas
sustancias y liberando gran cantidad de energı́a. Se ha comprobado experimentalmente que
la velocidad con que una sustancia radioactiva se descompone es directamente proporcional
a la cantidad de sustancia existente en dicho instante, es decir, satisface la ecuación
y ′ (t) = Ky(t),
donde y(t) representa la cantidad del isótopo radiactivo y K una constante que depende
de la sustancia considerada.
La ecuación anterior es un ejemplo de ecuación en variables separables ya que se puede
escribir en la forma
dy
= Kdt .
y
Integrando se obtiene
y(t) = CeKt .

Ejemplo 2.27 Se define la vida media tm de una sustancia radioactiva como el tiempo
necesario para que dicha sustancia se reduzca a la mitad. Si tenemos una cantidad inicial
N0 de una sustancia, calcular la vida media de la misma.
Para resolver este problema es necesario resolver el siguiente problema de valor inicial
 ′
y = Ky
y(0) = N0
La solución de este problema es
y(t) = N0 eKt .
Ası́, la vida media es el instante de tiempo tm tal que y(tm ) = N0 /2 .
Entonces, tm es la solución de la ecuación
N0
= N0 eKtm
2
esto es,
ln 2
tm = − .
K

48
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

2.7.4. Dinámica de poblaciones

Existen numerosos modelos sobre el crecimiento y decrecimiento de poblaciones. Uno
de los modelos más sencillos establece que la rapidez con la que crece una población es
proporcional a su tamaño. Se obtiene ası́ la ecuación diferencial

p′ (t) = kp(t) ,

siendo p(t) el número de individuos de la población en el instante de tiempo t. Resolviendo

este problema, se obtiene que la población en el instante t viene dada por la expresión

p(t) = Cek t C ∈ R .

Por ejemplo, las colonias de bacterias, en muestras en un laboratorio, siguen este

modelo; sin embargo, no es un modelo realista de la dinámica de una población general,
ya que
lı́m p(t) = ∞ .
t→∞

Un modelo más realista es la ecuación logı́stica que describe el crecimiento de muchas

poblaciones según el modelo

dy
= y(a − by), a, b > 0 .
dt

Esta ecuación sigue el modelo de ecuación de variables separables, bajo la condición

inicial y(0) = y0 > 0. Integrando esta ecuación se obtiene la siguiente solución.

ay0
y(t) = .
by0 + (a − by0 )e−at

Ejemplo 2.28 El número de células que componen un tumor es, inicialmente, 104 . El
crecimiento de dicho tumor responde a la ley logı́stica

1
y ′ (t) = ry(t)(1 − y(t)),
k

siendo y(t) el número de células en el tiempo t, midiendo el tiempo en dı́as, r = 0,2 y

k = 22 · 107 . Para este tipo de tumores la velocidad de crecimiento es máxima cuando
t = 50 dı́as. Calcular el tamaño de la población en dicho instante.
Para resolver este problema observamos que la población de células sigue la ecuación
r 0,2
logı́stica con a = r = 0,2 y b = = .
k 22 · 10−7
Sustituyendo estos valores en la solución de la ecuación logı́stica se obtiene

22 · 107
y(t) = .
1 + 21999e−0,2t

La población pasados 50 dı́as es

22 · 107 22 · 107 k
y(50) = −0,2·50
= −10
≈
1 + 21999e 1 + 21999e 2

49
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

2.7.5. Ley de enfriamiento de Newton

La Ley de enfriamiento de Newton afirma que la temperatura de un cuerpo T rodeado
por aire a temperatura T0 varı́a de modo que el ritmo de variación de su temperatura es
proporcional a la diferencia entre la temperatura del cuerpo y la temperatura ambiente,
es decir,
dT
= k(T − T0 ) ,
dt
donde T (t) es la temperatura del cuerpo en el instante t. Esta es una ecuación de variables
separable ya que se puede escribir en la forma
dT
= kdt .
T − T0
Integrando se obtiene

ln|T − T0 | = kt + c, c constante ⇐⇒ T − T0 = Aekt , A ∈ R .

Se obtiene finalmente
T (t) = T0 + Aekt , A ∈ R.

Ejemplo 2.29 Un objeto que inicialmente está a 120◦ C se pone en una habitación a
20◦ C. Al cabo de una hora, su temperatura es de 70◦ C. ¿Cuánto tiempo más tiene que
transcurrir para que baje a 40◦ C?
Según la Ley de Enfriamiento de Newton, la temperatura en el instante t viene dada
por la expresión
T (t) = T0 + Aekt , A ∈ R .
Sustituyendo los datos del problema, T0 = 20 y T (0) = 120 del obtiene el valor de la
constante A
T (0) = T0 + Aek·0 ⇐⇒ 120 = 20 + A ⇐⇒ A = 100 .
Sabemos que después de una hora la temperatura es 70◦ C, es decir,

T (1) = 70 ⇐⇒ 20 + 100ek = 70 ⇐⇒ k = − ln 2 .

Ası́, para este caso particular es

T (t) = 20 + 100e−t ln 2 .

Para calcular el tiempo necesario, en horas, para alcanzar los 40◦ C, sustituimos este
valor en la ecuación anterior, obteniendo

40 = 20 + 100e−t ln 2 ⇐⇒ ln 5 = t ln 2 ⇐⇒ t = ln 5/ ln 2 ≈ 2,3219 .

Es decir, al cabo de ln 5/ ln 2 horas del instante inicial se alcanzan los 40◦ C.

2.7.6. Circuitos eléctricos en serie

En el capı́tulo 1 hemos descrito la ecuación diferencial que permite calcular la intensi-
dad de corriente a lo largo de un circuito LRC.
Aunque la ecuación asociada a un circuito es, en general, una EDO de segundo orden
existen situaciones particulares donde se obtiene una EDO de primer orden. Algunas de
ellas son las siguientes:

50
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

En el caso de un circuito LR esta ecuación se reduce a la EDO lineal

Li′ (t) + Ri = E .

En el caso de un circuito RC la ED
1
Rq ′ + q = E(t) .
C

Ejemplo 2.30 Un condensador de capacidad C, con una carga Q, se descarga a través de

una resistencia R. Obtener la curva de descarga, es decir, la función q(t) que determina
la carga del condensador en cualquier instante t.
Para resolver este problema es necesario imponer la condición de que la caı́da de po-
tencial en todo el circuito es cero. Se obtiene ası́ la siguiente EDO lineal homogénea de
primer orden:
1
Rq ′ + q = 0
C
con la condición inicial q(0) = Q.
Resolviendo esta EDO se obtiene que
1
− t.
q(t) = Qe RC

2.7.7. Mezclas
En los problemas de mezclas aparecen sustancias, las cuales se mezclan dentro de un
recipiente de volumen V0 conocido.
Supongamos que inicialmente tenı́amos una cantidad de x0 kilogramos de una sustancia
diluida en una concentración de X0 /V0 Kg/m3 y que introducimos otras solución que
contiene una concentración b Kg/m3 de dicha sustancias que es introducida en el recipiente
a una velocidad de e m3 /sg.
Supongamos que además, sacamos parte de la solución que se produce dentro del
recipiente a una velocidad de f m3 /sg.
Si denotamos por y(t) la cantidad de sustancia dentro del recipiente por unidad de
tiempo, tenemos que la variación de dicha cantidad viene dada por la diferencia entre la
velocidad de entrada ve y la velocidad de salida vs de dicha sustancia.
Obsérvese que
ve = be Kg/sg
mientras que
vs = y(t)/V (t)
siendo V (t) el volumen de disolución en el recipiente por unidad de tiempo, es decir,

V (t) = V0 + et − f t .

El problema que se acaba de describir se corresponde con el siguiente PVI

 y
 y ′ = be − f

 V + et − f t
0


 y(0) = X0


51
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

Ejemplo 2.31 Una solución de salmuera fluye a razón constante de 6L/min hacia el
interior de un depósito que inicialmente contiene 5Kg de sal disueltos en 50L de disolu-
ción. La solución contenida en el depósito se mantiene bien agitada para que la mezcla
sea homogénea y fluye hacia el exterior con la misma rapidez. Si la concentración de sal
en la salmuera que entra en el depósito es de 0,5Kg/L, determinar la cantidad de sal
presente en el depósito después de t minutos. ¿Cuándo alcanzará la concentración de sal
en el depósito el valor de 0,3Kg/L?
Denotemos por y(t) los kg de sal dentro del depósito en el instante t.
Debido a que la velocidad de entrada y salida de la concentración es la misma, el
depósito siempre mantiene los 50L de disolución.
Teniendo en cuenta la descripción del problema que se ha realizado previamente sabe-
mos que esta situación se modeliza a partir del siguiente PVI

y(t) = 6 · 0,5 − 6x/50
y(0) = 5
Resolviendo esta ecuación de variables separables se obtiene

y(t) = 25 + Ce−3t/25 .

A partir de la condición inicial se obtiene que C = −20 y por tanto la cantidad de sal
presente en el depósito al cabo de t minutos es

y(t) = 25 − 20e−3t/25 .

Si denotamos por C(t) la concentración de sal en el instante t, se cumple

1 2 −3t/25
C(t) = y(t)/50 = − e .
2 5
25 ln 2
A partir de esta expresión, el valor de t para el cual C(t) = 0,3 es t = .
3

2.8. Ejercicios propuestos

1. Verificar que y = Cx + C 2 , siendo C una constante arbitraria, es una familia unipa-
ramétrica de soluciones de la ecuación diferencial y ′2 + xy ′ = y. Determinar el valor
K para que y = Kx2 sea una solución singular de la ecuación anterior.
2. Se considera el PVI √
y ′ = 4x −
 3
y−1
y(3) = 1

¿Se puede garantizar que existe solución única para este problema de valor inicial?
3. Comprobar que y(x) = 0 e y(x) = (x − 2)3 son soluciones del PVI,

y ′ = 3y 2/3


y(2) = 0

¿Contradice esto el teorema de existencia y unicidad de solución?

4. Demostrar que el siguiente PVI tiene solución única
 ′
y = x ln xy
y(1) = 2

52
Silvia Vilariño Fernández Curso 2025-2026 [Link]-GIE

5. Determinar si las siguientes ecuaciones diferenciales son lineales. En caso negativo

justifique su respuesta;

y ′ − (y − x)2 = 0.
√
y ′ + ex y = 2 x.
yy ′ + xy = x.
y ′ − 3exy = x3 .

6. Resolver, utilizando el método de separación de variables, el PVI

 ′
y = 2xy 2
y(2) = 1

7. Resolver, utilizando el método de separación de variables, el PVI

 ′
y = (2x + 3)(y 2 − 4)
y(0) = −1

8. Resolver la siguiente ecuación diferencial

y−x
y′ = .
y+x

9. Obtener la solución general de la EDO exacta

y
y′ = .
2y ln y + y − x

10. Obtener todas las funciones f (x) para que la ecuación diferencial

f (x)y ′ + y + x2 = 0

tenga un factor integrante µ(x) = x, e integrar la ecuación.

11. La ecuación diferencial

(y 4 − 2y 2 )dx + (3xy 3 − 4xy + y)dy = 0

tiene un factor integrante que es función de xy 2 . Obtenerlo y resolver la ecuación.

2
12. Resolver la ecuación lineal y ′ + 2xy = 2xe−x .

13. Resolver el siguiente PVI

xy ′ + 2y = sin x


y(π/2) = 0 .

14. Resolver la ecuación xy ′ + y = x4 y 3 .

15. Determinar las trayectorias ortogonales a las circunferencias centradas en el origen

x2 + y 2 = C .

16. Halla las trayectorias ortogonales a la familia de curvas y 2 + 2ax = a2 , donde a es

una constante arbitraria.

53
[Link]-GIE Curso 2025-2026 Silvia Vilariño Fernández

17. Hallar la familia de trayectorias que forman un ángulo π/4 con la familia de rectas
y = cx

18. Un depósito contiene 10 litros de salmuera con 2 gramos de sal disueltos. En un

cierto instante comienza a entrar en el depósito salmuera que contiene 1 gramo de
sal por litro, a razón de 3 litros por minuto y, al mismo tiempo, a salir del depósito
la disolución, supuesta homogénea en todo momento, a razón de 4 litros por minuto.
Hallar la cantidad de sal que contiene el depósito en un instante cualquiera t.

19. En un circuito eléctrico con autoinducción L y resistencia R en serie, la relación

entre voltaje e intensidad i es
Li′ + Ri = V.
Si la tensión y la resistencia varı́an uniformemente durante un minuto desde 0 hasta
120 voltios y desde 0 hasta 120 ohmios, respectivamente, la inductancia L = 1 henrio
y la corriente inicial I0 , obtener la función i(t) de la intensidad que recorre el circuito
respecto del tiempo en dicho minuto.

20. Compruebe que si y1 (x) es una solución de la ecuación no lineal homogénea de

primer orden y ′ + P (X)y 2 = 0, entonces Cy1 (X), C ∈ R no es una solución de dicha
ecuación.

21. Compruebe que si y1 (x) es una solución de la ecuación lineal homogénea de primer
orden y ′ + P (X)y = 0, entonces Cy1 (X), C ∈ R es una solución de dicha ecuación.

22. Compruebe que si y1 (x) es una solución de la ecuación lineal completa de primer
orden y ′ + P (X)y = Q(x), entonces Cy1 (X), C ∈ R no es una solución de dicha
ecuación.

54