15.05.2013 Views

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

Cálculo Numérico - Engenharia Civil UEM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CAPÍTULO 12. SOLUÇÃO NUMÉRICA DE EQUAÇÕES DIFERENCIAIS ORDINÁRIAS 391<br />

O método do trapézio é dado por (12.12). Assim:<br />

Portanto:<br />

α0 = −1 , β0 = 1 2 ,<br />

α1 = 1 , β1 = 1 2 .<br />

C0 = α0 + α1 ⇒ C0 = −1 + 1 = 0 ,<br />

C1 = α1 − (β0 + β1) ⇒ C1 = 1 − ( 1<br />

2<br />

C2 = α1<br />

2! − β1 ⇒ C2 = 1 1<br />

− = 0 ,<br />

2 2<br />

C3 = α1<br />

3!<br />

β1<br />

−<br />

(2)! ⇒ C3 = 1 1 1<br />

− = −<br />

6 4 12 .<br />

1<br />

+ ) = 0 ,<br />

2<br />

Logo, C0 = C1 = C2 = 0 e C3 = 0. Portanto a ordem do método do trapézio é q = 2 e a constante do<br />

erro é C3 = − 1 12 .<br />

Exercício<br />

12.3 - Determinar a ordem e a constante do erro para:<br />

a) regra do ponto médio,<br />

b) método de Simpson,<br />

c) método de Adams-Moulton,<br />

d) método de Adams-Bashforth,<br />

e) método 3 8<br />

de Simpson.<br />

12.3.4 Erro de Truncamento Local<br />

Agora, podemos definir formalmente o erro de truncamento local de um método linear de passo<br />

múltiplo.<br />

Definição 12.4 - Definimos Erro de Truncamento Local, em xn+k do método linear de passo<br />

múltiplo, definido por (12.7), por:<br />

Tn+k = L [y(xn); h] =<br />

onde y(x) é a solução exata do (p.v.i) (12.2).<br />

k<br />

j=0<br />

[αjy(xn+j) − hβjy ′ (xn+j)] ,

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

Saved successfully!

Ooh no, something went wrong!