Насчёт ТОТО я не очень в курсе.
Но когда речь о событиях разной вероятности, я пользуюсь своей прогой на Питоне.
Например, для подсчёта вероятностей встретить определённое количество скаттеров в каком-то слоте (барабаны часто разной длины и имеют разное количество скаттеров).

import pprint, collections

def binom_ext(p):
    '''p - массив вероятностей успеха в каждом испытании
    Возвращает для каждого числа успехов кортеж (вероятность, 1/вероятность)'''
    n = len(p)
    out = collections.defaultdict(int)
    for i in range(2**n):
        i = '{0:0{1}b}'.format(i, n)
        x = sum([ c=='1'  for c in i ])
        pp = 1
        for c,p_ in zip(i,p):  pp *= p_  if c=='1' else  1-p_
        if pp:  out[x] += pp
##        print(i, x, pp)
    print(sum(out.values()))
    out = { k:(v,1/v)  for k,v in out.items() }
    pprint.pprint(out)
    return out


binom_ext([3/29, 6/41, 3/42, 3/41, 3/42])   # Hall of Gods slot, scatters
binom_ext([0, 0, 9/42, 9/41, 9/42])         # Hall of Gods slot, bonus