clear all      % Finite-Elemente-Modell eines Kragtraegers ...
L  = 1000 ;    % ... der Laenge L ...
EI = 2e10 ;    % ... mit konstanter Biegesteifigkeit EI ...
q  = 2 ;       % ... unter konstamter Linienlast q, ...
n  = 10 ;      % ... eingeteilt in n Elemente

xy = zeros (n+1,1) ;       % Koordinaten der n+1 Knoten
for i=1:n+1
    xy(i,1) = (i-1)*L/n ;
end
km = zeros (n,2) ; ep = zeros (n,3) ;
for i=1:n
    km(i,:) = [i  i+1] ;   % Koinzidenzmatrix
    ep(i,:) = [EI q q] ;   % Elementparameter
end
kr = zeros (n+1,2) ;
kr(1,:) = [1 1] ;          % Knoten 1: Starre Einspannung
bk = zeros (n+1,2) ;       % Keine zusaetzlichen Knotenlasten

% Aufbau des linearen Gleichungssystems  Ksys * v = fq:
[succ Ksys fq]   = syswbc_m (xy , km , ep , kr , bk) ;
% Ksys ist eine symmetrische Bandmatrix, und ...

Ks = SymBand2Sparse (Ksys) ;        % Kss ist "Sparse matrix"
spy(Ks)                             % Muster zeichnen
GleichungsAnzahl = 2*n+2
Kondzahl = condest (Ks)             % Konditionszahl
v  = Ks \ fq ;                      % Matlab-Standard-Solver
format long
vEnd = v(2*n+1)                     % Maximale Durchbiegung
[Ksn fqn] = SkalGlsSparse (Ks,fq) ; % Zeilenskalierung
KondzahlSkal = condest (Ksn)        % Konditionszahl
v  = Ksn \ fqn ;                    % Matlab-Standard-Solver
vEndSkal = v(2*n+1)                 % Maximale Durchbiegung