answer := 1 - a!/((a-k)!*a^k); p := k -> evalf((1 - 365!/((365-k)!*365^k)),4): p10 := p(10); p20:=p(20); p23 := p(23); p40:=p(40); p50:=p(50);