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)