%Ce fichier recopie à peu de choses près le fichier hilb.m
%Il utilise une astuce : "Tony's trick"
function V=vdm(t,n)
t=t(:); %on met t en colonne si t est par erreur vecteur ligne
m = length(t);
j = [1:n+1];
V = t(:,ones(1,n+1)) .^( j(ones(m,1),:) -1);