11.08.2013 Views

Apunte

Apunte

Apunte

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.

180 8. RESOLUCIÓN DE ECUACIONES DIFERENCIALES ORDINARIAS<br />

para todo x, y ∈ IR y t ∈ [t0, T ]. Entonces, se tiene<br />

|x(T ) − xn| ≤ τmáx<br />

K (eK(T −t0) − 1), (8.20)<br />

donde τmáx = máx<br />

1≤i≤n {|τi|} con τi el error de truncamiento local del método en el paso i.<br />

Demostración. Sea ei = x(ti) − xi, el error global que se comete hasta el paso i.<br />

Consideramos la ecuación del método y la definición de τi (8.16),<br />

luego,<br />

xi+1 = xi + hΦ(ti, xi, h)<br />

x(ti+1) = x(ti) + hΦ(ti, x(ti), h) + hτi<br />

ei+1 = x(ti+1) − xi+1 = ei + h(Φ(ti, x(ti), h) − Φ(ti, xi, h)) + hτi.<br />

Usando que φ es Lipschitz se obtiene<br />

|ei+1| ≤ |ei| + h|Φ(ti, x(ti), h) − Φ(ti, xi, h)| + h|τi|<br />

≤ |ei| + hK|x(ti) − xi| + h|τi|<br />

≤ (1 + Kh)|ei| + h|τi|.<br />

Como |τi| ≤ τmáx, llamando A = 1 + Kh, tenemos<br />

|ei+1| ≤ A|ei| + hτmáx.<br />

Usando esta acotación para i = 0 y teniendo en cuenta que e0 = 0 se tiene<br />

|e1| ≤ hτmáx,<br />

|e2| ≤ A|e1| + hτmáx ≤ Ahτmáx + hτmáx = (1 + A)hτmáx,<br />

|e3| ≤ A|e2| + hτmáx ≤ A(1 + A)hτmáx + hτmáx<br />

= (A + A 2 )hτmáx + hτmáx<br />

= (1 + A + A 2 )hτmáx

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

Saved successfully!

Ooh no, something went wrong!