scipy имеет пакет для wishart:
import numpy as np
from scipy.stats import wishart
x = np.linspace(1e-5, 8, 100) # make an array
for y in x:
z = wishart.rvs(1, scale=y)
print(z)
output:
4.59465858669e-06
0.00403122342709
0.0268879506122
0.0100029090879
0.129477863995
0.372787021348
....
wishart.rvs - это функция wishart для случайной выборки вариаций.