I used ac=ac+z/n/2 when it should have been ac=ac+2*z/n and for RDPRP ac=ac+0.16*2*z/n.
Along with accumulated figures:
Code:
[5, 26, 0.03432083, 0.03432083]
[6, 98, 0.01985533, 0.05417616]
[7, 314, 0.02344581, 0.07762197]
[8, 1608, 0.02832678, 0.1059487]
[9, 15072, 0.1222199, 0.2281686]
[10, 101630, 0.08648315, 0.3146518]
For RDPRP:
Code:
[5, 26, 0.005491333, 0.005491333]
[6, 98, 0.003176853, 0.008668186]
[7, 314, 0.003751329, 0.01241951]
[8, 1608, 0.004532284, 0.01695180]
[9, 15072, 0.01955518, 0.03650698]
[10, 101630, 0.01383730, 0.05034429]
So pseudoprimes can be expected at about 14 and 19 digits (resp.). The prize money is within reach!