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

Ks = SymBand2Sparse (Ksys) ;          % Sparse matrix

v  = KonjugGrad (Ks , f) ;
LinkeObereEckeVonKs = full(Ks (1:6,1:6))
disp (['KondZahl von Ks:  ' num2str(condest(Ks))]) ;

[Ksq , fq , D] = SkalSymmGlsSparse (Ks , f) ;
disp (['KondZahl von Ksq: ' num2str(condest(Ksq))]) ;
vq = KonjugGrad (Ksq , fq) ;
v  = D*vq ;
v25_33 = v(25:33)