Эта задача связана с дискретным распределением Бернулли с плотностью
f(N,x,p) = C(N,x) * p^x * (1-p)^(N-x),
где в соответствии с условием задачи N=100, x=20..30, p=0.9
Искомая вероятность равна
P = sum(f(100,x,0.9), x=20..30) = 1.306316508 * 10^-46
Мат. ожидание M(x) = p*N = 90
Дисперсия D(x) = N*p*(1-p) = 9
В данном случае нормальное распределение, по которому проводились расчеты в предыдущих сообщениях, является приближенным. С его помощью, например, невозможно точно вычислить искомую вероятность, стремящуюся к нулю, как правильно заметил venja.