function ranges = intervs( f, range, nbInter)
% ranges = intervs( f, range, nbInter);
% tant donns f et un intervalle range = [x1, x2] la function divise
% l'intervalle en nbInter sous-intervalles gaux, et cherche les passages
% par zros de la fonction.
% ranges est un tableau 2xn (n <= nbinter) chaque colonne dfinissant
% un intervalle contenant au moins une racine (si f est continue).
%
ranges = linspace(range(1), range(2), nbInter+1);
ranges = [ranges(1:nbInter); ranges(2:nbInter+1)];
fv = feval(f, ranges);
ind = find(fv(1,:).*fv(2,:) < 0.0);
ranges = ranges(:,ind);
