% Implementering av halveringsmetoden for likninga f(x)=0 der
% f er kontinuerlig på [a,b] og f(a) og f(b) har ulike fortegn.

% Grenser
a=1;
b=2;

% Gir funksjonen som skal være null
NullFunk=@(x) x.^5-3*x.^3+x-1;

% Gir antall iterasjoner
N=20;

% Funksjonsverdier
Fa=NullFunk(a);
Fb=NullFunk(b);

% Starter for-løkke som kjøres N ganger
for i=1:N
  c=(a+b)/2;        % Midtpunktet
  Fc=NullFunk(c);   % Funksjonsverdi i midtpunktet
  if Fa*Fc<0
    b=c;            % Setter ny b til c
  else
    a=c;            % Setter ny a til c
  end
end

% Regner ut nytt midtpunkt og skriver svaret til skjerm
x=(a+b)/2