distlatlon(lat1, lon1, lat2, lon2 [, method] )
distlatlon function calculates the distance, in miles,
between two points, represented in latitude/longitude pairs in decimal
distgeocode, it uses the "Great Circle" method by
default, but can be overridden to use the faster, less accurate
Pythagorean method if "
pythagorean" is passed as the optional