import numpy as np X = np.random.uniform(size=(10,3)) n,m = X.shape # for generality X0 = np.ones((n,1)) Xnew = np.hstack((X,X0))