from CoaSim import *
import sys

if len(sys.argv) != 3:
    print 'Usage:', argv[0], 'beta noSeqs'
    sys.exit(2)

beta = float(sys.argv[1])
noSeqs = int(sys.argv[2])

lengths = []
for i in xrange(100000):
    arg = simulate([], noSeqs, beta=beta, keepEmptyIntervals=True)
    lengths.append(arg.intervals[0].tree.branchLength)
print sum(lengths)/len(lengths)