# Creating a matrix

X_1 = np.vstack([x1, x2])
X_2 = np.vstack([x3, x4])
X = np.hstack([X_1, X_2]).T
print(X.shape)

(1000, 2)