Creating digits of pi

[toml_12953]

(02-10-2018 12:44 PM)Mike (Stgt) Wrote:  Could someone of the interested (Decimal BASIC admirer, etc.) please confirm this Ramanujan–Sato of 1993 works? I fail with it, get Pi = 4.1748739, what is a bit too far off.
I assume a typo in the Wiki.

Ciao.....Mike

The quadratic convergence works:

Code:

OPTION ARITHMETIC DECIMAL_HIGH
LET a=SQR(2)
LET b=0
LET p=2+SQR(2)
LET n=1
DO UNTIL a=1
   LET an=a
   LET bn=b
   LET pn=p
   LET a=(SQR(an)+1/SQR(an))/2
   LET b=((1+bn)*SQR(an))/(an+bn)
   LET p=((1+a)*pn*b)/(1+b)
LOOP
PRINT p
END