syms r fa the x=r*sin(fa)*cos(the); y=r*sin(fa)*sin(the); z=r*cos(fa); ja=jacobian([x y z],[r,fa,the])