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"
[m,n]=size(Ks) ;
disp (['Anzahl der Gleichungen: ' , num2str(m)]) ;
xE = Ks \ fq ;                % Matlab-Standard-Solver

R = cholinc (Ks , '0') ;      % Incomplete Cholesky
spy (R) ;
% Preconditioned conjugate gradient:
xpcg = pcg (Ks , fq , 1e-6 , m*2 , R' , R) ;

mm = round (m/2) ;
x = [xE(mm-4:mm+4) xpcg(mm-4:mm+4)]
disp ('  El.-Verf.    pcg') ;