HP HP-15C - Storing and Recalling Matrix Elements

Section 4: Using Matrix Operations 1 41
^11 "^IS *" *^lo
X21 X22 ••• X2p
xnl Xn2 ••• xnp
and y =
V
y-
yn
The program does this for successively increasing values of n,
although the solution b = b(ra) is meaningful only when n ^p.
It is possible to factor the augmented n X (p + 1) matrix [X y] into
QTV, where Q is an orthogonal matrix,
V =
U g
0 q
0 0
(p rows)
(1 row),
(n—p — l rows)
and U is an upper-triangular matrix. If this factorization results
from including nrows rm = (xml, xm2,..., xmp, ym) form = 1, 2,..., n
in [X y], consider how to advance to n + 1 rows by appending row
rn + 1to[X y]:
X y
Q7 0
0 1
The zero rows of V are discarded.
Multiply the (p + 2) X (p + 1) matrix
V
U g
0 q
(p rows)
(1 row)
(1 row)