#Solutions A,B and E
#
ux:=simplify((diff(phi(x,y,z),x)+z*diff(omega(x,y,z),x)
+2*diff(psi(x,y,z),y))/(2*mu));
uy:=simplify((diff(phi(x,y,z),y)+z*diff(omega(x,y,z),y)
-2*diff(psi(x,y,z),x))/(2*mu));
uz:=simplify((diff(phi(x,y,z),z)+z*diff(omega(x,y,z),z)
-(3-4*nu)*omega(x,y,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));
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));
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));
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));
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));
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));