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