HP 67 Collection: October 2026

HP 67 Collection: October 2026

All examples are done in FIX 4. (DSP 4)

Earth's Gravity at a Specific Latitude

The true earth gravity force depends on several factors, such as latitude and altitude. This estimation takes the latitude (north/south), in degrees-minutes-seconds, and latitude, in kilometers, into account.

The constant 9.80665 m/s^2 is an accepted average, the true force varies.

g_s = 9.806 – 0.026 * cos(2 * φ)

g_h = g_s * (6378.1 / (6378.1 + h))

HP 67 Code:

LBL A: 001 31, 25, 11

DEG: 002 35, 41

H←: 003 31, 74

2 : 004 02

× : 005 71

COS: 006 31, 63

. : 007 83

0 : 008 00

2 : 009 02

6 : 010 06

× : 011 71

9 : 012 09

. : 013 83

8 : 014 08

0 : 015 00

6 : 016 06

x<>y: 017 35, 52

- : 018 51

x<>y: 019 35, 52

6 : 020 06

3 : 021 03

7 : 022 07

8 : 023 08

. : 024 83

1 : 025 01

+ : 026 61

LST x: 027 35, 82

x<>y: 028 35, 52

÷ : 029 81

× : 030 71

RTN : 031 35, 22

Inputs:

Y: h

X: φ

Examples:

Inputs: h = 0, φ = 2°33’: g ≈ 9.7801

Inputs: h = 0.762, φ = 52°17’: g ≈ 9.8114

Source 

Grainger Engineering Office of Marketing and Communications. (answer written by Rebecca H.) (2016, November 21). “How gravitational force varies at different locations on Earth.” Illinois. https://van.physics.illinois.edu/ask/listing/64061. Retrieved March 10, 2026.

Triangulation

With angles α and ß:

base = dist * (1/tan α + 1/tan ß)

dist = base ÷ (1/tan α + 1/tan ß)

Labels:

LBL A: Enter angles α and ß in the current angle mode. Returns 1/tan α + 1/tan ß.

LBL B: Base. When a number is entered prior to pressing [ B ], the number is stored as a base. Otherwise, the program solves for a base.

LBL C: Distance (of a line that bisects the base). When a number is entered prior to pressing [ C ], the number is stored as a base. Otherwise, the program solves for a base.

HP 67 Code:

LBL A: 001 31, 25, 11 (enter the angles)

TAN: 002 31, 64

1/x: 003 35, 62

x<>y: 004 35, 52

TAN: 005 31, 64

1/x: 006 35, 62

+: 007 61

STO 5: 008 33, 05

RTN: 009 35, 22

LBL B: 010 31, 25, 12 (base)

STO 3: 011 33, 03

F? 3: 012 35, 71, 03

RTN: 013 35, 22

RCL 4: 014 34, 04

RCL 5: 015 34, 05

× : 016 71

STO 3: 017 33, 03

RTN: 018 35, 22

LBL C: 019 31, 25, 13 (distance)

STO 4: 020 33, 04

F? 3: 021 35, 71, 03

RTN: 022 35, 22

RCL 3: 023 34, 03

RCL 5: 024 34, 05

÷ : 025 81

STO 4: 026 33, 04

RTN: 027 35, 22

Examples:

α = 30°, ß = 65°

Set degree mode. 30 [ ENTER ] 65 [ A ]. Result: 2.1984

Distance = 24, solve for base. 24 [ C ] [ B ]. Result: 52.7606

Base = 24, solve for distance. 24 [ B ] [ C ]. Result: 10.9172

Source:

“Triangulation (surveying)” Wikipedia. https://en.wikipedia.org/wiki/Triangulation_(surveying) (last edited October 24, 2025). Retried March 17, 2026

Time of Runoff on Turf

t = 2.5 * (L / (S * I^2))^(1/3)

t: time of runoff (rain off) in minutes

2.5: turf constant

L: distance from remote area in feet

S: slope of the terrain

I: rain intensity in inches per hour

Input Stack:

Z: L

Y: S

X: I

Output Stack:

X: t

HP 67 Code:

