T EMA 4: E CUACIONES DIFERENCIALES DE PRIMER

ORDEN

Dpto. de Matemática Aplicada y Ciencias de la Computación

Grado en Ingeniería de Tecnologías de Telecomunicación

25 de abril de 2024

Cálculo II EDOs de primer orden 0/29

C ONTENIDOS DEL TEMA

1 Introducción
2 Conceptos generales
3 EDOs de primer orden
4 Métodos elementales de resolución de EDOs
Ecuaciones de variables separables
Ecuaciones homogéneas
Ecuaciones exactas
Ecuaciones lineales
5 Modelado de procesos físicos y aplicaciones
6 Resolución numérica
Método de Euler básico
Método de Euler mejorado
Métodos de Euler
Método de Runge-Kutta
7 Otros métodos aproximados

Cálculo II EDOs de primer orden 1/29

 Introducción

M OTIVACIÓN

Existen infinidad de problemas en Ciencia e Ingeniería que pueden

ser expresados (modelizados) en términos de ecuaciones
diferenciales.
En ocasiones las ecuaciones diferenciales que debemos resolver
son relativamente simples y podremos encontrar soluciones exactas
al problema con más o menos dificultad recurriendo a métodos
analíticos.
En otras ocasiones, deberemos recurrir a esquemas numéricos
para determinar soluciones aproximadas.

Cálculo II EDOs de primer orden 2/29

 Introducción

A LGUNAS APLICACIONES

dy
La tasa de crecimiento de una población en función del tiempo: = αy
dt
dT
Velocidad de enfriamiento de una sustancia: = −α(T − T0 )
dt
Tasa de variación de volumen de agua en un depósito cilíndrico a una
dV
velocidad v: = −πr 2 v
dt
Conservación del momento lineal en la atmósfera:
dv
= −α∇p − ∇ϕ + F − 2Ω × v
dt
∂ρ
Conservación de la masa en la atmósfera: = −∇ · (ρv)
∂t

Cálculo II EDOs de primer orden 3/29

 Conceptos generales

Las ecuaciones diferenciales son ecuaciones que contienen las derivadas

