“Obvio” es la palabra más peligrosa del mundo en matemáticas.
– E. T. Bell

Introducción

Con esta entrada concluiremos el desarrollo de métodos de resolución de ecuaciones diferenciales de primer orden.

Presentaremos dos ecuaciones diferenciales no lineales más, conocidas como ecuación diferencial de Bernoulli y ecuación diferencial de Riccati en honor a sus formuladores Jacob Bernoulli y Jacopo Francesco Riccati, respectivamente.

Ecuación diferencial de Bernoulli

La ecuación diferencial de Bernoulli es una ecuación diferencial ordinaria de primer orden formulada por Jacob Bernoulli en el siglo XVll.

Si a la ecuación de Bernoulli la dividimos por la función $a_1(x)\ne 0$, obtenemos

$${\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{a_0(x)}{a_1(x)}}y={\displaystyle \frac{g(x)}{a_1(x)}}{y}^n$$

Definimos las siguientes funciones.

$$\begin{array}{}\text{(2)}& P(x)={\displaystyle \frac{a_0(x)}{a_1(x)}}yQ(x)={\displaystyle \frac{g(x)}{a_1(x)}}\end{array}$$

Entonces una ecuación de Bernoulli se puede reescribir como

$$\begin{array}{}\text{(3)}& {\displaystyle \frac{dy}{dx}}+P(x)y=Q(x)y^n\end{array}$$

La ecuación ($\text{3}$) es también una definición común de ecuación de Bernoulli.

Notemos que si $n=0$, la ecuación de Bernoulli se reduce a una ecuación diferencial lineal no homogénea.

$$\begin{array}{}\text{(4)}& {\displaystyle \frac{dy}{dx}}+P(x)y=Q(x)\end{array}$$

Y si $n=1$, la ecuación de Bernoulli se reduce a una ecuación diferencial lineal homogénea.

$$\begin{array}{rl}{\displaystyle \frac{dy}{dx}}+P(x)y& =Q(x)y\\ {\displaystyle \frac{dy}{dx}}+[P(x)-Q(x)]y& =0\end{array}$$

Si definimos

$$R(x)=P(x)-Q(x)$$

entonces

$$\begin{array}{}\text{(5)}& {\displaystyle \frac{dy}{dx}}+R(x)y=0\end{array}$$

Las ecuaciones ($\text{4}$) y ($\text{5}$) ya las sabemos resolver.

Nuestro objetivo será resolver la ecuación de Bernoulli para el caso en el que $n\ne 0$ y $n\ne 1$.

Una propiedad de las ecuaciones de Bernoulli es que la sustitución

$$\begin{array}{}\text{(6)}& u(x)=y^{1-n}\end{array}$$

la convierte en una ecuación lineal, de tal manera que podremos resolverla usando algún método de resolución visto para ecuaciones diferenciales lineales.

Consideremos la ecuación de Bernoulli en la forma ($\text{3}$).

$${\displaystyle \frac{dy}{dx}}+P(x)y=Q(x)y^n$$

Dividimos toda la ecuación por $y^n\ne 0$.

$$\begin{array}{}\text{(7)}& {\displaystyle \frac{1}{y^n}}{\displaystyle \frac{dy}{dx}}+P(x)y^{1-n}=Q(x)\end{array}$$

La derivada de la función ($\text{6}$) es

$${\displaystyle \frac{du}{dx}}=(1-n)y^{-n}{\displaystyle \frac{dy}{dx}}=(1-n){\displaystyle \frac{1}{y^n}}{\displaystyle \frac{dy}{dx}}$$

de donde,

$$\begin{array}{}\text{(8)}& {\displaystyle \frac{1}{y^n}}{\displaystyle \frac{dy}{dx}}={\displaystyle \frac{1}{1-n}}{\displaystyle \frac{du}{dx}}\end{array}$$

Sustituyamos ($\text{6}$) y ($\text{8}$) en la ecuación ($\text{7}$).

$$\begin{array}{}\text{(9)}& {\displaystyle \frac{1}{1-n}}{\displaystyle \frac{du}{dx}}+P(x)u=Q(x)\end{array}$$

Multipliquemos por $1-n$ en ambos lados de la ecuación.

$${\displaystyle \frac{du}{dx}}+(1-n)P(x)u=(1-n)Q(x)$$

Definimos las funciones

$$R(x)=(1-n)P(x)yS(x)=(1-n)Q(x)$$

