oldtests/gsl-intgamma/genvals-quadratures.py

   1 import random
   2 import math
   3
   4 nmu = 2
   5 mu = 1e-5
   6 ksigma = 1e7
   7
   8 for i in range(10000):
   9     R = abs(random.expovariate(1/mu))
  10     k = random.normalvariate(0,ksigma)
  11     n = math.ceil(abs(random.expovariate(1/nmu)))
  12     eta = abs(k/2/6**0.5)# abs(random.expovariate(8/k))
  13     print(n, k, R, eta)
  14
  15