simps(f, range, epsi)

          function int = simps(f, range, epsi)
          %simps  Methode de Simpson adaptative. version simple
          %         int = simps(f, range, epsi)
          %       calcul de l'intgrale de f sur l'intervalle [range(1), range(2)]
          %       par la mthode de simpson  la prcision epsi par dcoupage adaptatif.
          %       Entrees
          %            f           : Fonction  intgrer;
          %            range       : intervalle d'intgration
          %            {epsi}      : prcision dsire
          %       Sorties
          %            int         : rsultat de l'intgration

          int = []; 
          a = range(1); 
          b = range(2); 
          if nargin < 3 
            epsi = 0.5.*sqrt(eps)./(b-a); 
          else 
            epsi = 0.5.*epsi./(b-a); 
          end; 
          oldint = (feval(f,a) + 4.*feval(f,(a+b).*0.5) + feval(f,b)).*(b-a)./6.0;
          int = simpsau0(f, range, oldint, epsi);