
function k = kloub_vetknuti(E,I,L,A) %110_111

E = input ('enter the Young modulus'); %Modul pružnosti v tahu a v tlaku E [Pa]
I = input ('enter the moment inertia'); %% Moment setrvačnosti obdelnikového průřezu [m^4]
L = input ('enter the length of the beam'); %Délka prutu [m]
A = input ('enter the area of the beam'); %Plocha průřezu [m^2]

%Matice tuhosti prutu typ 111_101

k = [ E*A/L, 0, 0, -E*A/L, 0, 0;                             
      0, 3*E*I/L^3, 0, 0, -3*E*I/L^3, -3*E*I/L^2;
      0, 0, 0, 0, 0, 0;
      -E*A/L, 0, 0, E*A/L, 0, 0;
      0, -3*E*I/L^3, 0, 0, 3*E*I/L^3, 3*E*I/L^2;
      0, -3*E*I/L^2, 0, 0, 3*E*I/L^2, 3*E*I/L]; 

end