points = np.array([1, 2, 3, 4])
points[:, (None)] - points[(None), :]
Given a list of n points, how can I generate a matrix that contains the distance from every point to every other point using numpy?