Metode Numerik 2 - Universitas Indonesia
Metode Numerik 2 - Universitas Indonesia
Metode Numerik 2 - Universitas Indonesia
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
94<br />
Mencari wi<br />
:<br />
Untuk integrand f(x) dan s(x), yang merupakan polinomial orde 2N-1 berlaku<br />
integrasi numerik:<br />
b<br />
N<br />
∫ v(x)f(x)dx = ∑wif(x<br />
i)<br />
∫ v(x)s(x)dx = ∑<br />
i=<br />
1<br />
i=<br />
1<br />
a<br />
Integrasi numerik yang sama tentu berlaku juga untuk integrand polinomial orde<br />
lebih rendah, contohnya r(x), yang berorde N-1:<br />
b<br />
N<br />
∫ v(x)r(x)dx = ∑wir(x<br />
i)<br />
,<br />
N-1<br />
r(x) = ∑<br />
a<br />
i=<br />
1<br />
i=<br />
0<br />
b<br />
a<br />
N<br />
a x<br />
i<br />
w s(x )<br />
Dari penurunan rumus quadrature trapezoid, Simpson dll sebelum ini diketahui<br />
bahwa untuk mencari wi<br />
bisa digunakan r(x) sembarang polinomial orde N-1<br />
(koefisien ai<br />
tidak diperlukan). Karena itu, dipilih r(x) yang memudahkan:<br />
r(x) = l(x, xi<br />
) = ∏<br />
≠<br />
⎛ x − x ⎞ j ⎜ ⎟,<br />
⎜ xi<br />
x ⎟<br />
⎝<br />
−<br />
⎠<br />
j i<br />
j<br />
(i, j =<br />
b<br />
N<br />
Diperoleh: ∫ v(x)l(x, xj)dx<br />
= ∑wil(x<br />
i,<br />
xj)<br />
= wj<br />
a<br />
i=<br />
1<br />
1, ..., N)<br />
w<br />
j<br />
b<br />
= ∫ v(x)l(x, x<br />
a<br />
i<br />
i<br />
l(x , x ) = δ<br />
k<br />
j<br />
i<br />
i<br />
)dx<br />
ik<br />
(j =<br />
1, ..., N)