function A=MasseElem(d,vol)
  A=(ones(d+1)+eye(d+1))*vol/((d+1)*(d+2));
  %A=(ones(d+1))*vol/((d+1)*(d+2));
  %for i=1:d+1,A(i,i)=2*A(i,i);end
end
