Another solution in n-dimensional space might be:

import numpy as np
vector = np.array([1,2,3,4])
vector = vector[np.newaxis,:]
orthogonal_vertices = np.linalg.svd(vector)[-1]