Capı́tulo 2

Solución de las ecuaciones

diferenciales de primer orden

En el anterior capı́tulo hemos visto ya como solucionar de forma muy simple las
ecuaciones diferenciales de primer orden que son de variables separables. En este capı́tulo
aprenderemos cómo solucionar otras ecuaciones de primer orden: trataremos primero
las ecuaciones lineales, la cuales responden siempre al patrón y′ + p(x)y = g(x) y para
las cuales existe un método genérico para encontrar su solución de forma analı́tica. A
continuación nos adentraremos en el mundo de las ecuaciones no lineales, para las cuales
no siempre es posible encontrar la solución. Revisaremos los casos más comunes que uno
se puede encontrar: las ecuaciones homogéneas, que responden al patrón y′ = f (y/x), las
ecuaciones exactas o que se pueden transformar en exactas, y finalmente trabajaremos el
caso de las ecuaciones de Bernoulli [y′ + p(x)y = g(x)yn ] y Ricatti [y′ = q1 (x) + q2 (x)y +
q3 (x)y2 ], que aunque no lineales, todas ellas son ecuaciones diferenciales ordinarias de
primer orden. En la penúltima sección del capı́tulo, aprenderemos como mediante la
ayuda de Sage podemos dibujar el ası́ llamado campo de pendientes con el que es posible
conocer de forma cualitativa como deben comportarse las soluciones de una ecuación.
En esta última sección aprenderemos también a utilizar Sage en la resolución algebraica
de ecuaciones y qué se puede hacer cuando algebraicamente no es posible encontrar
una solución: como veremos, no se puede ir al tuntún cuando se usan ordenadores para
resolver numéricamente ecuaciones diferenciales, los conocimientos teóricos son muy
importantes. En la última sección se hallan diversos ejercicios resueltos para ayudar a la
comprensión de lo explicado durante el capı́tulo.

2.1. Ecuaciones diferenciales de primer orden lineales

La ecuación más general de este tipo es:

y′ + p(x)y = g(x), (2.1)

siendo p(x) y g(x) funciones dadas. En este caso, vamos a imponer ya la condición inicial
arbitraria y(x0 ) = y0 y vamos a resolver la ecuación con la condición inicial a la misma
32 Solución de las EDO de primer orden

vez. La técnica es bien conocida desde hace siglos (Newton y Leibniz ya sabı́an hacerlo).
Se trata de multiplicar a derecha e izquierda por
Z x 
µ (x) = exp ds p(s) . (2.2)
x0

⋆ Problema - 2.1. La ecuación resultante se puede escribir como:

y′ µ + y µ ′ = g µ , (2.3)

o
(yµ )′ = gµ , (2.4)
Si ahora integramos entre x0 y x y usamos que µ (x0 ) = 1, y(x0 ) = y0 , llegamos a la
solución: Z x 
1
y(x) = ds g(s)µ (s) + y0 (2.5)
µ (x) x0

que resuelve el problema de la ecuación diferencial y la condición inicial. Vemos que no

hay ninguna arbitrariedad en ningún paso y, por tanto, la solución encontrada es la única
que existe.

⋆ Problema - 2.2. Caso práctico. Sea el problema de condiciones iniciales:

y′ = −xy + 4x,

(2.6)
y(1) = 2.

Calculamos primero µ (x) = exp((x2 − 1)/2). La (única) solución es entonces

y(x) = 4 − 2 exp((1 − x2 )/2). (2.7)

⋆ Problema - 2.3. Encontrad la familia de soluciones de

y′ = −xy + 4x (2.8)

dependientes de un parámetro C. Imponed entonces la condición inicial y(1) = 2 y

recuperad la solución anterior.

2.2. Ecuaciones homogéneas

Hay que ir con cuidado porque la palabra “homogénea” quiere decir muchas cosas en
matemáticas (y también en el campo de las ecuaciones diferenciales). Hoy, “homogénea”
quiere decir que la ecuación diferencial tiene la forma:

y′ = f (y/x). (2.9)
2.3 Ecuaciones exactas 33

Por ejemplo y′ = sin(y/x). No siempre es tan fácil darse cuenta a primera vista de si una
ecuación diferencial es homogénea.

⋆ Problema - 2.4. Demostrar que

′ y2 + 2xy
y = 2 (2.10)
x + y2
es homogénea y encontrar la función f .

Las ecuaciones homogéneas se resuelven mediante el cambio de variables (que parece

hasta obvio) z = y/x.

⋆ Problema - 2.5. Esto lleva a:

f (z) − z
z′ = (2.11)
x
que es separable:
dz dx
= , (2.12)
f (z) − z x
y que se puede integrar sin más (bueno, hay que saber integrar, pero eso es cosa de
otro curso). Luego se deshace el cambio de variable.

⋆ Problema - 2.6. Demostrar que la familia de soluciones de la ecuación diferencial:

y2 + 2xy Cx2
y′ = es y(x) = .
x2 1 −Cx

2.3. Ecuaciones exactas

Puede suceder que nuestra ecuación no sea lineal, y no tengamos la suerte de que
sea homogénea. En estos casos, es conveniente comprobar si estamos ante una ecuación
diferencial exacta1 o que se pueda convertir fácilmente en exacta, dado que en dichos
casos es fácil encontrar la solución al problema general.
Una ecuación diferencial de primer orden, y′ = f (x, y), también se puede escribir
como
N(x, y) y′ + M(x, y) = 0, (2.13)
1
El nombre de ecuaciones diferenciales exactas, proviene de que dada una función F(x, y) llamamos
la diferencial exacta (o diferencial total) de la función F a dF = ∂x F dx + ∂y F dy, y por analogı́a
cuando tenemos una expresión cualquiera del estilo h(x, y) dx + g(x, y) dy, diremos que esta expresión
es o corresponde a una diferencial exacta si existe una función F tal que dicha expresión sea su dF. Es
decir, si existe una función F(x, y) tal que ∂x F = h(x, y) a la vez que ∂y F = g(x, y) lo cual nos permita
escribir dF = h(x, y) dx + g(x, y) dy.
34 Solución de las EDO de primer orden

