[VA] SRC #008 - 2021 is here !

[J-F Garnier]

(01-07-2021 09:37 PM)Valentin Albillo Wrote:  ... I've concocted this unoptimized 5-line (250-byte) program for the HP-71B, which exactly computes the 322 digits of the max. product P = 2*3⁶⁷³ in half a second using the excellent J-F Garnier's Emu71 emulator

Thanks :-)
Emu71/DOS is more than 20 years old but still going strong !

(01-07-2021 10:40 PM)Albert Chan Wrote:  

(01-07-2021 09:37 PM)Valentin Albillo Wrote:  

The funny thing about this integral is that its value remains the same if you replace 2.021 by any positive real value !!

Can you explain why ?

Yes, please Valentin; I was expecting to get an explanation or at least some indications.

I used this effect to actually compute the integral on the HP-32S with the value k=2 to replace the exponentiation operations with simple squaring, for a much faster computation.

Also I computed the integral in two parts, from 0 to pi/2 and pi/2 to pi, since the pi/2 value causes a 0/0 limit of the function. The HP algorithm is robust when this kind of limit is at the edges, but not always when it is within the integration range especially when using high target accuracy (I used accuracy= 1E-10).

J-F