syms F1 F2 F3 alfa1 alfa2 alfa3 F1x F1z F2x F2z F3x F3z x1 z1 x2 z2 x3 z3
F1x = F1 * cos(alfa1);
F1z = F1 * sin(alfa1);
F2x = F2 * cos(alfa2);
F2z = F2 * sin(alfa2);
F3x = F3 * cos(alfa3);
F3z = F3 * sin(alfa3);
Rx = F1x + F2x + F3x;
Rz = F1z + F2z + F3z;
R = [Rx sym(0) Rz];
My = x1*F1z - z1*F1x + x2*F2z - z2*F2x + x3*F3z - z3*F3x;
Mo = [sym(0) My sym(0)];
disp(R)
disp(Mo)
F1n=input('F1= ');
F2n=input('F2= ');
F3n=input('F3= ');
x1n=input('x1= ');
x2n=input('x2= ');
x3n=input('x3= ');
z1n=input('z1= ');
z2n=input('z2= ');
z3n=input('z3= ');
alfa1n=input('alfa1= ');
alfa2n=input('alfa2= ');
alfa3n=input('alfa3= ');
alfa1n=deg2rad(alfa1n);
alfa2n=deg2rad(alfa2n);
alfa3n=deg2rad(alfa3n);
R_n = eval(subs(R, [F1 F2 F3 alfa1 alfa2 alfa3], [F1n F2n F3n alfa1n alfa2n alfa3n]))
Mo_n = eval(subs(Mo, [F1 F2 F3 alfa1 alfa2 alfa3 x1 x2 x3 z1 z2 z3], [F1n F2n F3n alfa1n alfa2n alfa3n x1n x2n x3n z1n z2n z3n]))