% Berechnungsmodell (Aufgabe 19-3):
clear all
xy = [0 ; 0.25 ; 1 ; 1.75 ; 2] ;         % Koordinaten
km = [1 2 ; 2 3 ; 3 4 ; 4 5] ;           % Elemente/Knoten
ep = [1 0 0 2 ;   % Elementparameter, 1. Spalte: Biegesteifigkeit,
      1 0 0 2 ;   % 2. und 3. Spalte: Linienlastintensität
      1 0 0 2 ;   % ... an den Knoten des jeweiligen Elements,
      1 0 0 2 ] ; % 4. Spalte: Bettungsziffer
kr = [0 0 ;  0 0 ; 0 0 ;  0 0 ; 0 0] ;   % Verhinderte Verformungen
bk = [0 0 ; -1 0 ; 0 0 ; -1 0 ; 0 0] ;   % Einzellasten
 
% Finite-Elemente-Algorithmus:
[succ uv] = femalg_m (xy , km , ep , kr , bk) 