cog_faces_alu3d_hexa.m

function cog = cog_faces_alu3d_hexa(M);
%function cog = cog_faces_alu3d_hexa(M);
%
% function computing the cog of the boundary_faces of a hexaeder alu3d mesh
  
% Bernard Haasdonk 16.3.2006
  
  [XX,YY,ZZ] = coord_faces_alu3d_hexa(M);
  X = sum(XX)/4;
  Y = sum(YY)/4;
  Z = sum(ZZ)/4;
  cog = [X;Y;Z];
  
  

% TO BE ADJUSTED TO NEW SYNTAX
%| \docupdate