function M=MassePoidsElem(vol,w)
  % w (d+1)-by-1 
  w=w(:);
  d=size(w,1)-1;
  M=vol/((d+1)*(d+2)*(d+3))*(ones(d+1)+eye(d+1)).*(sum(w)+(w+w.'));
end