dy
o en una forma más mnemotécnica, aprovechando que y′ = también la podemos
dx
−2y
escribir como N(x, y) dy + M(x, y) dx = 0. Ası́ por ejemplo la ecuación y′ = se
3x
puede escribir en la forma (2.13) de varias maneras distintas:

3xy′ + 2y = 0, (2.14)
3x2 y2 y′ + 2xy3 = 0, (2.15)
3 ′ 2
y+ = 0. (2.16)
y x
Como probaremos más adelante, si conseguimos escribir cualquier ecuación diferencial
en la forma (2.13) de manera que se verifique que

∂N ∂M
= , (2.17)
∂x ∂y
entonces será una ecuación diferencial exacta.

⋆ Problema - 2.7. ¿ Resulta ser alguna de las formas anteriores de escribir la ecuación
−2y
y′ = una ecuación diferencial exacta?
3x

Como veremos, si conseguimos escribir una ecuación en una forma exacta, estare-
mos muy cerca de encontrar la familia de soluciones: tomemos del ejemplo anterior la
expresión 3x2 y2 y′ + 2xy3 = 0, observamos que N(x, y) = 3x2 y2 , M(x, y) = 2xy3 , y que
∂x N = ∂y M = 6xy2 por lo que sin duda estamos ante un caso de una ecuación de las que
hemos clasificado como exactas. Como probaremos más adelante, para cualquier ecuación
exacta existe una función Ψ(x, y), que en el caso que nos ocupa es Ψ(x, y) = x2 y3 , tal
que ∂∂Ψy = N y ∂∂Ψx = M, y por tanto la ecuación diferencial la puedo reescribir como2 :

dy
∂y Ψ + ∂x Ψ = 0.
dx
Fácilmente reconoceremos que el termino izquierdo no es más que el resultado de aplicar
la regla de la cadena cuando calculamos el diferencial total de Ψ(x, y) respecto a x:

dΨ(x, y)
=0
dx
y esta última ecuación si que la sabemos resolver muy fácilmente, puesto que necesari-
amente Ψ(x, y) = C donde C es una constante. Ası́ pues, Ψ(x, y) = C es una solución
implı́cita a nuestra ecuación diferencial. Aplicado a nuestro ejemplo, obtenemos x2 y3 = C,
y por tanto y(x) = (C/x2 )1/3 . Como vemos, en algunos casos, será fácil aislar y de la
expresión Ψ(x, y) = C y podremos obtener una solución explicita y(x) a la ecuación
diferencial, en otros casos, al no ser posible aislar y, [Link]. si Ψ(x, y) = C resulta ser
2
QUEDA PENDENT EXPLICAR
2.3 Ecuaciones exactas 35

una ecuación trascendente, entonces nos deberemos conformar con la solución implı́cita
Ψ(x, y) = C. En estos casos, para una C fijada, si uno desea obtener y(x), se puede ir
asignando valores a x y resolver numéricamente cuanto debe valer y para que se satisfaga
Ψ(x, y) = C.
Probemos ahora de manera rigurosa todo lo que hemos afirmado anteriormente,
empezando con el siguiente teorema que nos garantiza que si tenemos una ecuación
diferencial exacta, entonces, sin duda, existe una función Ψ como la que nosotros de-
seamos tener. Durante el proceso de demostración hallaremos una forma sistemática de
obtener la función Ψ(x, y).
∂M ∂N
Teorema 2.1. La ecuación diferencial N(x, y)dy + M(x, y)dx = 0, donde M, N, ,
∂y ∂x
∂N ∂M
son funciones continuas en una cierta región, es exacta, es decir = , si y sólo
∂x ∂y
si3 existe una función Ψ(x, y) tal que:
∂ Ψ(x, y) ∂ Ψ(x, y)
= M(x, y), = N(x, y). (2.18)
∂x ∂y
La demostración del teorema en el sentido ”⇒” es como sigue:
Demostración (⇒): Buscamos una función Ψ(x, y) que satisfaga ∂x Ψ = M y ∂y Ψ =
∂Ψ
N. Tomemos la primera igualdad = M e integremos con respecto a x (manteniendo
∂x
y constante): Z x
Ψ= dx M(x, y) + H(y) (2.19)

siendo H(y) una función arbitraria. Determinaremos esta función pidiendo que se cumpla
∂Ψ
la segunda condición = N. Esto lleva a
∂y
Z x
∂Ψ ∂ M dH(y)
N= = dx + (2.20)
∂y ∂y dy
de donde encontramos
Z x Z y  Z x 
dH(y) ∂M ∂M
=N− dx ⇒ H(y) = dy N − dx . (2.21)
dy ∂y ∂y
Para que eso tenga sentido, el lado derecho de esta última ecuación tiene que ser función
sólo de y. Es decir, necesitamos que
Z y  Z x 
∂ ∂M
dy N − dx = 0,
∂x ∂y
lo cual se puede probar usando el hecho de que tenemos una ecuación diferencial exacta,
∂N ∂M
es decir, = .
∂x ∂y
3
Observar el ”si y sólo si” que implica que si es cierta la premisa se verifica la conclusión pero que
también debe suceder al revés, si se verifica la conclusión debe ser cierta la premisa.
36 Solución de las EDO de primer orden

⋆ Problema - 2.8. Probar este último paso de la primera parte de la demostración

