import numpy as np

A = np.array([[1], [2], [3]])
B = np.array([[4], [5], [6]])

cross_product = np.zeros((3, 1))

for i in range(A.shape[0]):
    if i == A.shape[0] - 1:
        cross_product[i] = (A[i] * B[0]) - (B[i] * A[0])
    else:
        cross_product[i] = (A[i] * B[i + 1]) - (B[i] * A[i + 1])

print(cross_product)

        __
        __