from sympy import *
from IPython.display import display,Math
x,y,z = symbols("x y z", cls=Function)
r,theta,phi = symbols(r"r \theta \phi")
r_vector_symbol = symbols(r"\vec{r}")
x = r*sin(theta)*cos(phi)
y = r*sin(theta)*sin(phi)
z = r*cos(theta)
r_vector = Matrix([x,y,z])
dr_vectordr_symbol = symbols(r'\frac{\partial{\vec{r}}}{\partial{r}}')
dr_vectordtheta_symbol = symbols(r'\frac{\partial{\vec{r}}}{\partial{\theta}}')
dr_vectordphi_symbol = symbols(r'\frac{\partial{\vec{r}}}{\partial{\phi}}')
dr_vectordr = diff(r_vector,r)
dr_vectordtheta = diff(r_vector,theta)
dr_vectordphi = diff(r_vector,phi)
display(Math(f'{r_vector_symbol}={latex(r_vector)}'))
display(Math(f'{dr_vectordr_symbol}={latex(dr_vectordr)}'))
display(Math(f'{dr_vectordtheta_symbol}={latex(dr_vectordtheta)}'))
display(Math(f'{dr_vectordphi_symbol}={latex(dr_vectordphi)}'))