[VA] SRC #016 - Pi Day 2024 Special

[C.Ret]

(03-16-2024 11:13 AM)J-F Garnier Wrote:  Here is my keystroke sequence on the 71B, with my Math 2B to save a few keystrokes:

Nice, the method strangely match the method I used on my HP-15C.

You didn't indicate that you put your highly extended HP-71B in RADIANS trigonometric mode. Or isn't it necessary by using the top version of all the Math's module ?

Since, I only own the former Math 1A version, I get an erroneous result when omitting to set trigonometric mode right.


>RADIANS
>A=INTEGRAL(0,1,1E-8,ATAN(ATANH(IVAR))/IVAR)
            1.02576050997
>P=FNROOT(1,6,FVAR*LOG(GAMMA(1/FVAR)/GAMMA(1/FVAR+.5)/SQR(FVAR))-A)
            3.14159265196


And since I have my HP-71B in hands, I modestly present you my attempt for the second apparence:

It is a (too long to edit) one-liner:
10 WHILE K<N @ K=K+1 @ C=C+(PGCD(CEIL(N*RND),CEIL(N*RND))=1) @ END WHILE @ DISP SQR(6*N/C);N @ N=10*N @ RUN


>DESTROY ALL
>N=1
[RUN]
2.44948974278  1
2.73861278753  10
3.03821810125  100
3.12347523777  1000
3.14502731861  10000
3.14037146511  100000


The appearance isn't efficient at all, perhaps my code is wrong !?
I may have use 1+IP(N*RND) instead of the CEIL functions.