In[1] := f[r_,p_,l_]:={r Cos[p]Cos[l],r Cos[p]Sin[l], r Sin[p]} In[2] := jf[a_,b_,c_]:=D[f[r,p,l],{{r,p,l}}] /. {r->a,p->b,l->c} In[3] := MatrixForm[jf[r,p,l]] In[4] := Det[jf[r,p,l]] In[5] := Simplify[%]