%bisection_methods
syms x %deklarasi fungsi x
fx = x^2+x-3;
xl = 1; %tebakan bawah
xu = 2; %tebakan atas
xm = (xl+xu)/2;
fxl = subs (fx,x,xl); %mensubtitusikan xl ke fxl
fxu = subs (fx,x,xu); %mensubtitusikan xl ke fxu
if fxl*fxu > 0 %syarat program diteruskan
break
end
es = 0.0001; %eror batas
ea = 1; %eror komputasi
n = 100; %banyaknya pengulangan
N = [i xl xu xm ea]; %matriks penyimpan data
for i = 1 : n
xm = (xl+xu)/2;
fxl = subs (fx,x,xl); %mensubtitusikan xl ke fxl
fxu = subs (fx,x,xu); %mensubtitusikan xl ke fxu
fxm = subs (fx,x,xm);
if fxl*fxm < 0
ea = abs ((xm-xu)/xm);
xu = xm;
else
ea = abs ((xm-xl)/xm);
xl = xm;
end
if ea < es
break
end
N = [N;i xl xu xm ea];
end
N
0 komentar:
Posting Komentar