HHC 2019 will be in Reno, Nevada, Sept 21-22

[Thomas Okken]

Code for HP-42S or HP-41C:

Code:

00 { 72-Byte Prgm }
01▸LBL "QD"
02 2
03 STO 01
04 10↑X
05 STO 00
06▸LBL 00
07 0
08 RCL 00
09▸LBL 01
10 RCL ST X
11 10
12 STO÷ ST Z
13 MOD
14 ENTER
15 X↑2
16 ×
17 RCL ST Z
18 +
19 X<>Y
20 IP
21 X≠0?
22 GTO 01
23 R↓
24 RCL 00
25 X≠Y?
26 GTO 02
27 STO IND 01
28 1
29 STO+ 01
30▸LBL 02
31 1
32 STO+ 00
33 RCL 00
34 1ᴇ3
35 X>Y?
36 GTO 00
37 1/X
38 1
39 -
40 STO+ 01
41▸LBL 03
42 RCL IND 01
43 DSE 01
44 GTO 03
45 END

Nothing very clever or elegant, I'm afraid, but it gets the job done, taking 10 minutes and 30 seconds to return 153, 370, 371, and 407, on my HP-42S.