% ------ generate element matrices, assemble in global matrices 
nd = n_nodes*n_dof;
clear K M C
K = zeros(nd); M = zeros(nd); C = zeros(nd); 
nr_elem = 0;

loc = 'elements flex';
for ie = eFlex
  eo(ie,:) = [0 0 1];  
  [ke,me] = beam3d (Ex(ie,:),Ey(ie,:),Ez(ie,:),[0 0 1],Ep_flex);
  K = assem(Edof(ie,:),K,ke);
  M = assem(Edof(ie,:),M,me);
  nr_elem = nr_elem+1;
  end

loc = 'elements poutre';
for ie = ePoutre
  eo(ie,:) = [0 0 1];
  [ke,me] = beam3d (Ex(ie,:),Ey(ie,:),Ez(ie,:),eo(ie,:),Ep_poutre);
  K = assem(Edof(ie,:),K,ke);
  M = assem(Edof(ie,:),M,me);
  nr_elem = nr_elem+1;
  end

Check = nr_elem