de una o más variables dependientes respecto de una o más variables
independientes.
En las ecuaciones diferenciales, las incógnitas son funciones.
El orden n de una ecuación diferencial es el mayor orden de derivación que
aparece en la ecuación.
Tipos de ecuaciones diferenciales:
Ecuaciones diferenciales ordinarias (EDOs): relacionan la variable
independiente (por ejemplo, x) con alguna derivada de la función
incógnita (y(x), y ′ (x), ..., y (n (x)): F (x, y(x), y ′ (x), ..., y (n (x)) = 0
Ecuaciones en derivadas parciales (EDPs): tienen más de una variable
independiente.
Ejemplos:
dy
dx
+ 5y = ex , buscamos y = y(x), EDO, orden 1
d2 x dx
dt2
− dt
+ 6x = 0, buscamos x = x(t), EDO, orden 2
∂2u ∂2u
∂x2
+ ∂y 2
= 0, buscamos u = u(x, y), EDP, orden 2

Cálculo II EDOs de primer orden 4/29

 Conceptos generales

S OLUCIONES DE UNA ECUACIÓN DIFERENCIAL

S OLUCIÓN DE UNA ECUACIÓN DIFERENCIAL

Es una función y = y(x), definida en un intervalo I = (a, b) que tiene al menos n

derivadas continuas en I y que satisface la EDO en I.
F (x, y(x), y ′ (x), ..., y (n (x)) = 0

S OLUCIÓN GENERAL DE UNA ECUACIÓN DIFERENCIAL

Es una familia de funciones dependientes de n parámetros (constantes), de la cual

pueda extraerse la EDO por derivación y eliminación de los parámetros.

S OLUCIÓN PARTICULAR DE UNA ECUACIÓN DIFERENCIAL

Es cada una de las funciones que se obtiene de la solución general al dar valores a
los parámetros.

S OLUCIÓN SINGULAR DE UNA ECUACIÓN DIFERENCIAL

Es una solución que no puede extraerse de la solución general pero sí satisface la

EDO en I.

Cálculo II EDOs de primer orden 5/29

 Conceptos generales

E JEMPLO

Sea la ecuación diferencial: y ′ = 2

Se obtiene la solución general integrando:

dy
Z Z
y′ = 2 ⇒ = 2 ⇒ dy = 2dx ⇒ dy = 2dx
dx
La solución general es:
y = 2x + C

que es una familia uniparamétricas de rectas de pendiente 2.

Las soluciones particulares se hallan

dando valores a la constante C. Se
puede hallar la constante asociada a
la recta que pasa por un punto P .

Cálculo II EDOs de primer orden 6/29

 Conceptos generales

C ONDICIONES INICIALES PARA LA ECUACIÓN DIFERENCIAL

Es un conjunto de n datos que acompañan a la ecuación, formado por los valores

que toma la función incógnita y sus derivadas en x = x0 :
y(x0 ) = y0 , y ′ (x0 ) = y1 , y ′′ (x0 ) = y2 , ..., y (n−1 (x0 ) = yn−1

P ROBLEMA DE VALOR INICIAL DE ORDEN n

Es una ecuación diferencial de orden n acompañada de n condiciones iniciales:
 ′ (n
F (x, y(x), y (x), ..., y (x)) = 0
y(x0 ) = y0 , y ′ (x0 ) = y1 , y ′′ (x0 ) = y2 , ..., y (n−1 (x0 ) = yn−1

E CUACIÓN DIFERENCIAL LINEAL

Es toda ecuación lineal en las variables y, y ′ , y (n tomando la forma:

an (x)y (n (x) + ... + a1 (x)y ′ (x) + a0 (x)y(x) = b(x)
En una ecuación diferencial lineal la variable dependiente y sus derivadas sólo aparecen
elevadas a grado uno y además, cada coeficiente sólo depende de la variable independiente.

Cálculo II EDOs de primer orden 7/29

 EDOs de primer orden

E CUACIONES DIFERENCIALES ORDINARIAS DE PRIMER ORDEN

Forma implícita: F (x, y, y ′ ) = 0 Ejemplo: (y′ )2 − 2y′ + 4y − x = 0

Si la y ′ se puede despejar en la ecuación implícita, se obtiene la

forma explícita: y ′ = f (x, y) Ejemplo: y′ = y+2
1−x

dy
Forma diferencial, teniendo en cuenta que y ′ = dx :
M (x, y)dx + N (x, y)dy = 0 Ejemplo: (x − y)dx + (y2 − x)dy = 0
2

La solución general de una EDO de primer orden es una familia uniparamétrica de

funciones: y = ϕ(x, C), donde C es un parámetro real. Para cada valor de C se
obtendrá una solución particular.
 problema ′de valor inicial de primer
Un  orden se escribe de la forma:
F (x, y, y ) = 0 y ′ = f (x, y)
y(x0 ) = y0 y(x0 ) = y0
y se resuelve (se integra) encontrando la función y = ϕ(x, C) que verifica
y(x0 ) = y0

Cálculo II EDOs de primer orden 8/29

 EDOs de primer orden

S OLUCIÓN DE UNA EDO DE PRIMER ORDEN

Geométricamente, la familia solución general es el conjunto de curvas en las que

se cumple que la pendiente en un punto es el valor de la función f en ese punto.
De esas curvas, la que pasa por (x0 , y0 ) es la solución del problema de valor
inicial.
T EOREMA DE EXISTENCIA Y UNICIDAD

Dada la ecuación y ′ = f (x, y) con f definida en D = [a, b] × [c, d], a < x0 < b,
c < y0 < d, si f y fy′ son continuas en D, entonces existe una solución única de
la ecuación cumpliento que y(x0 ) = y0 . Esta solución única del problema de valor
inicial está definida en un intervalo centrado en x0 .

C AMBIO DE VARIABLE

En las EDOs de primer orden es frecuente utilizar un cambio de variable para

reducir la dificultad de la ecuación.

Cálculo II EDOs de primer orden 9/29

 EDOs de primer orden

E JERCICIO P ROPUESTO 2

Representación de las familias de curvas uniparamétricas (azul) y

aquella que pasa por el punto P (1, 2) (rojo).

Cálculo II EDOs de primer orden 10/29

 Métodos elementales de resolución de EDOs Ecuaciones de variables separables

E CUACIONES DE VARIABLES SEPARADAS

Una EDO de primer orden es de variables separables o separadas si se puede

escribir de la forma:
f1 (x)g1 (y)dx = f2 (x)g2 (y)dy
Las funciones que satisfacen esta ecuación deben cumplir alguna de las siguientes
ecuaciones:
f1 (x) g2 (y)
dx = dy f2 (x) = 0 g1 (y) = 0
f2 (x) g1 (y)

R ESOLUCIÓN
La solución general se obtiene por integración de la primera de las ecuaciones
anteriores: Z Z
f1 (x) g2 (y)
dx = dy
f2 (x) g1 (y)

y deberá analizarse si f2 (x) y g1 (y) = 0 proporcionan soluciones singulares.

Cálculo II EDOs de primer orden 11/29

 Métodos elementales de resolución de EDOs Ecuaciones homogéneas

E CUACIONES HOMOGÉNEAS

Una función f (x, y) es homogénea de grado n en sus argumentos si cumple:

f (tx, ty) = tn f (x, y)

Ejemplo 1:
f (x, y) = x2 y − 3xy 2 + 5y 3

f (tx, ty) = (tx)2 ty − 3tx(ty)2 + 5(ty)3 =

= t x y − 3t3 xy 2 + 5t3 y 3 = t3 (x2 y − 3xy 2 + 5y 3 ) = t3 f (x, y)
3 2

f(x,y) es homogénea de grado 3

Ejemplo 2:
x2 y − 3xy 2 + 5y 3
f (x, y) =
x3 − yx2
t (x y − 3xy + 5y 3 )
3 2 2
t3 (x2 y − 3xy 2 + 5y 3 )
f (tx, ty) = = =
t3 x3 − tyt2 x2 t3 (x3 − yx2 )
= t0 f (x, y) = f (x, y)
f(x,y) es homogénea de grado 0

Cálculo II EDOs de primer orden 12/29

 Métodos elementales de resolución de EDOs Ecuaciones homogéneas

E CUACIONES HOMOGÉNEAS

Un ecuación expresada de la forma

y ′ = f (x, y)

es homogénea si f (x, y) es una función homogénea de grado cero en sus

argumentos.
R ESOLUCIÓN
Las EDOs homogéneas se resuelven reduciéndolas a ecuaciones de
variables separables mediante el cambio de variable:
y
z=
x

Cálculo II EDOs de primer orden 13/29

 Métodos elementales de resolución de EDOs Ecuaciones exactas

E CUACIONES EXACTAS

La EDO:
M (x, y)dx + N (x, y)dy = 0
es una ecuación exacta si existe una función u(x, y) tal que:

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

Por lo tanto: M (x, y)dx + N (x, y)dy = 0 ⇒ du(x, y) = 0 con lo cual la solución
general de la ecuación será la familia uniparamétrica u(x, y) = C

C RITERIO DE ECUACIÓN EXACTA ( O TEST DE EXACTITUD )

Supongamos M (x, y) y N (x, y) funciones con derivadas parciales continuas. La

ecuación M (x, y)dx + N (x, y)dy = 0 es exacta si y solo si:

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

Ésta es la condición que verifica un campo vectorial conservativo. Buscar la solución

de la EDO equivale a determinar la función potencial u(x, y).

Cálculo II EDOs de primer orden 14/29

 Métodos elementales de resolución de EDOs Ecuaciones exactas

E CUACIONES EXACTAS ( CONT.)

R ESOLUCIÓN
∂u ∂u
Si la EDO es exacta, existe una función u(x, y) tal que du = dx + dy
∂x ∂y
∂u(x, y)
u′x (x, y) La solución general es
)
= = M (x, y)
∂x u(x, y) = C puesto que
∂u(x, y)
u′y (x, y) = = N (x, y) du = 0
∂y

1 Integrar M (x, y) respecto de x:

R
u′x (x, y) = M (x, y) ⇒ u(x, y) = M (x, y)dx + C(y)
2 Derivar u(x, y) respecto de y e igualar a N (x, y) para despejar C ′ (y):
u′y (x, y) = N (x, y) ⇒ C ′ (y) = ...
3 Integrar C ′ (y) respecto de y.

Este proceso es equivalente si comenzamos integrando N respecto de y para

luego derivar el resultado respecto de x.

Cálculo II EDOs de primer orden 15/29

 Métodos elementales de resolución de EDOs Ecuaciones lineales

E CUACIONES LINEALES

Una EDO lineal de primer orden (en y) es de la forma:

y ′ + p(x)y = q(x) Forma canónica

donde p(x), q(x) son funciones continuas en la región en la que se integra.

Si q(x) = 0 la ecuación es lineal homogénea.
Si q(x) ̸= 0 la ecuación es completa.

T EOREMA : MÉTODO DE SUPERPOSICIÓN

Si se denota por yp cualquier solución particular de la ecuación completa

′
y + p(x)y = q(x) y por yh la solución general de la homogénea asociada
y ′ + p(x)y = 0, entonces la solución general de la ecuación completa es:
yg = yp + yh

Por lo tanto, para resolver una ecuación lineal basta con encontrar la solución
general de su homogénea asociada y una solución particular de la completa.

Cálculo II EDOs de primer orden 16/29

 Métodos elementales de resolución de EDOs Ecuaciones lineales

E CUACIONES LINEALES ( CONT.)

R ESOLUCIÓN DE LA HOMOGÉNEA

Todas las EDOs lineales homogéneas son de variables separables:

y ′ = −p(x)y ⇒ y = Ceg(x) ó y = 0
R
donde g(x) = − p(x)dx

S OLUCIÓN PARTICULAR DE LA COMPLETA

Se puede utilizar el Método de variación de constantes, que consiste en buscar

una solución particular convirtiendo en función indeterminada la constante de la
familia uniparamétrica de soluciones de la homogénea asociada:

yh = Ceg(x) ⇒ yp = C(x)eg(x)

Para obtener C(x) se deriva y se impone que se verifique la ecuación diferencial

completa.

Cálculo II EDOs de primer orden 17/29

 Métodos elementales de resolución de EDOs Ecuaciones lineales

E CUACIONES LINEALES ( CONT.)

S OLUCIÓN PARTICULAR DE LA COMPLETA ( CONT.)

Para obtener C(x) se deriva y se impone que se verifique la ecuación diferencial

sustituyendo en la EDO completa: )
yp = C(x)eg(x)
⇒ yp′ + p(x)yp = C ′ (x)eg(x) = q(x)
yp′ = C ′ (x)eg(x) − C(x)eg(x) p(x)

Se despeja C ′ (x) y se obtiene C(x) por integración:

R
C ′ (x) = q(x)e−g(x) ⇒ C(x) = q(x)e−g(x) dx
R
Por lo tanto la solución general es: yg = yh + yp = Ceg(x) + eg(x) q(x)e−g(x) dx

Sacando factor común y agrupando la constante de integración con la que se

obtendrá al al calcular C(x), resulta:
R Z
− p(x)dx
yg = e q(x)ep(x)dx dx

Cálculo II EDOs de primer orden 18/29

 Métodos elementales de resolución de EDOs Ecuaciones lineales

E JEMPLO

Resolver la ecuación: y ′ = x + y
y′ = x + y ⇒ y′ − y = x Es la EDO completa.
EDO lineal, su solución general será de la forma: yg (x) = yh + yp
Solución de la EDO homogénea o incompleta:
dy dy
y′ − y = 0 ⇒ −y =0⇒ = dx EDO de variables separadas
R dy R dx y
y
= dx ⇒ log |y| = x + log C
Tomando exponenciales, la solución de la homogénea es:
yh (x) = ex+log C
= ex elog C
= Cex
Solución particular de la completa: aplicamos el método de variación de la(s)
constante(s). Para ello se toma como solución particular una de la forma de la
yh pero la constante varía con x, C(x):
yp (x) = C(x)ex ⇒ yp′ (x) = C ′ (x)ex + C(x)ex
Se sustituye en la ecuación completa: yp′ − yp = x
C ′ (x)ex + C(x)ex − C(x)ex = x ⇒ C ′ (x)ex = x

Cálculo II EDOs de primer orden 19/29

 Métodos elementales de resolución de EDOs Ecuaciones lineales

E JEMPLO (C ONT.)

Solución particular de la completa (cont.)

x
C ′ (x)ex = x ⇒ C ′ (x) = = xe−x
ex
Se resuelve integrando (por partes):
Z
→
n u=x du = dx
o
C(x) = xe−x dx = −x
=
dv = e dx → v = −e−x
Z
= −xe−x + e−x dx = −xe−x − e−x = −e−x (x + 1)

La solución particular de la completa:

yp (x) = C(x)ex = −e−x (x + 1)ex = −(x + 1)

Solución general de la completa:

yg (x) = yh + yp = Cex − (x + 1)

Cálculo II EDOs de primer orden 20/29

 Modelado de procesos físicos y aplicaciones

A PLICACIONES

Dado un problema físico, habrá que decidir cuál es la función incógnita (y), cuál es
la variable independiente (x) y expresar todas las magnitudes que intervienen en
función de x, y, y ′ .
1. D ISOLUCIÓN EN UN RECIPIENTE

Un recipiente con volumen inicial V0 litros consta de una disolución líquida o

gaseosa cuya cantidad inicial de soluto es S0 . En t = 0 se vierte disolución con
concentración inicial Ce (kg/l) a razón/velocidad de Qe (l/min) y se deja salir
disolución a razón de Qs (l/min). La incógnita es la variación de la masa de soluto
en función del tiempo (S(t)).
Balance de masa: masa acumulada = masa entrante - masa saliente
Al cabo de un tiempo ∆t el incremento de masa de soluto (S(t)) para un volumen
V (t) es:
∆S = S(t + ∆t) − S(t) = Ce Qe ∆t − VS(t) (t)
Qs ∆t
Si dividimos por ∆t y hacemos que ∆t → 0, resulta la EDO:
dS S
= Ce Qe − Qs
dt V

Cálculo II EDOs de primer orden 21/29

 Modelado de procesos físicos y aplicaciones

A PLICACIONES

1. D ISOLUCIÓN EN UN RECIPIENTE ( CONT.)

dS S
= Ce Qe − Qs
dt V
tasa acumulación = tasa de entrada - tasa de salida
dS Qs Qs
Qe = Qs ⇒ V (t) = V0 ⇒ = Ce Qe − S donde V0
= cte
dt V0
)
dS Qs
dt
= Ce Qe − V0
S
Problema de valor inicial (EDO de var. separadas)
S(t = 0) = S0
S(t)
Concentración de la disolución: C(t) = V0

Qe ̸= Qs ⇒ V (t) = V0 + (Qe − Qs )t depende del tiempo

Si Qe > Qs ⇒ se llega a desbordar el depósito
Si Qe < Qs ⇒ se llega a vaciar el depósito
)
dS Qs
dt
+ V0 +(Qe −Qs )t
S = Ce Qe
Problema de valor inicial (EDO lineal)
S(t = 0) = S0
S(t)
Concentración de la disolución: C(t) = V0 +(Qe −Qs )t

Cálculo II EDOs de primer orden 22/29

 Modelado de procesos físicos y aplicaciones

A PLICACIONES

2. R AZÓN DE CAMBIO DE UNA VARIABLE PROPORCIONAL AL VALOR DE LA VARIABLE

El ritmo de cambio de una población, la velocidad a la que una sustancia se

transforma en otra en una reacción química simple, los procesos de desintegración
radiactiva o el ritmo de cambio de la velocidad de un móvil debido a fuerzas de
rozamiento.
dy 
= ±ky
dt Problema de valor inicial
y(t = 0) = y0
k mide de la rapidez con la que se produce el cambio en la magnitud estudiada.

3. C AMBIO DE TEMPERATURA DE UN CUERPO

Experimentalmente se ha comprobado que la temperatura de un cuerpo cambia

proporcionalmente a la diferencia entre su temperatura y la del medio en que se
encuentra (ley de enfriamiento de Newton). Si en t = 0 la temperatura del cuerpo
es T (0) = T0 y la temperatura del medio es Tm :
dT 
= −k(T − Tm )
dt Problema de valor inicial
T (t = 0) = T0

Cálculo II EDOs de primer orden 23/29

 Resolución numérica Método de Euler básico

M ÉTODO DE E ULER BÁSICO

El método de las tangentes o método de Euler nos da una primera aproximación

de de los valores que toma la solución sobre un conjunto de puntos dado.

A PROXIMACIÓN LINEAL

Sea y = f (x) una función derivable en el punto x0 . La recta tangente en el punto

(x0 , f (x0 )) es una linealización (aproximación lineal) de la función en ese punto.

Dado que la ecuación de la recta

tangente en el punto (x0 , f (x0 )) tiene
por pendiente f ′ (x0 ), su ecuación es:
y − f (x0 ) = f ′ (x0 )(x − x0 )

Cálculo II EDOs de primer orden 24/29

 Resolución numérica Método de Euler básico

M ÉTODO DE E ULER BÁSICO

Supongamos que el problema de valor inicial que tiene solución única y = y(x) en
un intervalo centrado en x0 :
( La ecuación de recta tangente a la solución y = y(x)
y′ = f (x, y)
en el punto (x0 , y0 ), con pendiente igual a f (x0 , y0 ):
y(x0 ) = y0
y = y0 + (x − x0 )f (x0 , y0 )

Para el punto x1 = x0 + h el valor de la recta

tangente resulta
y1 = y0 + hf (x0 , y0 )
Nos situamos ahora en el punto (x1 , y1 ), con-
sideramos la recta con pendiente dada por
f (x1 , y1 ) y la evaluamos en
x2 = x1 + h, para obtener
y2 = y1 + hf (x1 , y1 )

n y = yn + hf (xn , yn )
n+1
n = 0, 1, 2, ...
xn+1 = xn + h

Cálculo II EDOs de primer orden 25/29

 Resolución numérica Método de Euler mejorado

M ÉTODO DE E ULER MEJORADO

El método de Euler mejorado es una

aproximación más fina pues toma el
valor medio de los valores de la fun-
ción en los puntos extremos, es decir,
tomar la media de las pendientes en
los puntos extremos, halladas con el
método de Euler básico.

h h
h i h i
y1 = y0 + f (x0 , y0 )+f (x0 +h, y0 +hf (x0 , y0 )) = y0 + f (x0 , y0 )+f (x1 , y1 )
2 2

hh i
yn+1 = yn + f (xn , yn ) + f (xn + h, yn + hf (xn , yn )) =
2
hh i
= yn + f (xn , yn ) + f (xn+1 , yn+1 ) n = 0, 1, 2, ...
2

Cálculo II EDOs de primer orden 26/29

 Resolución numérica Métodos de Euler

M ÉTODOS DE E ULER

Cálculo II EDOs de primer orden 27/29

 Resolución numérica Método de Runge-Kutta

M ÉTODO DE RUNGE -K UTTA

Se recurre a una media ponderada de los valores de f (x, y) en distintos puntos. La

elección de esos puntos y de los coeficientes de la combinación genera una gran
familia de métodos.
M ÉTODO DE RUNGE -K UTTA DE ORDEN 4
( h
yn+1 = yn + (K1 + 2K2 + 2K3 + K4 )
6
xn+1 = xn + h

K1 = f (xn , yn )
K2 = f (xn + h2 , yn + h2 K1 )
K3 = f (xn + h2 , yn + h2 K2 )
K4 = f (xn + h, yn + hK3 )

Cálculo II EDOs de primer orden 28/29

 Otros métodos aproximados

OTROS MÉTODOS APROXIMADOS

C AMPO DE DIRECCIONES

La ecuación y ′ = f (x, y) fija, en cada punto del plano donde f está definida, el
valor de la pendiente de la curva y = y(x).
La ecuación genera un campo de direcciones; la representación gráfica de este
campo mediante pequeños segmentos situados en una muestra de puntos del
plano es lo que se suele llamar campo de direcciones de la ecuación diferencial,
pues proporciona un ‘mapa’ aproximado de las soluciones.

I SOCLINAS
Es el lugar geométrico de los puntos del plano donde la pendiente de las curvas
solución es constante, siendo su ecuación f (x, y) = C.
Las pendientes de las soluciones en los puntos de corte con la isoclina son, por
tanto, y ′ = C.

Cálculo II EDOs de primer orden 29/29