METODE BAGI DUA
Sebelumya tentukan fungsi yang akan kita cari hampiran akarnya, taru di m-file misalnya :
% m-file fu.m
function y = f(x);
y = ((x+1)^2)* exp(x^2-2)-1
Fungsi fu itu nantinya akan dipanggil di program metode bagi duanya. Misalnya nama m-filenya bagidua.m (Mencoba menganalisis kodenya nih.. tolong koreksi ya bila salah…) :
% bagidua.m
a = 0.4; b = 1; tol = 0.0000000001; N = 15;
%keterangan: a, b = batas selang [a, b].
% tol = batas toleransi (lihat ***). N = banyak iterasi
hasil = [];
% hasil adalah matriks kosong yang
% disiapkan untuk menympan hasil perhitungan qta
for i = 1 : N, % sepanjang iterasi (yaitu dari 1 sampai N) lakukan :
x = (a + b)/2; % cari titik tengahnye…
hasil = [hasil; i a b x fu(x)];
% panggil fungsi fu dari fu.m
% isi matriks hasil dengan kolom-kolom i, a, b, fu(x)
% banyaknya baris matriks hasil adalah sebanyak N (jumlah iterasi)
if (fu(x) == 0| (b-a) < tol),
break; % jika syarat dah dipenuhi, hentikan iterasi (walau belum N gitu ga?)
end,
if (fu(a) * fu(x) < 0),
b = x;
else a = x;
end,
end
*** : Toleransi itu semacam bilangan yang kecil banget gitu ya untuk nentuin hampiran akar mau diambil hingga ketelitian berapa angka di belakang koma. Karena kadang pencarian akar ga bener2 bisa menemukan f(x) = 0… maybe..
Untuk yang metode regular falsi, tinggal mengganti baris x = (a + b)/2;
dengan x = (a * fu2(b) – b * fu2(a)) / (fu2(b) – fu2(a)) ; gitu ya..?
…
Allohua’lam
lab mat, 190908 3.35