En términos de estas funciones la ecuación ($\text{9}$) se puede escribir de la siguiente forma.

$$\begin{array}{}\text{(10)}& {\displaystyle \frac{du}{dx}}+R(x)u=S(x)\end{array}$$

Este resultado corresponde a una ecuación diferencial lineal de primer orden no homogénea y, por tanto, puede ser resuelta aplicando el algoritmo descrito para resolver ecuaciones diferenciales lineales.

Los pasos que se recomiendan seguir para resolver una ecuación diferencial de Bernoulli se presentan a continuación.

Método para resolver ecuaciones de Bernoulli

1. El primer paso es escribir a la ecuación de Bernoulli en la forma ($\text{3}$).

2. Dividimos toda la ecuación por $y^n$ y consideramos el cambio de variable $u=y^{1-n}$, así como la respectiva derivada $${\displaystyle \frac{du}{dx}}=(1-n){\displaystyle \frac{1}{y^n}}{\displaystyle \frac{dy}{dx}}$$

3. Sustituimos $$y^{1-n}=uy{\displaystyle \frac{1}{y^n}}{\displaystyle \frac{dy}{dx}}={\displaystyle \frac{1}{1-n}}{\displaystyle \frac{du}{dx}}$$ en la ecuación resultante del paso anterior y haciendo un poco de álgebra podremos reducir la ecuación de Bernoulli en una ecuación lineal de primer orden no homogénea.

4. Resolvemos la ecuación resultante usando el método de resolución de ecuaciones diferenciales lineales lo que nos permitirá obtener la función $u(x)$.

5. Regresamos a la variable original para obtener finalmente la solución $y(x)$.

Realicemos un ejemplo en el que apliquemos estos pasos.

Ejemplo: Resolver la ecuación de Bernoulli

$$3(1+x^2){\displaystyle \frac{dy}{dx}}=2xy(y^3-1)$$

Solución: El primer paso es escribir la ecuación de Bernoulli en la forma ($\text{3}$).

$$\begin{array}{rl}3(1+x^2){\displaystyle \frac{dy}{dx}}& =2xy(y^3-1)\\ {\displaystyle \frac{dy}{dx}}& ={\displaystyle \frac{2xy(y^3-1)}{3(1+x^2)}}\\ {\displaystyle \frac{dy}{dx}}& ={\displaystyle \frac{2x{y}^4}{3(1+x^2)}}-{\displaystyle \frac{2xy}{3(1+x^2)}}\\ {\displaystyle \frac{dy}{dx}}+\left({\displaystyle \frac{2x}{3(1+x^2)}}\right)y& =\left({\displaystyle \frac{2x}{3(1+x^2)}}\right)y^4\end{array}$$

La última relación muestra a la ecuación en la forma ($\text{3}$) con $n=4$, ahora dividamos toda la ecuación por $y^4$.

$$\begin{array}{}\text{(11)}& {\displaystyle \frac{1}{y^4}}{\displaystyle \frac{dy}{dx}}+\left({\displaystyle \frac{2x}{3(1+x^2)}}\right)y^{-3}={\displaystyle \frac{2x}{3(1+x^2)}}\end{array}$$

Consideremos la sustitución

$$u=y^{1-n}=y^{1-4}=y^{-3}={\displaystyle \frac{1}{y^3}}$$

y

$${\displaystyle \frac{du}{dx}}=-3y^{-4}{\displaystyle \frac{dy}{dx}}$$

De donde,

$${\displaystyle \frac{1}{y^4}}{\displaystyle \frac{dy}{dx}}=-{\displaystyle \frac{1}{3}}{\displaystyle \frac{du}{dx}}y{y}^{-3}=u$$

Sustituimos estos resultados en la ecuación ($\text{11}$).

$$\begin{array}{rl}-{\displaystyle \frac{1}{3}}{\displaystyle \frac{du}{dx}}+\left({\displaystyle \frac{2x}{3(1+x^2)}}\right)u& ={\displaystyle \frac{2x}{3(1+x^2)}}\\ \text{(12)}& {\displaystyle \frac{du}{dx}}+(-{\displaystyle \frac{2x}{1+x^2}})u& =-{\displaystyle \frac{2x}{1+x^2}}\end{array}$$

La última ecuación es una expresión en la forma ($\text{10}$). Con esto hemos logrado reducir la ecuación de Bernoulli en una ecuación diferencial lineal de primer orden no homogénea.

