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