function y = bisektion(f,a,b, delta)%%%%%%%
%%%%%%%%%  Wir suchen eine Approximation einer  Nullstelle von f, die zwischen a,b ist, mit
%%%%%%%%%  Fehler kleiner als delta
 if f(a)==0, y=a,    %%%%%%%%% Wenn f(a)=0 oder f(b)=0 dann ist die
 elseif f(b)==0, y=b,   %%%%%% Nullstelle a bzw. b
 elseif abs(b-a)<delta, y=(a+b)/2; %%%%%%%% |b-a| klein ist, dann ist der Fehler klein
 elseif f(a)*f((a+b)/2) <=0 y=bisektion(f,a,(a+b)/2, delta);%%%%%%% In diesem Fall wird Nullstelle in der Linke Hälfte gesucht
 elseif f(b)*f((a+b)/2) <=0 y=bisektion(f,(a+b)/2,b, delta);%%%%%%% In diesem Fall wird Nullstelle in der Linke Hälfte gesucht
end

%%%%%%%%% Das Programm funktioniert nur, wenn a<b, f reell und f(a)*f(b)<=0

%f=inline('t^2-2')

%f =

 %    Inline function:
 %    f(t) = t^2-2

%f(2)

%ans =

%     2


 
%bisektion(f,1,2,0.01)

%ans =

 %   1.4180

%format long
%bisektion(f,1,2,0.01)

%ans =

 %  1.417968750000000

%bisektion(f,1,2,0.00001)

%ans =

 %  1.414211273193359

%sqrt(2)

%ans =

 %  1.414213562373095

%bisektion(f,1,2,0.0000001)

%ans =

%   1.414213567972183