code: def cumulative_distribution(self, X):
          
          self.check_fit()
          return norm.cdf(X, loc=self.mean, scale=self.std)
docs: Cumulative distribution function for gaussian distribution.
  
          Arguments:
              X: `np.ndarray` of shape (n, 1).
  
          Returns:
              np.ndarray: Cumulative density for X.