(41C) Pythagorean Triples

[C.Ret]

(06-25-2022 03:40 PM)Thomas Klemm Wrote:  But it should also work with the HP-41C after the obvious transformations.

As usual, Thomas Klemm produces an excellent algorithm. I particularly like how the GCD determination is efficiently incorporated into the code.

So I grabbed my HP-41C and its trusty 82240A to make the obvious transformations needed to aim and print:

01 LBL"PYTH"
02 STO 04  XEQ 04
04 LBL 01
05  RCL 04  STO 05
07  LBL 00
08   DSE 05  XEQ 02  DSE 05  GTO 00
12  DSE 04  GTO 01
14 GTO 04

15 LBL 02
16  RCL 04  RCL 05
18  LBL 03
19   STO Z  MOD  X≠0?  GTO 03
23  10^X  X≠Y?  RTN
26  +  CLA
28  RCL 04  ARCL X  "┝,"  ST* Y  X^2  STO Z
34  RCL 05  ARCL X  ACA   ST* Z  X^2  ST+ T
40  -  FMT  ACX  RDN  ACX  RDN  ACX
47  LBL 04
48   ADV
49 END

   

Shorter than the original code, this one just displays triples in a different order.