Integral Numerik Sunkar E. Gautama, 2013 http://paradoks77.blogspot.

com

Integral numerik ialah metode untuk menghitung nilai integrasi suatu fungsi dalam suatu selang tanpa mempedulikan fungsi hasil integralnya dengan menggunakan metode numerik. Jadi, untuk menghitung nilai integral tentu = 2 2 + 5 3 terhadap secara numerik, tidak perlu mengetahui fungsi ( ). Dalam integral numerik, digunakan teorema integral tentu dari suatu fungsi terhadap dalam selang [, ] sama dengan luas daerah di bawah kurva yang dibatasi oleh sumbu X, garis = , dan = . Dengan demikian, digunakan suatu metode untuk menghitung luasan di bawah kurva tadi, dengan pengandaian luas daerah tadi setara dengan suatu segiempat dengan panjang = dan tinggi rata-rata ketinggian titik pada kurva, . Beberapa metode yang akan dibahas di bawah ini ialah metode trapesium, metode Simpson, dan metode Simpson 3/8, dengan perluasan integral komposit untuk kurva dengan loop yang tidak sederhana. 1. Metode Trapesium (trapezoidal rule) Metode trapesium ialah metode yang paling sederhana diantara tiga metode yang dibahas di sini. Metode trapesium mengandaikan daerah luasan yang ditinjau berbentuk trapesium dengan panjang = dan dipilih dua titik ketinggian yakni tinggi rusuk sejajar masing-masing 0 = () dan 1 = () yang berjarak = .

Pada trapesium, = panjang tinggi rata rata = Dengan =

0 dan 1 masing-masing ialah pembobotan untuk tiap ketinggian titik

yang diambil. Mengingat pada trapesium (dua sudut siku-siku) rusuk miringnya linier, maka pembobotannya seragam, yakni 0 = 1 = 1. Dengan demikian, luas daerah di bawah kurva ialah: 0 + 1 Luas = + (galat) 2 1

 + + (3 ) 2 1 = () = 0 + 1 + (3 ) 2 0 = ( ) =

2. Metode Simpson Jika pada metode trapesium hanya diambil dua titik ketinggian, maka pada metode simpson diambil tiga titik ketinggian yang berjarak sama, =
2

. 2

Tinggi rata-rata, =

0 0 + 1 1 + 2 2 . 0 + 1 + 2

Titik tinggi yang berada di tengah mendapatkan bobot

yang lebih tinggi, sebagaimana dapat dilihat bahwa ketinggian rata-rata cenderung mendekati ketinggian pada titik tengahnya. Menggunakan fungsi Lagrange orde-2: 1 2 0 2 0 1 2 = 0 + 1 + 2 0 1 0 2 1 0 1 2 2 0 2 1 Di mana =
, 0 2

= , 1 = + , 2 = = + 2. Dengan demikian
2 2

=
0

( ) =
0

2 () +

Di mana Rs ialah suku yang mengandung galat komputasi, O(h5). Dari kedua persamaan di atas, diperoleh rumus integral Simpson 2 = ( ) = 0 + 4 1 + 2 + (5 ) 3 0 Dengan demikian, diperoleh pembobotan 0 = 1, 1 = 4, 2 = 1.

3. Metode Simpson 3/8 Pada metode Simpson 3/8 dipilih empat titik ketinggian pada selang yang berjarak sama, =
3

, 3

sehingga 0 = , 1 = + , 2 = + 2, 3 = + 3 = .
0 0 + 1 1 + 2 2 + 3 3 0 + 1 + 2 + 3

Ketinggian rata-rata ialah = 3, 3 = 1.

, dengan 0 = 1, 1 = 3, 2 =

Luas daerah di bawah kurva:

=
3

( ) = + (galat) + (5 )

=
0

() =

3 0 + 3 1 + 3 2 + 3 8

4. Metode Bode Metode integral Bode menggunakan lima titik ketinggian yang berjarak sama, = tinggi rata-rata, 0 0 + 1 1 + 2 2 + 3 3 + 4 4 0 + 1 + 2 + 3 + 4 Dengan 0 = 7, 1 = 32, 2 = 12, 3 = 32, 4 = 7, sehingga =
4

dan

=
4

( ) = + (galat) + (7 )

=
0

( ) =

2 7 0 + 32 1 + 12 2 + 32 3 + 7 4 45

5. Integral Komposit Metode integral komposit yaitu metode menghitung integrasi numerik suatu fungsi dengan membaginya dalam selang-selang tertentu menjadi segmen-segmen luasan sebanyak N, selanjutnya segmen-segmen tadi dihitung luasnya menggunakan metode trapesium atau Simpson kemudian dijumlahkan untuk menghitung integrasi fungsi. Karena pembagian menjadi segmen-segmen, menggunakan integral komposit mengaruskan kita membuat vektor yang merepresentasikan sumbu X. Makin banyak segmen yang dibuat maka solusi numerik yang diperoleh akan makin mendekati solusi sebenarnya.
1 2

() =
0

( ) +
1

() + +
1

( )

Jika luasan tiap segmen dihitung menggunakan metode trapesium, diperoleh: = 0 + 1 + 1 + 2 + + 1 + 2 2 2 = 0 + 2 2

( ) + ( )
=1

