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