7.

3 Mtodo de Euler mejorado

El valor de y.x/ en x D x0 C h se aproxima entonces por

Derivando la ED con respecto a x:

ww

w.

Ejemplo 7.3.1 Considere el PVI y 0 D 2x C 1, con y.0/ D 1. Encuentre una aproximacin cuadrtica de la solucin y.x/ en x D 1. y 00 D 2:

Usando la condicin inicial y evaluando las derivadas: y.0/ D 1; y 0 .0/ D 1; y 00 .0/ D 2:

En consecuencia, la aproximacin cuadrtica y2 .x/ de la solucin y.x/ para x D x0 C h alrededor de x0 D 0 Q est dada por 1 y2 .0 C h/ D y.0/ C y 0 .0/h C y 00 .0/h2 D 1 C h C h2 : Q 2

M at

em

Como en general y.x0 C h/ y2 .x0 C h/, obtendremos errores al usar esta aproximacin los cuales pueden Q catalogarse como errores absoluto, relativo y porcentual; las frmulas para calcularlos son exactamente las mismas que se denieron en la seccin previa.

at

ic

a1

y.x0 C h/

1 y2 .x0 C h/ D y.x0 / C y 0 .x0 /h C y 00 .x0 /h2 : Q 2

.c om

Para determinar un valor aproximado de la solucin en x D 1 usamos h D 1: y2 .1/ D 1 C 1 C .1/2 D 3: Q Por otra parte la solucin analtica del PVI es y.x/ D x 2 C x C 1 por lo cual, al evaluar en x D 1, resulta yexacto D 3. En este caso y2 D yexacto . Q

El resultado del ejemplo anterior es de utilidad para describir una forma alternativa de interpretar la aproximacin cuadrtica. En efecto, la aproximacin lineal y1 .x/ est dada en general por Q y1 .x0 C h/ D y.x0 / C y 0 .x0 /h: Q As que en este caso particular: y1 .0 C h/ D y.0/ C y 0 .0/h D 1 C hI Q de donde, el cambio en la variable dependiente es y1 .0 C h/ Q Denotamos el cambio anterior por k1 , es decir, k1 D y1 .0 C h/ Q y.0/ D y 0 .0/h D f .0; 1/h D h: y.0/ D y 0 .0/h D h:

y1 .x0 C h/ D y.x0 / C y 0 .x0 /.x Q

Si ahora calculamos la aproximacin lineal de la solucin en este punto .h; 1 C h/ tenemos: y1 .h C h/ D y.h/ C y 0 .h/h D y.h/ C f .h; 1 C h/h D y.h/ C .2h C 1/h D y.h/ C 2h2 C h: Q

at

Denotando ahora el nuevo cambio en la variable dependiente como:

em
y;

y usando k1 D h, tenemos que

ww w.

k2 D y1 .h C h/ Q

at ic a1
k2 k1 C k2 : 2 1 .k1 C k2 / : 2

Nos preguntamos ahora si es posible evaluar la funcin en un segundo punto y reproducir la aproximacin cuadrtica. La respuesta es armativa; para mostrarlo observemos que el punto .x0 Ch; y0 Ck1 / D .h; 1Ch/ se encuentra sobre la recta denida por la aproximacin lineal, que est dada por x0/ D y0 C f .x0 ; y0 /h:

y.h/ D 2h2 C h; k1

: 2 En resumen, y.x0 / D y.0/ D 1, k1 D y 0 .0/h D h, y 00 .0/ D 2, por lo cual la aproximacin cuadrtica se puede escribir como: 1 k2 k1 k1 C k2 y2 .h/ D y.x0 / C y 0 .x0 /h C y 00 .x0 /h2 D 1 C h C h2 D 1 C k1 C Q D 1C : 2 2 2 En los ejemplos siguientes mostraremos que siempre ser posible determinar k1 & k2 de forma que y2 .h/ D y0 C Q

k2 D 2h2 C k1 ) h2 D

Es decir, la aproximacin cuadrtica est dada, al menos en este ejemplo, por el valor de la funcin f .x; y/ k1 C k2 en .x0 ; y0/ ms un cambio lineal promedio . 2 Ejemplo 7.3.2 Encuentre una aproximacin cuadrtica de la solucin de la ecuacin diferencial: y0 D x con y.1/ D 2;

