clear all

% FEM-Berechnungsmodell:
[xy km ep kr bk] = finmod_m ('Rahm17.dat') ;
% Systemsteifigkeitsmatrix und -belastungsvektor:
[succ Ksys fq]   = syswbc_m (xy , km , ep , kr , bk) ;

Ks = SymBand2Sparse (Ksys) ;       % "Sparse matrix"
spy (Ks) ;
[m,n]=size(Ks) ;
disp (['Anzahl der Gleichungen: ' , num2str(m)]) ;
xE = Ks \ fq ;                     % Matlab-Standard-Solver

[xPKG iter] = PKG (Ks , fq) ;      % Konjugierte Gradienten
disp (['PKG, ausgeführte Iterationen: ' , num2str(iter)]) ;

x = [xE(m-8:m) xPKG(m-8:m)]
disp ('  El.-Verf.     PKG') ;