Misalkan kita ingin menghitung integral fungsi () dalam selang (, ) dengan pembagian selang menjadi 100 upaselang (ingat, 100 upaselang berarti ada 100+1 titik). Berikut algoritmanya: 1) a batas bawah; // definisikan batas bawah integrasi 2) b batas atas; // definisikan batas atas integrasi 3) h = (b a)/100; // definisi lebar upaselang 4) x0 = a // membuat vektor X for i = 1, i <= 100 // bernilai awal a, naik dengan beda h, hingga b xi = x0 + h; yi = f(xi); // membuat vektor Y (hasil fungsi) end 5) for i = 0, i < 100 // perulangan dari i = 0 hingga 99 ui = (h/2)*(fi + fi+1) // menghitung luasan segmen end 6) L0 = 0 // menghitung luas total segmen for i = 1, i <= 100 Li = Li-1 + ui-1; end 7) I = L100 // diperoleh hasil integrasi = luas total Jika program/bahasa pemrograman yang digunakan mendukung fungsi barisan, maka poin (4) dapat diganti dengan x = a:h:b, atau dalam Matlab menyediakan fungsi x = linspace(a,b,101). Jika program menyediakan fungsi sumasi, maka poin (7) dan (8) dapat dihilangkan dan diganti menjadi I = sum(y) Patut diingat pula indeks yang digunakan, apakah berjalan dari 0 ataukah dari 1. Pada algoritma di atas, indeks berjalan dari nol, sehingga vektor X = (x0 x1 x2 x100), jika indeks berjalan dari satu, maka vektor X akan menjadi X = (x1 x2 x3 x101)

Lampiran: 1. Contoh program integral numerik dengan metode Simpson menggunakan bahasa C++: Fungsi: = 2 2 + 5 3.
//Integral Numerik Simpson by skaga #include<iostream> using namespace std; int main() { int i, kol=3; float x[kol], y[kol], a, b, deltax, h , ratay, L; char off; cout << "Masukkan batas bawah (a) : "; cin >> a; cout << "Masukkan batas atas (b) : "; cin >> b; deltax = b - a; h = deltax/2; x[0] = a, x[1] = a+h, x[2] = a+2*h; for (i=0; i<3; i++){ y[i]=2*x[i]*x[i] + 5*x[i] - 3; // masukkan fungsi } ratay = (y[0] + 4*y[1] + y[2])/6; L = deltax*ratay; cout << "integral y(x) = " << L << '\n'; cin >> off; return 0; }

2. Contoh program integral numerik metode Simpson 3/8 menggunakan Matlab: % Integral numerik dengan metode Simpson 3/8 a = input('masukkan batas bawah = '); b = input('masukkan batas atas = '); deltax = b - a; h = deltax/3; x = [a (a+h) (a+2*h) b]; y = input('masukkan fungsi y=f(x) : '); y = [1 3 3 1].*y; ratay = sum(y)/8; L = deltax*ratay; fprintf('integral fungsi y(x) dari x = %g hingga x = %g ialah %g \n', a, b, L);

3. Contoh perhitungan integral numerik metode trapesium komposit menggunakan MS Excel. = 2 sin ; = 2, = 3, = 0.1
i 0 1 2 3 4 5 6 7 8 9 10 sum x(i) 2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 y(i) 3.6371897 3.8067533 3.9131226 3.9447806 3.8906679 3.7404509 3.4847893 3.1155993 2.6263071 2.0120869 1.2700801 7.443943014 7.719875901 7.857903167 7.835448492 7.631118821 7.225240174 6.6003886 5.741906424 4.638393956 3.282166931 0.372197151 0.385993795 0.392895158 0.391772425 0.381555941 0.361262009 0.33001943 0.287095321 0.231919698 0.164108347 3.298819274 2 sin = 2 cos + 2 sin + 2 cos + y(i-1) + y(i) (h/2)*(x(i-1) + x(i))

Metode solusi

analitik

memberikan

Sehingga

3 2 2

sin = 3,30719

Dengan program Matlab pada lampiran 2, diperoleh hasil I = L = 3,3091.

4. Contoh program plot integral numerik fungsi Panas jenis zat padat sebagai fungsi temperatur berdasarkan model Debye diberikan oleh fungsi 3 / 4 = 9 3 1 2 0 Dengan temperatur Debye (konstanta), N jumlah atom, dan k tetapan Boltzmann. Plot fungsi CV terhadap T dibuat dengan memerikan vektor T, memasukkan fungsi dan mengintegralkannya, kemudian diplot. Berikut programnya menggunakan Matlab.
% Perhitungan Panas Jenis Zat Padat model Debye menggunakan Simpson 3/8 % @Sunkar Eka Gautama, 2011 clear; theta = input('nilai theta(Debye) = '); batasatas = input('batas atas temperatur = '); T = linspace(0,batasatas/theta,300); % T = T/theta t = 1./T; h = t/3; a = 0.*h; b = h; c = 2*h; d = 3*h; y = (3.*h./8).*(0 + 3.*(b.^4.*exp(b)./(exp(b)-1).^2)... + 3.*(c.^4.*exp(c)./(exp(c)-1).^2) + (d.^4.*exp(d)./(exp(d)-1).^2)); Y = 9.*T.^3.*y; plot(T,Y); title('Kurva Panas Jenis Zat Padat (CV) Model Debye'); xlabel('temperatur'); ylabel('CV (dalam R)');

Daftar Pustaka: Suarga, Fisika Komputasi: Solusi Problema Fisika dengan Matlab, Penerbit ANDI, Yogyakarta, 2005