clear all
clc

%% Zadání:
F_velikost=50; %(N) ... Velikost síly F
A=[4 3 2]; %(m) ... Působiště síly
% Směrové úhly nositelky f:
alfa_f=60;
beta_f=40;
%gamma_f<90°

%% Výpočet:
% 1) Dopočítáme úhel gamma_f
gamma_f=acosd(sqrt(1-cosd(alfa_f)^2-cosd(beta_f)^2))
f0=[cosd(alfa_f) cosd(beta_f) cosd(gamma_f)]; %jednotkový vektor ve směru nositelky f
% 2) Vektor síly F
F=F_velikost*f0
% 3) Moment síly k osám x,y,z
r=[A(1) A(2) A(3)];
M=cross(r,F)
% 4) Moment síly k ose o
%osa prochází počátkem a její směrové úhly jsou:
alfa_o=50;
beta_o=45;
%gamma_o>90°
gamma_o=acosd(-sqrt(1-cosd(alfa_o)^2-cosd(beta_o)^2))
o0=[cosd(alfa_o) cosd(beta_o) cosd(gamma_o)]; %jednotkový vektor ve směru osy o
Mo_velikost=dot(M,o0)
Mo=Mo_velikost*o0