LBL A: 001 31, 25, 11

x² : 002 32, 54

× : 003 71

÷ : 004 81

3 : 005 03

1/x: 006 35, 62

y^x: 007 35, 63

2 : 008 02

. : 009 83

5 : 010 05

× : 011 71

RTN: 012 35, 22

Examples:

Inputs: Z: 900 ft, Y: 0.11, X: 4 in/hr. Result: 19.9917 min

Inputs: Z: 650 ft, Y: 0.33, X: 0.06 in/hr. Result: 204.4744 min

Source

Texas Instruments. Slide Rule Calculator SR-50A Manual. 1975. Manual from datamath.org: http://www.datamath.net/Manuals/SR-50A_US.pdf pg. 59

Beamwidth of a Radio Telescope

This program calculates the beam width of a radio telescope in degree-minutes-seconds.

Θ = 67625 ÷ (D * f)

Θ = beam width

D = the diameter of a circular telescope in feet

f = observing frequency in MHz (megahertz)

Inputs:

Y: D

X: f

Outputs:

X: Θ

HP 67 Code:

LBL B: 001 31, 25, 12

× : 002 71

1/x : 003 35, 62

6 : 004 06

7 : 005 07

6 : 006 06

2 : 007 02

5 : 008 05

× : 009 71

→H.MS: 10 32, 74

RTN: 11 35, 22

Examples:

Inputs: X: 1980 MHz, Y: 106 ft. Output: Θ: 0.1920 (19’20”)

Inputs: X: 1665 MHz, Y: 84 ft. Output: Θ: 0.2901 (29’01”)

Source

Ball, John A. Algorithms for RPN Calculators John Wiley & Sons: New York, NY. 1978. ISBN 0-471-03078-8. pg. 259

Fan Laws

p2 = (rpm2 ÷ rpm1)^n * p1

(p1 and p2 stand for various parameters, see the input syntax below)

Inputs:

1^st law: (cfm: cubic feet per minute, rpm: revolutions per minute)

Calculate cfm2:

T: cfm1

Z: rpm2

Y: rpm1

X: 1

2^nd Law: (sp: static pressure, rpm: revolutions per minute)

Calculate sp2:

T: sp1

Z: rpm2

Y: rpm1

X: 2

3^rd Law: (hp: horsepower, rpm: revolutions per minute)

Calculate hp2:

T: hp1

Z: rpm2

Y: rpm1

X: 3

HP 67 Code:

LBL B: 001 31, 25, 12

R↓ : 002 35, 53

÷ : 003 81

R↑ : 004 35, 54

y^x: 005 35, 63

× : 006 71

RTN: 007 35, 22

Examples

1^st Law: T: cfm1 = 2000, Z: rpm2 = 1100, Y: rpm1 = 950, X: n = 1

Result: cfm2 ≈ 2315.7895

2^nd Law: T: sp1 = 0.5, Z: rpm2 = 1100, Y: rpm1 = 950, X: n = 2

Result: sp2 ≈ 0.6704

3^rd Law: T: hp1 = 1.5, Z: rpm2 = 1100, Y: rpm1 = 950, X: n = 3

Result: hp2 ≈ 2.3286

Source

Gladstone, John. Air Conditioning Testing and Balancing: A Field Practice Manual A Van Nostrand Reinhhold Book. August 1976. pp. 24-32

More HP 67 goodness to come in November,

Eddie

All original content copyright, © 2011-2026. 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.

Casio fx-50F/Radio Shack EC-4024 Program Collection: October 2026

Casio fx-50F/Radio Shack EC-4024 Program Collection: October 2026

Notes:

* The fx-50F/EC-4024 has a 29 step program space between two program spaces. Program steps can not be edited but can be cleared by [SHIFT] [ C ] (PCL) key sequence.

* Formula calls can not be included as a program step, but scientific constants can be. The values for the constants most likely are from the Japan Standards of 1978 (JISZ-8020-1978), so a little dated to 2026 standards.

* The Casio fx-10F is similar as fx-50F/EC-4024 except that the fx-10F does not have the scientific constants and has a slightly inclined keyboard.

* Mode changes can be included. SHIFT and Kin arithmetic calls are merged.

The programs presented today require inputs to be stored in registers 1-6 using [ SHIFT ] [ Kout ] (Kin) ((K)onstants In).

Angle Between 2 Lines (in Degrees)

Θ = arctan( abs((y – x) ÷ (1 + x * y)) )

abs(x) = √(x^2), key sequence: [ x² ] [ √ ]

x: slope of line 1 (Kin 1)

y: slope of line 2 (Kin 2)

Program Code (20 steps):

MODE 4 (degree)

(

(

Kout 2

-

Kout 1

)

÷

(

1

+

Kout 1

×

Kout 2

)

)

x²

√

TAN^-1

=

Maximum and Minimum of Two Numbers

The two numbers are stored in registers 1 and 2 (Kin 1, Kin 2).

max(x,y) = (x + y + abs(x – y)) ÷ 2 = min(x,y) + abs(x – y)

min(x,y) = (x + y - abs(x – y)) ÷ 2 = max(x,y) - abs(x – y)

Program Code (21 steps):

Kout 2

-

Kout 1

=

x²

√

Kin 3

(

Kout 1

+

Kout 2

+

Kout 3

)

÷

2

=

HLT (show maximum)

-

Kout 3

= (show minimum)

Terminal Velocity of a Spherical Object

