[VA] SRC #012a - Then and Now: Probability

[J-F Garnier]

(10-16-2022 03:07 PM)Albert Chan Wrote:  

(10-16-2022 02:10 PM)J-F Garnier Wrote:  HP-75 result is: 9.51234350244E-6

For machine that round-away-from-zero, I would scale by gcd(2,4,6) = 12 instead.
Or, scale by 1.2 for decimal machine, to keep denominator (1.2^S) from overflow, even with big S

< 40 FOR I=1 TO R @ FOR J=1 TO C @ W(I,J)=3/(3-(J=1)-(I=J)-(I=R)) @ NEXT J @ NEXT I
> 40 FOR I=1 TO R @ FOR J=1 TO C @ W(I,J)=0.6/(3-(J=1)-(I=J)-(I=R)) @ NEXT J @ NEXT I

< 180 DISP TIME-T0;R;S;T/6^S
> 180 DISP TIME-T0;R;S;T/1.2^S

With above changes, on emu71, I get: [..]
9.51234350204E-6

With these changes, I now get on the HP-75 (actually emu75):
9.51234350227E-6
which is indeed better, still not at the level of the best HP-71 performance but closer to my first HP-71 result (9.51234350213E-6).

J-F