Establecemos las siguientes funciones.

$$R(x)=-{\displaystyle \frac{2x}{1+x^2}}yS(x)=-{\displaystyle \frac{2x}{1+x^2}}$$

A partir de aquí aplicamos el método de resolución de ecuaciones diferenciales lineales.

La ecuación ya se encuentra en su forma canónica. Determinemos el factor integrante dado por

$$\begin{array}{}\text{(13)}& \mu (x)=e^{\int R(x)dx}\end{array}$$

Resolvamos la integral del exponente omitiendo la constante de integración.

$$\begin{array}{rl}\int R(x)dx& =-\int {\displaystyle \frac{2x}{1+x^2}}dx\\ & =-\ln |1+x^2|\end{array}$$

Por lo tanto,

$$\mu (x)=e^{-\ln |1+x^2|}={\displaystyle \frac{1}{1+x^2}}$$

Multipliquemos a la ecuación ($\text{12}$) por el factor integrante.

$${\displaystyle \frac{1}{1+x^2}}{\displaystyle \frac{du}{dx}}-{\displaystyle \frac{1}{1+x^2}}\left({\displaystyle \frac{2x}{1+x^2}}\right)u=-{\displaystyle \frac{1}{1+x^2}}\left({\displaystyle \frac{2x}{1+x^2}}\right)$$

Identificamos que el lado izquierdo de la ecuación es la derivada del producto del factor integrante $\mu (x)$ por la función $u(x)$, de esta manera

$${\displaystyle \frac{d}{dx}}\left({\displaystyle \frac{u}{1+x^2}}\right)=-{\displaystyle \frac{2x}{(1+x^2{)}^2}}$$

Integramos ambos lados de la ecuación con respecto a $x$. Por tratarse del último paso sí consideramos a la constante de integración.

$$\int {\displaystyle \frac{d}{dx}}\left({\displaystyle \frac{u}{1+x^2}}\right)dx=-\int {\displaystyle \frac{2x}{(1+x^2{)}^2}}dx$$

En el lado izquierdo aplicamos el teorema fundamental del cálculo y en el lado derecho consideramos la sustitución $a(x)=1+x^2$ para resolver la integral. El resultado que se obtiene es

$$\begin{array}{rl}{\displaystyle \frac{u}{1+x^2}}& ={\displaystyle \frac{1}{1+x^2}}+c\\ u& =1+(1+x^2)c\end{array}$$

Regresamos a la variable original $u=y^{-3}$.

$${\displaystyle \frac{1}{y^3}}=1+(1+x^2)c$$

Por lo tanto, la solución general (implícita) de la ecuación diferencial de Bernoulli

$$3(1+x^2){\displaystyle \frac{dy}{dx}}=2xy(y^3-1)$$

es

$$y^3(x)={\displaystyle \frac{1}{1+(1+x^2)c}}$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Ahora revisemos la ecuación de Riccati.

Ecuación diferencial de Riccati

La ecuación de Riccati es una ecuación diferencial ordinara no lineal de primer orden, inventada y desarrollada en el siglo XVlll por el matemático italiano Jacopo Francesco Riccati.

Resolver la ecuación de Riccati requiere del conocimiento previo de una solución particular de la ecuación, llamemos a dicha solución $\hat{y}(x)$. Si hacemos la sustitución

$$\begin{array}{}\text{(15)}& y(x)=\hat{y}(x)+u(x)\end{array}$$

La ecuación de Riccati adquiere la forma de una ecuación de Bernoulli, de tarea moral comprueba este hecho. Ya vimos que para resolver una ecuación de Bernoulli debemos reducirla a una ecuación lineal no homogénea, así que veamos directamente cómo reducir una ecuación de Riccati a una ecuación lineal no homogénea.

Sea $\hat{y}(x)$ una solución particular de la ecuación de Riccati y consideremos la sustitución

$$\begin{array}{}\text{(16)}& y(x)=\hat{y}(x)+{\displaystyle \frac{1}{u(x)}}\end{array}$$

Derivemos esta ecuación.

$$\begin{array}{}\text{(17)}& {\displaystyle \frac{dy}{dx}}={\displaystyle \frac{d\hat{y}}{dx}}-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\end{array}$$

Como $\hat{y}(x)$ es una solución de la ecuación de Riccati, entonces satisface la ecuación diferencial.

