% Lineares Gleichungssystem S. 296

clear all

ldL = 4 ;  % "Klein-l durch Gross-L"

A = [        1                 0                  1                   0         ;
             1                 1                 -1                   1         ;
     exp(ldL)*cos(ldL) exp(ldL)*sin(ldL) exp(-ldL)*cos(ldL)  exp(-ldL)*sin(ldL) ;
    -exp(ldL)*sin(ldL) exp(ldL)*cos(ldL) exp(-ldL)*sin(ldL) -exp(-ldL)*cos(ldL) ] ;

b = [-1 ; 0 ; -1 ; 0 ] ;   

% x = [C1 ; C2 ; C3 ; C4 ] * k / q0

format long ;
x = A \ b        % Ergebnis ---> Command window

C1 = x(1) ;
C2 = x(2) ;
C3 = x(3) ;
C4 = x(4) ;

for i=1:101
  zdL(i) = 0.01 * (i-1) ;
  arg    = 4*zdL(i);
  v(i)   = (1 +  exp(arg) * ( C1*cos(arg)+C2*sin(arg)) + exp(-arg) * (C3*cos(arg)+C4*sin(arg)))/1024 ;
  Mb(i)  =    - (exp(arg) * (-C1*sin(arg)+C2*cos(arg)) + exp(-arg) * (C3*sin(arg)-C4*cos(arg)))/32   ;
end
   
subplot (2,1,1) ; plot (zdL , v)  , grid on , title ('Biegelinie')
subplot (2,1,2) ; plot (zdL , Mb) , grid on , title ('Biegemoment')