del teorema.
Ası́ pues Ψ(x, y) vendrá dado por:
Z x Z y  Z x 
∂M
Ψ= dx M + dy N − dx (2.22)
∂y
y por construcción Ψ satisface ∂x Ψ = M y ∂y Ψ = N.
Llegados a este punto, hemos probado la primera parte del teorema, y además ahora
sabemos como encontrar de forma sistemática la función Ψ(x, y) lo cual nos será muy
útil una vez hayamos identificado una ecuación diferencial como exacta.
−2y
⋆ Problema - 2.9. Para la versión de y′ = que es exacta, encontrar la función
3x
Ψ(x, y).

La segunda parte de la demostración (en el sentido ⇐)es si cabe más simple:

Demostración (⇐): Imaginemos que existe una función Ψ(x, y) que satisface las
∂ Ψ(x, y) ∂ Ψ(x, y)
condiciones = M(x, y), y = N(x, y) . Entonces:
∂x ∂y
∂ 2Ψ ∗ ∂ 2Ψ
   
∂M ∂ ∂Ψ ∂ ∂Ψ ∂N
= = = = = (2.23)
∂y ∂y ∂x ∂ y∂ x ∂ x∂ y ∂ x ∂ y ∂x

⋆ Problema - 2.10. Justificar la igualdad marcada con *

−2y
⋆ Problema - 2.11. Para y′ = demostrar sin cambiar su forma que también
3x
Ψ(x, y) = x2 y3 y, por tanto, la solución es y(x) = Cx−2/3 . Este ejemplo se podı́a resolver
dándose cuenta de que la ecuación diferencial es de variables separadas. Comparad
los dos métodos.
Vemos como moraleja, que una misma ecuación diferencial puede, a veces, resolverse
de diversas maneras. Hay que estar pues atentos, y escoger aquella vı́a que nos resulte
menos tediosa.

2.4. Factor integrante

Si una ecuación diferencial escrita en la forma N(x, y)y′ + M(x, y) = 0 no es exac-
ta, puede ocurrir que exista una función µ (x, y) tal que al multiplicar la ecuación por
µ (x, y) se convierta en exacta. Esta función o factor µ (x, y) recibe el nombre de factor
integrante4
4
En principio una ecuación diferencial puede tener más de un factor integrante, es decir, para una
misma ecuación diferencial puede existir diversas funciones µ(x, y) tales que la conviertan en exacta. Se
puede probar, que únicamente en el caso de que µ dependa de una sola de las variables (µ = µ(x) o bien
µ = µ(y)), queda garantizado que no existe ningún otro factor integrante el cual dependa únicamente
de esa misma variable.
2.5 Ecuación de Bernoulli 37

⋆ Problema - 2.12. Demostrar que −x2 y′ + y2 + xy = 0 se convierte en exacta al

multiplicarla por µ (x, y) = xy12 . Resolver la ecuación.

El problema es pues ¿Cómo encontrar un factor integrante? Se debe verificar que

µ (x, y)N(x, y)y′ + µ (x, y)M(x, y) = 0 sea exacta, o sea que:

∂ (µ N) ∂ (µ M)
= (2.24)
∂x ∂y

Por consiguiente, lo que hay que hacer es resolver la anterior ecuación diferencial en
derivadas parciales para la función µ (x, y) ¡lo cual suele ser mucho más difı́cil que re-
solver la ecuación original! Sin embargo, hay un par de casos en los que sı́ se puede
encontrar fácilmente µ (x, y).

⋆ Problema - 2.13. Imaginemos que existan soluciones donde la función µ (x) no

dependa de y. En ese caso, se llega a:
 
1 ∂µ −1 ∂ M ∂N
=N − (2.25)
µ ∂x ∂y ∂x

Esta ecuación tendrá solución si el lado derecho depende sólo de la variable x.

⋆ Problema - 2.14. Encontrar un factor integrante para la ecuación (3xy + y2 )dx +

(x2 + xy)dy = 0 y resolver la ecuación. ¿Se os ocurre otra manera de resolverla?

⋆ Problema - 2.15. Encontrar las condiciones para que el factor integrante µ dependa
sólo de y.

2.5. Ecuación de Bernoulli

Tiene la forma
y′ + p(x)y = g(x)yn , (2.26)
con n 6= 1. Para n = 1 es fácil encontrar la solución,

⋆ Problema - 2.16. Resolver la ecuación de Bernoulli si n = 1.

Para otros valores de n, ya en su dı́a Bernoulli nos enseñó el camino a seguir: hacer
el cambio de variables z = y1−n .
38 Solución de las EDO de primer orden

⋆ Problema - 2.17. Esto lleva a la ecuación lineal:

z′ + p̂(x)z = ĝ(x), (2.27)

con p̂(x) = (1 − n)p(x), ĝ(x) = (1 − n)g(x), que se resuelve por la técnica del capı́tulo
2.1 para encontrar z(x). Luego se deshace el cambio de variable para encontrar y(x).

2.6. Ecuación de Ricatti

Es no lineal de la forma:

y′ = q1 (x) + q2 (x)y + q3 (x)y2 (2.28)

No existe una manera general de encontrar la solución de una ecuación de Ricatti váli-
da para funciones arbitrarias q1 (x), q2 (x), q3 (x). Sin embargo, es posible hacer las dos
afirmaciones siguientes:
1) Si conocemos una solución cualquiera de la ecuación de Ricatti y1 (x) entonces
la familia general de soluciones se encuentra mediante el cambio de variables y(x) =
1
y1 (x) + v(x) .

⋆ Problema - 2.18. Demostrar que v(x) satisface una ecuación lineal que se resolverá.
Aplicar este método para resolver el problema:
y
y′ = x2 + − y2 , y(x = 1) = 0 (2.29)
x
dándose cuanta de que y1 (x) = x es una solución de esta ecuación.

2) Una ecuación de Ricatti se puede convertir en una ecuación lineal de segundo

z′
orden, haciendo los cambios de variables a una función z(x) que satisface = −q3 y.
z

⋆ Problema - 2.19. Encontrar la ecuación diferencial lineal para la función z(x)

