Appendix
201
Irig1 REPEAT FOR GOING FROM 1 TO n
icol1 (irig1) = irig1
END REPEAT irig1
Irig1 REPEAT FOR GOING FROM 1 TO n-1
amax = vara (irig1, irig1)
imax = irig1
REPEAT FOR ICOL GOING TO BE A n +1 irig1
IF (ABS (vara (irig1, ICOL)). Gt.ABS (amax)) THEN
amax = vara (irig1, ICOL)
imax = ICOL
END IF
END REPEAT ICOL
VNORM = amax
icol1 (imax) = icol1 (irig1)
icol1 (irig1) = imax
vara (irig1, imax) = vara (irig1, icol1 (imax))
vara (irig1, icol1 (imax)) = amax
VNORM = vara (irig1, irig1)
varb (irig1) = varb (irig1) / VNORM
REPEAT FOR ICOL GOING TO BE A n irig1
vara (irig1, ICOL) = vara (irig1, ICOL) / VNORM
END REPEAT ICOL
REPEAT FOR irig2 GOING TO BE A n +1 irig1
varb (irig2) = varb (irig2) / VNORM-varb (irig1)
REPEAT FOR ICOL GOING TO BE A n irig2
vara(irig2,icol1(ICOL))=vara(irig2,icol1(ICOL))/VNORM
-vara(irig1,icol1(ICOL))
END REPEAT ICOL
END REPEAT irig2
END REPEAT irig1
REPEAT FOR irig1 GOING TO PASS n 1 -1
sum = varb (irig1)
SE (irig1.ne.n) THEN
REPEAT FOR ICOL GOING TO BE A n +1 irig1 OF STEP -1
sum = sum - vara (irig1, icol1 (ICOL))
END REPEAT ICOL
END IF
varX (irig1) = sum
END REPEAT irig1