distance.cdist is every pairwise distance. Rows are Nairobi, Mombasa, Kisumu as (x, y) in comments.

Goal

Print the 3×3 kilometre-ish distance matrix.

from scipy.spatial import distance
pts = np.array([[0.0, 0], [4.8, 0], [3.0, 2.0]])  # Nairobi, Mombasa, Kisumu toy km/100
print(distance.cdist(pts, pts).round(2))

from scipy.spatial import distance
a = np.array([[0.0, 0]])
b = np.array([[3.0, 4]])
print(distance.euclidean(a[0], b[0]))

from scipy.spatial import ConvexHull
pts = np.array([[0, 0], [1, 0], [0.5, 1], [0.2, 0.2]], dtype=float)
print(ConvexHull(pts).vertices)

from scipy.spatial import distance
print(distance.pdist(np.array([[0.0, 0], [1, 0], [0, 1]])).round(3))