code 36402;
procedure ORTHOG(N,LC,UC,X);value N,LC,UC;
integer N,LC,UC;array X;
begin integer I,J,K; real NORMX;
    for J:=LC step 1 until UC do 
    begin NORMX:=SQRT(TAMMAT(1,N,J,J,X,X));
        for I:=1 step 1 until N do X[I,J]:=X[I,J]/NORMX;
        for K:=J+1 step 1 until UC do 
        ELMCOL(1,N,K,J,X,X,-TAMMAT(1,N,K,J,X,X))
    end 
end ORTHOG;
        eop