Вероятность заключить контракт

//Функция расчета по формуле Пуассона
double Pt ( double lambda, double t, double k)
{
    double f = 1;
    //нахождение факториала
    for(int i = 1;i <= k; i++) f *= i;
    //вернуть вероятность - результат расчета по формуле Пуассона
    return pow(lambda*t,k)* exp(-lambda*t) / f;
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
    double  lambda = 3, t = 2, k = 4, P = 0;
    printf( 
"Среднее число заказов такси,поступающих на диспетчерский пункт в одну минуту,равно 3.\n\
Найти вероятность того, что за 2 мин поступит:\n\
а) четыре вызова;\n\
б) менее четырех вызовов;\n\
в) не менее четырех вызовов.\n\n");
    printf( 
"Решение:\n\
По условию, lambda=3 , t = 2, k = 4.\n\
Воспользуемся формулой Пуассона\n\
Pt(k)=((lambda*t)^k * e^(-lambda*t))/k!\n\n");
    printf( "а) Искомая вероятность того, что за 2 мин поступит четыре вызова :\n");
    printf( "Pt(4) = %lf\n\n",Pt(lambda,t,4.));
    printf(
"б) Событие «поступило менее четырех вызовов» произойдет,\n\
если наступит одно из следующих несовместных событий:\n\
1) поступило три вызова;\n\
2) поступило два вызова;\n\
3) поступил один вызов;\n\
4) не поступило ни одного вызова.\n\
Эти события несовместны, поэтому применима теорема сложения вероятностей\n\
несовместных событиий:\n");
    //цикл сложения и нахождения вероятностей
    for(int i = 0 ; i<4; i++)
        P += Pt(lambda,t,i);
    printf("Pt( k<4 ) = Рt(3)+ Рt(2) + Рt(1) + Рt(0) = %lf\n\n",P);
    printf(
"в) События «поступило менее четырех вызовов» и «поступило\n\
не менее четырех вызовов» противоположны, поэтому искомая вероятность того,\n\
что за 2 мин поступит не менее четырех вызовов:\n");
    printf( "1 - Pt( k<4 ) = %lf\n\n",1.-P);
    system("pause");
    return 0;
}