% Biegeproblem mit 2 Elementen und Linienlast
clear all

la =  1 ; lb = la*0.5 ; EIa = 1 ; EIb = EIa*0.25 ; 
q0 = -1 ;
% FEM-Berechnungsmodell:
xy = [0 ; la ; la+lb] ;
km = [1 2 ; 2 3] ;
ep = [EIa q0 q0 ; EIb 0 0] ;
kr = [1 1 ; 0 0 ; 1 0] ;
bk = [0 0 ; 0 0 ; 0 0] ;

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

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