SecantRoot1

function x = koreniSecice(f, x1, x2, itMax, errMax)
% Funkcija za racunanje korena funkcije f nad
% intervalom x1, x2
%   f funkcija ciji koren trazimo
%   x1 pocetak intervala
%   x2 kraj intervala
%   itMax broj iteracija
%   errMax maksimalna greska
%   x resenje
%primer
% >> x = koreniSecice(@nasaFunkcija, -1, 1, 100, 0.001)
    for it = 1:itMax
        y1 = feval(f, x1);
        y2 = feval(f, x2);

        a = (y1-y2) / (x1-x2);
        b = y1 - a*x1;
        x = -b/a;

        y = feval(f, x);
        if abs(y) < errMax
            return;
        end
        if y1 * y < 0
            x2 = x;
        else
            x1 = x;
        end
    end
end