code: def probability_density(self, X):
          
          self.check_fit()
          return norm.pdf(X, loc=self.mean, scale=self.std)
docs: Compute probability density.
  
          Arguments:
              X: `np.ndarray` of shape (n, 1).
  
          Returns:
              np.ndarray