function V=vdm1(t,n)

m = length(t);
V = zeros(m,n+1);    %on initialise la matrice V pour gagner du temps

for i=1:m
  for j=1:n+1
    V(i,j)=t(i)^(j-1);
  end
end