'Compute the variance of the dirichlet distribution.
Parameters
%(_dirichlet_doc_default_callparams)s
Returns
v : ndarray or scalar
Variance of the Dirichlet distribution.'

def var(self, alpha):

    alpha = _dirichlet_check_parameters(alpha)
    alpha0 = np.sum(alpha)
    out = ((alpha * (alpha0 - alpha)) / ((alpha0 * alpha0) * (alpha0 + 1)))
    return _squeeze_output(out)
