for(n in 1:99)cat(sprintf("n=%d, p=%.8f
",1+n,1-prod(1-1:n/365)))