% Differenzenverfahren, Biegeschwingungen
% (Beispiel 1, 5 Gleichungen)

clear all

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

nA  = 6    ;               % Anzahl der Abschnitte
h   = L/nA ;
n   = nA-1 ;

A = [7 -4  1  0  0 ;
    -4  6 -4  1  0 ;
     1 -4  6 -4  1 ;
     0  1 -4  6 -4 ;
     0  0  1 -4  5 ] ;

D = eig(A) ;
%for i=1:3
   Eigenfrequenzen = sqrt(D*EI/(rhoA*h^4))/(2*pi)
   %end