asociada a la ecuación de Ricatti del caso anterior.

2.7. Uso de Sage en la resolución de ecuaciones difer-

enciales de primer orden
2.7.1. Como extraer información cualitativa a partir de un cam-
po de pendientes.
En el caso de las ecuaciones diferenciales de primer orden, una forma cualitativa de
ver las soluciones y su comportamiento es mediante la representación gráfica del campo
de pendientes (o campo de direcciones): todas las ecuaciones de primer orden responden
2.7 Uso de Sage en la resolución de ecuaciones diferenciales de primer orden
39

al patrón y′ = f (x, y) y por tanto la función f (x, y) nos da información sobre el valor de
la pendiente de la solución y(x) que pasa por el punto (x, y). Ası́ pues, si a cada punto
(x, y) del espacio le asociamos un pequeño segmento rectilı́neo con pendiente igual al
valor de f (x, y), dado que dichos segmentos representaran la tangente en el punto (x, y)
a la solución y(x) que pase por ese punto, podremos entonces fácilmente deducir como
se comportarán la soluciones a la ecuación diferencial.
Dibujar un campo de pendientes, es muy fácil con Sage, supongamos que queremos
conocer cualitativamente que solución o soluciones tendrá la ecuación diferencial de
primer orden5 y′ = −2(1 − y/3)y en función de la condición inicial que impongamos
y(xo ) = yo . Para ello introduciremos en la lı́nea de comandos de Sage:

sage: var(’x,y’)
sage: plot_slope_field(-2*(1-y/3)*y, (x,-4,5), (y,-2,5))

Del dibujo resultante, ver figura 2.1, observamos que si nuestra solución debe pasar por
el punto (−2, −1,5) (es decir y(−2) = −1,5) entonces, a base de seguir los segmentos a
partir de dicho punto que representan las lı́neas tangentes a la solución y(x) vemos que
nuestra solución tenderá asintóticamente hacia cero a medida que x crece, y tenderá a
−∞ cuando x se haga negativo. Por el contrario, si nuestra condición inicial es tal que la
solución debe pasar por el punto (−2, 1,5) (es decir, y(−2) = 1,5), vemos que nuestra
solución debe ser tal que asintóticamente tienda a cero cuando x → ∞ y debe tender a
y = 3 cuando x → −∞. A su vez, si la condición inicial es y(−2) = 4,5 vemos que la
solución al problema de condiciones iniciales será tal que cuando x → ∞ el valor de y se
dispara hacia +∞ mientras que cuando x → −∞ la solución debe tender asintóticamente
hacia y = 3. También podemos observar que si la condición es y(xo ) = 0 cualquiera que
sea xo vemos que la solución al problema es y(x) = 0, mientras que si la condición inicial
es y(xo ) = 3 vemos que la solución será y(x) = 3. Observamos pues que la solución a
un problema de condiciones iniciales puede ser muy diferente dependiendo de cual sea
la condición inicial impuesta, aun cuando la ecuación diferencial sea exactamente la
misma. Otro aspecto a considerar es el hecho de que según el teorema de existencia y
unicidad sabemos que dada una condición inicial las soluciones a este problema deben
ser únicas. No obstante en la figura 2.1 podemos observar que dependiendo de cual sea
la condición inicial, las soluciones pueden ser válidas en todo el intervalo de las x o no:
cuando yo ∈ [0, 3] vemos que la solución y(x) abarca todos los valores de x, es decir y(x)
esta definida para todo x ∈ (−∞, +∞). Cuando el valor de yo ∈ / [0, 3] observamos que
las soluciones presentan siempre una ası́ntota hacia +∞ o hacia −∞ para un valor de x
finito y por tanto existe todo un rango de valores de x para el cual no es posible definir
una solución que cumpla con la condición inicial y(xo ) = yo . Por ejemplo, de la figura 2.1
vemos claramente que si la solución debe satisfacer la condición inicial de pasar por el
punto (2, −2), debido a la ası́ntota hacia −∞ no nos podemos plantear por ejemplo cual
es valor de dicha solución para x = −3, simplemente dicha solución nunca pasará por
5
Para aquellos más ansiosos de saber, esta ecuación recibe el nombre de ecuación logı́stica, y veremos
en detalle como caracterizar este tipo de ecuaciones más adelante, aquı́ simplemente la usamos como
ejemplo para ilustrar toda la información cualitativa que se puede obtener mediante un campo de
pendientes.
40 Solución de las EDO de primer orden

un punto que tenga por valor de x = −3. El teorema de existencia y unicidad nos dice
que existe una única solución pero no nos puede garantizar que el intervalo de validez de
las soluciones sea ∀x. Como vemos, hemos obtenido una información nada desdeñable
sobre toda una familia de problemas de condiciones iniciales con sólo pintar el campo de
pendientes asociado a la ecuación diferencial.

Figura 2.1: Campo de pendientes para la ecuación y′ = −2(1 − y/3)y.

