3 EDOs de orden 1

Despues de definir lo que son las EDOs, vamos a aprender tecnicas de resolucion de EDOs de primer orden.

Se escriben

  • ​

    En su forma explicita, y′=f(x,y)y^{\prime}=f(x,y)

  • ​

    En forma implicita, F(x,y,y′)=0F(x,y,y^{\prime})=0

  • ​

    En forma diferencial, M(x,y)dx+N(x,y)dy=0M(x,y)dx+N(x,y)dy=0

La funcion y(x)y(x) se supone diferenciable.

Ejemplo I.13.

La EDO

y′=x2+y2x−yy^{\prime}=\frac{x^{2}+y^{2}}{x-y}

esta escrita en forma explicita con f(x,y)=x2+y2x−yf(x,y)=\frac{x^{2}+y^{2}}{x-y}.

La forma implicita se obtiene agrupando todos los terminos en un mismo miembro

(x−y)y′−(x2+y2)=0(x-y)y^{\prime}-(x^{2}+y^{2})=0

con F(x,y,y′)=(x−y)y′−(x2+y2)=0F(x,y,y^{\prime})=(x-y)y^{\prime}-(x^{2}+y^{2})=0.

La forma diferencial se obtiene escribiendo la derivada y′=dydxy^{\prime}=\frac{dy}{dx}:

(x−y)dydx−(x2+y2)=0(x-y)\frac{dy}{dx}-(x^{2}+y^{2})=0
(x2+y2)dx−(x−y)dy=0(x^{2}+y^{2})\,dx-(x-y)\,dy=0

con

M(x,y)=x2+y2M(x,y)=x^{2}+y^{2}
N(x,y)=−(x−y)N(x,y)=-(x-y)

Una identidad util es la regla de la sustitucion. Consideramos las funciones continuas y diferenciables

  • ​

    HH funcion de yy: H(y)H(y)

  • ​

    hh su derivada: H′(y)=h(y)H^{\prime}(y)=h(y)

  • ​

    yy funcion de xx: y(x)y(x)

Aplicamos la regla de la cadena a HH:

ddxH(y(x))=dH(y)dydy(x)dx=H′(y)y′(x)=h(y)y′(x)\frac{d}{dx}H(y(x))=\frac{dH(y)}{dy}\frac{dy(x)}{dx}=H^{\prime}(y)y^{\prime}(x% )=h(y)y^{\prime}(x)

Entonces, si tenemos la integral ∫h(y)y′(x)dx\int h(y)y^{\prime}(x)dx podemos escribir

∫h(y)y′(x)dx=∫ddxH(y(x))dx=H(y)=∫h(y)dy\int h(y)y^{\prime}(x)dx=\int\frac{d}{dx}H(y(x))dx=H(y)=\int h(y)dy

Por lo tanto, obtenemos el resultado siguiente

∫h(y)y′(x)dx=∫h(y)dy\int h(y)y^{\prime}(x)dx=\int h(y)dy

Ahora consideramos la EDO de orden 11 en forma explicita

y′=g(x)h(y)y^{\prime}=\frac{g(x)}{h(y)}

donde y(x),g(x),h(y)y(x),g(x),h(y) son funciones continuas y diferenciables.

Aplicamos la sustitucion a la EDO

h(y)y′=g(x)h(y)y^{\prime}=g(x)

y, despues de integrar con respecto a xx,

∫h(y(x))y′(x)dx=∫g(x)dx+C\int h(y(x))y^{\prime}(x)dx=\int g(x)dx+C

Para obtener (por sustitucion del primer miembro)

∫h(y)dy=∫g(x)dx+C\int h(y)dy=\int g(x)dx+C

Conclusion: podemos integrar cada miembro de la EDO con respecto a su propia variable.

3.1 EDOs inmediatas

Las EDOs inmediatas se escriben y′=f(x)y^{\prime}=f(x). Son las EDOs mas sencillas, pues no aparece explicitamente la funcion yy. Tenemos

dydx=f(x)\frac{dy}{dx}=f(x)

Es el caso anterior con h(y)=1h(y)=1, con lo cual obtenemos

∫dy=∫f(x)dx+C⇒y=∫f(x)dx+C\int dy=\int f(x)dx+C\Rightarrow y=\int f(x)dx+C
Ejemplo I.14.
  • ​

    La EDO y′=cos(3x)+5y^{\prime}=\cos(3x)+5 es de orden 11, grado 11 y lineal. Ademas, viene escrita en forma explicita, con la particularidad de que es una EDO inmediata. Usando el resultado anterior,

    y=∫cos(3x)+5dx=∫cos(3x)dx+5∫dx=sin(3x)3+5x+Cy=\int\cos(3x)+5dx=\int\cos(3x)dx+5\int dx=\frac{\sin(3x)}{3}+5x+C

    por lo que ya hemos obtenido su solucion general (orden 1 - 1 parametro en la solucion general: OK).

  • ​

    La EDO y′=e2x−xy^{\prime}=e^{2x}-x. Resolvemos de nuevo calculando la integral.

    y=∫e2x−xdx=∫e2x−∫xdx=e2x2−x22+Cy=\int e^{2x}-xdx=\int e^{2x}-\int xdx=\frac{e^{2x}}{2}-\frac{x^{2}}{2}+C

    (orden 1 - 1 parametro en la solucion general: OK)