$$\begin{array}{}\text{(18)}& {\displaystyle \frac{d\hat{y}}{dx}}=q_0(x)+q_1(x)\hat{y}+q_2(x){\hat{y}}^2\end{array}$$

Sustituyendo ($\text{18}$) en ($\text{17}$) obtenemos la siguiente ecuación.

$$\begin{array}{}\text{(19)}& {\displaystyle \frac{dy}{dx}}=q_0(x)+q_1(x)\hat{y}+q_2(x){\hat{y}}^2-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\end{array}$$

Ahora podemos igualar la ecuación ($\text{19}$) con la ecuación de Riccati ($\text{14}$).

$$\begin{array}{rl}{q}_0(x)+q_1(x)y+q_2(x)y^2& =q_0(x)+q_1(x)\hat{y}+q_2(x){\hat{y}}^2-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\\ q_1(x)y+q_2(x)y^2& =q_1(x)\hat{y}+q_2(x){\hat{y}}^2-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\\ {\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}& =q_1(x)\hat{y}-q_1(x)y+q_2(x){\hat{y}}^2-q_2(x)y^2\\ {\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}& =q_1(x)(\hat{y}-y)+q_2(x)({\hat{y}}^2-y^2)\end{array}$$

En la última relación sustituimos la función ($\text{16}$).

$$\begin{array}{rl}{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}& =q_1(x)[\hat{y}-(\hat{y}+{\displaystyle \frac{1}{u}})]+q_2(x)[{\hat{y}}^2-{(\hat{y}+{\displaystyle \frac{1}{u}})}^2]\\ & =q_1(x)(\hat{y}-\hat{y}-{\displaystyle \frac{1}{u}})+q_2(x)({\hat{y}}^2-{\hat{y}}^2-2\hat{y}{\displaystyle \frac{1}{u}}-{\displaystyle \frac{1}{u^2}})\\ & =q_1(x)(-{\displaystyle \frac{1}{u}})+q_2(x)(-2{\displaystyle \frac{\hat{y}}{u}}-{\displaystyle \frac{1}{u^2}})\\ & =-{\displaystyle \frac{q_1(x)}{u}}-2q_2(x){\displaystyle \frac{\hat{y}}{u}}-{\displaystyle \frac{q_2(x)}{u^2}}\end{array}$$

Esto es,

$${\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}=-{\displaystyle \frac{q_1(x)}{u}}-2q_2(x){\displaystyle \frac{\hat{y}}{u}}-{\displaystyle \frac{q_2(x)}{u^2}}$$

Multipliquemos ambos lados de la ecuación por $u^2$.

$$\begin{array}{rl}{\displaystyle \frac{du}{dx}}& =-q_1(x)u-2q_2(x)\hat{y}u-q_2(x)\\ {\displaystyle \frac{du}{dx}}& =-(q_1(x)+2q_2(x)\hat{y})u-q_2(x)\end{array}$$

Vemos que

$$\begin{array}{}\text{(20)}& {\displaystyle \frac{du}{dx}}+(q_1(x)+2q_2(x)\hat{y})u=-q_2(x)\end{array}$$

Definamos las funciones

$$R(x)=q_1(x)+2q_2(x)\hat{y}(x)yS(x)=-q_2(x)$$

Por lo tanto, la ecuación ($\text{20}$) queda de la siguiente forma.

$$\begin{array}{}\text{(21)}& {\displaystyle \frac{du}{dx}}+R(x)u=S(x)\end{array}$$

Queda demostrado que la sustitución ($\text{16}$) convierte a la ecuación de Riccati en una ecuación diferencial lineal y, por tanto, puede ser resuelta con el método de resolución de ecuaciones lineales.

Como es usual, enunciemos la serie de pasos que se recomienda seguir para resolver las ecuaciones diferenciales de Riccati.

Método para resolver ecuaciones de Riccati

1. El primer paso es escribir a la ecuación de Riccati en la forma ($\text{14}$) y estar seguros de que conocemos previamente una solución particular $\hat{y}(x)$ de la ecuación.

2. Como queremos reducir la ecuación de Riccati en una ecuación lineal no homogénea consideramos la sustitución $$y(x)=\hat{y}(x)+{\displaystyle \frac{1}{u(x)}}$$ con $\hat{y}(x)$ la solución particular dada.
   
   Si se deseara reducirla a una ecuación de Bernoulli se hace la sustitución $$y(x)=\hat{y}(x)+u(x)$$

