% Modes et vecteurs propres
% on obtient:  n_modes = nombre de modes calculés
%              freq = fréquences propres (Hz)
%              Egv = vecteurs propres (eigenvectors)

ex1_mass;  % on ajoute d'abord les masses ponctuelles

b = (n_fix-1)+[1:6];   % on fixe le noeud 1
[L,Egv] = eigen (K,M,bc);
freq = sqrt(L)/(2*pi);
disp ('fréq (Hz) =')
disp (freq(1:4))
n_modes = length (freq)

for i = [1:3]
  plt_mode (Coord,Ex,Ey,Ez,Edof,real(Egv),real(freq),i,0.01);
end