5   REM LEAST SQUARES METHOD WITH DIB
6   REM EXAMPLE TAKEN FROM THE www.hpmuseum.org FORUM
10  READ M
20  READ N
30  DIM A(M,N), C(M), B(N,N), X(N), Z(N,N), T(N,M), D(N)
40  PRINT"READING MAT A (COEFFICIENT)" 
50  MAT READ A
60  MAT PRINT A
70  PRINT 
80  PRINT"READING MAT C (KNOWNS)"
90  MAT READ C
100 MAT PRINT C
110 PRINT
120 PRINT
130 PRINT"CALCULATING..."
140 MAT T=TRN(A)
150 MAT B=T*A
160 PRINT"NEW COEFFICIENT MATRIX B"
170 MAT PRINT B
180 PRINT
190 MAT D = T*C
200 PRINT"NEW KNOWNS MATRIX D"
210 MAT PRINT D
220 PRINT
230 MAT Z=INV(B)
240 MAT X=Z*D
250 PRINT"RESULTS"
260 MAT PRINT X
400 DATA 9,3,.95,1.09,1.01,1.09,1.94,1.95,2.02,1.91,2.92,3,2.93,2.98
410 DATA 4.01,3.95,3.97,3.94,1,0,0,-1,1,0,-1,0,1,-1,0,0,0,1,0,-1,0,1
420 DATA -1,-0,0,0,1,-1,1,2,3,4,0,0,0,0,0,0
500 END


