code: static double havDistance(double lat1, double lat2, double dLng) {
          return hav(lat1 - lat2) + hav(dLng) * cos(lat1) * cos(lat2);
      }
docs: Returns hav() of distance from (lat1, lng1) to (lat2, lng2) on the unit sphere.