function A=Masse(nq,me,vols)
  % vols=FEM.volumes(q,me);
  nme=length(vols);
  d=size(me,1)-1;
  A=sparse(nq,nq);
  for k=1:nme
    I=me(:,k);
    A(I,I)=A(I,I)+simplex.MasseElem(d,vols(k));
  end
end
