I don't think you need to simulate this, the number of "deep work" slots should be (λD+1)(1-d)^(λD) where d=(Δ+θ)/D and D is the length of your workday.

I also invoke the Knuth Defence ("I have only proved it correct, not tried it"), and Cunningham's Law.