% Berechnungsmodell (Aufgabe 19-10):
clear all
xy = [0 ; 100 ; 200 ; 300 ; 400 ; 500 ; 1200] ;
km = [1 2 ; 2 3 ; 3 4 ; 4 5 ; 5 6 ; 6 7] ;           
ep = [6e9  0   0 12  ;   % 1. Spalte: Biegesteifigkeit,
      6e9  0   0 12  ;   % 2. und 3. Spalte: Linienlastintensitaet
      6e9  0   0 12  ;   % ... an den Knoten des jeweiligen Elements,
      6e9  0   0 12  ;   % 4. Spalte: Bettungsziffer 
      6e9  0   0 12  ;
      8e10 0 -10  0] ;        
kr = [1 1 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 1   0] ;  % Randbedingungen  
bk = [0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 8e5] ;  % Kraefte/Momente 
sc = [0 0 ; 0 0 ; 0 0 ; 0 0 ; 0 0 ; 800 0 ; 0 0] ;  % Federn

% Finite-Elemente-Algorithmus:
[succ uv] = femalg_m (xy , km , ep , kr , bk , sc) 