Syntax:

psi = angdist(phi1,theta1,phi2,theta2)

This function calculates the angular distance between two points on a sphere. Given (latitude, longitude) for two points and , the angle between them is: