エカつきのブログ = Eka tsuki no blog

clc;
clear;
%f(x)=exp(x)-9*x
a=0;
b=1;
epsilon = 0.0000001;
iterasi=0;
disp('METODE BAGI DUA untuk persamaan f(x)=exp(x)-9*x');
disp('a=0, b=1, epsilon=0.0000001');
disp('Iterasi a b c f(a) f(b) f(c) selang baru lebar');
while abs(b-a) > epsilon;
c=(a+b)/2;
fc=exp(c)-9*c;
fa=exp(a)-9*a;
fb=exp(b)-9*b;
fprintf(' %3g %3.4f %5.4f %5.4f %10.7f %10.7f %10.7f',iterasi,a,b,c,fa,fb,fc);
if fa*fc < 0;
b=c;
selang = '[a,c]';
else;
a=c;
selang = '[c,b]';
end;
fprintf(' %8s %10.7f\n',selang,abs(b-a));
iterasi=iterasi+1;
end;
akar = c;
fprintf('Akarnya adalah = %8.7f\n', akar);
fprintf('Jumlah Iterasi = %g\n', iterasi);