import numpy as np
import matplotlib.pyplot as plt
from sklearn import datasets
import cvxpy as cvx
X ,y = datasets.make_blobs(n_samples=300, n_features=2,centers=2,cluster_std=[2.0,2.0],random_state=123)
N = X.shape[0]
d = X.shape[1]
y = np.vstack(y)
print(X)
I = np.vstack(np.ones((N,1)))
print(I)
X = np.hstack((I,X))
print(X)