Este proceso de integracion para EDOs inmediatas de orden 11 se puede generalizar para EDOs de orden nn de la forma

y(n)=f(x)y^{(n)}=f(x)

como

y(n−1)(x)=∫f(x)dx+Cy^{(n-1)}(x)=\int f(x)dx+C

Podemos continuar un paso mas y obtener

y(n−2)(x)=∫(∫f(x)dx)dx+C1x+C2y^{(n-2)}(x)=\int\left(\int f(x)dx\right)dx+C_{1}x+C_{2}
y(n−3)(x)=∫∫∫f(x)dx+C12x2+C2x+C3y^{(n-3)}(x)=\int\int\int f(x)dx+\frac{C_{1}}{2}x^{2}+C_{2}x+C_{3}

Finalmente,

y(x)=∫⋯∫∫f(x)(dx)n+K1xn−1+K2xn−2+⋯+Kn−1x+Kny(x)=\int\cdots\int\int f(x)\,(dx)^{n}+K_{1}x^{n-1}+K_{2}x^{n-2}+\cdots+K_{n-1% }x+K_{n}

Asi, ya sabemos como calcular EDOs inmediatas de cualquier orden.

Ejemplo I.15 (Problema del objeto).

Tras plantear el problema, llegamos al sistema diferencial siguiente, definido por la EDO de orden 22 y 22 condiciones iniciales (problema de Cauchy (de valores iniciales)

{z′′(t)=−10z(0)=0z′(0)=40\begin{cases}z^{\prime\prime}(t)=-10\\ z(0)=0\\ z^{\prime}(0)=40\end{cases}

Con la variable independiente tt, la variable dependiente zz, la EDO es de orden 2 y grado 1, es lineal y en forma explícita.

Observamos que no aparece explícitamente la función zz en la EDO (es de la forma z′′(t)=f(t)z^{\prime\prime}(t)=f(t)): es una EDO inmediata pero de orden 2, por lo que hay que aplicar el proceso de resolución de EDOs inmediatas 2 veces (integrar 2 veces).

Aplicamos la formula, de modo que

z(t)=∫(∫−10dt)dt+C1t=∫−10tdt=−10∫tdt+C1t=−5t2+C1t+C2z(t)=\int\left(\int-10\;dt\right)\;dt+C_{1}t=\int-10tdt=-10\int tdt+C_{1}t=-5t% ^{2}+C_{1}t+C_{2}

es la solucion general de la EDO (orden 2 y 2 parametros, por lo que es coherente). Impongamos ahora las condiciones iniciales:

z(0)=−10⋅02+C1⋅0+C2=0⇒C2=0z(0)=-10\cdot\frac{0}{2}+C_{1}\cdot 0+C_{2}=0\Rightarrow C_{2}=0
z′(0)=−10⋅0+C1=40⇒C1=40z^{\prime}(0)=-10\cdot 0+C_{1}=40\Rightarrow C_{1}=40

y obtenemos que la solucion particular del problema es

z(t)=−5t2+40tz(t)=-5t^{2}+40t

Las respuestas a las preguntas iniciales son:

  • ​

    Posición del objeto a los 22 segundos:

    z(t=2)=−5⋅22+40⋅2=60z(t=2)=-5\cdot 2^{2}+40\cdot 2=60

    El objeto se encontrará a 60 m60\text{ m} de altura.

  • ​

    Altura maxima del objeto: punto estacionario de z(t)z(t), es decir, z′(t)=0z^{\prime}(t)=0 (velocidad 0) ⇒−10t+40=0⇒t=4\Rightarrow-10t+40=0\Rightarrow t=4.

    Entonces z(t=4)=−5⋅42+40⋅4=80z(t=4)=-5\cdot 4^{2}+40\cdot 4=80 m es la altura maxima.

3.2 EDOs de variables separables

Las EDOS de variables separables se escriben en forma diferencial de la siguiente forma

F(x)G(y)dx+H(x)P(y)dy=0F(x)G(y)dx+H(x)P(y)dy=0

Podemos separar las variables xx e yy dividiendo por H(x)G(y)H(x)G(y):

F(x)H(x)dx+P(y)G(y)dy=0⇒P(y)G(y)dy=−F(x)H(x)dx\frac{F(x)}{H(x)}dx+\frac{P(y)}{G(y)}dy=0\Rightarrow\frac{P(y)}{G(y)}dy=-\frac% {F(x)}{H(x)}dx

Procedemos como antes, integrando cada miembro con respecto a su propia variable:

∫P(y)G(y)dy=−∫F(x)H(x)dx+C\int\frac{P(y)}{G(y)}dy=-\int\frac{F(x)}{H(x)}dx+C

Esto es correcto, pero hemos excluido las posibles soluciones que existen a partir de H(x)G(y)=0H(x)G(y)=0. Debemos comprobar si existen soluciones de la EDO en las que este producto es nulo. Hay que sustituir las soluciones de H(x)G(y)=0H(x)G(y)=0 y comprobar si son soluciones de la EDO o no. Si lo son, se añaden a la solucion general.

Ejemplo I.16.
  • ​

    Consideremos la EDO por variables separables 2ydx+3xdy=02ydx+3xdy=0. Dividimos entre H(x)G(y)=3x⋅2y=6xyH(x)G(y)=3x\cdot 2y=6xy, de forma que se tiene

    13xdx+12ydy=0⇒12ydy=−13xdx\frac{1}{3x}dx+\frac{1}{2y}dy=0\Rightarrow\frac{1}{2y}dy=-\frac{1}{3x}dx

    Integramos cada miembro con respecto a su propia variable:

    ∫12ydy=−∫13xdx+C⇒12ln(y)=−13ln(x)+C\int\frac{1}{2y}dy=-\int\frac{1}{3x}dx+C\Rightarrow\frac{1}{2}\ln(y)=-\frac{1}% {3}\ln(x)+C

    Entonces

    ln(y1/2)=ln(x−1/3)+C⇒ln(y1/2)−ln(x−1/3)=C⇒\ln(y^{1/2})=\ln(x^{-1/3})+C\Rightarrow\ln(y^{1/2})-\ln(x^{-1/3})=C\Rightarrow
    ln(y1/2x−1/3)=Cln(y1/2x1/3)=C\ln\left(\frac{y^{1/2}}{x^{-1/3}}\right)=C\quad\ln(y^{1/2}x^{1/3})=C
    y1/2x1/3=eC=Knotación (K>0)y^{1/2}x^{1/3}=e^{C}=K\quad\text{notación }(K>0)
    y=(Kx1/3)2y=\left(\frac{K}{x^{1/3}}\right)^{2}
    y=K2x−2/3notacióny(x)=Kx−2/3y=K^{2}x^{-2/3}\quad\text{notación}\quad\boxed{y(x)=Kx^{-2/3}}

    Asi, ya hemos obtenido la solución general de la ED  (orden 1,1 parámetro eficiente)\text{ (orden }1,1\text{ parámetro eficiente)}.

    Hemos dividido entre 6xy6xy

    6xy=0⟹{x=0y=06xy=0\implies\begin{cases}x=0\\ y=0\end{cases}

    Hemos podido excluir las soluciones de la EDO que verifican x=0x=0 e y=0y=0.

    • •
      ​

      Caso y=0⟶y′=0y=0\longrightarrow y^{\prime}=0

      En la EDO (2y+3xy′=02y+3xy^{\prime}=0):

      2⋅0+3x⋅0\displaystyle 2\cdot 0+3x\cdot 0
      =0\displaystyle{}=0
      0+0\displaystyle 0+0
      =0\displaystyle{}=0
      0\displaystyle 0
      =0se verifica la EDO.\displaystyle{}=0\quad\text{se verifica la EDO.}

      Por tanto, y=0y=0 es otra solución de la EDO, que no puede incluirse en la solución general, porque K=eC>0K=e^{C}>0 (≠0\neq 0).

      (Es solución singular).

      Sin embargo, la incluimos forzando la posibilidad k=0k=0 en la solucion general (sabiendo que tiene otro origen).

    • •
      ​

      Caso x=0x=0: En la EDO (2y+3xy′=0)(2y+3xy^{\prime}=0):

      2y+3⋅0⋅y′\displaystyle 2y+3\cdot 0\cdot y^{\prime}
      =0\displaystyle{}=0
      2y+0\displaystyle 2y+0
      =0\displaystyle{}=0
      2y\displaystyle 2y
      =0\displaystyle{}=0
      y\displaystyle y
      =0⟹recaemos en el caso anterior y=0\displaystyle{}=0\implies\text{recaemos en el caso anterior }y=0

      Comprobación de la solución general:

      EDO: 2y+3xy′=02y+3xy^{\prime}=0

      y(x)=kx−2/3y(x)=kx^{-2/3}
      y′(x)=−23kx−5/3y^{\prime}(x)=\frac{-2}{3}kx^{-5/3}
      ⇒\displaystyle\Rightarrow
      2kx−2/3+3x−23kx−5/3=\displaystyle{}2kx^{-2/3}+3x\frac{-2}{3}kx^{-5/3}=
      =2kx−2/3−2kx−2/3=0¡OK!\displaystyle{}=2kx^{-2/3}-2kx^{-2/3}=0\quad\text{¡OK!}
Ejemplo I.17.

Vamos a resolver la EDO por variables separables

exdx−(1+ex)ydy=0e^{x}dx-(1+e^{x})ydy=0

Dividiendo entre 1+ex1+e^{x}, nos queda

ex1+exdx−ydy=0⇒ex1+ex=ydy\frac{e^{x}}{1+e^{x}}dx-ydy=0\Rightarrow\frac{e^{x}}{1+e^{x}}=ydy

Integramos cada miembro con respecto a su propia variable:

∫ex1+exdx+C=∫ydy⇒ln(1+ex)+C=y22\int\frac{e^{x}}{1+e^{x}}dx+C=\int ydy\Rightarrow\ln(1+e^{x})+C=\frac{y^{2}}{2}

Por tanto,

y=(2(ln(1+ex)+C))12y=(2(\ln(1+e^{x})+C))^{\frac{1}{2}}

y ya hemos obtenido la solucion general de la EDO (orden 1 y 1 parametro, por lo que es coherente).

Comentario: en realidad ln|1+ex⏟>0|=ln(1+ex)\ln|\underbrace{1+e^{x}}_{>0}|=\ln(1+e^{x}).

Hemos dividido entre 1+ex≠01+e^{x}\neq 0, por lo que no hemos excluido ninguna solucion.

Ejemplo I.18.

Resolvamos la EDO por variables separables

(x−4)y4dx−x3(y2−3)dy=0(x-4)y^{4}dx-x^{3}(y^{2}-3)dy=0

Dividiendo entre y4x3y^{4}x^{3}, obtenemos

x−4x3dx−y2−3y4dy=0⇒x−4x3dx=y2−3y4\frac{x-4}{x^{3}}dx-\frac{y^{2}-3}{y^{4}}dy=0\Rightarrow\frac{x-4}{x^{3}}dx=% \frac{y^{2}-3}{y^{4}}

Integramos cada miembro con respecto a su propia variable

∫x−4x3dx+C=∫y2−3y4dy\int\frac{x-4}{x^{3}}dx+C=\int\frac{y^{2}-3}{y^{4}}dy

Resolvemos cada integral:

∫x−4x3dx=∫x−2dx−4∫1x3dx=−1x+2x2\int\frac{x-4}{x^{3}}dx=\int x^{-2}dx-4\int\frac{1}{x^{3}}dx=-\frac{1}{x}+% \frac{2}{x^{2}}

y

∫y2−3y4dy=∫1y2dy−3∫1y4dy=−1y+1y3=1−y2y3\int\frac{y^{2}-3}{y^{4}}dy=\int\frac{1}{y^{2}}dy-3\int\frac{1}{y^{4}}dy=-% \frac{1}{y}+\frac{1}{y^{3}}=\frac{1-y^{2}}{y^{3}}

por lo que nos queda

−1x+2x2+1y−1y3=C-\frac{1}{x}+\frac{2}{x^{2}}+\frac{1}{y}-\frac{1}{y^{3}}=C

solucion general de la EDO expresada de forma implicita (orden 1 y 1 parametro, por lo que es coherente).

Hemos dividido entre x3y4x^{3}y^{4}

x3y4=0⟶{x=0y=0x^{3}y^{4}=0\longrightarrow\begin{cases}x=0\\ y=0\end{cases}

Hemos podido excluir las soluciones de la EDO que verifican x=0x=0 e y=0y=0.

x=0ey=0x=0\quad\text{e}\quad y=0

En el caso y=0⇒y′=0y=0\Rightarrow y^{\prime}=0 y en la EDO (x−4)y4−x3(y2−3)y′=0(x-4)y^{4}-x^{3}(y^{2}-3)y^{\prime}=0 se tiene:

(x−4)⋅04−x3(02−3)⋅0=0(x-4)\cdot 0^{4}-x^{3}(0^{2}-3)\cdot 0=0
0+0=00+0=0
0=0sí verifica la EDO0=0\quad\text{sí verifica la EDO}

Por tanto, y=0y=0 es otra solución de la EDO.

Por ultimo, en el caso x=0x=0,

(0−4)y4−03(y2−3)y′\displaystyle(0-4)y^{4}-0^{3}(y^{2}-3)y^{\prime}
=0\displaystyle{}=0
−4y4−0\displaystyle-4y^{4}-0
=0\displaystyle{}=0
−4y4\displaystyle-4y^{4}
=0\displaystyle{}=0

Caemos en el caso anterior y=0y=0.