en el punto x D 1:2. Posteriormente demuestre que dicha aproximacin se puede reescribir como y2 D y0 C Q

.c

om

Derivando la ecuacin diferencial con respecto a x: y 00 D 1 y0 D 1 x C y:

Usando la condicin inicial, obtenemos: y.1/ D 2; y 0 .1/ D 1 y.1/ D 1 2D 1; y 00 .1/ D 1 1 C y.1/ D 2:

Para determinar un valor aproximado de la solucin en el punto x D 1:2, desarrollamos y.x/ en polinomio de Taylor de orden dos en x0 D 1: 1 y2 .1 C h/ D y.1/ C y 0 .1/h C y 00 .1/h2 D 2 Q 2 Con h D 0:2 tenemos que y2 .1:2/ D 2 Q 0:2 C .0:2/2 D 1:84 : h C h2 :

Por otra parte, como la solucin analtica es y D x 1 C 2e 1 x (ver ejemplo ??), al evaluar en x D 1:2 se tiene que yexacto D 1:2 1 C 2e 1 1:2 D 0:2 C 2e 0:2 1:8375 : La aproximacin cuadrtica proporciona una cifra decimal exacta con un error porcentual de: EP D 100

Notemos que, en este caso, f .x; y/ D x k1 D y 0 .1/h D f .1; 2/h D .1

y. Si denotamos k1 & k2 como h/ h D h C 2h2 D k1 C 2h2 ;

2/h D h;

k2 D f .1 C h; 2 C k1 /h D 1 C h entonces tendremos: Por lo que

at em at
h2 D k2 2 y0 D x y ; xCy con y2 D y0 C Q yD xC 0:9 C

.2 C k1 / h D 1 C h k1

ic a1
: k2 2 1 .k1 C k2 / : 2 2x 2 C 23; y/

Como se quera demostrar.

Ejemplo 7.3.3 Encuentre una aproximacin cuadrtica de la solucin en x D 0:9 de la ecuacin diferencial y.1/ D 4:

Muestre despus que

En este caso, la solucin analtica, ver ejemplo ??, est dada por

que evaluada en x D 0:9 produce el valor exacto: yexacto D 2.0:9/2 C 23 4:0619 :

Para determinar un valor aproximado de la solucin, derivamos la ecuacin diferencial; as obtenemos : y 00 D .1 y 0 /.x C y/ .1 C y 0 /.x .x C y/2 D 2y 2xy 0 : .x C y/2

ww

w.

y2 .1 C h/ D 2 Q

h C h2 D 2 C k1 C

k1

.c om
.2

yexacto y2 Q 1:8375 1:84 % D 100 % D 0:1361% : yexacto 1:8375

D 2C

k1 C k2 : 2

De la condicin inicial y.1/ D 4, tenemos: 1 4 3 D D 0:6; 1C4 5 8 C 1:2 9:2 y 00 .1/ D D D 0:368 : 2 .1 C 4/ 25 y 0 .1/ D Para determinar un valor aproximado de la solucin en el punto x D 0:9, consideremos la aproximacin cuadrtica y2 .x/ de la solucin alrededor de x D 1: Q 1 y2 .1 C h/ D y.1/ C y 0 .1/h C y 00 .1/h2 D 4 Q 2 En este caso tenemos h D 0:1, por lo cual y2 .0:9/ D 4 Q 0:6. 0:1/ C 0:184. 0:1/2 D 4:0618 : 0:6h C 0:184h2 : (7.1)

En este ejemplo, la aproximacin cuadrtica proporciona cuatro cifras decimales exactas de la solucin con un error porcentual de

em

denotamos, como en los ejemplos previos,

at

ste es un error muy pequeo que mejora el calculado por el mtodo de Euler, en el ejemplo ?? donde EP D 0:04554. Por otra parte, como x y f .x; y/ D ; xCy

ww

w.

y tambin

at

k1 D y 0 .1/h D f .1; 4/h D

ic

