Autor Tema: Solución aproximada de la regla de Euler

0 Usuarios y 1 Visitante están viendo este tema.

07 Diciembre, 2018, 07:21 am
Respuesta #10

manooooh

  • $$\Large \color{#9c57a6}\pi\,\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 4,791
  • País: ar
  • Karma: +1/-0
  • Sexo: Masculino
Hola

O sea, que la solución a la primera parte del ejercicio sería:

\( (x+\sqrt[ ]{x})y'=y \)              ecuación original

\( \displaystyle\frac{1}{y}y'=\displaystyle\frac{1}{x+\sqrt[ ]{x}} \)           reescribir en forma de ODE

\( \displaystyle\int_{}^{}\displaystyle\frac{1}{y}dy=\displaystyle\int_{}^{}\displaystyle\frac{1}{x+\sqrt[ ]{x}}dx \)          integrar cada lado de la ecuación

\( ln(y)+c2=2.ln(\sqrt[ ]{x}+1)+c1 \)

\( ln(y)=2.ln(\sqrt[ ]{x}+1)+c1 \)                          combinar constantes


Aplicar las condiciones iniciales:

\( ln(4)=2.ln(\sqrt[ ]{1}+1)+c1 \)                       \( y(1)=4 \)


Aislar c1:

\( ln(4)=2.ln(\sqrt[ ]{1}+1)+c1 \)

\( 2.ln(\sqrt[ ]{1}+1)+c1=ln(4) \)

\( 2.ln(2)+c1=ln(4) \)

\( 2.ln(2)+c1-2.ln(2)=ln(4)-2.ln(2) \)                 restar \( 2.ln(2) \) de ambos lados

\( 2.ln(2)-2.ln(2)=c1 \)

\( c1=0 \)


Sustituir \( c1=0 \) en \( y=2.ln(\sqrt[ ]{x}+1)+c1 \):

\( ln(y)=2.ln(\sqrt[ ]{x}+1)+c1 \)

\( ln(y)=2.ln(\sqrt[ ]{x}+1) \)                     \( c1=0 \)


Aislar y:

\( ln(y)=2.ln(\sqrt[ ]{x}+1) \)

\( ln(y)=ln((\sqrt[ ]{x}+1)^2) \)

\( y=(\sqrt[ ]{x}+1)^2 \)                         los logaritmos tienen la misma base


Comprobación: comprobar las soluciones aplicando las condiciones iniciales

\( y(1)=4 \)                 condición inicial dada en el enunciado

\( y=(\sqrt[ ]{x}+1)^2 \)

\( 4=(\sqrt[ ]{1}+1)^2 \)

\( 4=4 \)


Conclusión: la solución final para \( (x+\sqrt[ ]{x})y'=y, y(1)=4 \) es:

\( y=(\sqrt[ ]{x}+1)^2 \)

Bien. Siendo un poco más molesto, toda función se caracteriza por su dominio, codominio y regla de correspondencia.

Creo que ya solucioné la segunda parte del ejercicio.

No necesitamos el número de pasos.

Sí se necesita. Como dije, si no conocemos el número de pasos \( n \) no podemos saber cada cuánto nos tenemos que mover (desplazamiento a izquierda o a derecha).

Como ya tenemos la solución general: \( y(x)=(\sqrt[ ]{x}+1)^2 \), y nos piden la solución en x=1.5, solo hay que sustituir x=1.5 en la solución general:

Un pequeño matiz: al ser un PVI la solución general es la función con su constante arbitraria, mientras que en la solución particular se considera la constante como ya calculada.

\( y(x)=(\sqrt[ ]{x}+1)^2 \)             solución general

\( y(1.5)=(sqrt(1.5)+1)^2 \)       sustituir x=1.5 en la solución general

\( y(1.5)=4.9494\approx{5} \)    solución aproximada para x=1.5

Bien, pero ahí no has aplicado del todo el método de Euler, sino que sustituiste directamente.

La Wikipedia lo dice así:

\( \begin{array}{|c|c|c|}
\hline\dfrac{dy}{dx}&x&y\\\hline &x_0&y_0\\\hline f(x_0,y_0)=\dfrac{y_0}{x_0+\sqrt{x_0}}&x_1=x_0+h&y_1=y_0+h\cdot f(x_0,y_0)\\\hline f(x_1,y_1)=\dfrac{y_1}{x_1+\sqrt{x_1}}&x_2=x_1+h&y_2=y_1+h\cdot f(x_1,y_1)\\\hline f(x_2,y_2)=\dfrac{y_2}{x_2+\sqrt{x_2}}&x_3=x_2+h&y_3=y_2+h\cdot f(x_2,y_2)\\\hline
\end{array} \)

Consideremos

\( (x+\sqrt{x})y'=y\implies y'=\dfrac{dy}{dx}=f(x,y)=\dfrac y{x+\sqrt x} \)

con \( n=1 \) como proponés (nos desplazamos una unidad a derecha, podríamos haber elegido \( n=-1 \)), entonces

\( h=\dfrac{x_1-x_0}n=\dfrac{1.5-1}1=0.5 \).

Ahora hay que armar la tabla:

\( \begin{array}{|c|c|c|}
\hline\dfrac{dy}{dx}&x&y\\\hline &1&4\\\hline f(1,4)=\dfrac{4}{1+\sqrt{1}}=2&x_1=1+0.5=1.5&y_1=4+0.5\cdot f(1,4)=\boxed5\\\hline f(1.5,5)=\dfrac{5}{1.5+\sqrt{1.5}}=1.83&x_2=1.5+0.5=2&y_2=5+0.5\cdot f(1.5,5)=5.92\\\hline
\end{array} \)

Como el número de pasos es "chico", podemos aproximar mejor el resultado aumentando la cantidad de pasos, como por ejemplo \( n=5 \). Luego \( h=\dfrac{1.5-1}5=0.1 \). Ahora la tabla:

\( \begin{array}{|c|c|c|}
\hline\dfrac{dy}{dx}&x&y\\\hline &1&4\\\hline f(1,4)=2&x_1=1.1&y_1=4.2\\\hline f(1.1,4.2)=1.95&x_2=1.2&y_2=4.4\\\hline f(1.2,4.4)=1.97&x_3=1.3&y_3=4.6\\\hline f(1.3,4.6)=1.89&x_4=1.4&y_4=4.79\\\hline f(1.4,4.79)=1.85&x_5=1.5&y_5=\boxed{4.98\approx4.9494\dots}\\\hline
\end{array} \)

Es decir, a medida que aumentamos el número de pasos estamos más cerca del resultado.

Pero ya te digo, no conocemos \( n \).

Saludos

P.D. De todas maneras, la ecuación diferencial tiene solución exacta. Es raro que te pidan utilizar el método de Euler; eso se usa para EDOs no resolubles, aunque por ahí podrían habértelo mandado como práctica :laugh:.

09 Diciembre, 2018, 02:08 am
Respuesta #11

mgranadosgg

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 131
  • País: co
  • Karma: +0/-0
  • Sexo: Masculino
Hola. Gracias de antemano.

Es un ejercicio de un examen. Quizá por eso despiste un poco.

Tal vez la sustitución que yo propongo funciona por casualidad, pero cuando vi la solución se me ocurrió que como tiene forma parecida a la de una recta pasando los primeros valores, podría servir simplemente la sustitución del valor. Por otro lado, también veo que si la solución del ejercicio hubiera sido más pequeña seguramente no habría funcionado la simple sustitución.

Supongo, a propósito de lo que comentas, que en la respuesta se podría argumentar que como el enunciado no dice nada nosotros decidimos "n", haciendo, por ejemplo, 5 pasos y realizando así el método de Euler completamente, quedando acreditado que lo sabemos utilizar.

Quizá habría alguna forma de explicar esa sustitución como un añadido tras haber resuelto el ejercicio correctamente o como una posible solución alternativa.

En todo caso, muchas gracias de nuevo.

09 Diciembre, 2018, 06:23 am
Respuesta #12

mgranadosgg

  • $$\Large \color{#5e8d56}\pi\,\pi\,\pi$$
  • Mensajes: 131
  • País: co
  • Karma: +0/-0
  • Sexo: Masculino
Hola. Gracias de antemano.

propongo 2 soluciones para el método de Euler: como en el enunciado no especifican nada, entonces sugiero que

 solución 1.- o bien, el número de pasos sea n=1.
 solución 2.- o bien, el tamaño de paso sea h=0.1, por ejemplo.

solución 1
-----------
El problema sería: Encontrar \( y_(1.5) \) para \( y'=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \), cuando \( y_(1)=4, n=1 \).

El método de Euler establece que

 \( y_n+1=y_ n+h*f(x,y) \), cuando \( x_n+1=x_n+h \).

Tenemos que \( h=0.5, x_0=1, y_0=4, f(x,y)=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \).
Entonces sería,

Paso 1:

 \( x_1=x_0+h=1+0.5=1.5 \)

 \( y_(x_1)=y_(1.5)=y_1=y_0+h*f(x_0,y_0) \)
                                 \( = 4+h*f(1,4) \)
                                 \( = 4+0.5*(2) \)
                                 \( = 5 \)

Por tanto, \( y_(1.5)=5 \)


solución 2
-----------
El problema sería: Encontrar \( y_(1.5) \) para \( y'=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \), cuando \( y_(1)=4, h=0.1 \) (h=0.1, por ejemplo).
El método de Euler establece que

 \( y_n+1=y_ n+h*f(x,y) \), cuando \( x_n+1=x_n+h \).

Tenemos que \( h=0.1, x_0=1, y_0=4, f(x,y)=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \).

Entonces sería,

Paso 1:

\( x_1=x_0+h=1+0.1=1.1 \)
\( y_(x_1)=y*(1.1)=y_ 1=y_0+h*f(x_0,y_0) \)
                                 \( = 4+h*f(1,4) \)
                                 \( = 4+0.1*(2) \)
                                 \( = 4.2 \)

Paso 2:

\( x_2=x_1+h=1.1+0.1=1.2 \)
\( y_(x_2)=y*(1.2)=y_2=y_1+h*f(x_1,y_1) \)
                                 \( = 4.2+h*f(1.1,4.2) \)
                                 \( = 4.2+0.1*(2) \)
                                 \( = 4.3954 \)

Paso 3:

\( x_3=x_2+h=1.2+0.1=1.3 \)
\( y_(x_3)=y*(1.3)=y_3=y_2+h*f(x_2,y_2) \)
                                 \( = 4.3954+h*f(1.2,4.3954) \)
                                 \( = 4.3954+0.1*(1.9148) \)
                                 \( = 4.5869 \)

Paso 4:

\( x_4=x_3+h=1.3+0.1=1.4 \)
\( y_(x_4)=y*(1.4)=y_4=y_3+h*f(x_3,y_3) \)
                                 \( = 4.5869+h*f(1.3,4.5869) \)
                                 \( = 4.5869+0.1*(1.879) \)
                                 \( = 4.7749 \)

Paso 5:

\( x_5=x_4+h=1.4+0.1=1.5 \)
\( y_(x_5)=y*(1.5)=y_5=y_4+h*f(x_4,y_4) \)
                                 \( = 4.7749+h*f(1.4,4.7749) \)
                                 \( = 4.7749+0.1*(1.8484) \)
                                 \( = 4.9597 \)

Por tanto, \( y_(1.5)=4.9597 \)

No sé cuál sería la solución más correcta en un examen, la solución 1 o la solución 2. ¿Cuál preferirían?

09 Diciembre, 2018, 07:34 am
Respuesta #13

manooooh

  • $$\Large \color{#9c57a6}\pi\,\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 4,791
  • País: ar
  • Karma: +1/-0
  • Sexo: Masculino
Hola

Es un ejercicio de un examen. Quizá por eso despiste un poco.

No hay ningún despiste; si te piden hallar una solución aproximada sin el paso es como calcular un límite sin saber a qué tiende la variable (salvo que en este último podemos "intuir" en dónde nos piden calcular, ya que casi siempre son indeterminaciones, pero en este peor, no sabemos cuántas iteraciones tenemos que hacer).

Tal vez la sustitución que yo propongo funciona por casualidad, pero cuando vi la solución se me ocurrió que como tiene forma parecida a la de una recta pasando los primeros valores, podría servir simplemente la sustitución del valor. Por otro lado, también veo que si la solución del ejercicio hubiera sido más pequeña seguramente no habría funcionado la simple sustitución.

¿A qué sustitución te referís, a reemplazar \( x=1.5 \) en la solución particular \( y_P(x)=(\sqrt x+1)^2 \)? Pues eso es lo que deberías hacer si no te piden calcular el valor por el método de Euler; sustituir directamente y dejarlo expresado como \( (\sqrt{1.5}+1)^2 \).

Supongo, a propósito de lo que comentas, que en la respuesta se podría argumentar que como el enunciado no dice nada nosotros decidimos "n", haciendo, por ejemplo, 5 pasos y realizando así el método de Euler completamente, quedando acreditado que lo sabemos utilizar.

Sí claro; pero podrías haber hecho uso de cualquier otro paso. Fijate que con \( n=5 \) la \( x_i \) crece a razón de \( 0.1 \); o sea que, si empezamos en \( x=1 \) (del PVI), luego de \( 5 \) iteraciones tendremos el valor esperado de \( y(x) \) cuando \( x=1.5 \). Vendría a ser el "mejor", porque si aumentás \( n \) vas a tener que realizar muchas más iteraciones.

Quizá habría alguna forma de explicar esa sustitución como un añadido tras haber resuelto el ejercicio correctamente o como una posible solución alternativa.

Si es como creí, no es ninguna novedad; o lo hacés por sustitución directa o por el método de Euler.

Propongo 2 soluciones para el método de Euler: como en el enunciado no especifican nada, entonces sugiero que

 solución 1.- o bien, el número de pasos sea n=1.
 solución 2.- o bien, el tamaño de paso sea h=0.1, por ejemplo.

solución 1
-----------
El problema sería: Encontrar \( y_(1.5) \) para \( y'=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \), cuando \( y_(1)=4, n=1 \).

El método de Euler establece que

 \( y_n+1=y_ n+h*f(x,y) \), cuando \( x_n+1=x_n+h \).

Tenemos que \( h=0.5, x_0=1, y_0=4, f(x,y)=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \).
Entonces sería,

Paso 1:

 \( x_1=x_0+h=1+0.5=1.5 \)

 \( y_(x_1)=y_(1.5)=y_1=y_0+h*f(x_0,y_0) \)
                                 \( = 4+h*f(1,4) \)
                                 \( = 4+0.5*(2) \)
                                 \( = 5 \)

Por tanto, \( y_(1.5)=5 \)


solución 2
-----------
El problema sería: Encontrar \( y_(1.5) \) para \( y'=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \), cuando \( y_(1)=4, h=0.1 \) (h=0.1, por ejemplo).
El método de Euler establece que

 \( y_n+1=y_ n+h*f(x,y) \), cuando \( x_n+1=x_n+h \).

Tenemos que \( h=0.1, x_0=1, y_0=4, f(x,y)=\displaystyle\frac{y}{x+\sqrt[ ]{x}} \).

Entonces sería,

Paso 1:

\( x_1=x_0+h=1+0.1=1.1 \)
\( y_(x_1)=y*(1.1)=y_ 1=y_0+h*f(x_0,y_0) \)
                                 \( = 4+h*f(1,4) \)
                                 \( = 4+0.1*(2) \)
                                 \( = 4.2 \)

Paso 2:

\( x_2=x_1+h=1.1+0.1=1.2 \)
\( y_(x_2)=y*(1.2)=y_2=y_1+h*f(x_1,y_1) \)
                                 \( = 4.2+h*f(1.1,4.2) \)
                                 \( = 4.2+0.1*(2) \)
                                 \( = 4.3954 \)

Paso 3:

\( x_3=x_2+h=1.2+0.1=1.3 \)
\( y_(x_3)=y*(1.3)=y_3=y_2+h*f(x_2,y_2) \)
                                 \( = 4.3954+h*f(1.2,4.3954) \)
                                 \( = 4.3954+0.1*(1.9148) \)
                                 \( = 4.5869 \)

Paso 4:

\( x_4=x_3+h=1.3+0.1=1.4 \)
\( y_(x_4)=y*(1.4)=y_4=y_3+h*f(x_3,y_3) \)
                                 \( = 4.5869+h*f(1.3,4.5869) \)
                                 \( = 4.5869+0.1*(1.879) \)
                                 \( = 4.7749 \)

Paso 5:

\( x_5=x_4+h=1.4+0.1=1.5 \)
\( y_(x_5)=y*(1.5)=y_5=y_4+h*f(x_4,y_4) \)
                                 \( = 4.7749+h*f(1.4,4.7749) \)
                                 \( = 4.7749+0.1*(1.8484) \)
                                 \( = 4.9597 \)

Por tanto, \( y_(1.5)=4.9597 \)

¡Es lo que escribí en la respuesta #10, decimales más decimales menos!

No sé cuál sería la solución más correcta en un examen, la solución 1 o la solución 2. ¿Cuál preferirían?

Pues si tenemos que resolver un enunciado incompleto es subjetivo ;). Para no "terminar demasiado rápido", yo elegiría con \( n=5 \), aunque es una pérdida de tiempo cuando podemos hallar la respuesta sustituyendo directamente en la solución particular. Si el enunciado no está completo (como este caso) debés preguntarle a tu profesor.

Saludos