function nicla=bisekcija(f,a,b,eps)
%BISEKCIJA vrne niclo funkcije po metodi bisekcije
%nicla=BISEKCIJA(f,a,b,eps)
%f je funkcija, katere niclo iscemo
%a in b sta krajisci intervala
%eps je toleranca

while abs(b-a)>eps
  c= a+(b-a)/2; %bolje kot c=(a+b)/2!!!
  if sign(feval(f,a))==sign(feval(f,c))
    a=c;
  else
    b=c;
  end
end
nicla=a+(b-a)/2;