function res = regulafalsi(xim1,xi,tol)
    i = 0;                                              % # iterace
    g = 1;                                              % f(koren)
    while(g > tol)
        i = i + 1;
        xip1 = (funkce(xi)*xim1-funkce(xim1)*xi)/(funkce(xi)-funkce(xim1));
        if(funkce(xip1)*funkce(xim1) > 0)
            xi = xip1;
            g = abs(funkce(xi));
            res = xi;
        else
            xim1 = xip1;
            g = abs(funkce(xim1));
            res = xim1;
        end
        fprintf('%3d%17.8f%17.8f%17.8f%17.8f%17.8f%17.8f\n', ...
            i,xim1,xi,xip1,funkce(xim1),funkce(xi),funkce(xip1));
    end
end