a1

1 4 h D 0:6h 1C4

k2 D f .1 C h; 4 C k1 /h D f .1 C h; 4

0:6h/h D

3 C 1:6h 5 C 0:4h

3 C 1:6h 5 hD h D 5 C 0:4h 5

.c om

EP D 100

yexacto y2 Q 4:0619 4:0618 % D 100 % D 0:00025% : yexacto 4:0618

0:6 C 0:32h 1 C 0:08h

h:

Si usamos ahora la frmula siguiente, que representa el desarrollo de una serie geomtrica, 1 D1 1Cx obtenemos: k2 D . 0:6 C 0:32h/.1 de donde: 0:184h2 Por lo que, usando (7.1): y2 .h/ D 4 Q 0:6h C 0:184h2 4 C k1 C k2 2 k1 D 4C k1 C k2 : 2 0:08h C 0:0064h2 C /h D 0:6h C 0:368h2 C k2 2 k1 k1 C 0:368h2; x C x2 ; cuando j x j < 1;

Ejemplo 7.3.4 Considere la ecuacin diferencial y0 D x y; con y.1/ D 2:

Encuentre una aproximacin cuadrtica de la solucin en x D 1:1; posteriormente use esta aproximacin para determinar la solucin en x D 1:2. H Sabemos que y 00 D 1 y0 D 1 x C y:

La aproximacin cuadrtica de la solucin alrededor de x D 1; con y.1/ D 2; y 0 .1/ D 1 es Considerando h D 0:1: y2 .1:1/ D 2 Q 2 D 1 & y 00 .1/ D 1 . 1/ D 2;

1 y2 .1 C h/ D y.1/ C y 0 .1/h C y 00 .1/h2 D 2 Q 2 0:1 C .0:1/2 D 1:91 :

h C h2 :

Si ahora calculamos la funcin y.x/ y sus dos primeras derivadas en x D 1:1: y.1:1/ D 1:91; y 0 .1:1/ D 1:1 y.1:1/ D 1:1 1:91 D 0:81; y 00 .1:1/ D 1 1:1 C y.1:1/ D 1:81 :

y2 .1:2/ D 1:91 Q

Comparando con la solucin exacta obtenida al evaluar la solucin y D x

at e

0:81.0:1/ C 0:905.0:1/2 D 1:838 : 1 C 2e 1

at

Evaluando en h D 0:1, obtenemos la aproximacin pedida:

ic

1 y2 .1:1 C h/ D y.1:1/ C y 0 .1:1/h C y 00 .1:1/h2 D 1:91 Q 2

a1

.c
0:81h C 0:905h2 :

Por lo cual la aproximacin cuadrtica en x D 1:1 est dada por

om

en x D 1:2, tenemos:

EP D 100

que, en efecto, es menor que el cometido al calcular y2 .1:2/ con una sola aproximacin cuadrtica, como Q puede vericarse en el ejemplo 7.3.2, donde EP D 0:1382%. Ejemplo 7.3.5 Encuentre una aproximacin cuadrtica de la solucin en x D 0:9 del PVI y0 D x y ; xCy con y.1/ D 4:

Primero determine una aproximacin de la solucin en x D 0:95 y, posteriormente, use su resultado para obtener la aproximacin pedida. H Procedamos como en los tres ejemplos anteriores, primero observemos que la solucin analtica est dada por y D x C 2x 2 C 23; y evaluada en x D 0:9 produce el valor exacto yexacto D 0:9 C 2.0:9/2 C 23 4:0619 :

ww w.

En este caso el error porcentual cometido en la aproximacin est dado por yexacto y2 Q 1:8375 1:838 % D 100 % D 0:0272%: yexacto 1:8375

yexacto D 1:2

1 C 2e 1

1:2

D 0:2 C 2e

0:2

1:8375 :

Para determinar la solucin aproximada, necesitamos la derivada de la ecuacin diferencial, esto es, y 00 D De la condicin inicial, tenemos que y.1/ D 4; y 0 .1/ D 1 4 3 D D 0:6 1C4 5 & y 00 .1/ D 8 C 1:2 9:2 D D 0:368 : .1 C 4/2 25 2y 2xy 0 : .x C y/2

