clc %czyszczenie ekranu
clear all %usunięcie wszystkich zmiennych z pamięci

hy = 0.005;
hz = 0.05;

Iz = hz*hy^3 / 12;
E = 2.1e11; %2.1*10^11
l =  2;

K1 = MacierzSztywnosciBelka(E, Iz, l/2);
K2 = MacierzSztywnosciBelka(E, Iz, l/2);

%globalna macierz sztywności

K = zeros(6,6);

K = [K1(1,1) K1(1,2) K1(1,3) K1(1,4) 0 0;
     K1(2,1) K1(2,2) K1(2,3) K1(2,4) 0 0;
     K1(3,1) K1(3,2) K1(3,3)+K2(1,1) K1(3,4)+K2(1,2) K2(1,3) K2(1,4);
     K1(4,1) K1(4,2) K1(4,3)+K2(2,1) K1(4,4)+K2(2,2) K2(2,3) K2(2,4);
           0       0 K2(3,1) K2(3,2) K2(3,4) K2(3,4);
           0       0 K2(4,1) K2(4,2) K2(4,3) K2(4,4)]