Contoh Program - METODE POSISI PALSU - METODE NUMERIK

Dibawah adalah sintax metode numerik dengan metode posisi palsu (regulasi falsi). Yang kebetulan software yang penulis gunakan adalah program QtOctave di linux, jadi yang belum ada install seperti post ini.

Catatan : sintax yang di blok hijau, ganti sesuai soal yang di dapat . wkwkwk

Langsung aja di copi om di kode editor nya!!

clear;

function y=f1(x)

y=x^2-x-3;

return;

end

disp('===============================')

disp('*** METODE POSISI PALSU (REGULAFALSI) ***')

disp(' ** CREATED BY : SUDAHLAHHH!! **')

disp('===================================')

iterasi=1;

a=input('Masukkan interval (a0) : ');

aawal=a;

b=input('Masukkan interval (b0) : ');

bawal=b;

epsilon=input('Masukkan Nilai Epsilon : ');

clama=2*b-a;

c=b-(f1(b)*(b-a))/(f1(b)-f1(a));

error=abs((c-clama)/c);

if error <= epsilon

disp('Nilai akar adalah : %2.8f',c);

else

disp('===============================================================================') 

disp('Iterasi a b w f(w) cek selang interval baru [a,b]') 

disp('===============================================================================')

while error > epsilon

cekselang=f1(a)*f1(c);

c=b-(f1(b)*(b-a))/(f1(b)-f1(a));

error=abs((c-clama)/c);

if cekselang < 0

bbaru=c;

abaru=a;

else

abaru=c;

bbaru=b;

end

fprintf('%2d %15.8f %14.8f %13.8f %13.8f %13.8f [%3.8f , %3.8f]\n',iterasi,a,b,c,cekselang,error,abaru,bbaru);

clama=c;

a=abaru;

b=bbaru;

iterasi=iterasi+1;

end

end

disp('===============================================================================') 

fprintf('Hampiran akar pada interval [%1.8f , %1.8f] dengan epsilon = %f adalah : %.8f \n\n',aawal,bawal,epsilon,c);