a=-2 ; %en verdi til venstre for et nullpunkt
b=3 ; %en verdi til høyre for et nullpunkt
N=30 ;  %antall iterasjoner

f =@(x) 3*x -2 ; 
%exp(x) -3*x ;
% x^2 -1
%x^2 -2 
%x^3 -x 
% x /abs(x)   
if f(a)*f(b) >= 0 
    disp('funksjonsverdiene i a og b har samme fortegn')
else

for n=1:N
    c=(a+b)/2; %punktet midt mellom a og b
   if f(c)  == 0
       disp(['et nullpunkt er ' num2str(c)])
       a=c ;
   elseif  f(c) * f(a) > 0 
       a=c ;
   else
       b=c ;
   end
end

(a+b)/2
end



