This program prints the % accuracy at each level. The idea is to bracket the % to form an estimate of the next unknown count of primes at 10^n. Over 100% the true value is less. Under the true value is more. A ratio may come close to the true value (by manipulating between 9.1151 and 10). 10 point 2 20 'A141306, Enoch Haga, 28 Jun 2008 30 X=2:Y=2:goto 60 40 B=odd(X):if B=1 and A=0 then Y=Y+2-2 50 A=even(X):if A=1 and B=0 then Y=Y+4-2 60 C=C+Y:if C>=1 then 70 else 80 70 N=(log(9.1151^X)*10^(X-1))/C 80 print A;B;X;C;N; 90 if X<=23 then X=X+1 100 read D 110 data 25,168,1229,9592,78498,664579,5761455,50847534,455052511 120 data 4118054813,37607912018,346065536839,3204941750802 130 data 29844570422669,279238341033925,2623557157654233,24739954287740860 140 data 234057667276344607,2220819602560918840,21127269486018731928 150 data 201467286689315906290,1925320391606803968923 160 R=N/D*100:print R;"% accuracy" 170 goto 40