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) 
xE = Ks \ fq ;                     % Matlab-Standard-Solver

xKG = KonjugGrad     (Ks , fq) ;   % Konjugierte Gradienten

xbicg = bicg         (Ks , fq) ;   % Biconjugate gradient
xbicgstab = bicgstab (Ks , fq) ;   % Biconjugate gradient stabilized
xcgs = cgs           (Ks , fq) ;   % Conjugate gradient squared
xgmres = gmres       (Ks , fq) ;   % Generalized minimum residual
xlsqr = lsqr         (Ks , fq) ;   % LSQR (Conjugate Gradients on the Normal Equations)
xminres = minres     (Ks , fq) ;   % Minimum residual
xpcg = pcg           (Ks , fq) ;   % Preconditioned conjugate gradient
xqmr = qmr           (Ks , fq) ;   % Quasiminimal residual
xsymmlq = symmlq     (Ks , fq) ;   % Symmetric LQ

x = [xE(m-8:m) xKG(m-8:m) xbicg(m-8:m) xbicgstab(m-8:m) xcgs(m-8:m) xgmres(m-8:m) ...
     xlsqr(m-8:m) xminres(m-8:m) xpcg(m-8:m) xqmr(m-8:m) xsymmlq(m-8:m)]
disp ('  El.-Verf. KonjugGrad   bicg    bicgstab     cgs      gmres      lsqr    minres      pcg       qmr     symmlq') ;