3. Debido a que $\hat{y}(x)$ es solución de la ecuación de Riccati, el siguiente paso es derivar la sustitución $y=\hat{y}+{\displaystyle \frac{1}{u}}$ y en el resultado sustituir ${\displaystyle \frac{d\hat{y}}{dx}}$ por la ecuación de Riccati para la solución particular, esto es

$${\displaystyle \frac{dy}{dx}}={\displaystyle \frac{d\hat{y}}{dx}}-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}=[q_1(x)+q_2(x)\hat{y}+q_3(x){\hat{y}}^2]-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}$$

4. Igualamos la ecuación anterior con la ecuación de Riccati original en la forma ($\text{14}$) y hacemos la sustitución $$y(x)=\hat{y}(x)+{\displaystyle \frac{1}{u(x)}}$$

5. Hecho lo anterior y haciendo un poco de álgebra podremos reducir la ecuación de Riccati en una ecuación lineal de primer orden y así aplicar el método de resolución para este tipo de ecuaciones.

6. Una vez obtenida la función $u(x)$ la sustituimos en $y(x)$ para obtener la solución deseada.

Realicemos un ejemplo para poner en practica este método.

Ejemplo: Resolver la ecuación de Riccati

$${\displaystyle \frac{dy}{dx}}=-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{y}{x}}+y^2$$

dada la solución particular $\hat{y}={\displaystyle \frac{2}{x}}$.

Solución: La ecuación diferencial prácticamente se encuentra en la forma de la ecuación ($\text{14}$), sólo para que sea claro escribimos

$${\displaystyle \frac{dy}{dx}}=(-{\displaystyle \frac{4}{x^2}})+(-{\displaystyle \frac{1}{x}})y+y^2$$

Comencemos por verificar que la solución particular dada efectivamente satisface la ecuación de Riccati. Por un lado,

$${\displaystyle \frac{d\hat{y}}{dx}}=-{\displaystyle \frac{2}{x^2}}$$

Por otro lado,

$$\begin{array}{rl}-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{\hat{y}}{x}}+{\hat{y}}^2& =-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{1}{x}}\left({\displaystyle \frac{2}{x}}\right)+{\left({\displaystyle \frac{2}{x}}\right)}^2\\ & =-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{4}{x^2}}\\ & =-{\displaystyle \frac{2}{x^2}}\end{array}$$

En efecto,

$${\displaystyle \frac{d\hat{y}}{dx}}=-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{\hat{y}}{x}}+{\hat{y}}^2=-{\displaystyle \frac{2}{x^2}}$$

El siguiente paso es hacer la sustitución ($\text{16}$).

$$y(x)=\hat{y}(x)+{\displaystyle \frac{1}{u(x)}}={\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{u}}$$

De acuerdo a ($\text{19}$), tenemos

$${\displaystyle \frac{dy}{dx}}=-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{1}{x}}\left({\displaystyle \frac{2}{x}}\right)+{\left({\displaystyle \frac{2}{x}}\right)}^2-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}$$

Igualemos este resultado con la ecuación de Riccati original.

$$\begin{array}{rl}-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{y}{x}}+y^2& =-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\\ -{\displaystyle \frac{y}{x}}+y^2& ={\displaystyle \frac{2}{x^2}}-{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}\\ {\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}& ={\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{y}{x}}-y^2\end{array}$$

En la última ecuación sustituimos $y={\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{u}}$.

$$\begin{array}{rl}{\displaystyle \frac{1}{u^2}}{\displaystyle \frac{du}{dx}}& ={\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{1}{x}}({\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{u}})-{({\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{u}})}^2\\ & ={\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{2}{x^2}}+{\displaystyle \frac{1}{xu}}-({\displaystyle \frac{4}{x^2}}+{\displaystyle \frac{4}{xu}}+{\displaystyle \frac{1}{u^2}})\\ & ={\displaystyle \frac{4}{x^2}}+{\displaystyle \frac{1}{xu}}-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{4}{xu}}-{\displaystyle \frac{1}{u^2}}\\ & =-{\displaystyle \frac{3}{xu}}-{\displaystyle \frac{1}{u^2}}\end{array}$$

De donde,

$${\displaystyle \frac{du}{dx}}+{\displaystyle \frac{3}{x}}u=-1$$

Esta expresión tiene la forma de una ecuación diferencial lineal ($\text{21}$), de donde podemos determinar que

