#loading library
import numpy as np

# Create a random vector of coordinates
vec = np.random.rand(5, 2)

# Calculate the distances between all pairs of points using broadcasting
distances = np.sqrt(((vec[:, np.newaxis] - vec) ** 2).sum(axis=2))

#printing result
print(distances)