<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def probability_density(self, X):
 """Compute probability density. Arguments: X: `np.ndarray` of shape (n, 1). Returns: np.ndarray """

        self.check_fit()
        return norm.pdf(X, loc=self.mean, scale=self.std)