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

Ks = SymBand2Sparse (Ksys) ;         % "Sparse matrix"
subplot (2,1,1) ; spy(Ks) ;          % Muster zeichnen
R = chol (Ks) ;
subplot (2,1,2) ; spy(R) ;           % Muster zeichnen