--------
Instantiate a distribution with the desired parameters:

>>> from scipy import stats
>>> X = stats.Normal(mu=1., sigma=2.)

Evaluate the variance:

>>> X.mean()
1.0
>>> X.mean() == X.moment(order=1, kind='raw') == X.mu
True

r