code: def add_intercept_term(self, x):
          
          nr_x,nr_f = x.shape
          intercept = np.ones([nr_x,1])
          x = np.hstack((intercept,x))
          return x
docs: Adds a column of ones to estimate the intercept term for 
          separation boundary