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.
Substitute eq182b (the unit vector in spherical coordinates) in eq33, noting that and rearranging, we have