MatLab / Dihotomia

a = 1;
b = 2;
n = 1;
eps = 0.000001;
k = 0;
while (abs(b - a) > eps)
fa = fun(a);
fb = fun(b);
if (fa*fb < 0)
c = b - (b - a)/2;
fc = fun(c);
if (fc == 0)
disp(c)
disp(n)
k = 1;
break
end
if (fa*fc < 0)
b = c;
else
a = c;
end
else
if (fa == 0)
disp(a)
disp(n)
k = 1;
break
else
if (fb == 0)
disp(b)
disp(n)
k = 1;
break
else
disp('Roots do not exist')
disp(n)
k = 1;
break
end
end
end
n = n+1;
end
if (k == 0)
disp(c)
disp(n-1);
end