Xk = np.atleast_2d(X).T