VIGILADA MINEDUCACIÓN

ESCUELA Departamento de Matemáticas

COLOMBIANA Guı́a de trabajo 1: Campo de pendientes
DE INGENIERÍA
JULIO GARAVITO Competencias: PA2, R2, RM2, CM2, SP2,C2
2024 - 1
UNIVERSIDAD

Objetivo: Realizar las gráficas de soluciones de ecuaciones diferenciales ordinarias de primer orden y de proble-
mas de valor inicial a partir del campo de pendientes, usando un asistente computacional y usrlas para resolver
problemas de aplicación.

Curvas solución sin solución: Campo de pendientes

dy
En términos geométricos, resolver el problema de valor inicial de primer orden = f ( x, y), con y( x0 ) = y0 ,
dx
significa encontrar la función y( x ) = φ( x ) cuya pendiente de recta tangente a la curva de φ( x ) en el punto ( x0 , y0 )
es m = f ( x0 , y0 ) y que pasa por el punto( x0 , y0 ). De aquı́ que realizar el campo de pendientes o de direcciones de
la ecuación diferencial signifique realizar la gráfica de un pequeño segmento de recta con pendiente m en cada
punto ( x0 , y0 ).

dy
Las siguientes instrucciones permiten realizar el campo de pendientes de la ecuación diferencial = f ( x, y),
dx
usando Geogebra
Instrucción
CampoDirecciones(f(x,y))
CampoDirecciones(f(x,y),n)

En la segunda instrucción el valor de n indica el tamaño del campo

dy x − 1 dy
Ejemplo: Realizar el campo de pendientes de las ecuaciones diferenciales = y = x2 + y.
dx y dx
Ejemplo

campo1=CampoDirecciones((x-1)/y)
campo2=CampoDirecciones(x^(2)+y,20)

Si ahora lo que se quiere es realizar la gráfica de una curva solución para un problema de valor inicial,
dy

dx = f ( x, y)
y( a, b)
se usa la instrucción:
Instrucción
ResuelveEDO(f(x,y),(a,b))

dy 4x2 − 2y
Ejemplo: Realizar el campo de pendientes de la ecuación diferencial = y la solución que pasan
dx x
por el punto (-2,3) y la que pasa por el punto (1,2).
 Ejemplo

campo3=CampoDirecciones((4x^2-2y)/x)
ResuelveEDO((4x^2-2y)/x,(-2,3)
ResuelveEDO((4x^2-2y)/x,(1,2))

Si al usar la instrucción no se indica el punto si no que se deja ( a, b) se obtiene la solución de la ecuación dife-
rencial y usando deslizadores para a y b se realizan las gráficas para diferentes puntos

Ejemplo

ResuelveEDO((4x^2-2y)/x)
ResuelveEDO((4x^2-2y)/x,(a,b))

Si se quiere hacer la gráfica de la familia de funciones, en GeoGebra se determina la solución de la ecuación

en un punto variable ( a, b) y luego se da click en el resultado obtenido, se selecciona propiedades del objeto y se
elige mostrar rastro, al movimiento del deslizador se mantienen las gráficas de las soluciones.
Para realizar la gráfica de la solución a un problema de valor inicial especı́fico, simplemente se resuelve la
ecuación y se define especı́ficamente el punto, para realizar la gráfica del punto se selecciona punto en el menú
dy
y se definen las coordenadas de él. Ası́, para resolver el problema de valor inicial = f ( x, y) con y( a) = b
dx
entonces se escribe:

Instrucción
ResuelveEDO(f(x,y),(a,b))
A=(a,b)

Ecuaciones Autónomas
dy
Una ecuación autónoma de primer orden = f ( x, y) es una ecuación diferencial para la cual la función
dx
f ( x, y) = g(y) es una función que solo depende de la incógnita y y no de su variable independiente, esto es
dy
= g ( y ).
dx
Los puntos crı́ticos c para una ecuación diferencial autónoma son los valores reales que hacen que g(y) = 0,
dy
esto significarı́a que = g(c) = 0 y por tanto la función constante y( x ) = c serı́a una solución constante de la
dx
ecuación diferencial autónoma. A dicha solución se le denomina solución de equilibrio o solución estacionaria. Los
puntos crı́ticos de la ecuación autónoma proporcionan las únicas soluciones constantes de la ecuación.

El diagrama de fase para la ecuación diferencial autónoma de primer orden permite determinar el comporta-
miento de las soluciones de la ecuación diferencial. Si c > 0 entonces la función solución y( x ) es creciente y si
c < 0 la función solución es decreciente. A la recta vertical sobre la cual se analizan las condiciones de y( x ) se le
llama recta de fase.

El punto crı́tico y = c se dirá que es asintóticamente estable (atractor) si todas las posibles soluciones y( x )
presentan un comportamiento asintótico hacia c, esto es, lı́m y(t) = c ; si por el contrario las curvas de las posi-
t→∞
bles soluciones se ”alejan”de y = c, entonces el punto crı́tico se dice que es inestable (repulsor) y finalmente, si
para valores ”un poco”mayores de c la gráfica de la posible solución se aleja (se acerca) de c y para valores ün
poco”menores de c la gráfica se acerca (se aleja) entonces el punto crı́tico se denomina semiestable.

Usando WolframAlpha se obtienen las raı́ces de un polinomio usando Roots[P(m)==0,m], y para resolver una
desigualdad se usa Reduce[g(c) > 0, c] o Reduce[g(c) < 0, c]

dy
Ejemplo: Realizar el campo de pendientes de la ecuación diferencial autónoma = y2 (4 − y2 ).
dx
Para hallar las soluciones de equilibrio de la ecuación autónoma y clasificar su estabilidad se hallan los puntos
crı́ticos:

Ejemplo con WolframAlpha

Roots[c^2*(4-c^2)==0,c];
Reduce[c^2*(4-c^2)>0,c];
Reduce[c^2*(4-c^2)<0,c];

Con los resultados obtenidos se concluye que las soluciones de equilibrio son y = −2, y = 0 y y = 2, y
c2 (4 − c2 ) > 0 cuando c ∈ (−2, 0) ∪ (0, 2), lo que significa que las funciones solución posibles son crecientes en
dichos intervalos, y serán decrecientes para c ∈ (−∞, −2) ∪ (2, ∞).

Analizando el cambio de signo se tiene que y = 2 es asintóticamente estable, en el caso de y = 0, es semiestable

y finalmente el punto crı́tico y = −2 es inestable, con esta información se realiza el diagrama de fase. Al graficar
el campo de pendientes se puede visualizar esta conclusión.

Ejercicios
1. Usando un asistente computacional, realice el campo de pendientes junto con la gráfica de la familia de
soluciones de la ecuación diferencial. Además, realice la gráfica de la solución particular para el problema
de valor inicial dado, si no tiene condición, elija usted una.

a) y′ = x, con y(0) = −3
b) y′ = x + y, con y(−2) = 2
c) y′ = x2 − y
d) ( x 2 + 1)y′ + 3xy = 6x
e) y′ = xey

