% Test "m" file	
a = 2.1111111111
b = 3.121212121
c = 1
epsilon = 10^(-25)
N = log((b-a)/epsilon)/log(2) -1
for i=1:N
   i
   b = (a+b)/2
end
a
b
