function [x,X] = bisekcija(f,a,b,N)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Vhodni podatki:
% f ........ funkcija, ki ji iscemo niclo
% a, b ..... krajisci intervala, na katerem iscemo niclo
% N ........ stevilo korakov bisekcije
%-------------------------------------
%Izhodni podatki:
% x ........ priblizek za niclo
% X ........ vektor dimenzije Nx1 vmesnih priblizkov
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fa = f(a);
fb = f(b);
if sign(fa) == sign(fb)
    error('Funkcija enako predznacena v a in b')
end

X = zeros(N,1);
s = b-a;
for k = 1:N
    s = s/2;
    x = a + s;
    fx = f(x);
    if sign(fx) == sign(fa)
        a = x;
        fa = fx;
    end
    X(k) = x;

end