$$R(x)={\displaystyle \frac{3}{x}}yS(x)=-1$$

La ecuación de Riccati ha sido reducida a una ecuación lineal no homogénea, ahora apliquemos el método de resolución de ecuaciones diferenciales lineales.

Calculemos el factor integrante $\mu (x)=e^{\int R(x)dx}$.

$$\int R(x)dx=\int {\displaystyle \frac{3}{x}}dx=3\ln |x|$$

El factor integrante es

$$\mu (x)=e^{3\ln |x|}=x^3$$

Multipliquemos la ecuación diferencial por el factor integrante.

$$\begin{array}{rl}{x}^3{\displaystyle \frac{du}{dx}}+x^3\left({\displaystyle \frac{3}{x}}\right)u& =-x^3\\ x^3{\displaystyle \frac{du}{dx}}+3x^{2}u& =-x^3\end{array}$$

Identificamos que el lado izquierdo de la ecuación corresponde a la derivada del producto entre el factor integrante $\mu (x)$ y la función $u(x)$, entonces

$${\displaystyle \frac{d}{dx}}\left(x^{3}u\right)=-x^3$$

Integramos ambos lados de la ecuación con respecto a $x$.

$$\begin{array}{rl}\int {\displaystyle \frac{d}{dx}}\left(x^{3}u\right)dx& =\int -x^{3}dx\\ x^{3}u& =-{\displaystyle \frac{x^4}{4}}+c\\ u(x)& =-{\displaystyle \frac{x}{4}}+{\displaystyle \frac{c}{x^3}}\end{array}$$

Ya determinamos el valor de $u(x)$, ahora sólo lo sustituimos en la función $y={\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{u}}$.

Por lo tanto, la solución general de la ecuación de Bernoulli

$${\displaystyle \frac{dy}{dx}}=-{\displaystyle \frac{4}{x^2}}-{\displaystyle \frac{y}{x}}+y^2$$

es

$$y(x)={\displaystyle \frac{2}{x}}+{\displaystyle \frac{1}{{\displaystyle \frac{c}{x^3}}-{\displaystyle \frac{x}{4}}}}={\displaystyle \frac{2}{x}}+{\displaystyle \frac{4x^3}{4c-x^4}}$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Hemos concluido con el estudio de las ecuaciones diferenciales de primer orden.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

1. Resolver las siguientes ecuaciones de Bernoulli.

• ${\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{1}{x}}y={\displaystyle \frac{2}{3}}{x}^4{y}^4$

• $3x{\displaystyle \frac{dy}{dx}}-2y=x^3{y}^{-2}$

• $x^2{\displaystyle \frac{dy}{dx}}-2xy=3y^4$ con la condición inicial $y(1)={\displaystyle \frac{1}{2}}$

2. Resolver las siguientes ecuaciones de Riccati.

• $x^3{\displaystyle \frac{dy}{dx}}=x^4{y}^2-2x^{2}y-1$ con solución particular $\hat{y}={\displaystyle \frac{1}{x^2}}$

• ${\displaystyle \frac{dy}{dx}}=x{y}^2+y+{\displaystyle \frac{1}{x^2}}$ con solución particular $\hat{y}=-{\displaystyle \frac{1}{x}}$

3. Demostrar que la sustitución $$y(x)=\hat{y}(x)+u(x)$$ convierte a una ecuación de Riccati en una ecuación de Bernoulli. $\hat{y}(x)$ es una solución particular de la ecuación de Riccati.

Más adelante…

Con esta entrada concluimos el estudio de las ecuaciones diferenciales de primer orden, a lo largo de la unidad vimos una descripción cualitativa y posteriormente una descripción analítica en la que desarrollamos varios métodos para resolver ecuaciones diferenciales de primer orden tanto lineales como no lineales.

Antes de pasar a la siguiente unidad y comenzar con el estudio de las ecuaciones diferenciales de segundo orden, es importante hacer un estudio con mayor detalle sobre el teorema de existencia y unicidad ya que es este teorema el que justifica toda la teoría que hemos desarrollado a lo largo de la unidad.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales l
• Entrada anterior del curso: Ecuaciones diferenciales exactas
• Siguiente entrada del curso: Teorema de Existencia y Unicidad – Ecuación Integral, Funciones Lipschitzianas y Lema de Gronwall

• Video relacionado al tema: Ecuaciones de Bernoulli y Riccati

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»