Práctica 1: Estudio cualitativo de ecuaciones
diferenciales de primer orden
1. Objetivos
Obtener el campo de direcciones de una ecuación diferencial ordinaria
(EDO) de primer orden, tanto por el método manual de las isoclinas,
como utilizando el programa Maxima.
Dibujar curvas solución de una EDO de primer orden a partir del
campo de direcciones.
Analizar el crecimiento o decrecimiento de las soluciones de EDO’s de
primer orden.
Calcular soluciones de algunas EDO’s y de problemas de valor inicial
relativos a ellas con Maxima.
2. Campos de pendientes y curvas solución
Desde el punto de vista geométrico, una curva solución de la EDO de
primer orden:
y ′ (x) = f (x, y) (1)
es una curva en el plano tal que la pendiente de la recta tangente en cada
punto (x, y) es m = f (x, y). Esto sugiere un método gráfico para construir
soluciones aproximadas de (1): en puntos representativos (x, y) del plano,
dibujamos segmentos de pendiente m = f (x, y). El conjunto de esos segmen-
tos se llama campo de pendientes ó campo de direcciones de la ecuación (1).
Los vectores (1, f (x, y)) tienen la pendiente deseada, pero tienen diferentes
longitudes en los distintos puntos (x, y), por lo que suelen normalizarse.
Siguiendo el campo de pendientes, podemos dibujar curvas soluciones de
la ecuación diferencial, ya que cada curva solución es tangente en cada uno de
sus puntos al segmento del campo de pendientes en ese punto. Ası́, el campo
de direcciones nos permite realizar un estudio cualitativo de las soluciones
directamente desde la ecuación diferencial, sin necesidad de obtener una
fórmula para ellas.
2.1. Método de las isoclinas
Describiremos en primer lugar un método manual para dibujar el campo
de pendientes de (1), basado en las isoclinas. Una isoclina de la ecuación
1
diferencial (1) es una curva en el plano de la forma:
f (x, y) = C, C constante,
es decir, es una curva a lo largo de la cual la pendiente de las rectas tangentes
es la misma. Por tanto, para obtener el campo de pendientes de (1), se
pueden seguir los siguientes pasos:
1. Se dibujan varias isoclinas f(x,y)=C.
2. En puntos representativos de cada isoclina se dibujan vectores con la
misma pendiente: m=C, los cuales formarán el campo de pendientes.
En particular, si la ecuación diferencial es autónoma, es decir, de la forma
y ′ = f (y), las isoclinas son f (y) = C, luego son rectas horizontales. Y si es
de la forma y ′ = f (x), las isoclinas son f (x) = C, luego rectas verticales.
En la isoclina de nivel 0: f (x, y) = 0, la pendiente de la solución es
0. Esta isoclina suele separar las isoclinas f (x, y) = C, C > 0, donde las
soluciones tienen pendiente positiva y por tanto son funciones crecientes, de
las isoclinas f (x, y) = C, C < 0, en las que las soluciones tienen pendiente
negativa y en consecuencia son funciones decrecientes.
Si la ecuación diferencial es de la forma F (x, y, y ′ ) = 0, las isoclinas son
las curvas de la familia F (x, y, C) = 0.
Para dibujar con Maxima las isoclinas de (1) para x ∈ [a, b], y ∈ [c, d],
utilizaremos la orden:
wxcontour plot(f, [x,a,b], [y,c,d])$
2.2. Campos de pendientes con Maxima
El paquete plotdf de Maxima permite dibujar campos de pendientes de
EDO’s de primer orden y de sistemas de dos EDO’s autónomas de primer
orden. Puesto que se trata de un paquete adicional, para poderlo utilizar
primero lo cargaremos con la orden
load(plotdf )$
Para dibujar el campo de pendientes de la EDO de primer orden (1)
utilizaremos la orden plotdf, dándole como argumento la expresión de la
función f :
plotdf(f )$
En esta orden plotdf la variable independiente debe ser siempre x y la
variable dependiente y. Ninguna de estas dos variables podrá tener asociado
un valor numérico.
Por ejemplo, para dibujar el campo de direcciones de la ecuación
y ′ = e−x + y,
2
podemos escribir:
(%i1) load(plotdf)$
(%i2) plotdf(exp(-x)+y)$
y aparece una ventana gráfica con el campo de pendientes buscado. Sobre él
se pueden dibujar de forma interactiva soluciones particulares: basta hacer
click en un punto del gráfico para que se dibuje la solución que pasa por
él. Como puedes ver, la ventana gráfica que ha aparecido dispone de varios
menús cuya utilidad puedes ir explorando.
La orden plotdf admite varias opciones, cada una de las cuales es una
lista entre corchetes con dos ó más elementos. Destacamos aquı́ las siguientes:
[x, a, b] especifica que el rango de valores de x es el intervalo [a, b].
[y, c, d] especifica que el rango de valores de y es el intervalo [c, d].
[trajectory at, x0 , y0 ] dibuja sobre el campo de pendientes la solución
que verifica y(x0 ) = y0 .
[xfun, “expr1; . . . ; exprn”] dibuja sobre el campo de pendientes las fun-
ciones de x especificadas en las expresiones entre las comillas.
Por ejemplo, escribiendo
(%i3) plotdf(exp(-x)+y, [trajectory_at,2,-1])$
se obtiene la gráfica de la Figura 1.
-4
-8
-8 -4 0 4 8
Figura 1: campo de pendientes de y ′ = e−x + y y solución verificando y(2) =
−1
3
El menú Ayuda de Maxima te puede ayudar para conocer la utilidad y
sintaxis de este comando ó de cualquier otro.
3. Soluciones de EDO’s con Maxima
Maxima dispone de varias funciones para resolver analı́ticamente ciertos
tipos de ecuaciones diferenciales ordinarias. La función ode2 resuelve EDO’s
de primer y segundo orden. Su sintaxis es:
ode2(ecu, vard, vari);
donde ecu es la EDO a resolver, vard la variable dependiente y vari la varia-
ble independiente. Si encuentra la solución, la devuelve en forma explı́cita
ó implı́cita.
Para expresar y ′ (x) = dx
dy
en ecu se utiliza ’diff(y,x), siendo obligatorio
el uso del apóstrofo al comienzo para evitar el cálculo de la derivada (que
en nuestro caso darı́a cero al no haber definido y como función de x). Ası́,
para resolver la ecuación (de variables separadas)
dy
x3 (y − 1) + (x − 1)y 3 = 0,
dx
podemos escribir:
(%i4) ecu:x^3*(y-1)*’diff(y,x)+(x-1)*y^3=0;
(%o4) x^3*(y-1)*(’diff(y,x,1))+(x-1)*y^3=0
(%i5) ode2(ecu,y,x);
(%o5) (2*y-1)/(2*y^2)=%c-(2*x-1)/(2*x^2)
En la última salida vemos que Maxima utiliza el sı́mbolo %c para denotar
una constante arbitraria.
Si lo que se quiere es resolver un problema de valor inicial para una EDO
de primer orden, y la condición inicial es y(x0 ) = y0 , utilizaremos la orden
ic1(sol, x = x0 , y = y0 );
donde sol es la solución general obtenida con ode2.
Ası́, para encontrar la solución de la ecuación anterior que verifica y(2) =
−3, escribiremos:
(%i6) ic1(%,x=2,y=-3);
(%o6) (2*y-1)/(2*y^2)=-(x^2+72*x-36)/(72*x^2)
Poniendo % como primer argumento de la función ic1 estamos indicando
que dicho argumento es la última salida obtenida. También podrı́amos haber
puesto %o5, ya que la solución general se habı́a obtenido en la salida (%o5).
4
4. Problemas
1. Dada la ecuación diferencial y ′ = −y, se pide:
a) Obtener las isoclinas y desde ellas dibujar manualmente el campo
de pendientes y algunas curvas solución.
b) ¿A qué tiende cada solución y(x) cuando x → +∞?
c) Dibujar con Maxima el campo de pendientes para x, y ∈ [−4, 4].
d ) Calcular con Maxima las soluciones que satisfacen y(−4) = 4,
y(−4) = −4.
e) Dibujar sobre el campo de pendientes dichas soluciones.
2. Para la ecuación diferencial y ′ = sen(x − y), se pide:
a) ¿Cómo son las isoclinas de esta ecuación diferencial?
b) Calcular con Maxima las soluciones que satisfacen y(0) = −3/4,
y(0) = −π/2.
c) Dibujar dichas soluciones sobre el campo de pendientes para x ∈
[−π, π], y ∈ [−2π, 2π].
3. Sea la ecuación diferencial y 2 + x2 y ′ = 0. Se pide:
a) Dibujar el campo de pendientes para x, y ∈ [−3, 3].
b) Estudiar el crecimiento o decrecimiento de sus soluciones.
c) ¿Es única la solución de esta ecuación que satisface la condición
inicial y(0) = 0?
d ) ¿Existe alguna solución que satisfaga y(0) = b, con b ̸= 0?
e) Si a ̸= 0 y b arbitrario, ¿cuántas soluciones satisfacen y(a) = b?
4. Se considera la ecuación y ′ = x − y. Se pide:
a) Dibujar el campo de pendientes para x ∈ [−4, 4], y ∈ [−8, 8].
b) Representar en la misma figura las isoclinas f (x, y) = C, para
C = −3, −2, −1, 0, 1, 2, 3.
c) Representar en la misma figura las soluciones que satisfacen y(−4) =
6, y(−4) = −8.
d ) Estudiar el crecimiento o decrecimiento de las soluciones de la
ecuación dada.
5. Para la ecuación y ′ = y(1 − y), se pide:
a) Dibujar el campo de direcciones para x ∈ [−4, 4], y ∈ [−5, 5].
b) Obtener sus isoclinas.
c) Estudiar el crecimiento y decrecimiento de sus soluciones.
