Solucionador de EDO de Bernoulli

Solucionador de EDO de Bernoulli

El Solucionador de EDO de Bernoulli aborda una de las ecuaciones diferenciales de primer orden no lineales más famosas — la ecuación de Bernoulli y' + P(x)y = Q(x)yⁿ — y convierte la derivación clásica de los libros de texto en un recorrido interactivo paso a paso. Linealiza la ecuación mediante la sustitución v = y¹⁻ⁿ, construye el factor integrante μ(x) y superpone la curva resultante de forma cerrada sobre una solución numérica RK4 y un campo de pendientes para que puedas ver cada detalle a la vez.

¿Qué es una ecuación diferencial de Bernoulli?

Introducida por Jacob Bernoulli en 1695, una ecuación de Bernoulli es una EDO de primer orden de la forma

Cuando n = 0 la ecuación ya es lineal; cuando n = 1 es separable. Para cualquier otro n real la ecuación es no lineal, pero la sustitución clásica v = y¹⁻ⁿ la convierte en una EDO lineal en v, que puede resolverse con el truco estándar del factor integrante.

El método de Bernoulli en seis pasos

Partiendo de y' + P(x)y = Q(x)yⁿ:

1. Divide por yⁿ: \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \).
2. Sustituye v = y¹⁻ⁿ: observa que \( v' = (1-n)y^{-n}y' \), por lo que \( y^{-n}y' = v'/(1-n) \).
3. Linealiza: \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) — una EDO lineal de primer orden en v.
4. Factor integrante: \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \), así que \( (\mu v)' = \mu(1-n)Q(x) \).
5. Resuelve para v(x): \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. Sustitución inversa: \( y(x) = v(x)^{1/(1-n)} \).

Cuando las integrales involucradas son elementales se obtiene una forma cerrada limpia; cuando no lo son, la calculadora las evalúa numéricamente usando la regla de Simpson para trazar la curva de solución.

Casos especiales manejados automáticamente

Ejemplo resuelto — n = 2, tipo logística

Considera y' + y/x = x·y² con la condición inicial y(1) = 1. Aquí P(x) = 1/x, Q(x) = x, y n = 2, por lo que 1 − n = −1.

1. Sustituye v = y⁻¹ = 1/y. Entonces v' = −y⁻²y' y la ecuación se convierte en v' − (1/x)v = −x.
2. Factor integrante: μ(x) = exp(∫−1/x dx) = 1/x.
3. (μ·v)' = μ·(−x) = −1. Integrando: (1/x)·v = −x + C, es decir, v = −x² + Cx.
4. Aplica la CI: en x = 1, v = 1/1 = 1, por lo que 1 = −1 + C ⇒ C = 2. Por lo tanto v(x) = −x² + 2x.
5. Sustitución inversa: y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x)).

La solución de forma cerrada y = 1/(x(2−x)) tiene asíntotas verticales en x = 0 y x = 2 — exactamente el tipo de cosas que un campo de pendientes hace evidentes de un vistazo.

Cómo usar esta calculadora

1. Completa el constructor de ecuaciones. Escribe P(x) y Q(x) en los espacios azules, y el exponente n en el pequeño cuadro de superíndice. El diseño refleja la forma estándar y' + P(x)y = Q(x)yⁿ.
2. Establece la condición inicial (x₀, y₀) y el rango de graficación [x mín, x máx]. El rango debe contener x₀.
3. Haz clic en Resolver. La calculadora detecta si te encuentras en un caso especial (n = 0 o n = 1) y muestra la derivación correspondiente. De lo contrario, ejecuta la sustitución completa de Bernoulli de seis pasos con ecuaciones renderizadas por MathJax.
4. Analiza el gráfico. La curva naranja es la solución numérica RK4. La curva discontinua azul es la forma cerrada evaluada a través del factor integrante. El campo de flechas muestra y' en todas partes, permitiéndote visualizar otras soluciones también.
5. Copia un CSV de los puntos de muestra si deseas importar la trayectoria en otro programa.

Consejos, dificultades y casos límite