Por lo cual, la aproximacin cuadrtica de la solucin alrededor de x D 1 es 1 y2 .1 C h/ D y.1/ C y 0 .1/h C y 00 .1/h2 D 4 Q 2 Si usamos h D 0:05, tenemos: y2 .0:95/ D 4 Q 0:6. 0:05/ C 0:184. 0:05/2 D 4:0305 : 0:6h C 0:184h2 :

Calculando y.x/ y sus dos primeras derivadas en x D 0:95, tenemos:

Considerando ahora la aproximacin cuadrtica de la solucin alrededor de x D 0:95: 1 y2 .0:95 C h/ D y.0:95/ C y 0 .0:95/h C y 00 .0:95/h2 D 4:0305 Q 2 Evaluando esta expresin cuando h D 0:05, obtenemos nalmente: y2 .0:9/ D 4:0305 Q 0:6185. 0:05/ C 0:1862. 0:05/2 D 4:0619 :

at

em

y.0:95/ D 4:0305; x y 0:95 4:0305 y 0 .x/ D ) y 0 .0:95/ D D 0:6185; xCy 0:95 C 4:0305 2y 2xy 0 2.4:0305/ 2.0:95/. 0:6185/ y 00 .x/ D ) y 00 .0:95/ D D 0:3723 : .x C y/2 .0:95 C 4:0305/2

at

ic a1

.c om
0:6185h C 0:1862h2:

Resultado que es exactamente igual en cuatro cifras decimales al resultado exacto, claramente el error porcentual es cero por ciento.

En los primeros tres ejemplos de esta seccin mostramos que la aproximacin cuadrtica se puede reescribir usando dos trminos k1 y k2 que se obtienen evaluando la funcin f .x; y/ en puntos adecuados. En los ltimos dos ejemplos, mostramos cmo utilizar una doble aproximacin cuadrtica para estimar la solucin de la ecuacin diferencial en un punto dado. Desde luego que podemos estimar la solucin con mayor nmero de aproximaciones cuadrticas, mejorando sustancialmente los resultados que arroja el mtodo de Euler. Para ilustrar este nuevo mtodo, consideremos los cambios en la variable dependiente y, que se obtienen al calcular con el mtodo de Euler las dos primeras aproximaciones; estos cambios son: k1 Dy1 k2 Dy2 y0 D h f .x0 ; y0/ ; y1 D h f .x1 ; y1/ D h f .x0 C h; y0 C k1 / :

ww

w.

y Error absoluto

y1 y.x1 / k1 y0

k2

x x0 h x1 h x2

Observe que k1 & k2 han conservado sus ndices. En general, para obtener una aproximacin yi de y.xi / D y.x0 C ih/, seguimos el llamado mtodo de Euler mejorado que establecemos a continuacin: Mtodo de Euler mejorado La solucin numrica de la ecuacin diferencial y 0 D f .x; y/, con la condicin inicial y.x0 / D y0 y con tamao de paso h est formada por los puntos.xi C1; yi C1 / que se obtienen mediante las frmulas de recurrencia xi C1 D xi C hI k1 D h f .xi ; yi /I k2 D h f .xi C h; yi C k1 /I 1 yi C1 D yi C .k1 C k2 / ; 2 con i D 0; 1; 2; 3 : : : y los segmentos rectilneos entre cada par de puntos consecutivos. A continuacin, mostraremos que este mtodo produce, en general, mejores resultados que el mtodo de Euler. En efecto, sustituyendo k1 y k2 en la yi C1 : yi C1 D yi C h f .xi ; yi / C f .xi C h; yi C k1 /: 2 (7.3)

ww

w.

y1 D h f .x1 ; y1 / I y2 D h f .x2 ; y2 / D h f .x1 C h; y1 C k1 / I 1 y2 Dy1 C .k1 C k2 / : 2

at e

k1 Dy2 k2 Dy3

at

Siguiendo este proceso, podemos calcular ahora una aproximacin y2 de y.x2 / D y.x0 C 2h/ utilizando:

ic a

y1 D y0 C

