fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
44<br />
2.6.2 Integração pela Regra do Trapézio<br />
A fórmula <strong>de</strong> integração pela regra do trapézio geralmente é mais acurada que a regra<br />
do retângulo (KREYSZIG, 1999). A Fig. 2.8 mostra que a regra do trapézio é obtida tomando<br />
a mesma subdivisão da regra do retângulo e aproximando φ por segmentos <strong>de</strong> linha com<br />
pontos finais nos seus extremos.<br />
y = φ<br />
( x)<br />
K<br />
0<br />
K<br />
x1<br />
x2<br />
x N − 1<br />
L<br />
x<br />
Figura 2.8: Regra do trapézio<br />
A área sob a curva <strong>de</strong> ( φ ) entre [ 0 , L]<br />
é aproximada por ( )<br />
N trapezói<strong>de</strong>s <strong>de</strong> área:<br />
1<br />
2<br />
[ φ ( 0) φ( x )] ∆x<br />
, [ φ ( x ) φ( x )] ∆x<br />
, K, [ ( ) + φ( L)<br />
] ∆x<br />
+<br />
1<br />
1<br />
2<br />
1<br />
+<br />
2<br />
1<br />
2<br />
φ x N −1<br />
(2.23)<br />
Tomando sua soma, obtém-se a regra do trapézio:<br />
on<strong>de</strong><br />
L<br />
⎡1<br />
1 ⎤<br />
φ = ∫φ( x) dx ≈ ∆x⎢<br />
φ( 0) + φ( x1<br />
) + φ( x2<br />
) + K + φ( x N −1<br />
) + φ( L)<br />
⎣2<br />
2 ⎥<br />
(2.24)<br />
⎦<br />
0<br />
∆ x assume o mesmo valor empregado na regra do retângulo.<br />
2.7 MÉTODO DIRETO DE SOLUÇÃO<br />
Para resolver o sistema <strong>de</strong> equações algébricas representado pela Eq. (2.4), ou seja,<br />
para obter os valores das soluções numéricas para a variável <strong>de</strong> interesse ( φ ), utilizou-se o