Как можно реализовать данную формулу?

Формула показывает ненормированную вероятность P(i) - вероятность того, что система находится в состоянии i, причем P (0) = 1, где I () - индикаторная функция, то есть она принимает значение 0, когда внутреннее выражение неверно и 1 в обратном случае. b_k и r_k - некоторые постоянные параметры, которые подаются на вход. Это формула реализует первый шаг алгоритма.

Также есть иллюстрирующий пример алгоритма.

Полностью весь алгоритм привден в книге Степанова С. Н. "Теория телетрафика: концепции, модели, приложения" на стр. 700-701. Буду рад любым предложениям и советам по реализации.

По-моему мнению должно быть что-то такое, но я не совсем понимаю как именно внедрить туда индикаторную функцию именно рекурсивно.

def prob(i, rs, bs): 
    if i <= 0:
        return 1
    return (1 / i) * sum(b * prop(i - abs(r)) for (r, b) in zip(rs, bs))