HP 32SII: Random Number Utilities
Note: The following programs should work with the original HP 32S. The checksums listed are for the HP 32SII only.
HP 32SII Program: Fill a Set of Registers with Increasing Integers
This program fills up the alpha registers with consecutive increasing integers. Up to 25 slots are available, A- Y, with the lowest being stored in A. Register Z is used to store the difference.
The difference from both numbers should not be more than 24.
F01 LBL F
F02 x<>y
F03 1
F04 -
F05 STO Z
F06 -
F07 STO i
I01 LBL I
I02 RCL i
I03 IP
I04 RCL+ Z
I05 STO(i)
I06 DSE i
I07 GTO I
I08 RTN
F: 10.5 bytes, 32SII checksum 9FB9
I: 12.0 bytes, 32SII checksum E189
Total: 22.5 bytes
Examples
18 to 27: 18 ENTER 27 XEQ F
Results:
A = 18, B = 19, C = 20, D = 21, E = 22,
F = 23, G = 24, H = 25, I = 26, J = 27
HP 32SII Program: Pick from Registers A through Y
Enter the number of registers to pick from into register Z. For example, if you have registers A through J loaded, store 10 into register Z.
Note: When used for indirect addressing, the absolute value of integer of register i is used. This will allows us to save a few steps.
G01 RANDOM
G02 RCL× Z
G03 1
G04 +
G05 STO i
G06 RCL(i)
G07 RTN
G: 12.0 bytes, 32SII checksum FE18
Example:
Use the data from the last problem.
10 STO Z,
XEQ G: 22
XEQ G: 21
(results will vary)
Source for the next two programs:
Chamberlain, Gary. "HP-29C Random Number Generators" PPC Journal V6 N7 October 1979
HP 32SII Program: Uniform Distribution
This program picks random numbers between two parameters.
U01 LBL U
U02 STO C
U03 -
U04 STO D
W01 LBL W
W02 RANDOM
W03 RCL× D
W04 RCL+ C
W05 STOP
W06 GTO W
U: 6.0 bytes, 32SII checksum 8022
W: 9.0 bytes, 32SII checksum 88EE
Total: 15.0 bytes
Example:
Random numbers between -1 and 1:
XEQ U: -0.74076 (press R/S for more random numbers)
0.28000
-0.80090
0.20653
0.50614
HP 32SII Program: Gaussian Distribution of Random Numbers
Generate random numbers that fit on a normal distribution with mean M and deviation D. Two random numbers are generated, press x<>y to see both results.
X01 LBL X
X02 INPUT M
X03 INPUT D
X04 RAD (radians mode)
Y01 LBL Y
Y02 RANDOM
Y03 2
Y04 ×
Y05 π
Y06 ×
Y07 RANDOM
Y08 LN
Y09 -2
Y10 ×
Y11 SQRT
Y12 θ,r→y,x
Y13 RCL× D
Y14 RCL+ M
Y15 x<>y
Y16 RCL× D
Y17 RCL+ M
Y18 STOP
Y19 GTO Y
Example:
Mean = 3
Deviation: 1.5
XEQ X
M? 3 R/S
D? 1.5 R/S
1.70434 x<>y 2.89313 R/S
2.22287 x<>y 2.50113 R/S
1.67644 x<>y 2.87034 ...
(results will vary)
Happy programming,
Eddie
All original content copyright, © 2011-2022. Edward Shore. Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited. This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author.
