% Aufgabe 32-7 (n Elemente)
clear all

EI   = 3000 ;
rhoA = 3    ;
L    = 1    ;

n = 20 ;                                               % Anzahl der Elemente  
for i=1:n+1                    
    xy(i,1)=(i-1)*L/n ;                                % Knotenkoordinaten  
end
for i=1:n
    km(i,:) = [ i  i+1] ;                              % Koinzidenzmatrix
    ep(i,:) = [EI rhoA] ;                              % Elementparameter
end
kr = zeros (n+1,2) ;                                   % Matrix der Randbedingungen
kr(1,:) = [1 1] ;                                      % Einspannung am Knoten 1
mid = round (n/2) + 1 ;                                % Nummer des Knotens in der Mitte
kr(mid,1) = 1 ;                                        % Lager in der Mitte
 
% Finite-Elemente-Algorithmus:
[nc omega x] = femeig_m (3 , xy , km , ep , kr) ;      % 3 Eigenkreisfrequenzen
if nc > 0
    f = omega/(2*pi)                                   % Frequenzen -> Command Window
    draweigvecs (xy , km , kr , x) ;                   % Schwingungsformen zeichnen
end