bisektionsSchritt.m

function [a1, b1] = bisektionsSchritt(f, a, b)
% macht einen Schritt des Bisektionsverfahrens für die Funktion f
% ausgehend vom Intervall [a, b].
% Liefert das neue Intervall zurück.

sa = sign(f(a));

xm = (a + b)/2;
sm = sign(f(xm)); 
if sm == sa
  a = xm;
else
  b = xm; 
end

if a < b
  a1 = a;
  b1 = b;
else
  a1 = b;
  b1 = a;
end