function [z, iters] = regfalsi(a,b, epsilon,max_iter,fnc)
x0 = a;
xk = x0;
x = xk;
fxk = eval(fnc);
x = b;
fb = eval(fnc);
xkol = x;
xk = xk - (fxk/(fb-fxk))*(b-xk);
z=[xk];
i = 0;
while abs(xkol - xk) > epsilon && i < max_iter,
    i = i + 1;
    xkol = xk;
    x = xk;
    fxk = eval(fnc);
    xk = xk - (fxk / (fb - fxk)) * (b - xk);
    z = [z; xk];
end
iters = i;