"CODE" 34606;
"PROCEDURE" HSH3COL(LA,LB,U,I,A1,A2,A3,A,B);
"VALUE"LA,LB,U,I,A1,A2,A3;"INTEGER"LA,LB,I,U;"REAL"A1,A2,A3;"ARRAY"A,B;
"IF" A2^=0 "OR" A3^=0 "THEN"
"BEGIN" "REAL" R,T,C;"ARRAY" V[I:I+2];
    R:="IF" A1<0 "THEN" -SQRT(A1*A1+A2*A2+A3*A3)
    "ELSE" SQRT(A1*A1+A2*A2+A3*A3);
    T:=A1+R;C:=-T/R;V[I]:=1;V[I+1]:=A2/T;V[I+2]:=A3/T;
    HSHVECMAT(I,I+2,LA,U,C,V,A);HSHVECMAT(I,I+2,LB,U,C,V,B)
"END" HSH3COL

        "EOP"