Closure in RPL

[Albert Chan]

My mistake. It take 125 primes to reach 12 digits accuracy.

XCas> M(n) := Mp(makelist(ithprime,n+1,1,-1));
XCas> Mp(p) := min(abs(proot(p-p[1:]))); // n-degree prime-gap polynomial min abs root

XCas> M(20)        → 0.809455733629
XCas> M(30)        → 0.807108146818
XCas> M(40)        → 0.806549211933
XCas> M(50)        → 0.806513445127
XCas> M(123)      → 0.806513599262
XCas> M(124)      → 0.806513599261
XCas> M(125)      → 0.806513599261
XCas> M(126)      → 0.806513599261