'Draw random samples from a multivariate normal distribution.
Parameters
%(_doc_default_callparams)s
size : integer, optional
Number of samples to draw (default 1).
Notes
%(_doc_callparams_note)s
Returns
rvs : ndarray or scalar
Random variates of size (`size`, `N`), where `N` is the
dimension of the random variable.'

def rvs(self, mean=None, cov=1, size=1):

    (dim, mean, cov) = _process_parameters(None, mean, cov)
    out = np.random.multivariate_normal(mean, cov, size)
    return _squeeze_output(out)
