15.05.2013 Views

Salinas, P. (2010). - Repositorio Digital - Instituto Politécnico Nacional

Salinas, P. (2010). - Repositorio Digital - Instituto Politécnico Nacional

Salinas, P. (2010). - Repositorio Digital - Instituto Politécnico Nacional

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Esto suponiendo que pudimos encontrar esa constante M 1 que funcione como cota para<br />

la segunda derivada de la solución, como la encontramos en el primer paso. Siendo este el<br />

caso entonces el error total después del segundo paso está acotado por el error total<br />

después del primer paso, más el error por truncamiento asociado al segundo paso, y un<br />

tercer término, f t 1 , yt 1 <br />

f t 1 , y 1 <br />

t<br />

que estamos usando la pendiente 1, 1<br />

pendiente f t , yt <br />

en el punto , <br />

solución exacta.<br />

1 1<br />

que se toma en cuenta por el hecho de<br />

f t y en el punto t1, y 1en<br />

lugar de la<br />

t y t que es la que corresponde a la curva de la<br />

1 1<br />

En este momento, los autores, en la búsqueda de una cota relativamente conveniente<br />

f t , y t y f t , y considerando el valor<br />

para 2<br />

e , acotan la diferencia entre <br />

máximo de la derivada parcial<br />

, <br />

f<br />

y<br />

1 1<br />

1 1<br />

a lo largo del segmento de recta entre , <br />

t y t y<br />

1 1<br />

t1 y 1 ya que estos dos valores de f provienen de puntos a lo largo de la recta t t1<br />

Nuevamente es la continuidad de esa parcial la garantía de la existencia de la acotación<br />

Por lo tanto, se obtiene la estimación<br />

<br />

<br />

f t , y t f t , y M e<br />

1 1 1 1 2 1<br />

<br />

2 2<br />

<br />

<br />

t t<br />

e 2 e 1 M 2 e 1 t M 1 1 M 2 t e 1 M 1<br />

2 2<br />

Los autores afirman que el análisis del error total en el paso k-ésimo es esencialmente el<br />

mismo que el recién hecho para el segundo paso, y la acotación queda expresada de la<br />

siguiente manera.<br />

.<br />

Si 1 M es una cota para la segunda derivada a lo largo de la gráfica de la solución, y si 2 M<br />

f<br />

es una cota para la derivada parcial<br />

y<br />

, entonces el error e k está acotado por<br />

<br />

e 1 M t e M<br />

k 2 k11<br />

t<br />

t <br />

donde esta fórmula recursiva incluye el término M 1<br />

2<br />

<br />

que toma en cuenta el error por<br />

truncamiento en el k-ésimo paso, y el término 1M2 t e<br />

k1<br />

representa la magnificación<br />

del error en el paso k-1 debido al hecho de que estamos calculando el lado derecho de la<br />

t , y en vez de en el valor exacto del<br />

ecuación diferencial en la aproximación de Euler k k<br />

punto t , yt <br />

de la solución.<br />

k k<br />

2<br />

2<br />

2<br />

78

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!