function A=Rigidite(nq,me,vols,G)
  % G=FEM.GradBaCoS(q,me); 
  % vols=FEM.volumes(q,me);
  nme=length(vols);
  A=sparse(nq,nq);
  for k=1:nme
    I=me(:,k);
    A(I,I)=A(I,I)+simplex.RigiditeElem(vols(k),squeeze(G(k,:,:)));
  end
end
