% Eigenschwingungen eines Rahmens (FEM-Modell)
clear all
xy = [0 0 ; 0 2 ; 0 4 ; 0 6 ; 4 0 ; 4 2 ; 4 4 ; 4 6] ;
km = [1 2 ; 2 3 ; 3 4 ; 2 6 ; 3 7 ; 4 8 ; 5 6 ; 6 7 ; 7 8] ;
for i=1:9  ep(i,:) = [3.5e6 2.9e9 100] ; end
kr = [1 1 1 ; 0 0 0 ; 0 0 0 ; 0 0 0 ; 1 1 1 ; 0 0 0 ; 0 0 0 ; 0 0 0] ;
[succ K M] = eigwbc_m (xy , km , ep , kr) ;
KQ = SymBand2Sparse (K) ; MQ = SymBand2Sparse (M) ;

[Z D] = eigs (KQ , MQ  , 5 , 'sa') ;

Eigenkreisfrequenzen = (sqrt(diag(D)))'
Eigenvektoren = Z
draweigvecs (xy , km , kr , [zeros(3,5) ; Z( 1: 9,1:5) ; ...
                             zeros(3,5) ; Z(10:18,1:5)] , 1 , 3)