#	   Solutions A,B,E and T in cylindrical polar coordinates
#
#The commented lines contain solution T only, the uncommented
#include all four solutions.
#
#ur:=simplify(z*diff(chi(r,theta,z),r)/(2*mu));
#ut:=simplify(z*diff(chi(r,theta,z),theta)/(2*mu*r));
#uz:=simplify((chi(r,theta,z)+z*diff(chi(r,theta,z),z))/(2*mu));
#
#srr:=simplify(z*diff(chi(r,theta,z),r,r)-2*diff(chi(r,theta,z),z));
#stt:=simplify(z*diff(chi(r,theta,z),r)/r-z*diff(chi(r,theta,z),theta,theta)/r^2
#-2*diff(chi(r,theta,z),z));
#szz:=simplify(z*diff(chi(r,theta,z),z,z));
#srz:=simplify(z*diff(chi(r,theta,z),r,z)+diff(chi(r,theta,z),r));
#stz:=simplify(z*diff(chi(r,theta,z),theta,z)/r+diff(chi(r,theta,z),theta)/r);
#srt:=simplify(z*diff(chi(r,theta,z),r,theta)/r-z*diff(chi(r,theta,z),theta)/r^2);
#
ur:=simplify((diff(phi(r,theta,z),r)+z*diff(omega(r,theta,z),r)
+2*diff(psi(r,theta,z),theta)/r+z*diff(chi(r,theta,z),r))/(2*mu));
ut:=simplify((diff(phi(r,theta,z),theta)/r+z*diff(omega(r,theta,z),theta)/r
-2*diff(psi(r,theta,z),r)+z*diff(chi(r,theta,z),theta)/r)/(2*mu));
uz:=simplify((diff(phi(r,theta,z),z)+z*diff(omega(r,theta,z),z)
-(3-4*nu)*omega(r,theta,z)+chi(r,theta,z)+z*diff(chi(r,theta,z),z))/(2*mu));
#
srr:=simplify(diff(phi(r,theta,z),r,r)+z*diff(omega(r,theta,z),r,r)
-2*nu*diff(omega(r,theta,z),z)
+2*diff(psi(r,theta,z),r,theta)/r-2*diff(psi(r,theta,z),theta)/r^2
+z*diff(chi(r,theta,z),r,r)-2*diff(chi(r,theta,z),z));
stt:=simplify(diff(phi(r,theta,z),r)/r+diff(phi(r,theta,z),theta,theta)/r^2
+z*diff(omega(r,theta,z),r)/r+z*diff(omega(r,theta,z),theta,theta)/r^2
-2*nu*diff(omega(r,theta,z),z)
-2*diff(psi(r,theta,z),r,theta)/r+2*diff(psi(r,theta,z),theta)/r^2
+z*diff(chi(r,theta,z),r)/r-z*diff(chi(r,theta,z),theta,theta)/r^2
-2*diff(chi(r,theta,z),z));
szz:=simplify(diff(phi(r,theta,z),z,z)+z*diff(omega(r,theta,z),z,z)
-2*(1-nu)*diff(omega(r,theta,z),z)
+z*diff(chi(r,theta,z),z,z));
srz:=simplify(diff(phi(r,theta,z),r,z)+z*diff(omega(r,theta,z),r,z)
-(1-2*nu)*diff(omega(r,theta,z),r)
+diff(psi(r,theta,z),theta,z)/r
+z*diff(chi(r,theta,z),r,z)+diff(chi(r,theta,z),r));
stz:=simplify(diff(phi(r,theta,z),theta,z)/r
+z*diff(omega(r,theta,z),theta,z)/r
-(1-2*nu)*diff(omega(r,theta,z),theta)/r-diff(psi(r,theta,z),r,z)
+z*diff(chi(r,theta,z),theta,z)/r+diff(chi(r,theta,z),theta)/r);
srt:=simplify(diff(phi(r,theta,z),r,theta)/r-diff(phi(r,theta,z),theta)/r^2
+z*diff(omega(r,theta,z),r,theta)/r
-z*diff(omega(r,theta,z),theta)/r^2+diff(psi(r,theta,z),r)/r
-diff(psi(r,theta,z),r,r)+diff(psi(r,theta,z),theta,theta)/r^2
+z*diff(chi(r,theta,z),r,theta)/r-z*diff(chi(r,theta,z),theta)/r^2);
#
T:=simplify(((1-nu)/(mu*alpha*(1+nu)))*diff(chi(r,theta,z),z));
qz:=-simplify(((1-nu)*K/(mu*alpha*(1+nu)))*diff(chi(r,theta,z),z,z));
qr:=-simplify(((1-nu)*K/(mu*alpha*(1+nu)))*diff(chi(r,theta,z),z,r));