Matlab schema logica !

syms F1x F2x F3x F1y F2y F3y F1z F2z F3z real
F1 = sqrt(F1x^2 + F1y^2 + F1z^2);
F2 = sqrt(F2x^2 + F2y^2 + F2z^2);
F3 = sqrt(F3x^2 + F3y^2 + F3z^2);
alfa1 = acos(F1x/F1);
beta1 = acos(F1y/F1);
gama1 = acos(F1z/F1);
alfa2 = acos(F2x/F2);
beta2 = acos(F2y/F2);
gama2 = acos(F2z/F2);
alfa3 = acos(F3x/F3);
beta3 = acos(F3y/F3);
gama3 = acos(F3z/F3);
F1F = [F1x F1y F1z].';
F2F = [F2x F2y F2z].';
F3F = [F3x F3y F3z].';
E1 = -F1F -F2F - F3F;
E2 = cross(cross(F1F, F2F),F3F);
display(F1)
display(F2)
display(F3)
display(alfa1)
display(beta1)
display(gama1)
display(alfa2)
display(beta2)
display(gama2)
display(alfa3)
display(beta3)
display(gama3)
display(E1)
display(E2)
F1xn = 0;
F1yn = 1;
F1zn = 2;
F2xn = 3;
F2yn = 10;
F2zn = 5;
F3xn = -23;
F3yn = 23;
F3zn = 8;
F1n = eval(subs(F1, [F1x F1y F1z], [F1xn F1yn F1zn]));
display(F1n)
F2n = eval(subs(F2, [F2x F2y F2z], [F2xn F2yn F2zn]));
display(F2n)
F3n = eval(subs(F3, [F3x F3y F3z], [F3xn F3yn F3zn]));
display(F3n)
alfa1n = rad2deg(eval(subs(alfa1, [F1x F1y F1z], [F1xn F1yn F1zn])));
display(alfa1n)
beta1n = rad2deg(eval(subs(beta1, [F1x F1y F1z], [F1xn F1yn F1zn])));
display(beta1n)
gama1n = rad2deg(eval(subs(gama1, [F1x F1y F1z], [F1xn F1yn F1zn])));
display(gama1n)
alfa2n = rad2deg(eval(subs(alfa2, [F2x F2y F2z], [F2xn F2yn F2zn])));
display(alfa2n)
beta2n = rad2deg(eval(subs(beta2, [F2x F2y F2z], [F2xn F2yn F2zn])));
display(beta2n)
gama2n = rad2deg(eval(subs(gama2, [F2x F2y F2z], [F2xn F2yn F2zn])));
display(gama2n)
alfa3n = rad2deg(eval(subs(alfa3, [F3x F3y F3z], [F3xn F3yn F3zn])));
display(alfa3n)
beta3n = rad2deg(eval(subs(beta3, [F3x F3y F3z], [F3xn F3yn F3zn])));
display(beta3n)
gama3n = rad2deg(eval(subs(gama3, [F3x F3y F3z], [F3xn F3yn F3zn])));
display(gama3n)
E1n = eval(subs(E1, [F1x F1y F1z F2x F2y F2z F3x F3y F3z],...  
	[F1xn F1yn F1zn F2xn F2yn F2zn F3xn F3yn F3zn]));
display(E1n)
E2n = eval(subs(E2, [F1x F1y F1z F2x F2y F2z F3x F3y F3z],...	
	[F1xn F1yn F1zn F2xn F2yn F2zn F3xn F3yn F3zn]));
display(E2n)

Edited by MarianG, 20 May 2019 - 19:50.
CODE