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);