% Moment síly v prostoru
A = [3,2,4];
F_vel = 40;
alfaf = 80*pi/180;
betaf = 60*pi/180;
alfao = 60*pi/180;
betao = 120*pi/180;

% Řešení4
gamaf = acos(sqrt(1 - (cos(alfaf))^2- (cos(betaf))^2));
Fx = F_vel*cos(alfaf);
Fy = F_vel*cos(betaf);
Fz = F_vel*cos(gamaf);
F = [Fx Fy Fz];
% Moment k počátku
M0 = cross(A,F);
% Moment k ose
gamao = acos(-sqrt(1 - (cos(alfao))^2- (cos(betao))^2));
oo = [cos(alfao) cos(betao) cos(gamao)];
Mo_vel = dot(M0,oo);
Mo = Mo_vel*oo;
% Jednotkový vektor momentu M0
M0_vel = sqrt(M0(1)^2 + M0(2)^2 + M0(3)^2);
m0 = M0/M0_vel;