%% Lecture 7:            Bisektionsverfahren
%                        Skriptfile

% Fehlertoleranz: halbe Intervallänge
 tol = 1e-4;
 
 a=0;             % Definition des Ausgangsintervalles
 b=pi/2;
 
 %fa = exp(-a)-sin(a);  % Funktionswerte an Intervallenden
 x = a;
 func
 fa = y;
 fb = exp(-b)-sin(b);
 
 if fa*fb <=0   % Intervallbedingung erfüllt?
     while (b-a)/2 > tol
         m = (a+b)/2;
         fm = exp(-m)-sin(m);
         if fm*fa <=0
             b=m;
             fb=fm;
         else
             a=m;
             fa = fm;
         end
     end % end of while
     
     m = (a+b)/2;
     fm = exp(-m)-sin(m);
     disp('Nullstelle bei x= '), m
     disp('Funktionswert f(x) = '), fm
 else
     disp('unzulässsiges Intervall')
 end
 
 