Crecimiento exponencial
Campo de pendientes Solución (Euler, h=0,1) Condición inicial (0, y₀)
Explora ODEs de forma visual e interactiva: campo de direcciones, sistema predador-presa animado y aplicaciones reales en ingeniería, biología y física.
100% visual · Sin instalación · Método de Euler y Runge-Kutta 4
De las poblaciones al circuito eléctrico: matemáticas que modelan el cambio
Una ecuación diferencial ordinaria (ODE) relaciona una función desconocida y(x) con sus derivadas. En lugar de decir "y vale esto", dice "la tasa de cambio de y depende de y y x de esta forma". Son la herramienta fundamental para modelar cualquier sistema que cambia en el tiempo: poblaciones, temperatura, corriente eléctrica, posición de un planeta.
En cada punto del plano (x, y) la ODE dy/dx = f(x, y) nos da una pendiente. Al dibujar todas las pendientes simultáneamente obtenemos el campo de direcciones: una "textura" del plano que muestra hacia dónde fluiría cualquier solución. La curva solución es la que en todo punto es tangente a las flechas del campo.
En los años 1920, Vito Volterra modeló las oscilaciones de peces en el Adriático y Alfred Lotka lo hizo independientemente para reacciones químicas. El sistema dR/dt = αR − βRL / dL/dt = δRL − γL captura la coexistencia oscilatoria entre presa y depredador: cuando los conejos abundan, los lobos prosperan; cuando los lobos son muchos, los conejos escasean; entonces los lobos también disminuyen, y el ciclo comienza de nuevo.
Cuando no existe solución analítica (la mayoría de casos reales), usamos métodos numéricos. El método de Euler avanza paso a paso: y(t+h) ≈ y(t) + h·f(t, y(t)). Simple pero poco preciso. El Runge-Kutta 4 evalúa la pendiente en cuatro puntos del intervalo y promedia ponderadamente — el error es O(h⁴) vs O(h) de Euler.
| Tipo | Forma | Método de resolución | Ejemplo | Aplicación |
|---|---|---|---|---|
| Separable | dy/dx = f(x)g(y) | Separar variables e integrar | dy/dx = ky | Crecimiento/decaimiento exponencial |
| Lineal 1er orden | dy/dx + p(x)y = q(x) | Factor integrante μ = e^∫p(x)dx | dy/dt + y/RC = V/RC | Circuito RC |
| Autónoma | dy/dx = f(y) | Análisis de puntos de equilibrio | dy/dt = ky(1−y/K) | Logística (poblaciones) |
| Sistema 2D | dx/dt=f(x,y), dy/dt=g(x,y) | Eigenvalores de la matriz Jacobiana | Lotka-Volterra | Predador-presa, circuitos LC |
El modelo Lotka-Volterra predice oscilaciones de poblaciones de presa y depredador en ecosistemas reales.
Ejemplo: fluctuaciones históricas de liebres y linces en Canadá siguen el patrón Lotka-Volterra.
Circuitos RC, RL y RLC son EDO de primer y segundo orden; el comportamiento transitorio se calcula con ellas.
Ejemplo: el tiempo de carga de un condensador es τ = RC — una EDO de primer orden.
La ley de enfriamiento de Newton predice cuándo un cuerpo alcanza el equilibrio térmico con su entorno.
Ejemplo: calcular cuánto tarda un café en enfriarse a temperatura óptima de consumo.
La cinética de eliminación de fármacos (vida media) sigue una EDO exponencial: dC/dt = −kC.
Ejemplo: la vida media de un antibiótico determina la frecuencia de dosificación correcta.
Las ecuaciones diferenciales ordinarias (EDO) tienen una sola variable independiente, normalmente el tiempo t o una dimensión espacial. Las ecuaciones en derivadas parciales (EDP) tienen varias variables independientes simultáneas (ej: temperatura en función de x, y y t).
En cada punto (x, y) del plano, la EDO dy/dx = f(x, y) indica una pendiente. El campo de direcciones dibuja esa pendiente como una flecha pequeña en cada punto. La solución particular es la curva que sigue las flechas partiendo de una condición inicial.
El método de Euler aproxima la solución avanzando paso a paso: y(t+h) ≈ y(t) + h·f(t, y(t)). Nunca es exacto (salvo que la función sea constante), pero el error se reduce al reducir h. Para h suficientemente pequeño, la aproximación es útil en la práctica.
Con parámetros positivos (α, β, δ, γ > 0) y condiciones iniciales positivas, el sistema siempre oscila de forma periódica — nunca diverge ni converge a un punto fijo. Las oscilaciones son conservativas: el sistema conserva una cantidad llamada "integral primera" de Lotka-Volterra.
Un punto de equilibrio es estable si pequeñas perturbaciones hacen que el sistema vuelva a él (atractor). Es inestable si las perturbaciones se amplifican (repulsor). Para EDO autónomas dy/dt = f(y), un equilibrio y* es estable si f'(y*) < 0 e inestable si f'(y*) > 0.
Euler usa solo la pendiente al inicio del intervalo. Runge-Kutta 4 evalúa la pendiente en cuatro puntos del intervalo [t, t+h] y calcula un promedio ponderado. Esto elimina los errores de curvatura que acumula Euler, logrando un error proporcional a h⁴ en lugar de h.
¿Es separable? ¿Lineal de primer orden? ¿Autónoma? La clasificación determina el método de resolución correcto. Busca si puedes escribirla como dy/dx = g(x)·h(y) (separable) o dy/dx + p(x)y = q(x) (lineal).
Pasa dy y todo lo que dependa de y a un lado, y dx con todo lo que dependa de x al otro. Luego integra ambos lados por separado: ∫(1/h(y)) dy = ∫g(x) dx.
Calcula μ(x) = e^∫p(x)dx. Multiplica toda la ecuación por μ(x). El lado izquierdo se convierte en d/dx[μ(x)·y], que se integra directamente: μ(x)·y = ∫μ(x)·q(x) dx.
La solución general contiene una constante C. Sustituye la condición inicial y(x₀) = y₀ para determinar C. Así obtienes la solución particular que describe el fenómeno concreto.
Diferencia la solución obtenida y sustitúyela en la EDO original. Si ambos lados son iguales, la solución es correcta. Si la EDO describe un fenómeno físico, comprueba también que las unidades son coherentes.
Antes de resolver analíticamente, dibuja el campo de direcciones. La solución debe seguir las flechas — si no lo hace, hay un error.
Los puntos donde dy/dt = 0 son "atractores" o "repulsores". Identifícalos antes de simular para entender el comportamiento a largo plazo.
El método de Euler tiene error O(h); Runge-Kutta 4 tiene error O(h⁴). Para la misma precisión, puedes usar un paso h mayor con RK4.
Las oscilaciones de Lotka-Volterra son siempre cerradas (conservan energía): el diagrama de fase es una curva cerrada alrededor del punto de equilibrio.
Si la EDO modela un fenómeno físico, las unidades de ambos lados deben coincidir. Un error de unidades delata un modelo incorrecto.