    function [a1,b1]=inthalv(f,a,b,tol)
    % Intervallhalveringsmetoden för f(x)=0, 
    % med startintervall [a,b]
    % Halvera tills intervallet blir mindre än tol
    a1=a; b1=b; fb=feval(f,b);
    while b1-a1 > tol,
      disp([a1 b1])
      x=(a1+b1)/2;
      if feval(f,x)*fb> 0,
        b1=x;
      else
        a1=x;
      end
    end