No obstante, hay que ir con cuidado, por ejemplo si nos concentramos en el problema
y′ = 2y1/2 , y dibujamos el campo de pendientes introduciendo
sage: var(’x,y’)
sage: plot_slope_field(2*yˆ(1/2), (x,-4,5), (y,-2,5))
obtendremos un dibujo como el de la figura 2.2, pero esto sólo es parte de la realidad
puesto que Sage sólo nos da la raı́z positiva de y1/2 , pero hay que considerar también la
raı́z negativa, es decir, realmente tenemos ±|y1/2 |. Esto nos debe llevar a pensar que si
en un mismo punto puede haber dos pendientes diferentes, la solución no será en general
única. De hecho, si recordamos el enunciado del teorema de existencia y unicidad, vemos
que aquı́ no se puede aplicar dicho teorema y por tanto no hay ninguna garantı́a de que
la solución sea única.
Para pintar el campo de pendientes completo correspondiente a f (x, y) = 2y1/2 ,
podemos superponer los gráficos para cada una de las raı́ces de la siguiente manera 6 :
sage: var(’x,y’)
sage: dibujo1=plot_slope_field(2*yˆ(1/2), (x,-4,5), (y,-2,5))
sage: dibujo2=plot_slope_field(-2*yˆ(1/2), (x,-4,5), (y,-2,5))
sage: mitexto=text(’Campo para $f(x,y)=\pm 2yˆ{1/2}$’,(0,-1))
sage: show(dibujo1+dibujo2+mitexto)
6
Como vemos en este ejemplo los gráficos o dibujos son objetos que también se pueden guardar en
variables, y posteriormente los podemos mostrar combinados con la instrucción show(). Observar que
la función para crear un objeto de texto se compone de dos partes, la cadena de texto que queremos
escribir entre comillas, y a continuación las coordenadas (x, y) en donde debe colocarse el centro de la
cadena de texto. También conviene observar para escribir fórmulas matemáticas en un texto, se puede
utilizar en Sage la sintaxis que se utiliza en LaTex. Para más detalles sobre LaTex consultar por ejemplo
[Link]
2.7 Uso de Sage en la resolución de ecuaciones diferenciales de primer orden
41

El resultado será similar al de la figura 2.3, y en estos casos no es posible obtener

demasiada información cualitativa del campo de pendientes: por ejemplo vemos que no
habrá soluciones reales que den valores negativos para y, pero serı́a un error concluir
del campo de pendientes que sólo existen dos soluciones para cualquier condición inicial
dada y(xo ) = yo . De hecho se puede ver que en muchos casos hay ¡infinitas soluciones
posibles!,

⋆ Problema - 2.20. Discutir el número de soluciones que admite y′ = 2y1/2 según

cual sea la condición inicial impuesta y(xo ) = yo . En particular probar que todas las
siguientes familias de funciones son de hecho soluciones posibles a nuestro problema
de condiciones iniciales:
 p
0p x < xo ± p(yo ),
y(x) =
(x − xo ± (yo ))2 x ≤ xo ± (yo )
(x − p)2

 x p < p,
y(x) = p0 p≤ p x ≤ xo − (yo ),
(xo − (yo ))2 xo − (yo ) < x


 p p
 (x − xo + (yo ))2 p x < xo + (yo ),
y(x) = 0 xo + (yo ) ≤ x ≤ p,
(x − p)2 p< x


Figura 2.2: Campo de pendientes (incompleto) para la ecuación y′ =.

2.7.2. Extraer soluciones analı́ticas de una ecuación de primer

orden
El programa Sage también nos permite obtener en muchos casos soluciones analı́ticas,
ya sea implı́citas o explı́citas, para ecuaciones diferenciales de primer y segundo grado me-
diante la función desolve(). Por ejemplo, si queremos obtener las solución a la ecuación
y′ = (x + y)/x, escribiremos:
42 Solución de las EDO de primer orden

Figura 2.3: Campo de pendientes (completo) para la ecuación y′ =.

sage: var(’x’)
sage: y=function(’y’,x)
sage: ecuacion=(diff(y,x)==(x+y)/x)
sage: desolve(ecuacion,y)

que producirá el output (c + log(x)) ∗ x, el cual se puede comprobar fácilmente que

verifica la ecuación. Observar que en este caso, dado que no conocemos la función y(x)
a priori, debemos definir y explı́citamente como una función de x mediante el comando
y = f unction(′ y′ , x). En algunos casos, como en el caso de las ecuaciones de Ricatti,
dentro de la función dsolve debemos añadir el parámetro contrib ode=True el cual
esta como opción por defecto desactivado, y sirve para decirle al programa que siga
adelante en el caso de ecuaciones diferenciales que potencialmente pueden ser muy
lentas de resolver. A modo de ejemplo, supongamos por ejemplo que queremos obtener
las soluciones de la siguiente ecuación de Ricatti: y′ = 1 + x2 − 2xy + y2 , entonces
escribiremos

sage: var(’x’)
sage: y=function(’y’,x)
sage: ecuacion=(diff(y,x)==1+xˆ2-2*x*y+yˆ2)
sage: desolve(ecuacion,y,contrib_ode=True)

el resultado será

[[x == c + 1/sqrt(t - 1), y(x) == x - sqrt(t - 1)],

[x == c - 1/sqrt(t - 1), y(x) == x + sqrt(t - 1)]]

el cual si se simplifica nos indica que la solución es: y(x) = x + 1/(c − x). La función
dsolve permite calcular la solución para condiciones de contorno especı́ficas entre otras
muchas cosas. Los parámetros asociados a dsolve se pueden consultar introduciendo en
la lı́nea de comandos dsolve?. Dicha función no es la única de que dispone Sage para la
resolución de ecuaciones diferenciales, el comando dsolvers? nos mostrará información
sobre las diversas funciones disponibles para este fin.
2.7 Uso de Sage en la resolución de ecuaciones diferenciales de primer orden
43

2.7.3. Extraer soluciones numéricas de una ecuación de primer