1 .k1 C k2 / : 2

1.c

om

El signicado geomtrico de ambas cantidades se ilustra en la gura anterior. Observemos que k1 es el cambio en la altura cuando x se incrementa una cantidad h partiendo de x0 , mientras que k2 es el cambio en la variable y cuando x pasa de x0 C h a x0 C 2h. Observe ahora que el promedio de ambas cantidades produce un cambio en la altura que aproxima mejor el cambio que sufre la solucin y.x/. Es decir, una mejor aproximacin a y.x1 / que la obtenida por el mtodo de Euler est dada por (7.2)

. /

Desarrollando (*) en serie de Taylor, la funcin f .x; y/ alrededor del punto .xi ; yi /: f .xi C h; yi C k1 / D f .xi ; yi / C hfx .xi ; yi / C k1 fy .xi ; yi /; luego: yi C1 D yi C h f .xi ; yi / C f .xi ; yi / C h fx .xi ; yi / C k1 fy .xi ; yi / C 2 h2 k1 fx .xi ; yi / C fy .xi ; yi / C D yi C h f .xi ; yi / C 2 h D

Usando k1 D h f .xi ; yi / y simplicando obtenemos: yi C1 D yi C h f .xi ; yi / C O bien yi C1 D yi C h y 0 .xi ; yi / C h2 00 y .xi ; yi / C 2 (7.4) h2 fx .xi ; yi / C f .xi ; yi / fy .xi ; yi / C 2

Para obtener la ltima igualdad utilizamos y 0 D f .x; y/. En consecuencia: y 00 D

Ejemplo 7.3.6 Considere la ecuacin diferencial y 0 D x 2 xy 2 , con la condicin inicial y.1/ D 2. Utilice el mtodo de Euler mejorado para estimar y.1:2/; considere h D 0:1. H Como el tamao de paso es h D 0:1, debemos repetir el proceso de Euler mejorado 2 veces, con i D 1, para calcular k1 , k2 & y1 . Si consideramos que .x0 ; y0 / D .1; 2/ y que f .x; y/ D x 2 xy 2 , obtenemos: k1 D h f .x0 ; y0 / D 0:1f .1; 2/ D 0:1 . 3/ D 0:3I k2 D h f .x0 C h; y0 C k1 / D 0:1f .1:1; 1:7/ D 0:1 .1:1/2 1:1.1:7/2 D 0:1969I 1 y1 D y0 C .k1 C k2 / D 2 C 0:5 . 0:3 0:1969/ D 1:7516 : 2 Repetimos el proceso de (7.3), con i D 2 y obtenemos en este caso: k1 D h f .x1; y1 / D 0:1f .1:1; 1:7516/ D 0:2165I k2 D h f .x1 C h; y1 C k1 / D 0:1f .1:2; 1:5351/ D 0:1388I 1 y2 D y1 C .k1 C k2 / D 1:7516 C 0:5 . 0:2165 0:1388/ D 1:5740 : 2 En conclusin, y2 D 1:5740 es una aproximacin de y.1:2/.

ww w.

Claramente la relacin (7.4) muestra que la aproximacin del mtodo de Euler mejorado es equivalente a un desarrollo de Taylor de orden dos; en consecuencia el error cometido en cada paso es proporcional a h3 , lo cual es una mejora sustantiva con respecto al mtodo de Euler. Sin embargo, tal mejora se obtiene mediante el aumento en el nmero de evaluaciones de la funcin f .x; y/ que pasa de uno a dos en cada aproximacin.

at e

at ic

df @f @f dy @f @f D C D Cf : dx @x @y dx @x @y

a1 .c

om

