function [c] = bisection2(f, a, b, eps_step)

format long

while (b - a >= eps_step)
    c = (a + b)/2;
    if ( f(c) == 0 )
       break;
    elseif ( f(a)*f(c) < 0 )
       b = c;
    else
       a = c;
    end
end
