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

I just learn to be good one...^_~,

clc;
clear;
%f(x) = exp(x) - 9*x
a=0;
b=1;
epsilon= 0.0000001;
iterasi=0;
disp('METODE REGULA FALSI 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
fa=exp(a)-9*a;
fb=exp(b)-9*b;
c=b-(fb*(b-a)/(fb-fa));
fc=exp(c)-9*c;
fprintf(' %3g %5.4f %5.4f %5.4f %10.7f %10.7f %10.7f ',iterasi,a,b,c,fa,fb,fc);
if abs(fc)<epsilon;
a=c;
b=c;
elseif fa*fc<0
b=c;
selang = '[a,c]';
else
a=c;
selang='[c,b]';
end;
fprintf(' %8s %9.8f\n',selang,abs(b-a));
iterasi=iterasi+1;
end;
akar = c;
fprintf('Akarnya adalah = %8.7f\n', akar);
fprintf('Jumlah Iterasi = %g\n', iterasi);

0 komentar:

Post a Comment