Note that we have a constant probability, fixed number of tests, independent patients and only two outcomes, so we can write a binomial distribution X~Bin(4000,0.0005) where X is the random variable representing the number of patients who test positive.Use a Poisson approximation with L=np=2, so P(X=5) = ~(e^(-2)2^5)/5! = 0.036089