% biseccion(f, a, b, e, N)
% dada f continua, que se supone que cambia de signo entre a y b
% devuelve un numero r:
% si |f(c)|< e antes de N pasos, entonces c = r
% si no, c = -inf
%
% esta funcion NO COMPRUEBA que f cambie de signo
function [c n] = biseccion(f, a, b, e, N)
  % vale empezar en 1 o en 0, no pasa nada
  n = 1;
  c = (b + a)/2;
  % valor absoluto: abs
  while (abs(f(c)) > e && n < N)
    if (f(c)*f(a) < 0)
      b = c;
    else
      a = c;
    end
    c = (b + a)/2;
    n = n + 1;
  end
  % OJO: lo contrario de n<N es n>=N
  if n >= N
    c = -inf;
  end
end
