stat110 pset3 code

// by Raphael Rouvinov
// for stat 110 pset 3

#include <stdio.h>

int main(void)
{
    double sum = 0;
    
    for(int j = 1; j <= 48; j++)
    {
        double num = 3.0 / (52 - j);
        num *= 4.0 / (53 - j);
        
        for(int i = 1; i <= j-1; i++)
        {
            num *= (49.0-i)/(53.0-i);
        }
        
        sum += num;
    }
    
    printf("sum: %f\n", sum);
}