• n no entero requiere y > 0. El solucionador marca combinaciones como n = 1/2 con y₀ ≤ 0, donde yⁿ sería complejo.
• y₀ = 0 suele ser singular. Cualquier ecuación de Bernoulli con Q ≠ 0 y n > 0 tiene la solución trivial y ≡ 0, que normalmente no es la rama que buscas.
• Evita explosiones de P(x) cerca de x₀. Expresiones como 1/x requieren x₀ ≠ 0; el solucionador valida esto antes de ejecutarse.
• Exponentes grandes (|n| > 20) son rechazados para evitar desbordamientos. En la práctica, las ecuaciones de Bernoulli con una n tan grande casi nunca aparecen en problemas reales.
• Asíntotas verticales. Si RK4 diverge, intenta estrechar el rango de x hacia el lado de x₀ donde la solución se mantiene finita.

Dónde aparecen las ecuaciones de Bernoulli

• Dinámica de poblaciones — la ecuación logística y' = ry(1 − y/K) es una ecuación de Bernoulli disfrazada (n = 2 después de reordenar).
• Cinética química — las reacciones autocatalíticas suelen obedecer y' ∝ y − y².
• Circuitos eléctricos — ciertos circuitos RL con resistores no lineales producen la forma de Bernoulli.
• Mecánica de fluidos — ecuaciones de capa límite tras la reducción de similitud.
• Modelos epidémicos — la fracción de susceptibles del modelo SIR puede reducirse a la forma de Bernoulli.
• Crecimiento económico — el modelo de Solow–Swan con tasa de ahorro constante es de Bernoulli con n = α.

Preguntas frecuentes

¿Qué es una ecuación diferencial de Bernoulli?

Una ecuación de Bernoulli es una EDO de primer orden de la forma y' + P(x)y = Q(x)yⁿ, donde P y Q son funciones continuas y n es cualquier número real. Es un ejemplo clásico de una EDO no lineal que puede transformarse en una lineal mediante la sustitución v = y¹⁻ⁿ.

¿Cómo funciona la sustitución v = y¹⁻ⁿ?

Multiplica la ecuación original por y⁻ⁿ para que cada término y se convierta en y¹⁻ⁿ o y⁻ⁿy'. Establecer v = y¹⁻ⁿ da v' = (1−n)y⁻ⁿy'. La sustitución transforma la ecuación de Bernoulli en v' + (1−n)P(x)v = (1−n)Q(x), que es lineal en v y resoluble con un factor integrante.

¿Qué ocurre cuando n = 0 o n = 1?

Cuando n = 0 la ecuación ya es lineal de primer orden, por lo que no se requiere sustitución. Cuando n = 1 la receta de Bernoulli dividiría por 1 − n = 0, por lo que lo manejamos por separado: la ecuación se reduce a y' = (Q(x) − P(x))·y, que es separable con la solución de forma cerrada y = y₀·exp(∫(Q−P) dx).

¿Se pueden resolver siempre las ecuaciones de Bernoulli en forma cerrada?

En principio sí, pero las integrales resultantes que involucran el factor integrante pueden no tener antiderivadas elementales. Cuando esto sucede, la calculadora las evalúa numéricamente con la regla de Simpson y traza la curva de solución. El método en sí siempre reduce una EDO de Bernoulli a cuadraturas.

¿Por qué las y negativas y las n no enteras causan problemas?

Si n no es un entero, yⁿ se define como exp(n·ln y) y solo es real para y > 0. Introducir una y negativa produciría un número complejo. El solucionador advierte esta situación y solicita y₀ > 0 o un exponente entero para que la solución siga siendo de valor real.

¿Qué muestra el campo de pendientes?

El campo de pendientes es una cuadrícula de pequeños segmentos tangentes cuyo ángulo es igual a y' en ese punto (x, y). Cualquier curva de solución está obligada a seguir estas tangentes, por lo que el campo de pendientes permite ver la forma cualitativa de todas las soluciones a la vez, con la condición inicial resaltando la curva particular.

Lecturas adicionales

• Ecuación diferencial de Bernoulli — Wikipedia
• Factor integrante — Wikipedia
• Función logística — Wikipedia
• Campo de pendientes — Wikipedia