% FEMSET-Berechnung: Stabdreischlag
clear all

xy = [0 240 ; 0 0 ; 0 -320 ; 450  0  ] ;   % Kn.-Koord.
FK = [0  0  ; 0 0 ; 0   0  ;  0 -3000] ;   % Kn.-Lasten
Lg = [1  1  ; 1 1 ; 1   1  ;  0   0  ] ;   % u-verhindert
%
KM = [ 1  4 ;  2    4    ;  3  4] ; 
EP = [5e6 0 ; 8e6 240e-5 ; 2e6 0] ; 
  
% Finite-Elemente-Algorithmus:
[succ uv] = femalg_m (xy , KM , EP , Lg , FK) 

if succ == 1                                           % Nach erfolgreicher Rechnung: 
    drawfw2duv (xy , KM , Lg , FK , uv) ;              % Verformtes System zeichnen
    sk = Stabkraefte2D (xy , KM , EP , uv) ;           % Berechnung der Stabkraefte ...
    for i = 1:3
       disp (['Stab ',num2str(i), ': Stabkraft F = ',num2str(sk(i)),' N']) ;
    end
    figure ;                                           % ... und zeichnen des Systems ..
    drawfw2dsk (xy , KM , Lg , FK , sk , 'e') ;        % ... in ein neues Graphik-Fenster
end