Ejemplo 7.3.7 Considere el PVI y 0 D x y, con y.0/ D 1. Determine una solucin numrica en el intervalo 0; 1 utilizando primero h D 0:1 y despus h D 0:05. Compare los resultados de los dos casos anteriores con la solucin numrica obtenida mediante el mtodo de Euler y la solucin analtica de la ecuacin diferencial. H Para h D 0:1, necesitaremos repetir diez veces el proceso de Euler mejorado. Basta hacerlo un par de veces para ilustrar el procedimiento. Consideremos el caso i D 1 del proceso (7.3) para calcular k1 , k2 & y1 ; necesitamos adems .x0 ; y0 / D .0; 1/ y f .x; y/ D x y. Obtenemos en este caso: k1 D h f .x0 ; y0 / D 0:1f .0; 1/ D 0:1 . 1/ D 0:1I k2 D h f .x0 C h; y0 C k1 / D 0:1f .0:1; 0:9/ D 0:1 .0:1 0:9/ D 0:08I 1 y1 D y0 C .k1 C k2 / D 1 C 0:5 . 0:1 0:08/ D 0:91 : 2 Nuevamente repetimos el proceso con i D 2; obtenemos ahora: k1 D h f .x1 ; y1/ D 0:1f .0:1; 0:91/ D 0:081I k2 D h f .x1 C h; y1 C k1 / D 0:1f .0:2; 0:829/ D 0:0629I 1 y2 D y1 C .k1 C k2 / D 0:91 C 0:5 . 0:081 0:0629/ D 0:8381 : 2

at ic

Repitiendo el proceso otras 8 veces obtenemos los resultados que se muestran en la tabla siguiente. [Hemos incluido adems los resultados correspondientes a la solucin analtica y.x/ D x 1 C 2e x y los errores absoluto y porcentual.]

a1 .c
y.exacta/ 1 0.9097 0.8375 0.7816 0.7406 0.7131 0.6976 0.6932 0.6987 0.7131 0.7358

om

0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000

0.9100 0.8381 0.7825 0.7417 0.7142 0.6989 0.6945 0.7000 0.7145 0.7372

ww w.

0.0000

1.0000

at e

xi

yi

k1

k2

EA 0 0.0003 0.0006 0.0009 0.0011 0.0011 0.0013 0.0013 0.0013 0.0014 0.0014

EP (%) 0 0.03 0.07 0.12 0.15 0.15 0.19 0.19 0.19 0.2 0.19

-0.1000 -0.0810 -0.0638 -0.0483 -0.0342 -0.0214 -0.0099 0.0005 0.0100 0.0186 0.0263

-0.0800 -0.0629 -0.0474 -0.0334 -0.0208 -0.0093 0.0011 0.0105 0.0190 0.0267 0.0337

