Xs := -1/9*3^(1/2)*b*(2*2^(1/2)*cos(beta)*sin(beta)*x1^2-2*2^(1/2)*cos(beta)*
sin(beta)*x1*x2-2*2^(1/2)*cos(beta)*sin(beta)*x1*x3-2^(1/2)*cos(beta)*sin(beta
)*x2^2+4*2^(1/2)*cos(beta)*sin(beta)*x2*x3-2^(1/2)*cos(beta)*sin(beta)*x3^2+8*
cos(beta)^2*x1^2+10*cos(beta)^2*x1*x2+10*cos(beta)^2*x1*x3-4*cos(beta)^2*x2^2-\
2*cos(beta)^2*x2*x3-4*cos(beta)^2*x3^2-6*x0*x2+6*x0*x3-4*x1^2-8*x1*x2-8*x1*x3+
2*x2^2-2*x2*x3+2*x3^2)/cos(beta)/sin(beta);
Ys := 1/9*3^(1/2)*b*(2^(1/2)*cos(beta)*sin(beta)*x1^2+2*2^(1/2)*cos(beta)*sin(
beta)*x1*x2-4*2^(1/2)*cos(beta)*sin(beta)*x1*x3-2*2^(1/2)*cos(beta)*sin(beta)*
x2^2+2*2^(1/2)*cos(beta)*sin(beta)*x2*x3+2^(1/2)*cos(beta)*sin(beta)*x3^2+4*
cos(beta)^2*x1^2-10*cos(beta)^2*x1*x2+2*cos(beta)^2*x1*x3-8*cos(beta)^2*x2^2-\
10*cos(beta)^2*x2*x3+4*cos(beta)^2*x3^2-6*x0*x1+6*x0*x3-2*x1^2+8*x1*x2+2*x1*x3
+4*x2^2+8*x2*x3-2*x3^2)/cos(beta)/sin(beta);
Zs := 1/9*3^(1/2)*b*(2^(1/2)*cos(beta)*sin(beta)*x1^2-4*2^(1/2)*cos(beta)*sin(
beta)*x1*x2+2*2^(1/2)*cos(beta)*sin(beta)*x1*x3+2^(1/2)*cos(beta)*sin(beta)*x2
^2+2*2^(1/2)*cos(beta)*sin(beta)*x2*x3-2*2^(1/2)*cos(beta)*sin(beta)*x3^2+4*
cos(beta)^2*x1^2+2*cos(beta)^2*x1*x2-10*cos(beta)^2*x1*x3+4*cos(beta)^2*x2^2-\
10*cos(beta)^2*x2*x3-8*cos(beta)^2*x3^2+6*x0*x1-6*x0*x2-2*x1^2+2*x1*x2+8*x1*x3
-2*x2^2+8*x2*x3+4*x3^2)/cos(beta)/sin(beta);