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

Kf = SymBand2Quad (Ksys) ;          % "Voll besetzte Matrix"
t1 = cputime ;
xV = Kf \ fq ;
disp (['Zeit fuer volle Matrix:  ' , num2str(cputime-t1)]) ;

Ks = SymBand2Sparse (Ksys) ;        % "Sparse matrix"
spy(Ks)                             % Muster zeichnen
t1 = cputime ;
xS = Ks \ fq ;
disp (['Zeit fuer Sparse Matrix:  ' , num2str(cputime-t1)]) ;

t1 = cputime ;
xI = PKG (Ks , fq) ; 
disp (['Zeit fuer PKG-Verfahren: ' , num2str(cputime-t1)]) ;

x = [xV(end-7:end) xS(end-7:end) xI(end-7:end)]
disp ('    Voll     Sparse   Iteration') ;