% Aufgabe 32-9 (beliebig feine Diskretisierung)
clear all

EI   = 5000  ;
rhoA = 4.8   ;
a    = 0.6   ;
L    = 1.4   ;
m    = 4     ;
cT   = 60000 ;

n = 14 ;                                               % 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 0] ;                                    % Lager am Knoten 1
kr(n+1,:) = [1 1] ;                                    % Einspannung rechts
nMasse = round (n*a/L + 1.1) ;                         % Nummer des Knotens mit der Masse
mk = zeros (n+1,2) ;                                   % Matrix der Knotenmassen
mk(nMasse,1) = m ;                                     % Masse m am rechten Knoten
sc = zeros (n+1,2) ;                                   % Matrix der Federsteifigkeiten
sc(1,2) = cT ;
 
% Finite-Elemente-Algorithmus:
[nc omega x] = femeig_m (3 , xy , km , ep , kr , mk , sc) ; % 3 Eigenkreisfrequenzen
if nc > 0
    f = omega/(2*pi)                                   % Frequenzen -> Command Window
    draweigvecs (xy , km , kr , x) ;                   % Schwingungsformen zeichnen
end