orden
Hemos visto que para las ecuaciones diferenciales de primer orden que son lineales
siempre es posible obtener la solución de forma algebraica. Pero en el caso de las ecua-
ciones diferenciales no lineales de primer orden, no siempre es posible. ¿Qué podemos
hacer en estos casos?
A parte de emplear el teorema de existencia y unicidad para saber si al menos tenemos
garantizado que existe una única solución, y a parte de estudiar el campo de pendientes
que nos puede dar información cualitativa sobre las soluciones, tenemos una tercera
herramienta que consiste en el cálculo numérico de las soluciones. Evidentemente, no
sólo se pueden aplicar los métodos numéricos a ecuaciones de primer orden, sino que
el uso de métodos numéricos se extiende a todo tipo de ecuaciones diferenciales, tanto
ordinarias como en derivadas parciales, sean o no lineales. Como es de imaginar, este
es un campo de conocimiento muy vasto que puede llevar perfectamente varios años el
dominarlo. Ası́ pues, aquı́ sólo ofreceremos una pincelada muy básica para abrir boca y
dejaremos el resto para cursos posteriores .
Dado el problema y′ = f (x, y), el método numérico más básico proviene de recordar
la definición de la derivada de una función
y(x + h) − y(x)
f (x, y(x)) = y′ ≈
h
que podemos escribir como
y(x + h) = y(x) + h f (x, y(x)).
Observemos que la anterior fórmula nos permite conocer el valor de la función en la
posición x + h a partir del conocimiento de la posición actual, x, y del valor de la función
en la posición actual, y(x). Ası́ pues, si conocemos la condición inicial, y(xo ) = yo ,
podemos calcular el valor de y para otras posiciones avanzando (o retrocediendo) a
pequeños intervalos h. Por ejemplo si tenemos y′ = ax procederı́amos como sigue:
Itera →
#0 y(xo ) (2.30)
Itera →
#1 y(xo + 1 h) = y(xo ) + h ax (2.31)
Itera →
#2 y(xo + 2 h) = y(xo + h) + h ax (2.32)
Itera →
#3 y(xo + 3 h) = y(xo + 2 h) + h ax (2.33)
.... (2.34)
Itera #n → y(xo + n h) = y(xo + (n − 1) h) + h ax (2.35)
y ası́ sucesivamente. Este método se conoce como el método de Euler , y evidentemente,
hacer esto a mano es aburrido y cansado, pero por suerte Sage puede hacerlo por nosotros
muy fácilmente mediante la función eulers method( f (x, y), xo , yo , h, xstop ). Supongamos
para ilustrarlo que queremos resolver el problema7 y′ = −2∗x hasta llegar a xstop = 0 para
7
Este problema corresponde al caso de un objeto que a tiempo ’cero’ se encuentra en la posición
100 y que esta frenando con una aceleración de −2, donde x representa el tiempo, y y la posición del
objeto.
44 Solución de las EDO de primer orden

la condición inicial y(0) = 100, como paso tomamos h = 0,1. Ası́ pues introduciremos
los siguientes comandos8 :

sage: var(’x’)
sage: y=function(’y’,x)
sage: f(x,y)=-2*x
sage: eulers_method(f,0,100,0.1,1)

lo cual producirá una tabla con los valores de x, y en cada iteración y una tercera
columna que nos da el resultado de h f (x, y) para cada iteración. Si ahora queremos ver
el resultado no en una tabla sino en un gráfico, nada más sencillo, el último comando lo
reemplazamos por

sage: listadepuntos=eulers_method(f,0,100,0.1,1,method=’none’)
sage: P=list_plot(listadepuntos)
sage: show(P)

que producirá un gráfico como el mostrado en la figura 2.4 izquierda. Observar que en
la función eulers method() hemos añadido un sexto parámetro method =′ none′ para
evitar una salida en forma de tabla, y que los datos salgan formateados como una lista
de valores (x, y). Para observar como se estructura la lista basta con teclear el comando
eulers method( f , 0, 100, 0,1, 1, method =′ none′ ) sin asignarlo a la variable puntos. La
función list plot() se encarga de producir un dibujo a partir de una lista de valores (x, y).
Si en vez de dibujar los puntos solos, queremos dibujar una lı́nea que una los puntos
entre si, bastará con añadir:

sage: L=line(listadepuntos)
sage: show(P+L)

donde la función line() se encarga de dibujar a partir de las coordenadas (x, y) de la

lista una lı́nea continua uniendo los puntos. Recordemos que al introducir line? nos
aparecerá ayuda sobre los parámetros adicionales de dicha funcion, como [Link]. como
cambiar el color de la lı́nea a dibujar, etc. También conviene observar de nuevo como con
Sage podemos combinar o superponer dos dibujos o objetos diferentes con el operando
+ en la instrucción show(P + L).
Desafortunadamente no todo el monte es orégano, el método de Euler es muy sencillo,
pero tiene un problema: para predecir más o menos correctamente el valor de y(x), cuanto
más lejos se encuentre x del valor de la condición inicial xo , más pequeño deberá ser h,
hasta llegar un punto tal que el valor de h a utilizar es tan pequeño que el numero de
iteraciones es enorme, e incluso con la ayuda de un potente ordenador, podrı́amos tardar
años en ejecutar todos los pasos.
8
Dependiendo de la configuración que uno tenga de Sage es posible que haga falta cargar previamente
la función euler method() mediante el comando: from [Link] import eulers method
2.7 Uso de Sage en la resolución de ecuaciones diferenciales de primer orden
45

Figura 2.4: Solución numérica al problema y′ = −2x con y(0) = 100 y paso h = 0,1.

⋆ Problema - 2.21. Experimentar que pasa, en el problema de condiciones iniciales

anterior, cuando h = 0,5, 0,1, 0,05, 0,01, 0,001. Dibujar conjuntamente la solución
numérica y(x) para cada uno de los valores de h y comparar con la solución analı́tica
del problema. ¿Si vamos a valores de x mayores, como ahora x = 100 que observa-
mos?, ¿como podemos reducir el error?, ¿que ocurre al disminuir h hacia valores muy
pequeños?

El problema del método de Euler radica en que los errores que cometemos al aproxi-
mar y′ por (y(x + h) − y(x))/h con un valor de h que no es infinitamente pequeño se van
acumulando, no se cancelan entre ellos, de forma que cuantas más iteraciones hagamos,
más nos distanciaremos de la solución correcta. Por esta razón, a lo largo de los últimos
siglos se han ido ideando métodos en los que los errores que se comenten al hacer las
aproximaciones tiendan a cancelarse los unos con los otros. Pero dichos métodos son
harina de otro costal y nosotros dejaremos aquı́ el tema de la resolución numérica de las
ecuaciones diferenciales.
46 Solución de las EDO de primer orden