vt = √(2 * g * mass ÷ (ρ * cd * exposed area of the object) (units: m/s)

≈ √(34.05656 * mass ÷ (π * radius^2))

The following constants are assumed:

g = 9.80665 m/s² (Earth’s gravitational constant)

ρ ≈ 1.225 kg/m³ (air pressure)

cd = 0.47 (coefficient drag of a sphere)

2 * g ÷ (ρ * cd) ≈ 34.06565 m^4/(s^2 kg)

Program Code (19 steps):

Kout 1

×

3

4

.

0

6

5

6

5

÷

(

π

×

Kout 2

x²

)

=

√

=

Ideal Thermal Voltage and Shockley Diode Equation

Thermal Voltage:

vt = k * T ÷ q * 10^3 (mV)

Shockley Diode Equation:

Id = Is * (exp(vD ÷ (n * vt)) – 1)

Ideal: n = 1

Inputs:

Is: scale current in milliamps (mA) Kin 1

Vd: diode voltage in millivolts (vA) Kin 2

T: temperature in degrees Kelvin (K) Kin 3

Outputs:

vt: Thermal Voltage (mA),

Id: diode current (mV)

Constant used (see comments on the scientific constants above):

k: Boltzmann Constant ([SHIFT] [ 8 ]): 1.380662 * 10^-23 J/K

(q) e: Electron Charge ([SHIFT] [ 4 ]): 1.6021892 * 10^-19 C

Program Code (21 steps):

Kout 3

×

k (SHIFT 8)

÷

e (SHIFT 4)

×

3

10^x

=

HLT (see vt: thermal voltage)

1/x

×

Kout 2

=

e^x

-

1

=

×

Kout 1

=

Source:

“Shockley Diode Equation” Wikipedia. Last Edited November 13, 2025.

https://en.wikipedia.org/wiki/Shockley_diode_equation

Retrieved June 19, 2026.

Eddie

All original content copyright, © 2011-2026. 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.

Earth's Radius by Latitude

Earth's Radius by Latitude

Introduction: Calculating the Earth’s Radius

In quick, general calculations, we assume that the shape of the Earth approximates the shape of sphere, with the radius to generally be approximated as 6371 km. However, the shape of the Earth is closer to an oblate ellipsoid. The Earth is flattened at the poles and bulges the greatest at the Equator.

There are several calculations to estimate the Earth’s radius given the latitude (degrees North or South from the Equator), this blog is working with the WGS 84 formula (see Source):

R = √[ ((a² * cos(φ))² + (b² * sin(φ))²) ÷ ((a * cos(φ))² + (b * sin(φ))²) ]

Which is simplified from:

R = √[ (a^4 + b^4 * tan²(φ)) ÷ (a² + b² * tan²(φ))

where:

φ: latitude (in degrees)

a: semi-major axis: 6378.137 km

b: semi-minor axis: 6356.7523142 km

R: radius of kilometers

HP 15C Code: Earth’s Radius by Latitude

Step; Key; Key Code

001: LBL E; 42, 21, 15

002: DEG; 43, 7

003: →H; 43, 2

004: TAN; 25

005: x²; 43, 11

006: STO 1; 44, 1

007: 6; 6

008: 3; 3

009: 7; 7

010: 8; 8

011: . ; 48

012: 1; 1

013: 3; 3

014: 7; 7

015: STO 2; 44, 2

016: 4; 4

017: 14; y^x

018: 6; 6

019: 3; 3

020: 5; 5

021: 6; 6

022: . ; 48

023: 7; 7

024: 5; 5

025: 2; 2

026: 3; 3

027: 1; 1

028: 4; 4

029: 2; 2

030: STO 3; 44, 3

031: 4; 4

032: y^x; 14

033: RCL 1; 45, 1

034: ×; 20

035: +; 40

036: RCL 2; 45, 2

037: x²; 43, 11

038: RCL 3; 45, 3

039: x²; 43, 11

040: RCL 1; 45, 1

041: ×; 20

042: +; 40

043: ÷; 10

044: √; 11

045: RTN; 43, 32

TI-60 Code: Earth’s Radius by Latitude

Step; Key; Key Code

Set degrees mode before beginning.

00: DMS-DD; 39

01: TAN; 34

02: x²; 96

03: STO; 61

04: 1; 01

05: ( ; 53

06: 6; 06

07: 3; 03

08: 7; 07

09: 8; 08

10: . ; 93

11: 1; 01

12: 3; 03

13: 7; 07

14: STO; 61

15: 2; 02

16: y^x; 45

17: 4; 04

18: +; 85

19: 6; 06

20: 3; 03

21: 5; 05

22: 6; 06

23: . ; 93

24: 7; 07

25: 5; 05

26: 2; 02

27: 3; 03

28: 1; 01

29: 4; 04

30: 2; 02

31: STO; 61

32: 3; 03

33: y^x; 45

34: 4; 04

35: ×; 25

36: RCL; 71

37: 1; 01

38: ) ; 54

39: ÷; 55

40: ( ; 53

41: RCL; 71

42: 2; 02

43: x²; 96

44: +; 85

45: RCL; 71

46: 3; 03

47: x²; 96

48: ×; 65

49: RCL; 71

50: 1; 01

51: ) ; 54

52: = ; 95

53: √; 86

54: R/S; 13

55: RST; 22

Examples

Latitude: 80° 00’; Radius (km): 6357.402412

Latitude: 57° 24’; Radius (km): 6362.996788

Latitude: 43° 40’; Radius (km): 6367.986902

Latitude: 9° 15’; Radius (km): 6377.588959

Latitude: 2° 56’; Radius (km): 6378.081467

Source

Planetcalc. “Earth Radius by Latitude (WGS 84)” Timur. 2021. https://planetcalc.com/7721/ Retrieved March 22, 2026.

Eddie

All original content copyright, © 2011-2026. 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.

TI-60 and HP 65: Distance by Stadia Tacheometry

TI-60 and HP 65: Distance by Stadia Tacheometry

Introduction

The stadia calculation measures the distance from the level or theodolite (measuring device) to a graduated staff. Two readings are taken from the graduated staff, an upper reading (UR) and a lower reading (LR). The word tacheometry means "swift calculation". The distance between the theodolite and the staff is calculated as:

dist = k * s

k = is a factor taken of a radius of focal distance over image distance. Typically, k is set to 100, and will be assumed that k = 100 for these programs.

s = the difference between the upper reading and lower reading on the graduated staff. (UL - RL). The readings are assumed to be in meters.

If the theodolite is tilted at angle α, known as the vertical angle, then the distance becomes:

dist = k * s * cos² α

α is usually given in degrees-minutes-seconds and must be converted to decimal degrees.

TI-60 Program: Stadia Tracheotomy

Store before running: R1: UR (m), R2: LR (m), R3: α (D.MMSS)

Code:

00: 1 ; 01

01: 0 ; 00

02: 0 ; 00

03: × ; 65

04: ( ; 53

05: RCL; 71

06: 1 ; 01

07: - ; 75

08: RCL; 71

09: 2 ; 02

10: ) ; 54

11: × ; 65

12: RCL; 71

13: 3 ; 03

14: DMS-DD; 39

15: cos; 33

16: x²; 96

17: = ; 95

18: R/S; 13

19: RST; 22

HP 65 Program: Stadia Tracheotomy

Input Stack:

Z: vertical angle in degrees-minutes-second (V.MS)

Y: upper reading (UR)

X: lower reading (LR)

Code:

23: LBL

14: D

51: -

02: 2

32: f^-1

08: LOG (10^x)

71: ×

35 07: x<>y

32: f^-1

03: →D.MS (→DD)

31: f

05: COS

32: f^-1

09: √ (x²)

71: ×

24: RTN

Examples

Example 1:

V: 1°05' (Z stack, R3)

UR: 2.014 m (Y stack, R1)

LR: 1.668 m (X stack, R2)

Result: distance ≈ 34.5876 m

Example 2 (from HP 35):

V: 4°17' (Z stack, R3)

UR: 3.144 m (Y stack, R1)

LR: 1.761 m (X stack, R2)

Result: distance ≈ 137.5285 m

Example 3 (from HP 35):

V: -7°21' (Z stack, R3)

UR: 2.817 m (Y stack, R1)

LR: 0.731 m (X stack, R2)

Result: distance ≈ 205.1860 m

Sources

CivilFerba "Measure the distance by stadia method" Video posted on YouTube on December 8, 2018. https://www.youtube.com/watch?v=oon5ayl9DYs Retrieved March 22, 2026

HP-35 Surveying. Hewlett Packard. Buchs, Switzerland. February 1973. pp. 19-20

fx-FD10 Pro. User's Guide. Casio. Tokyo, Japan. 2014. pg. α-16 (alpha-16)

Eddie

All original content copyright, © 2011-2026. 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.

HP 32SII: Volume of a Square Frustum (Three Approaches in RPN)

HP 32SII: Volume of a Square Frustum (Three Approaches in RPN)

The volume of a square frustum is:

v = height ÷ 3 × (s1^2 + s2^2 + √(s1^2 × s2^2))

which is based off the formula for the volume of general frustum:

v = height ÷ 3 × (base_area_1 + base_area_2 + √(base_area_1 × base_area_2))

Here are three approaches, which were all done with a HP 32SII calculator. The beauty of RPN is on display here. The three programs are labeled A, C, and D.

s1, s2: side of the square base h: height

Stack: Z: s1 Y: s2 X: h

Program A: (24 steps, 36 bytes) 

V = h/3 * (s1^2 + s2^2 + √(s1^2 * s2^2))

LBL A

STO H

R↓

x²

x<>y

x²

ENTER

ENTER

CLx

-

R↓

x<>y

R↓

×

√

R↓

+

R↑

+

RCL H

×

3

÷

RTN

Program C: (17 steps, 25.5 bytes)

Note: √(s1^2 * s2^2) = s1 * s2

V = h/3 * (s1^2 + s2^2 + (s1 * s2))

LBL C

ENTER

3

÷

R↓

×

LST x

x<>y

R↓

x²

x<>y

x²

+

R↑

+

×

RTN

Program D: (15 steps, 22.5 bytes) 

Note: s1^2 + s2^2 + s1 * s2 = s1^2 + s2^2 + 2 * s1 * s2 – s1 *s2 = (s1 + s2)^2 - s1*s2 

V = h/3 * ((s1 + s2)^2 - s1 * s2)

LBL D

3

÷

R↓

×

LST x

x<>y

R↓

+

x²

R↑

-

×

RTN

Examples

Example 1:

s1: 24

s2: 30

h: 5.8

Result: 4245.6

Example 2:

s1: 2.54

s2: 12.15

h: 10

Result: ≈ 616.4503

Eddie

All original content copyright, © 2011-2026. 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.