Solucionador de EDO de Primer Orden

Solucionador de EDO de Primer Orden

El Solucionador de EDO de Primer Orden toma una ecuación diferencial ordinaria en la forma dy/dx = f(x, y), clasifica automáticamente su estructura (separable, lineal, autónoma, exacta o general) y produce tanto una solución simbólica de forma cerrada donde sea posible como una solución numérica de alta precisión en todos los casos. Una visualización en vivo del campo de pendientes con la curva de solución superpuesta hace que el significado geométrico de la ecuación sea inmediatamente obvio: las soluciones son exactamente las curvas tangentes a cada flecha.

¿Qué es una EDO de primer orden?

Una ecuación diferencial ordinaria de primer orden involucra únicamente a una función desconocida y(x) y su primera derivada y'(x). La forma explícita estándar es:

Combinada con una condición inicial y(x₀) = y₀, esto define un problema de valor inicial (PVI). El teorema de Picard-Lindelöf garantiza una solución única en algún entorno de x₀ siempre que f sea continua de Lipschitz en y cerca de (x₀, y₀). Geométricamente, el PVI busca la curva única que pasa por (x₀, y₀) cuya pendiente en cada punto coincide con f en ese punto; exactamente la curva tangente al campo de pendientes.

Seis clases que el solucionador reconoce

Método de forma cerrada: Lineal con coeficientes constantes

Cuando el lado derecho se simplifica a dy/dx = a·y + b con constantes a y b, el factor integrante μ(x) = e^(-a·x) da una solución exacta. La solución general es:

Aplicar la condición inicial y(x₀) = y₀ determina la constante C y produce la única solución particular. Esta clase única cubre una enorme cantidad de problemas de libros de texto:

• Crecimiento exponencial — dy/dx = k·y, solución particular y(t) = y₀·e^(k·t).
• Decaimiento exponencial — dy/dx = -k·y, vida media ln 2 / k.
• Ley de enfriamiento de Newton — dy/dx = -k·(y - T_amb), la temperatura del cuerpo se relaja exponencialmente hacia la ambiente.
• Carga de un circuito RC — dV/dt = (1/RC)·(V_in - V), el voltaje del capacitor se aproxima a la fuente.
• Aclaramiento de fármacos — farmacocinética de primer orden con tasa de eliminación k.

Lectura de un campo de pendientes

En cada punto de la cuadrícula (x, y), la herramienta dibuja un segmento de línea corto cuya pendiente es igual a f(x, y). Tres observaciones útiles:

• Equilibrios son puntos donde f(x, y) = 0; el campo de pendientes es horizontal. Para ecuaciones autónomas, estos son puntos fijos y* que satisfacen f(y*) = 0; las trayectorias cercanas se acercan (estable) o se alejan (inestable) de y*.
• Isoclinas son curvas donde f(x, y) es igual a una constante c, por lo que todas las flechas a lo largo de la curva tienen la misma pendiente c.
• Las curvas de solución nunca se cruzan (cuando f es Lipschitz); es visualmente obvio porque dos curvas que se cruzan necesitarían pendientes diferentes en la intersección.

Método numérico: Runge-Kutta clásico (RK4)

Dado (x_n, y_n), el siguiente valor se calcula promediando cuatro estimaciones de pendiente:

RK4 tiene un error de truncamiento local de O(h⁵) y un error global de O(h⁴), lo que proporciona una precisión de aproximadamente seis dígitos con el conteo de pasos predeterminado para ecuaciones no rígidas. El solucionador integra hacia afuera desde el punto inicial en ambas direcciones de x y se detiene limpiamente si la magnitud de y excede 10¹⁵, algo típico en soluciones que explotan en tiempo finito, como dy/dx = y².

Cómo usar esta calculadora

1. Ingrese el lado derecho en el campo dy/dx = .... Use x e y como variables, * para multiplicación, ^ o ** para potencias, y funciones estándar como sin, cos, exp, log, sqrt. Se reconocen las constantes pi y e.
2. Especifique la condición inicial (x₀, y₀); la curva de solución única pasará por este punto.
3. Elija el rango de x sobre el cual graficar el campo de pendientes y la curva de solución. El rango de y se ajusta automáticamente a partir de la solución integrada.
4. Haga clic en Resolver y Visualizar. El clasificador se ejecuta primero; si su ecuación coincide con un patrón de forma cerrada (lineal con coeficientes constantes), obtendrá la respuesta simbólica. El campo de pendientes y la curva de solución siempre se renderizan.
5. Active o desactive el campo de pendientes para enfocarse en la curva de solución, o repita la animación del trazado de la curva para ver cómo progresa la integración desde el punto inicial.

