The angle between two vectors in spherical coordinates is depicted in the diagram below:
We have
where and are the unit vectors of and , respectively.
Substituting eq182b (the unit vector in spherical coordinates) in eq33, noting that and rearranging gives