% Aufgabe 32-9 (2 Elemente):
clear all

xy = [ 0  ;  0.6  ;  1.4 ] ;  % Koordinaten der drei Knoten 1, 2 und 3
km = [  1  2   ;   2  3  ] ;  % Koinzidenzmatrix: Elem. 1 -> Knoten 1 u. 2, Elem. 2 -> Knoten 2 u. 3
ep = [5000 4.8 ; 5000 4.8] ;  % Elementparameter: EI und rhoA fuer beide Elemente
kr = [1 0     ; 0 0 ; 1 1] ;  % Verhinderte Verformungen: Knoten 1: v, Knoten 3: v und phi
mk = [0 0     ; 4 0 ; 0 0] ;  % Einzelmassen an den Knoten: m am Knoten 2
sc = [0 60000 ; 0 0 ; 0 0] ;  % Federn: Drehfeder am Knoten 1
 
[nc omega] = femeig_m (1 , xy , km , ep , kr , mk , sc) ;  % Berechnung einer Eigenkreisfrequenz
Eigenfrequenz = omega/(2*pi)                               % Eigenfrequenz -> Command Window