import numpy as np
# Define X and Y arrays
X = np.array([1, 2, 3])
Y = np.array([4, 5, 6])
# Compute the Cauchy matrix using broadcasting
Cauchy = 1 / (X[:, np.newaxis] - Y)
print("Cauchy matrix:")
print(Cauchy)