clc,clear
disp('metode biseksi')
disp('untuk mencari akar persamaan')
disp('========================')
disp('pers nonlinear')
f=inline('x^2-3','x')
ezplot(f)
grid on
a=input('a = ')
b=input('b = ')
if f(a)*f(b)>0
disp('interval tdak mengapit akar')
disp('ulangi lagi masukkan interval')
elseif f(a)*f(b)<0 p=""> e=input('masukkan batas toleransi maks abs(f(xt)) = ')
disp('=================================')
disp('')
xa(1)=a
xb(1)=b
xt(1)=(xa+xb)/2
g(1)=nan
i=1
while abs(f(xt(i)))>e
if f(xa(i)*f(xt(i)))<=0
xa(i+1)=xa(i)
xb(i+1)=xt(i)
else
xa(i+1)=xt(i)
xb(i+1)=xb(i)
end
xt(i+1)=(xa(i+1)+xb(i+1))/2
g(i+1)=xt(i+1)-xt(i)
i=i+1
end
disp('hasil komputasi')
disp('iterasi xa xb xt f(xa) f(xb) f(xt) galat')
disp([(1:i)' xa' xb' xt' f(xa)' f(xb)' f(xt)' g'])
disp('==========================')
fprintf('akar persamaanya adalah=%3.4f\n',xt(i))
end
disp('metode biseksi')
disp('untuk mencari akar persamaan')
disp('========================')
disp('pers nonlinear')
f=inline('x^2-3','x')
ezplot(f)
grid on
a=input('a = ')
b=input('b = ')
if f(a)*f(b)>0
disp('interval tdak mengapit akar')
disp('ulangi lagi masukkan interval')
elseif f(a)*f(b)<0 p=""> e=input('masukkan batas toleransi maks abs(f(xt)) = ')
disp('=================================')
disp('')
xa(1)=a
xb(1)=b
xt(1)=(xa+xb)/2
g(1)=nan
i=1
while abs(f(xt(i)))>e
if f(xa(i)*f(xt(i)))<=0
xa(i+1)=xa(i)
xb(i+1)=xt(i)
else
xa(i+1)=xt(i)
xb(i+1)=xb(i)
end
xt(i+1)=(xa(i+1)+xb(i+1))/2
g(i+1)=xt(i+1)-xt(i)
i=i+1
end
disp('hasil komputasi')
disp('iterasi xa xb xt f(xa) f(xb) f(xt) galat')
disp([(1:i)' xa' xb' xt' f(xa)' f(xb)' f(xt)' g'])
disp('==========================')
fprintf('akar persamaanya adalah=%3.4f\n',xt(i))
end