/* ASG1605 */ output file = ASG1605.OUT reset; format 8,4;
let x1 = 1 1 1 1 1; let x2 = 2 4 3 5 2; let y = 14 17 8 16 3;
X = x1~x2; Q = X'X; dq = det(Q); QI = invpd(Q);
A = QI*X'; N = X*A; I = eye(5); M = I - N;
trn = sumc(diag(N)); trm = sumc(diag(M)); b = A*y; yh = N*y; e = M*y;
"Q =     "         Q; ?; "det Q  =   "       dq; ?;
"Q inverse =   "  QI; ?;
"N =    "          N; ?; "  M    =   "        M; ?;
"tr(N) =    "    trn;    "tr(M)  =   "   trm; ?;
"b' =    "        b'; ?;
"yhat' = "       yh'; ?;"  e'  = "       e'; end;