% Berechnungsmodell (Biegeproblem mit Federn und Gelenk):
clear all
xy = [0 ; 200 ; 900 ; 900 ; 1090 ; 1330  ; 1590 ; 1870] ;
km = [1 2 ; 2 3 ; 4 5 ; 5 6 ; 6 7 ; 7 8] ;
ep = [6.0e9   0 0 ;   % 1. Spalte: Biegesteifigkeit
      1.2e10 -3 0 ;   % 2. und 3. Spalte: Linienlastintensität
      1.2e10  0 0 ;   % ... an den Knoten des jeweiligen Elements
      1.2e10  0 0 ;  
      1.0e9   0 0 ; 
      1.0e9   0 0] ; 
kr = [1   0    ; 0 0 ; 0 0 ; 0 0 ;  0  0 ;   0   0 ; 1 0 ;  0 0] ;
bk = [0 200000 ; 0 0 ; 0 0 ; 0 0 ;  0  0 ; -2000 0 ; 0 0 ;  0 0] ; 
sc = [0   0    ; 0 0 ; 0 0 ; 0 0 ; 100 0 ;   0   0 ; 0 0 ; 50 0] ;
li = [0   0    ; 0 0 ; 1 0 ; 1 0 ;  0  0 ;   0   0 ; 0 0 ;  0 0] ; 
% Knoten 1     ;  2  ;  3  ;  4  ;    5  ;     6   ;  7  ;   8

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