% Allgemeines Matrizeneigenwertproblem auf Seite 400

clear all

nA = 8 ;
n  = nA - 1 ;

A = zeros(n,n) ;
B = zeros(n,n) ;

for i=3:n-2
    A(i,i-2:i+2) = [1 -4 6 -4 1] ;
    B(i,i-1:i+1) = [-i+0.5 i*2 -i-0.5] ;
end

A(1,1:3) = [5 -4 1]    ;
A(2,1:4) = [-4 6 -4 1] ;
A(n-1,n-3:n) = [1 -4 6 -4] ;
A(n,n-2:n)   = [1 -4 5]    ;

B(1,1:2) = [2 -1.5]    ;
B(2,1:3) = [-1.5 4 -2.5] ;
B(n-1,n-2:n) = [-(n-1)+0.5 (n-1)*2 -(n-1)-0.5] ;
B(n,n-1:n)   = [-n+0.5 n*2] ;

d = eig (A,B) ;             % ... berechnet alle Eigenwerte
Eigenwerte = d.^(1/3) * nA  % ... berechnet die kritischen Längen
                            % ("Punkt-Operation" .^(1/3) berechnet 3. Wurzel elementweise)

