Update docstring for multivariate_normal().

The change is to make a stronger statement that the covariance matrix must
be positive-semidefinite and that if it is not, then the results are not
guaranteed across NumPy versions.
Merge pull request #5266 from chebee7i/multivariate_normal

Possible regression in `multivariate_normal()`