code: def sample(self, n_samples):
          
          if self.tau > 1 or self.tau < -1:
              raise ValueError("The range for correlation measure is [-1,1].")
  
          v = np.random.uniform(0, 1, n_samples)
          c = np.random.uniform(0, 1, n_samples)
  
          u = self.percent_point(c, v)
          return np.column_stack((u, v))
docs: Generate specified `n_samples` of new data from model. `v~U[0,1],v~C^-1(u|v)`
  
          Args:
              n_samples: `int`, amount of samples to create.
  
          Returns:
              np.ndarray: Array of length `n_samples` with generated data from the model.