The Museum of HP Calculators

3 x 3 Matrix Operations for the HP-95C

This program is Copyright © 1977 by Hewlett-Packard Company and is used here by permission. It was originally printed in the HP-95C Applications book. This program was entered by André Wilhelmus and the curator added comments from the HP manual.

This program is supplied without representation or warranty of any kind. André Wilhelmus, Hewlett-Packard and The Museum of HP Calculators therefore assume no responsibility and shall have no liability, consequential or otherwise, of any kind arising from the use of this program material or any part thereof.

Program Listing

CODE KEYS
A-000    lbl A f LBL A
A-001    35  2 STO 2 Input
A-002       12 R↓
A-003    35  1 STO 1
A-004       12 R↓
A-005    35 .1 STO . 1
A-006       64 R/S
A-007    35  4 STO 4
A-008       12 R↓
A-009    35  3 STO 3
A-010       12 R↓
A-011    35 .2 STO . 2
A-012       64 R/S
A-013    35  6 STO 6
A-014       12 R↓
A-015    35  5 STO 5
A-016       12 R↓
A-017    35 .3 STO . 3 - - - - - - - - - - - - - - - - - -
A-018        1 1
A-019    35 .0 STO . 0
A-020    45 .1 RCL . 1
A-021       22 CHS mi1 ← − ai/a1, i = 2,3
A-022 35 24  1 STO ÷ 1
A-023 35 24  2 STO ÷ 2
A-024    45  1 RCL 1
A-025    45 .2 RCL . 2
A-026       39 ×
A-027 35 59  3 STO + 3 b2 ← b2 + m21b1
A-028    45  1 RCL 1
A-029    45 .3 RCL . 3
A-030       39 ×
A-031 35 59  5 STO + 5 c2 ← c2 + m21c1
A-032    45  2 RCL 2
A-033    45 .2 RCL . 2
A-034       39 ×
A-035 35 59  4 STO + 4 b3 ← b3 + m31b1
A-036    45  2 RCL 2
A-037    45 .3 RCL . 3
A-038       39 ×
A-039 35 59  6 STO + 6 c3 ← c3 + m31c1
A-040    45  3 RCL 3
A-041    62 42 g x≠0 - - - - - - - - - - - - - - - - - -
A-042    63  0 GTO 0 If b2 = 0, swap:
A-043        1 1
A-044       22 CHS
A-045    35 .0 STO . 0 b2 ⇄ b3, c2 ⇄ c3
A-046    45  5 RCL 5
A-047    45  6 RCL 6 Set pivot to -1.
A-048    35  5 STO 5
A-049       12 R↓
A-050    35  6 STO 6
A-051    45  4 RCL 4
A-052    45  3 RCL 3
A-053    35  4 STO 4
A-054       12 R↓
A-055    35  3 STO 3 - - - - - - - - - - - - - - - - - -
A-056 61 63  0 f LBL 0
A-057       22 CHS
A-058 35 24  4 STO ÷ 4 m32 ← − b3/b2
A-059    45  4 RCL 4
A-060    45  5 RCL 5
A-061       39 ×
A-062 35 59  6 STO + 6
A-063    45 .1 RCL . 1 c3 ← c3 + m32c2
A-064    45  3 RCL 3
A-065       39 ×
A-066    45  6 RCL 6 Determinant.
A-067       39 ×
A-068    45 .0 RCL . 0
A-069       39 ×
b-000    lbl b f LBL B - - - - - - - - - - - - - - - - - -
b-001        1 1 Inverse.
b-002       21 ENTER ↑
b-003        0 0
b-004       21 ENTER ↑
b-005        C C
b-006        0 0
b-007       21 ENTER ↑
b-008        1 1
b-009       21 ENTER ↑
b-010        0 0
b-011        C C
b-012        0 0
b-013       21 ENTER ↑
b-014       21 ENTER ↑
b-015        1 1
b-016        C C
C-000    lbl C f LBL C - - - - - - - - - - - - - - - - - -
C-001    35  9 STO 9 Solve Ax = d.
C-002       12 R↓
C-003    35  8 STO 8
C-004       12 R↓
C-005    35  7 STO 7
C-006    45  1 RCL 1
C-007    45  7 RCL 7 Ld
C-008       39 ×
C-009 35 59  8 STO + 8 d2 ← d2 + m21d1
C-010    45  2 RCL 2
C-011    45  7 RCL 7
C-012       39 ×
C-013 35 59  9 STO + 9 d3 ← d3 + m31d1
C-014    45 .0 RCL . 0
C-015    62 44 g x>0
C-016    63  0 GTO 0 If Pivot = -1, swap:
C-017    45  8 RCL 8
C-018    45  9 RCL 9 d2 ⇄ d3
C-019    35  8 STO 8
C-020       12 R↓
C-021    35  9 STO 9
C-022 61 63  0 f LBL 0 - - - - - - - - - - - - - - - - - -
C-023    45  4 RCL 4
C-024    45  8 RCL 8 d3 ← d3 + m32d2
C-025       39 ×
C-026 35 59  9 STO + 9
C-027    45  6 RCL 6
C-028 35 24  9 STO ÷ 9 z = d3/c3
C-029    45  9 RCL 9
C-030    45  5 RCL 5
C-031       39 ×
C-032 35 49  8 STO − 8
C-033    45  3 RCL 3
C-034 35 24  8 STO ÷ 8 y = (d2 - c2 z)/b2
C-035    45  9 RCL 9
C-036    45 .3 RCL . 3
C-037       39 ×
C-038    45  8 RCL 8 x = (d1 - b1 y − c1 z)/a1
C-039    45 .2 RCL . 2
C-040       39 ×
C-041       59 +
C-042 35 49  7 STO − 7
C-043    45 .1 RCL . 1
C-044 35 24  7 STO ÷ 7
C-045    45  7 RCL 7
C-046       14 PRINT x Ouput x, y, z.
C-047    45  8 RCL 8
C-048       14 PRINT x
C-049    45  9 RCL 9
C-050       14 PRINT x
C-051    61 14 f SPACE

Go back to the HP-95C software library
Go back to the general software library
Go back to the main exhibit hall