vnd(x)

              function v = vnd(x)

              x = x(:);              % x doit tre en colonne
              v = ones(size(x));     % v est au dpart une colonne de 1
              for i = 1:length(x)-1  %  chaque tape, on ajoute devant
                v = [ v(:,1).*x v];  % la premire colonne de v 
              end;                   % multiplie terme  terme par x.