a=1 ;  %venstre endepunkt for intervallet
b=2;   %høyre endepunkt for intervallet

N=20  ;  %antall itterasjoner


f=@(x) x.^5 - 3.*x -1  ;   %Funksjonen vi skal finne nullpunkt til

v=linspace (a,b, 100);
plot(v,f(v), 'r ' ) ;
hold on 
format long
for n=1:N
    c=(a+b)/2 ;
    plot(c,0,'linewidth', 5 )
    
    if f(c)==0 ;
     fprintf(   'Et nullpunkt er %d \n' , c )
    end
        
    if f(a)*f(c) > 0 ; 
        a=c ; 
        
    else b=c ;
    end
    
end
   hold off 
     
fprintf( [num2str(a, ceil(N/3)+1 ) ' < nullpunkt < \n' num2str(b, ceil(N/3) +1)  '\n'  ])
%fprintf(  'Bredden på dette intervallet er %d ' , b -a)

%Vi har forsøkt å skrive ut resultatet på en mer 
%brukervennlig måte.  De to grensene er nå rett
%over hverandre slik at det er lettere å sammenligne dem
%Vi har og forsøkt å tilpasse antall siffer til nøyaktigheten til
%utregningen.  ceil tar det minste heltallet større enn eller lik sitt
%argument. 

 %fprintf(   'Nullpunktet er større enn %d \n' , a )
  %fprintf(  'Nullpunktet er mindre enn eller lik  %d \n' , b )
% fprintf(  'Bredden på dette intervallet er %d ' , b -a)