clc; clear all;
format short g;
%SILY NA HRIDELI

%Zadane hodnoty
m_hridel=63;
m_kolo=635;
m_dvojkolo=10;
f_Ac=0.3;
f_Bc=0.3;
r_Ac=0.125;
r_Bc=0.0625;
l_1=1.6;
l_2=0.4;
b_1=0.44;
b_2=0.8;
b_3=0.06;
h_1=0.325;
F_Cz=111.64;%Dosadit hodnotu potrebnou pro rozhybani klece ze scriptu "Vodozdviz_sily_v_kleci.m"
g=9.81;

F_ghridel=m_hridel*g;
F_gkolo=m_kolo*g;
F_gdvojkolo=m_dvojkolo*g;


    %M_hnaci   F_Az     F_Ay         F_Bz       F_By   
 A=[ 0,        -1,        0,        1,            0;           %Fz=0
     0,         0,         0,        0,          -(l_1+l_2);     %Mz_A=0
     0,         0,         1,        0,           -1;            %Fy=0
     0,         0,         0,        -(l_1+l_2),   0;            %My_A=0
    -1,         0.4*r_Ac*(f_Ac/sqrt(1+f_Ac^2)),        0.96*r_Ac*(f_Ac/sqrt(1+f_Ac^2)),        0.4*r_Bc*(f_Bc/sqrt(1+f_Bc^2)),        0.96*r_Bc*(f_Bc/sqrt(1+f_Bc^2))     %Mz_A=0 (v z-ose zadna sila nepusobi, proto je vynechana)
    ];              
 
 B=[-F_Cz;
    F_gkolo*b_1/2+F_ghridel*b_2+2*F_gdvojkolo*(l_1+l_2/2);
    F_gkolo+F_ghridel+2*F_gdvojkolo;
    F_Cz*(l_1+b_3);
    -F_Cz*h_1];

 x=A\B 
 
   %Hledana hodnota pro rozhybani systemu je prvni vypsana - tedy M_hnaci