(*  SOLUTIONS A B E AND T IN SPHERICAL POLAR COORDINATES  *)

uR=Simplify[(D[ϕ,R]+R*Cos[β]*D[ω,R]-(3-4*ν)*Cos[β]*ω+(2/R)*D[ψ,θ]+R*Cos[β]*D[χ,R]+χ*Cos[β])/(2*μ)]
uθ=Simplify[((1/R*Sin[β])*D[ϕ,θ]+Cot[β]*D[ω,θ]-2*Sin[β]*D[ψ,R]-(2*Cos[β]/R)*D[ψ,β]+Cot[β]*D[χ,θ])/(2*μ)]
uβ=Simplify[((1/R)*D[ϕ,β]+Cos[β]*D[ω,β]+(3-4*ν)*Sin[β]*ω+(2*Cot[β]/R)*D[ψ,θ]+Cos[β]*D[χ,β]-χ*Sin[β])/(2*μ)]

σRR=Simplify[D[ϕ,R,R]+R*Cos[β]*D[ω,R,R]-2*(1-ν)*Cos[β]*D[ω,R]+(2*ν/R)*Sin[β]*D[ω,β]+(2/R)*D[ψ,R,θ]-(2/R^2)*D[ψ,θ]+R*Cos[β]*D[χ,R,R]+2*Sin[β]*D[χ,β]/R]
σRθ=Simplify[(1/(R*Sin[β]))*D[ϕ,R,θ]-(1/(R^2*Sin[β]))*D[ϕ,θ]+Cot[β]*D[ω,R,θ]-(2*(1-ν)*Cot[β]/R)*D[ω,θ]+(1/(R^2*Sin[β]))*D[ψ,θ,θ]-Sin[β]*D[ψ,R,R]-(Cos[β]/R)*D[ψ,R,β]+(2*Cos[β]/R^2)*D[ψ,β]+(Sin[β]/R)*D[ψ,R]+Cot[β]*D[χ,R,θ]]
σθθ=Simplify[(1/R)*D[ϕ,R]+(Cot[β]/R^2)*D[ϕ,β]+(1/(R*Sin[β])^2)*D[ϕ,θ,θ]+(1-2*ν)*Cos[β]*D[ω,R]+(2*ν/R)*Sin[β]*D[ω,β]+(Cos[β]^2/(R*Sin[β]))*D[ω,β]+(Cot[β]/(R*Sin[β]))*D[ω,θ,θ]+(2/(R*Sin[β])^2)*D[ψ,θ]-(2/R)*D[ψ,R,θ]-(2*Cot[β]/R^2)*D[ψ,β,θ]+(Cot[β]/(R*Sin[β]))*D[χ,θ,θ]-Cos[β]*D[χ,R]+((1+Sin[β]^2)/(R*Sin[β]))*D[χ,β]]
σβR=Simplify[(1/R)*D[ϕ,β,R]-(1/R^2)*D[ϕ,β]+(1-2*ν)*Sin[β]*D[ω,R]+Cos[β]*D[ω,β,R]-(2*(1-ν)*Cos[β]/R)*D[ω,β]+(1/R^2)*D[ψ,θ,β]+(Cot[β]/R)*D[ψ,θ,R]-(2*Cot[β]/R^2)*D[ψ,θ]+Cos[β]*D[χ,R,β]-Sin[β]*D[χ,R]]
σθβ=Simplify[(1/(R^2*Sin[β]))*D[ϕ,θ,β]-(Cot[β]/(R^2*Sin[β]))*D[ϕ,θ]+(Cot[β]/R)*D[ω,β,θ]+(2*(1-ν)/R)*D[ω,θ]-(1/(R^2*Sin[β]))*D[ω,θ]+(Cot[β]/(R^2*Sin[β]))*D[ψ,θ,β]-(Sin[β]/R)*D[ψ,R,β]-(Cos[β]/R^2)*D[ψ,β,β]+(1/(R^2*Sin[β]))*D[ψ,β]+(Cot[β]/R)*D[χ,β,θ]-D[χ,θ]/(R*Sin[β]^2)]
σββ=Simplify[(1/R)*D[ϕ,R]+(1/R^2)*D[ϕ,β,β]+(Cos[β]/R)*D[ω,β,β]+(1-2*ν)*Cos[β]*D[ω,R]+(2*(1-ν)*Sin[β]/R)*D[ω,β]+(2*Cot[β]/R^2)*D[ψ,θ,β]-2*(Cot[β]/R)^2*D[ψ,θ]+Cos[β]*D[χ,β,β]/R-Cos[β]*D[χ,R]]

T=((1-ν)/(μ*α*(1+ν)))*(Cos[β]*D[χ,R]-Sin[β]*D[χ,β]/R)
qR=((1-ν)*K/(μ*α*(1+ν)))*(Sin[β]*D[χ,R,β]/R-Cos[β]*D[χ,R,R]-Sin[β]*D[χ,β]/R^2)
qθ=((1-ν)*K/(μ*α*(1+ν)))*(D[χ,β,θ]/R^2-(Cot[β]/R)*D[χ,R,θ])
qβ=((1-ν)*K/(μ*α*(1+ν)))*(Sin[β]*D[χ,R]/R-Cos[β]*D[χ,R,β]/R+Sin[β]*D[χ,β,β]/R^2+Cos[β]*D[χ,β]/R^2)