2.8. Problemas Resueltos

⋆ Problema Resuelto - 2.1. Encontrar la familia de soluciones para la ecuación:

(x2 + y2 + x)dx + xydy = 0.
Esta es una EDO de primer orden. En general no sabemos resolver ecuaciones no
lineales, pero si son de primer orden cabe la esperanza de que bien sea de variables
separables (no lo es), o sea de alguno de los tipos de ecuaciones vistas anteriormente.
En este caso lo más razonable es comprobar si es exacta o si no lo es, intentar
transformarla en una exacta: M = x2 + y2 + x, N = xy. ∂y M = 2y 6= ∂x N = y. Por tanto
no es exacta, intentemos encontrar un factor integrante µ . Empezemos por lo más
sencillo, suponer que µ = µ (x), esto nos lleva a

∂y M − ∂x N 2y − y 1
= = ,
N xy x
y por tanto el factor integrante obedece la siguiente ecuación diferencial,
1 dµ 1
= ,
µ x x

que es de variables separables, y nos conduce a µ (x) = x. Ası́ pues la ecuación

original se puede multiplicar por x para obtener la siguiente EDO que es exacta:
(x3 + y2 + x)dx+ x2 ydy = 0 = dΨ(x, y), y por tanto se debe verificar que Ψ(x, y) =

x4 x3 x2 y2
4 + 3 + 2 con lo que llegamos a que la solución implı́cita a la ecuación diferen-
cial es
x 4 x 3 x 2 y2
+ + =C
4 3 2
donde C es una constante. Es este caso en particular seria posible dar la solución como
y(x), pero hay que ser consciente que en muchas ocasiones la ecuación algebraica a
resolver seria trascendente.
2.8 Problemas Resueltos 47

⋆ Problema Resuelto - 2.2. Encontrar la familia de soluciones para la ecuación:

y′ = − ax+by
bx+cy .
En este caso es fácil observar que si multiplicamos por 1/x, obtenemos

a + by/x
y′ = −
b + cy/y

que es una ecuación homogénea, en el sentido de que y′ = f (x/y). Ası́ pues un simple
cambio z = y/x nos conduce a la ecuación

−a − 2bz − c2
xz′ = .
b + cz
Que si nos fijamos bien, resulta ser de variables separables, es decir:
b + cz 1
dz = − dx
a + 2bz + cz2 x

que nos conduce a ln(a + 2bz + cz2 )/2 = −ln(x) + K siendo K una constante. Y por
tanto la solución implı́cita es: ax2 + 2bxy + cy2 = K. Nuevamente, en este caso serı́a
posible expresar la solución de forma implı́cita como y = y(x).

⋆ Problema Resuelto - 2.3. Encontrar la familia de soluciones de la siguiente

ecuación diferencial:
yln(y)dx + (x − ln(y))dy = 0.
Aunque normalmente tendemos a asociar y = y(x), hemos de tener la mente abierta
a otras posibilidades. Fijémonos aquı́ lo complicado que puede resultarnos resolver la
ecuación si nos obcecamos en resolver en función de y, y lo fácil que puede ser si
suponemos en cambio que x = x(y). En este último caso podemos escribir la ecuación
como:
dx 1 1
+ x= ,
dy yln(y) y
que es una ecuación lineal de primer orden que se puede resolver fácilmente tal y como
hemos explicado en la primera sección de este capı́tulo:
dy
Z
µ (y) = exp( ) = ln(y).
yln(y)

Y por tanto x(y) = ln(y)/2 +C/ln(y).

48 Solución de las EDO de primer orden

⋆ Problema Resuelto - 2.4. Encontrar la familia de soluciones de la siguiente

ecuación diferencial : y′ = x/(x2 y + y3 ).
Una vez más si pensamos en términos de x = x(y), podemos reescribir la anterior
ecuación como
dx 1
= xy + y3 ,
dy x
que resulta ser una ecuación de Bernouilli con n = −1. De acuerdo con lo explicado
anteriormente en este capı́tulo, el cambio de variables z = x1−n = x2 nos conduce a la
EDO lineal de primer orden,
dz
− 2zy = 2y3 .
dy
2 2
Ası́ pues, utilizando 2.2, se obtiene µ (y) = e−y . Y por tanto z = x2 = −(1 + y2 ) +Cey
siendo C una constante proveniente de la integración.

⋆ Problema Resuelto - 2.5. Encontrar la solución general de x2 y′ = y2 + 2xy y

razonar en que intervalo de valores de x está garantizada la existencia de una y
sólo una solución cuando la condición inicial es: y(x = 1) = 2. Si la condición inicial
fuera y(x = 0) = 0, ¿que podemos decir acerca de la existencia y unicidad? ¿Son las
predicciones consistentes con la solución general hallada?
Se puede observar fácilmente que esta EDO no lineal de primer orden es homogénea,
en el sentido y′ = f (y/x). Uno puede reescribir la ecuación como y′ = (y/x)2 + 2(y/x).
Propongamos pues el cambio y = zx siendo y′ = z′ x + z que nos conduce a z′ x − z = z2 .
Que resulta ser una ecuación de Bernouilli con n = 2. Si efectuamos el cambio
l = z−1 , obtenemos la ecuación l ′ + l/x = −1/x que es una EDO de primer orden
lineal a la que aplicamos 2.2 para obtener finalmente la solución y = x2 /(−x + C)
siendo C una constante de integración.

Para examinar las consecuencias del teorema de existencia y unicidad, véase

