% Aufgabe 18-14
clear all

% FEM-Berechnungsmodell:
xy = [0 ; 1 ; 2 ; 3] ;
km = [1 2 ; 2 3 ; 3 4] ;
ep = [1 0 0 ; 2 0 0 ; 1 -2 -1] ;
kr = [1 1 ; 0 0 ; 0 0 ; 1 0] ;
bk = [0 0 ; -1 0 ; 0 0 ; 0 0] ;

% Elementsteifigkeitsbeziehungen:
[Ke1 fe1] = elemat_m (1 , 2 , xy , km , ep)
[Ke2 fe2] = elemat_m (2 , 2 , xy , km , ep)
[Ke3 fe3] = elemat_m (3 , 2 , xy , km , ep)

[succ Ksys f] = system_m (xy , km , ep , kr , bk) ;
K = SymBand2Quad (Ksys) , f     % Systemsteifigkeitsbeziehung

[succ Ksys fq] = syswbc_m (xy , km , ep , kr , bk) ;
Kq = SymBand2Quad (Ksys) , fq   % Systemsteifigkeitsbeziehung
               % (verhinderte Verschiebungen beruecksichtigt)

% Verschiebungsberechnung:
[succ vphi] = femalg_m (xy , km , ep , kr , bk)