clear all
% Allgemeines symmetrisches Matrizen-
% eigenwertproblem (A-lambda*B)*x=0 ...
A = [  24  0 -12   6   0   0   0 ;
        0  8  -6   2   0   0   0 ;
      -12 -6  24   0 -12   6   0 ;
        6  2   0   8  -6   2   0 ;
        0  0 -12  -6  24   0   6 ;
        0  0   6   2   0   8   2 ;
        0  0   0   0   6   2   4 ] ;
B = [ 312  0  54 -13   0   0   0 ;
        0  8  13  -3   0   0   0 ;
       54 13 312   0  54 -13   0 ;
      -13 -3   0   8  13  -3   0 ;
        0  0  54  13 312   0 -13 ;
        0  0 -13  -3   0   8  -3 ;
        0  0   0   0 -13  -3   4 ] ;

% ... wird reduziert auf spezielles    
% Matrizeneigenwertproblem (Aq-lambda*E)*x=0
% mit unsymmetrischer Matrix Aq:
Aq = inv(B) * A

% ... wird reduziert auf spezielles    
% Matrizeneigenwertproblem (As-lambda*E)*x=0
% mit symmetrischer Matrix As:
R  = inv(chol(B)) ;
As = R'*A*R

lamAllg  = eig(A,B) ;
lamSpUns = eig(Aq)  ;
lamSpSym = eig(As)  ;
Eigenwerte = [lamAllg lamSpUns lamSpSym]