Consideremos ahora h D 0:05; debemos repetir el proceso 20 veces para obtener la solucin numrica. Aplicando las ecuaciones (7.3), con i D 1: k1 D h f .x0 ; y0 / D 0:05f .0; 1/ D 0:05 . 1/ D 0:05I k2 D h f .x0 C h; y0 C k1 / D 0:05f .0:05; 0:95/ D 0:05 .0:05 1 y1 D y0 C .k1 C k2 / D 1 C 0:5 . 0:05 0:045/ D 0:9525 : 2 Repitiendo otra vez el proceso con i D 2 obtenemos: k1 D h f .x1 ; y1 / D 0:05f .0:05; 0:9525/ D 0:0451I k2 D h f .x1 C h; y1 C k1 / D 0:05f .0:1; 0:9074/ D 0:0404I 1 y2 D y1 C .k1 C k2 / D 0:91 C 0:5 . 0:0451 0:0404/ D 0:9098 : 2 Los resultados para el intervalo 0; 1 se muestran en la tabla siguiente, donde se han incluido los clculos de Euler (para h D 0:05/ y exactos. Note que la solucin con el mtodo de Euler mejorado con h D 0:1 reproduce mejor los resultados exactos que la solucin numrica obtenida con el mtodo de Euler con h D 0:05. Observe adems que las soluciones numricas calculadas con el mtodo de Euler mejorado estn ligeramente arriba de la solucin analtica. 0:95/ D 0:045I

Euler mejorado xi 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 yi .h D 0:1/ 1 0.91

a1

Euler

.c
Exacto yexacto 1 0.9097 0.8375 0.7816 0.7406 0.7131 0.6976 0.6932 0.6987 0.7131 0.7358 yi .h D 0:05/ 1

em

at ic
0.782

yi .h D 0:05/ 1.0000 0.9050 0.8290 0.7702 0.7268 0.6975 0.6807 0.6754 0.6803 0.6944 0.7170

0.8381

at

0.9098 0.8377

ww w.

0.7825 0.7417 0.7142 0.6989 0.6945 0.7 0.7145 0.7372

0.7411 0.7135 0.6982 0.6938 0.6993 0.7138 0.7364

Las grcas de las curvas obtenidas se muestran en la siguiente gura. El origen de coordenadas en la grca se ha colocado en .0; 0:6/.

om

-0.1 0.5

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

En conclusin, al igual que en el mtodo de Euler podemos reducir los errores de aproximacin y de propagacin haciendo ms pequeo el tamao de paso h, lo que implica un mayor esfuerzo de clculo y, en consecuencia, para determinar rpidamente la solucin de la ecuacin considerada es necesario utilizar alguna herramienta computacional. En los siguientes dos ejemplos ilustramos el clculo de la solucin numrica utilizando Excel y en Mathematica hemos implementado el pseudocdigo asociado a este mtodo de Euler mejorado.

1. Proporcionar f; x0; y0; h; n. 2. Imprimir x0; y0 . 3. Desde i D 1 hasta i D n. a. Calcular: k1 D h k2 D h f .x0 ; y0 /I f .x0 C h; y0 C k1 / I k1 C k2 y1 D y0 C : 2 b. Hacer y0 D y1 ; x0 D x0 C h; c. Imprimir x0 ; y0 . 4. Terminar.

Ejemplo 7.3.8 Use el mtodo de Euler mejorado en una hoja de clculo de Excel para determinar un valor aproximado de y.1/, si y.x/ es la solucin de la ecuacin diferencial y 0 D x 2 y C y; con la condicin inicial y.0/ D 1. Considere que el tamao de paso es h D 0:1. H Utilizamos las siguientes instrucciones en una hoja de clculo de Excel para resolver el ejemplo: El mtodo de Euler mejorado en Excel. 1. En las celdas A1, A2, A3 se escriben las etiquetas: "x0=, y0=, h=". 2. En las celdas B1, B2, B3 se escriben "=0, =1, =0.1", respectivamente. 3. En las celdas A5, B5, C5, D5, E5 se escriben las etiquetas: "i; xi ; yi ; k1 ; k2". 4. Se escriben en las celdas A6-A16 los nmeros "0; 1; 2; : : : ; 10". 5. En las celdas B6 y C6 se escriben, respectivamente: "=B1, =B2".

ww w.

at

em

at ic a1

Pseudocdigo del mtodo de Euler mejorado

.c

om

0.7

1.0

0.8

0.9

1.0

h D 0:1 Euler mejorado h D 0:05 Euler mejorado h D 0:05 Euler

6. En la celda D6 se escribe "=$B$3*(B62*C6+C6)". Observe que, en este paso, se evala la funcin f .x; y/ D x 2 y C y en el punto .x0 ; y0 / y se multiplica por h; a esta expresin la llamamos k1 . 7. En la celda E6 se escribe "=$B$3*((B6+$B$3)2*(C6+D6)+(C6+D6))". Observe que en este paso, se evala la funcin f .x; y/ D x 2y C y en el punto .x0 C h; y0 C k1 /y se multiplica por h; a esta expresin la llamamos k2 . 8. En la celda B7 se escribe "D B6 C $B$3" . 9. En la celda C7 escribimos ahora "=C6+(D6+E6)/2 ". Observe que, en este paso, estamos usando el mtodo de Euler mejorado. 10. Se seleccionan las celdas D6-E6 y se copian en D7-E7. 11. Se seleccionan las celdas B7-E7 y se arrastran hasta llegar a las celdas B16-E16. 12. Se graca la solucin utilizando el asistente de grcos con la opcin de XY-Dispersin. En la tabla siguiente se muestran los resultados numricos obtenidos.

i 0 1 2 3 4 5 6 7 8 9 10

xi 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

yi 1 1.10555 1.224675224 1.362044594 1.523885457

k1

k2

0.127366223 0.148462861 0.176770713 0.214820228 0.266212101 0.336175319 0.432447248 0.566687581 0.756794044

ic a1
y.0/ D 1;

0.11166055

at
3y;

em

at

1.718561824 1.957441917 2.256210192 2.636873463 3.130870612 3.783970221

ww w.

Ejemplo 7.3.9 Resolver el PVI y0 D x utilizando el mtodo de Euler mejorado repitiendo el proceso n D 10 veces con h D 0:3 implementando el mtodo en Mathematica. H Aplicamos el pseudocdigo del mtodo para elaborar el siguiente programa en el paquete Mathematica, hemos incluido comentarios para que dicho programa quede ms claro.

.c o

0.126589897 0.147372518 0.175218865 0.212582021 0.262939959 0.331324449 0.425151224 0.555547049 0.739511639 1.003508903

0.1

0.1111

El mtodo de Euler mejorado en Mathematica f[x_,y_]:=x-3y; x0=0; y0=1; h=0.3; n=10; lista={{x0,y0}}; Do[ k1=h*f[x0,y0]; k2=h*f[x0+h,y0+k1]; y1=y0+(k1+k2)/2; y0=y1; x0=x0+h; AppendTo[lista,{x0,y0}], {i,1,n}]; ListPlot[lista] (* (* (* (* (* (* (* (* (* (* (* (* (* (* Definir f *) Abscisa del punto inicial *) Ordenada del punto inicial *) Incremento en el paso *) Total de pasos a realizar *) Definir lista con punto inicial *) Calcular k1 *) Calcular k2 *) Determinar y1 *) Intercambiar y0 con y1 *) Incrementar x0 *) Incluir punto en la lista *) Terminar el proceso *) Graficar los puntos obtenidos *)

Despus de correr este programa, obtenemos los resultados que se muestran en la tabla siguiente; hemos incluido los valores exactos, valores que se obtiene evaluando la funcin

x 0

at
y2 Q 1 0.55 0.37225 0.331986 0.361153 0.425382 0.507318 0.598196 0.693589 0.791262 0.890087

em

M at

0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.

ww

w.

ic a
yexacto 1. 0.440633 0.272554 0.263562 0.319249 0.401232 0.493907 0.590929 0.689718 0.789226 0.889026

en los puntos x D 0; 0:3; 0:6; : : : ; 3.

1.c

y.x/ D

1 3x 9

1 C 10e

3x

om
;

Ejercicios 7.3.1 Mtodo de Euler mejorado. Soluciones en la pgina 15 Determine una aproximacin cuadrtica de la solucin y.x/ de cada uno de los siguientes PVI utilizando el h proporcionado y calcule el error porcentual. En los casos que se requiera aplique dos veces el proceso de aproximacin cuadrtica para obtener una estimacin de la solucin. 1. y 0 D 2x 2. y 0 D x 3. y 0 D y 4. y 0 D x 2 y, xy, x C 5, y, con con con con y.0/ D 3 en x D 0:2, para h D 0:2.

Use el mtodo de Euler mejorado para determinar una aproximacin numrica de la solucin en el punto indicado de cada una de los siguientes PVI; utilice el tamao de paso proporcionado y utilice tambin tambin redondeo a cuatro cifras decimales en todos sus clculos. 5. y 0 D x 2 y, con 1, con y.1/ D 2; calcule y.1:5/, para h D 0:1. con y.2/ D 1; calcule y.2:5/, para h D 0:1.
ww w.

6. y 0 D x C 2y 7. y 0 D x2 C 4 , y

y.0/ D 1; calcule y.0:25/, para h D 0:05.

8. y 0 D xy,

con

y.2/ D 1; calcule y.3/, para h D 0:2.

at

em

y.0/ D 3 en x D 0:4, para h D 0:2.

at

ic

y.1/ D 1 en x D 1:2, para h D 0:1.

a1

.c om

y.1/ D 2 en x D 1:1, para h D 0:1.

Ejercicios 7.3.1 Mtodo de Euler mejorado. Pgina 14

1. 2:5. 2. 1:9. 3. 2:0841. 4. 2:0324.

a1 at at ic M

.c

om

5. 1:8592. 6. 3:7297. 7. 1:7352. 8. 11:3026.