% ------ topologie -----------------------------------------------
[n_nodes,n_dof,n_elem,n_nel,Dof,Edof] = topol (Coord,Elem);
[Ex,Ey,Ez] = coordxtr(Edof,Coord,Dof,n_nel);

% Check elements
Number_of_elements = size(Elem,1)
Check = length([ eFlex ePoutre])
  
Coord_xy(:,1) = Coord(:,1);	Coord_xy(:,2) = Coord(:,2);
figure; femdraw2 (Coord_xy,Ex,Ey,[1 4 2 5 0.]); 
ylabel('y'); grid;