# Ecuación diferencial de Bernoulli

La ecuación diferencial de Bernoulli es una ecuación diferencial ordinaria de primer orden, formulada por Jacob Bernoulli. Esta ecuación fue transformada, por Gottfried Leibniz en 1693 y por Johann Bernoulli en 1697, en una ecuación diferencial lineal de primer orden mediante el cambio de variable $y^{1-\alpha }=v$, esta ecuación es de la forma
${\frac {dy}{dx}}+P(x)y=Q(x)y^{\alpha }$
donde $P(x)$ y $Q(x)$ son funciones continuas en un intervalo abierto $(a,b)\subseteq \mathbb {R}$ con $\alpha \in \mathbb {R}$.

## Solución
### Caso general ($\alpha \neq 0,1$)

Dividimos la ecuación diferencial entre $y^{\alpha }$ y obtenemos
    ${\frac {1}{y^{\alpha }}}{\frac {dy}{dx}}+P(x)\;{\frac {y}{y^{\alpha }}}=Q(x)$
o, equivalentemente
    ${\frac {1}{y^{\alpha }}}{\frac {dy}{dx}}+P(x)y^{1-\alpha }=Q(x)$
Definiendo $z=y^{1-\alpha }$ obtenemos las igualdades
    ${\begin{aligned}{\frac {dz}{dx}}&=(1-\alpha )y^{-\alpha }{\frac {dy}{dx}}={\frac {1-\alpha }{y^{\alpha }}}{\frac {dy}{dx}}\end{aligned}}$
o
    ${\frac {1}{1-\alpha }}{\frac {dz}{dx}}={\frac {1}{y^{\alpha }}}{\frac {dy}{dx}}$
Reemplazando en la ecuación diferencial
    ${\frac {1}{1-\alpha }}{\frac {dz}{dx}}+P(x)z=Q(x)$
    ${\frac {dz}{dx}}+(1-\alpha )P(x)z=(1-\alpha )Q(x)$
Ecuación que resulta ser una ecuación diferencial lineal cuya solución está dada por
    $z=e^{-(1-\alpha )\int P(x)dx}\left(\int e^{(1-\alpha )\int P(x)dx}(1-\alpha )Q(x)dx+C\right)$
donde $C\in \mathbb {R}$ es una constante arbitraria, como $z=y^{1-\alpha }$ entonces
    $y^{1-\alpha }=e^{-(1-\alpha )\int P(x)dx}\left((1-\alpha )\int e^{(1-\alpha )\int P(x)dx}Q(x)dx+C\right)$
Finalmente
    $y={\sqrt[{1-\alpha }]{e^{-(1-\alpha )\int P(x)dx}\left((1-\alpha )\int e^{(1-\alpha )\int P(x)dx}Q(x)dx+C\right)}}$
    $y=e^{-\int P(x)dx}\left[{\sqrt[{1-\alpha }]{(1-\alpha )\int e^{(1-\alpha )\int P(x)dx}Q(x)dx+C}}\right]$

### Casos particulares

Cuando $\alpha =0$ entonces la ecuación
    ${\frac {dy}{dx}}+P(x)y=Q(x)y^{\alpha }$
se reduce a la ecuación
    ${\frac {dy}{dx}}+P(x)y=Q(x)$
cuya solución está dada por
    $y=e^{-\int P(x)dx}\left(\int Q(x)e^{\int P(x)dx}dx+C\right)$
Cuando $\alpha =1$ entonces la ecuación
    ${\frac {dy}{dx}}+P(x)y=Q(x)y^{\alpha }$
se reduce a
    ${\frac {dy}{dx}}+P(x)y=Q(x)y$
que puede resolverse mediante variables separables, dicha solución está dada por
    $\ln y=\int (Q(x)-P(x))dx+C$

## Ejemplo

Para resolver la ecuación:
   (*)$\qquad xy'+y=x^{4}y^{3}$
Se hace el cambio de variable $z=y^{-2}\;$, que introducido en (*) da simplemente:
   (**)$y^{2}={\frac {1}{z}}\Rightarrow 2yy'=-{\frac {1}{z^{2}}}z'$
Multiplicando la ecuación anterior por el factor: ${\frac {2y}{x}};$ se llega a:
   $\qquad 2yy'+{\frac {2}{x}}y^{2}=2x^{3}y^{4}$
Si se sustituye (**) en la última expresión y operando:
   $-{\frac {z'}{z^{2}}}+{\frac {2}{x}}{\frac {1}{z}}={\frac {2x^{3}}{z^{2}}}\quad \Rightarrow \quad z'-{\frac {2z}{x}}=-2x^{3}$
Que es una ecuación diferencial lineal que puede resolverse fácilmente. Primeramente se calcula el factor integrante típico de la ecuación de Bernouilli:
   $e^{\int P(x)dx}=e^{\int -{\frac {2}{x}}dx}=e^{-2\ln(x)}={\frac {1}{x^{2}}}$
Y se resuelve ahora la ecuación:
   $\left({\frac {z}{x^{2}}}\right)'=-2x^{3}{\frac {1}{x^{2}}}=-2x\qquad {\frac {z}{x^{2}}}=\int {-2xdx}=-2\int {xdx}=-2{\frac {x^{2}}{2}}+C_{1}=-x^{2}+C_{1}$
Deshaciendo ahora el cambio de variable:
   ${\frac {z}{x^{2}}}=-x^{2}+C_{1}\quad \Rightarrow \quad z=C_{1}x^{2}-x^{4}$
Teniendo en cuenta que el cambio que hicimos fue $z=y^{-2}\;$:
   ${\frac {1}{y(x)^{2}}}=C_{1}x^{2}-x^{4}\quad \Rightarrow \quad y(x)={\frac {\pm 1}{\sqrt {C_{1}x^{2}-x^{4}}}}$