Ejemplo resuelto: Ley de enfriamiento de Newton

Una taza de café a 80 °C se enfría en una habitación a 20 °C. La tasa de transferencia de calor es proporcional a la diferencia de temperatura:

Esta es lineal con coeficientes constantes (a = -0.1, b = 2). La forma cerrada es:

Después de 30 minutos: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. La vista del campo de pendientes hace que el comportamiento límite sea obvio: cada curva de solución, independientemente de la temperatura inicial, tiende asintóticamente a la línea horizontal T = 20.

Aplicaciones comunes

• Dinámica de poblaciones — modelos de crecimiento exponencial, logístico y efecto Allee.
• Farmacocinética — absorción y eliminación de fármacos, cálculos de vida media.
• Transferencia de calor — ley de enfriamiento de Newton, modelos de capacitancia agrupada.
• Circuitos RC y RL — transitorios eléctricos lineales de primer orden.
• Decaimiento radiactivo — cadenas de decaimiento de un solo isótopo.
• Tanques de mezcla — concentración de un soluto bajo flujo de entrada / salida.
• Objeto en caída con resistencia — análisis de velocidad terminal dv/dt = g - kv.

Preguntas frecuentes

¿Qué es una ecuación diferencial ordinaria de primer orden?

Una ecuación diferencial ordinaria (EDO) de primer orden es una ecuación de la forma dy/dx = f(x, y) que involucra a la función desconocida y(x) y su primera derivada. Resolver la EDO significa encontrar la función y(x) cuya derivada coincida con el lado derecho. Con una condición inicial y(x₀) = y₀, la solución es única bajo suposiciones de regularidad moderadas (teorema de Picard-Lindelöf).

¿Qué es un campo de pendientes?

Un campo de pendientes (o campo de direcciones) traza un pequeño segmento de línea en cada punto de la cuadrícula (x, y) cuya pendiente es igual a f(x, y). Las curvas de solución de la EDO son exactamente las curvas tangentes a estos segmentos en cada punto. El campo de pendientes proporciona una intuición visual instantánea del comportamiento global de las soluciones sin resolver la ecuación simbólicamente.

¿Qué clases de EDO de primer orden resuelve esta herramienta?

La herramienta clasifica automáticamente la ecuación en una de las siguientes: integrable (depende solo de x, resuelta por integración directa), lineal con coeficientes constantes y' = a·y + b (se proporciona la forma cerrada completa), autónoma (depende solo de y), separable (se factoriza como g(x)·h(y)), lineal con coeficientes variables (P(x)·y + Q(x)) o general. Para cada clase, se genera una solución numérica de Runge-Kutta de alta precisión y una visualización del campo de pendientes.

¿Qué método numérico se utiliza?

Se aplica el método clásico de Runge-Kutta de cuarto orden (RK4) con 300 subpasos en cada dirección desde el punto inicial. El RK4 tiene un error de truncamiento local de O(h⁵) y es el estándar para EDO no rígidas a esta escala. El solucionador detecta divergencias (desbordamiento o NaN) y detiene la integración limpiamente para que la gráfica siga siendo válida.

¿Qué es el método del factor integrante para EDO lineales?

Para una EDO lineal de primer orden y' + P(x)·y = Q(x), se multiplican ambos lados por el factor integrante μ(x) = e^∫P(x) dx. El lado izquierdo se convierte en la derivada exacta d/dx[μ·y], por lo que y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Cuando P y Q son constantes, esto se reduce a la forma cerrada y = -b/a + C·e^(a·x), que la herramienta devuelve automáticamente.

¿Puede esta herramienta manejar ecuaciones rígidas o sistemas de EDO?

Este solucionador está destinado a EDO escalares de primer orden no rígidas. Los problemas muy rígidos (donde la solución tiene múltiples escalas de tiempo que difieren en muchos órdenes de magnitud) pueden requerir un método implícito como Euler hacia atrás o Rosenbrock; los sistemas acoplados requieren un solucionador con valores vectoriales. Para esos casos, use un paquete dedicado como solve_ivp de SciPy o solucionadores especializados para EDO rígidas.

Lectura adicional

• Ecuación diferencial ordinaria — Wikipedia
• Campo de pendientes — Wikipedia
• Métodos de Runge-Kutta — Wikipedia
• Factor integrante — Wikipedia
• Teorema de Picard-Lindelöf — Wikipedia