x= rho Sin[phi]Cos[theta];
y=rho Sin[phi]Sin[theta];
z=rho Cos[phi];
T={x,y,z}

Jac=D[T,{{rho,phi,theta}}]
MatrixForm[Jac]
Det[Jac]
Simplify[%]