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

Kf = SymBand2Quad   (Ksys) ;       % "Voll besetzte Matrix"

t1 = cputime ;
x  = Kf \ fq ;
disp (['Zeit für Standard-Solver (Full):   ' num2str(cputime-t1)]) ;

t1 = cputime ;
opts.SYM = true ; opts.POSDEF = true ;
x  = linsolve (Kf , fq , opts) ;
disp (['Zeit für linsolve (opts):          ' num2str(cputime-t1)]) ;

Kf(1,1) = -Kf(1,1) ;   % ... und damit ist es vorbei mit der 
                       %     positiven Definitheit
t1 = cputime ;
x  = Kf \ fq ;
disp (['Zeit für Standard-Solver (Full):   ' num2str(cputime-t1)]) ;

t1 = cputime ;
opts.SYM = true ; opts.POSDEF = true ;
x  = linsolve (Kf , fq , opts) ;