04.05.2013 Views

Metode Numerik 2 - Universitas Indonesia

Metode Numerik 2 - Universitas Indonesia

Metode Numerik 2 - Universitas Indonesia

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.

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)

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

Saved successfully!

Ooh no, something went wrong!