Программа для островной модели без мутации

Необходима программа подсчитывающая вероятности островной модели с L островов и N аллей, без мутации.

Для 2-х островов и 2-х аллелей составлял таблицу с формулами в Exel, дальше уже не получилось. Формулы для общего случая в книжке на странице 23-26.

Требуется: чтобы выводил элементы переходной вероятности матрицы цепи Маркова, а дальше генерация, процесс которой описан в скриншоте (могу скинуть для 2-х островов и 2-х аллей как делал), после графики.

Реализация: R, Octave, Python. Могу докинуть кучей материала, который нашел. Знаю, что это не сильно запарно, если знать тему :)