1.4, conviene observar que f (x, y) = (y/x)2 + 2(y/x) y que ∂y f = 2y/x2 + 2/x.
Potencialmente tendremos ”problemasçuando tengamos las condiciones iniciales
(xo , yo ) = (0, 0) puesto que entonces ∂y f es discontinua. Para el caso y(x = 1) = 2
es decir xo = 1 e yo = 2, obtenemos C = 3/2. Por tanto en este caso la solución
es y = x2 /(3/2 − x) que diverge cuando x = 3/2. Por tanto tenemos garantizado
que existe una y sólo una solución para el intervalo x ∈ (0, 3/2). Recordar que el
teorema afirma para nuestro caso que tendremos una solución única en un intervalo
(x1 , x2 ) ⊂ (0, ∞).

En el caso de las condiciones iniciales y(x = 0) = 0, el teorema de existencia y

unicidad no es aplicable puesto que falla la condición que impone el teorema sobre la
continuidad de ∂y f . De hecho, al exigir que y(0) = 0 en la solución general encontrada,
obtenemos que C puede tomar cualquier valor real o complejo, lo cual implica que
existen infinitas soluciones que verifican la ecuación diferencial y que obedecen que
y = 0 cuando x = 0.
2.8 Problemas Resueltos 49

⋆ Problema Resuelto - 2.6. Probar que −2x + y cos(x) + 5exy = 0 es una solución
implı́cita de la ecuación diferencial y′ (cos(x) + 5xexy ) + y(sin(x) + 5exy ) = 2.
La complicación que surge en este ejercicio es que no puedo obtener analı́ticamente
una función y(x) que pueda substituir directamente en la ecuación. En este caso lo que
haremos será usar la diferenciación implı́cita, es decir, aplicaremos d/dx a la solución
implı́cita. Si lo hacemos comprobaremos fácilmente que obtenemos como resultado
una ecuación diferencial igual a la del enunciado del problema.a
a
Para los más rigurosos y exigentes: ¿Quien nos asegura que existe una función y(x) tal que se
puede definir vı́a −2x + y cos(x) + 5exy = 0 y además sea diferenciable? El teorema de la función
implı́cita es la clave, véase [Link]. el libro Cálculo Vectorial de Marsden-Tromba.
50 Solución de las EDO de primer orden

2.9. Ejercicios
1. Resolver y decir en que intervalo es válida la solución de:
 ′
tx + x = t + 1
(a)
x(1) = 0.

tx′ + x = t + 1

(b)
x(−1) = 2.
Dibujad gráficas de las dos soluciones y comprobad la diferencia que hay entre las
dos.

2. Resolver y decir en que intervalo es válida la solución de

 ′
x cost + x sint = 1,
x(0) = 1.

¿es única la solución?

3. Encontrar la familia de soluciones de las siguientes ecuaciones diferenciales:

1. xy′ − y + 3x2 y − x2 = 0.
2
2. (2y + 1)dx + x x−y dy = 0, probar el cambio z = x2 .
3. y′ + y = 1+e 1
x.

4. y + 2xy + xy4 = 0, probar el cambio y = z−1/3 .

′

5. y′ + y = y2 (cos x − sin x), probar el cambio y = z−1 .

6. y′ + 1 = 4e−y sin x, probar el cambio y = ln(z).
7. y′ − 2y cot(2x) = 1 − 2x cot(2x) − 2 csc(2x).
9. y′ = x2 y+y
x 2
3 , probar el cambio z = x .

10. (cos y + y cos x)dx + (sin x − x sin y)dy = 0.

4. Encontrar la familia de soluciones de las siguientes ecuaciones diferenciales:

1. (4x3 y3 − 2xy)dx + (3x4 y2 − x2 )dy = 0.
2. (2x3 y2 + 4x2 y + 2xy2 + xy4 + 2y)dx + 2(y3 + x2 y + x)dy = 0.
3. y2 dx + (x2 − xy − y2 )dy = 0.

5. Resolver el problema de valores iniciales

1/2
xy′ − y = x2 − y2 , con y(1) = 0.
es única la solución? Se contradice el teorema de existencia y unicidad de la
solución?

6. Encontrar la familia de soluciones de las siguientes ecuaciones diferenciales:

1. y′ + 3y = e3x y2 .
2. e−x y′ = 1 − e2x + 2ex y − y2 , sabiendo que y1 = ex es una solución.
2.9 Ejercicios 51

7. Encontrar la familia de soluciones de las siguientes ecuaciones diferenciales:

1. (x4 + y4 )dx − xy3 dy = 0.
2. (2x3 + 3y)dx + (3x + y − 1)dy = 0.
3. xy′ + y − y2 e2x = 0.

8. Encontrar las condiciones para que Mdx + Ndy = 0 admita un factor integrante del
tipo µ (x, y) = µ (x). Usar dicho resultado para encontrar la familia de soluciones
de la siguiente ecuación diferencial: (x2 + y2 + x)dx + xydy = 0.

9. Encontrar las condiciones para que Mdx + Ndy admita un factor integrante del
tipo µ (x, y) = µ (xy).

10. Demostar que si Mdx + Ndy = 0 es una ecuación homogénea, un factor integrante
viene dado por:
1
µ (x, y) = .
xM + yN
11. Demostar que si conocemos dos soluciones particulares y1 (x), y2 (x) de la ecuación
de Ricatti, es posible encontrar la solución completa haciendo uso del cambio:
y1 − zy2
y= .
1−z

12. (a) Encontrar la solución general de la siguiente ecuación diferencial:

x2 y′ = y2 + 2xy.
(b) Enunciar el teorema de existencia y unicidad que podemos aplicar a la anterior
ecuación diferencial, y razonar que intervalo de valores de x nos garantizan la
existencia de una y solo una solución si la condición inicial es y(x = 1) = 2. ¿ Es
la predicción del teorema consistente con la solución encontrada a partir de (a)?
(c) Si la condición inicial fuera y(x = 0) = 0, que información nos proporcionaria el
teorema de existencia y unicidad? ¿Es consistente esta información con la solución
encontrada a partir del apartado (a)?