
 a := 1;
 evalf(a);
 b := 2;
 c := 3;
 E := a*z^2+b*z+c;
 niy := .4;
 G := E/(2*(1+niy));
 Alfa11 := 1/E; Alfa12 := -(1*niy)/E; Alfa13 := -(1*niy)/E;
 Alfa22 := 1/E; Alfa23 := -(1*niy)/E; Alfa33 := 1/E;
 Alfa44 := 1/G; Alfa55 := 1/G; Alfa66 := 1/G;
 Liambda := 2*niy*G/(1-2*niy);
 Liambda11 := 2*G+Liambda;
 Liambda12 := Liambda; Liambda13 := Liambda;
 Liambda22 := 2*G+Liambda;
 Liambda23 := Liambda; Liambda33 := 2*G+Liambda;
 Liambda44 := G; Liambda55 := G; Liambda66 := G;
 A := Matrix(3, 3, [[Alfa11, Alfa12, Alfa13], [Alfa12, Alfa22, Alfa23], [Alfa13, Alfa23, Alfa33]]);
 B := evalm(1/A);
 Xalfa := -1/G;
 Yalfa := -1/G;
 Zalfa := -1/(2*G+Liambda);
 a1 := -(2*G+Liambda)/G;
 b1 := -1;
 c1 := -(G+Liambda)/G;
 d1 := -1-Liambda/G;
 a2 := -1;
 b2 := -(2*G+Liambda)/G;
 c2 := -(G+Liambda)/G;
 d2 := -1-Liambda/G;
 a3 := -(G+Liambda)/(2*G+Liambda);
 b3 := -(G+Liambda)/(2*G+Liambda);
 c3 := -G/(2*G+Liambda);
 d3 := G/(2*G+Liambda);
 