2. Usando el diagrama de fase, realizar el análisis de los puntos crı́ticos para cada una de las siguientes ecua-
ciones autónomas.

a) y′ = y(3 − y)(y − 2)
b) y′ = y2 − y3
c) y′ = 10 + 3y − y2
d) y′ = y5 − 4y3
e) y′ = (y − 2)4

3. Use el análisis para las ecuaciones autónomas (diagrama de fase) para responder las preguntas sobre el
problema planteado.
Sea P(t) la población de cierta especie en un parque natural, con t tiempo en años. La ecuación diferencial
dP P(1500 − P
= describe la tasa de cambio de la población de la especie.
dt 3200
 a) Si la población inicial de dicha especie es de 1200 elementos, ¿qué ocurre con la población después de
mucho tiempo?. Realice un bosquejo de la posible solución sin resolver la ecuación diferencial.
b) Si ahora la población inicial es de 1800 elementos, ¿qué ocurre con la población después de mucho
tiempo?

Si en ciertas temporadas del año se permite cazar a una razón de s unidades por año:
c) Escriba la ecuación resultante que muestra la tasa de cambio de la población al cazar s número de
unidades anuales.
d) Realice el diagrama de fase de la tasa de cambio de la población si se permite cazar s = 0, 50, 100,
125, 150, 175, 200 elementos.
e) Explique el significado de los puntos de equilibrio.
f ) ¿Hasta qué valor de s se debe permitir cazar para que la especie no desaparezca?

dP
4. Se sabe que la tasa de cambio de la población de cierta especie (en miles) está dada por la ecuación =
dt
3P − 2P2 . Considere t en años.

a) Use el diagrama de fase de la ecuación para determinar el comportamiento de la población de la espe-

cie.
b) Si la población inicial es de 2000, ¿qué ocurre con la población después de mucho tiempo?. Explique.
c) Si ahora se considera que la especie inició con cien especı́menes, ¿qué ocurre con la población después
de mucho tiempo?
d) Mediante análisis se concluye que la tasa de nacimientos es de 150 por cada trimestre y que mueren s
ejemplares en ese mismo periodo, ¿cuál es la ecuación diferencial que describe la tasa de cambio de la
población de la especie anualmente?
e) Analice el comportamiento de las posibles soluciones teniendo en cuenta esa tasa de muertes por tri-
mestre.