syms v1x v2x v3x v1y v2y v3y v1z v2z v3z real
v1 = sqrt(v1x^2 + v1y^2 + v1z^2);
v2 = sqrt(v2x^2 + v2y^2 + v2z^2);
v3 = sqrt(v3x^2 + v3y^2 + v3z^2);
alfa1 = acos(v1x/v1);
beta1 = acos(v1y/v1);
gama1 = acos(v1z/v1);
alfa2 = acos(v2x/v2);
beta2 = acos(v2y/v2);
gama2 = acos(v2z/v2);
alfa3 = acos(v3x/v3);
beta3 = acos(v3y/v3);
gama3 = acos(v3z/v3);
v1v = [v1x v1y v1z].';
v2v = [v2x v2y v2z].';
v3v = [v3x v3y v3z].';
E1 = (v1v + v2v + v3v).^2;
E2 = v1v*dot(v2v, v3v);
display(v1)
display(v2)
display(v3)
display(alfa1)
display(beta1)
display(gama1)
display(alfa2)
display(beta2)
display(gama2)
display(alfa3)
display(beta3)
display(gama3)
display(E1)
display(E2)
v1xn = 1;
v2xn = 0;
v3xn = 0;
v1yn = 0;
v2yn = 1;
v3yn = 0;
v1zn = 0;
v2zn = 0;
v3zn = 1;
v1n = eval(subs(v1, [v1x v1y v1z], [v1xn v1yn v1zn]));
display(v1n)
v2n = eval(subs(v2, [v2x v2y v2z], [v2xn v2yn v2zn]));
display(v2n)
v3n = eval(subs(v3, [v3x v3y v3z], [v3xn v3yn v3zn]));
display(v3n)
alfa1n = rad2deg(eval(subs(alfa1, [v1x v1y v1z], [v1xn v1yn v1zn])));
display(alfa1n)
beta1n = rad2deg(eval(subs(beta1, [v1x v1y v1z], [v1xn v1yn v1zn])));
display(beta1n)
gama1n = rad2deg(eval(subs(gama1, [v1x v1y v1z], [v1xn v1yn v1zn])));
display(gama1n)
alfa2n = rad2deg(eval(subs(alfa2, [v2x v2y v2z], [v2xn v2yn v2zn])));
display(alfa2n)
beta2n = rad2deg(eval(subs(beta2, [v2x v2y v2z], [v2xn v2yn v2zn])));
display(beta2n)
gama2n = rad2deg(eval(subs(gama2, [v2x v2y v2z], [v2xn v2yn v2zn])));
display(gama2n)
alfa3n = rad2deg(eval(subs(alfa3, [v3x v3y v3z], [v3xn v3yn v3zn])));
display(alfa3n)
beta3n = rad2deg(eval(subs(beta3, [v3x v3y v3z], [v3xn v3yn v3zn])));
display(beta3n)
gama3n = rad2deg(eval(subs(gama3, [v3x v3y v3z], [v3xn v3yn v3zn])));
display(gama3n)
E1n = eval(subs(E1, [v1x v1y v1z v2x v2y v2z v3x v3y v3z], [v1xn v1yn v1zn v2xn v2yn v2zn v3xn v3yn v3zn]));
display(E1n)
E2n = eval(subs(E2, [v1x v1y v1z v2x v2y v2z v3x v3y v3z], [v1xn v1yn v1zn v2xn v2yn v2zn v3xn v3yn v3zn]));
display(E2n)