#	   Solutions A,B,E and T in Cartesian coordinates
#
#The commented lines contain solution T only, the uncommented
#include all four solutions.
#
#ux:=simplify(z*diff(chi(x,y,z),x)/(2*mu));
#uy:=simplify(z*diff(chi(x,y,z),y)/(2*mu));
#uz:=simplify((chi(x,y,z)+z*diff(chi(x,y,z),z))/(2*mu));
#
#sxx:=simplify(z*diff(chi(x,y,z),x,x)-2*diff(chi(x,y,z),z));
#sxy:=simplify(z*diff(chi(x,y,z),x,y));
#syy:=simplify(z*diff(chi(x,y,z),y,y)-2*diff(chi(x,y,z),z));
#sxz:=simplify(z*diff(chi(x,y,z),x,z)+diff(chi(x,y,z),x));
#syz:=simplify(z*diff(chi(x,y,z),y,z)+diff(chi(x,y,z),y));
#szz:=simplify(z*diff(chi(x,y,z),z,z));
#
T:=simplify(((1-nu)/(mu*alpha*(1+nu)))*diff(chi(x,y,z),z));
qx:=-simplify(((1-nu)*K/(mu*alpha*(1+nu)))*diff(chi(x,y,z),z,x));
qy:=-simplify(((1-nu)*K/(mu*alpha*(1+nu)))*diff(chi(x,y,z),z,y));
qz:=-simplify(((1-nu)*K/(mu*alpha*(1+nu)))*diff(chi(x,y,z),z,z));
#
ux:=simplify((diff(phi(x,y,z),x)+z*diff(omega(x,y,z),x)
+2*diff(psi(x,y,z),y)+z*diff(chi(x,y,z),x))/(2*mu));
uy:=simplify((diff(phi(x,y,z),y)+z*diff(omega(x,y,z),y)
-2*diff(psi(x,y,z),x)+z*diff(chi(x,y,z),y))/(2*mu));
uz:=simplify((diff(phi(x,y,z),z)+z*diff(omega(x,y,z),z)
-(3-4*nu)*omega(x,y,z)+chi(x,y,z)+z*diff(chi(x,y,z),z))/(2*mu));
#
sxx:=simplify(diff(phi(x,y,z),x,x)+z*diff(omega(x,y,z),x,x)
-2*nu*diff(omega(x,y,z),z)+2*diff(psi(x,y,z),x,y)+z*diff(chi(x,y,z),x,x)-2*diff(chi(x,y,z),z));
sxy:=simplify(diff(phi(x,y,z),x,y)
+z*diff(omega(x,y,z),x,y)+diff(psi(x,y,z),y,y)-diff(psi(x,y,z),x,x)+z*diff(chi(x,y,z),x,y));
szz:=simplify(diff(phi(x,y,z),z,z)+z*diff(omega(x,y,z),z,z)
-2*(1-nu)*diff(omega(x,y,z),z)+z*diff(chi(x,y,z),z,z));
syy:=simplify(diff(phi(x,y,z),y,y)+z*diff(omega(x,y,z),y,y)
-2*nu*diff(omega(x,y,z),z)-2*diff(psi(x,y,z),x,y)+z*diff(chi(x,y,z),y,y)-2*diff(chi(x,y,z),z));
sxz:=simplify(diff(phi(x,y,z),x,z)+z*diff(omega(x,y,z),x,z)
-(1-2*nu)*diff(omega(x,y,z),x)+diff(psi(x,y,z),y,z)+z*diff(chi(x,y,z),x,z)+diff(chi(x,y,z),x));
syz:=simplify(diff(phi(x,y,z),y,z)
+z*diff(omega(x,y,z),y,z)-(1-2*nu)*diff(omega(x,y,z),y)
-diff(psi(x,y,z),x,z)+z*diff(chi(x,y,z),y,z)+diff(chi(x,y,z),y));