OSTI Preprint With 41C Code

[Howard Owen]

While searching for something else, I ran across this 1981 preprint from LLNL. It contains an HP41C program in an appendix. There are several other documents on the site that mention the HP41C, but this was the only one with code.

The DOE's Office of Scientific and Technical Information looks like it might be a useful source of HP calculator information. I only searched for HP41C, but other searches might be fruitful too.

[Peter Murphy]

I was a technical editor at LLNL when that was written! Too bad I wasn't into HP calculators then.

The pdf version of the paper is searchable, but the listing (in case you need it) looks nearly unreadable in places. A paper copy is available for purchase, however.

[Thomas Klemm]

(09-16-2015 10:26 PM)Peter Murphy Wrote:  the listing (in case you need it) looks nearly unreadable in places.

Code:

58 RCL 01
59 12.00
60 /
61 FRC
62 12.00
63 *
64 "INCHES="
65 ARCL X
66 AVIEW

Of course this could be much simpler:

Code:

RCL 01
12
MOD

But it appears that memory consumption wasn't an issue yet.

Thanks for sharing the document. How would that measuring be done today? Probably with a laser?

Cheers
Thomas

[Howard Owen]

Here's a transcription from the PDF. I made what I think are reasonable interpretations of the dodgy characters in the listing. Sample inputs and program, outputs appear after the listing. I have no idea if the sample inputs are reasonable for the application. The comments in the PDF make clear that the program is correcting measurement results for thermal expansion of various metals. That, combined with the idea that the environment involves nuclear power led me to use high values for the temperature. Since the temperature range goes over the melting point of aluminum, you can assume that this run was for iron or steel tape.

Code:


  01◆LBL "TAPE"
 02 FIX 3
 03 "BASE DEG C ?"
 04 PROMPT
 05 STO 00
 06◆LBL 01
 07 "TAPE FEET?"
 08 PROMPT
 09 12.00
 10 *
 11 STO 01
 12 "TAPE INCHES?"
 13 PROMPT
 14 ST+ 01
 15 "DIAL INCHES?"
 16 PROMPT
 17 STO 02
 18 "TEMP DEG C ?"
 19 PROMPT
 20 RCL 00
 21 -
 22 STO 03
 23 STO 04
 24 STO 05
 25 17.3 E-6
 26 ST* 05
 27 23.5 E-6
 28 ST* 04
 29 11.6 E-6
 30 ST* 03
 31 8.50
 59 12.00
 60 /
 61 FRC
 62 12.00
 63 *
 64 "INCHES="
 65 ARCL X
 66 AVIEW
 67 STOP
 68 GTO 01
 69◆LBL 09
 70 *
 71 ST+ 01
 72 .END.
IN=600, 6, 5, 2, 920
FEET=6.000
INCHES=5.380

[Sylvain Cote]

Hello Howard

IMHO you made two errors in your transcription ...

25 1.73 E-6 (should be) 17.3 E-6 (see figure 3)
27 33.5 E-6 (should be) 23.5 E-6 (see figure 3)

Thank you for sharing.

Best regards,

Sylvain

[Howard Owen]

(09-18-2015 08:18 AM)Sylvain Cote Wrote:  IMHO you made two errors in your transcription ...

Both good catches Sylvain. The change on line 25 is clearly correct. The 2 vs the 3 in line 27 is a tougher call. I think you are correct because 23 vs. 33 in the mantissa seems more in line with the other values.

I've corrected the listing. Thanks!

[Sylvain Cote]

(09-18-2015 01:19 PM)Howard Owen Wrote:  Both good catches Sylvain. The change on line 25 is clearly correct. The 2 vs the 3 in line 27 is a tougher call. I think you are correct because 23 vs. 33 in the mantissa seems more in line with the other values.

Howard,

Look at page 22 (Figure 3) of the document, the values are clearly indicated.

Sylvain