code 34605;
procedure HSH2COL(LA,LB,U,I,A1,A2,A,B);value LA,LB,U,I,A1,A2;
integer LA,LB,U,I;real A1,A2;array A,B;
if A2^=0 then 
begin real R,T,C;array V[I:I+1];
    R:=if A1<0 then -SQRT(A1*A1+A2*A2) else SQRT(A1*A1+A2*A2);
    T:=A1+R;C:=-T/R;V[I]:=1;V[I+1]:=A2/T;
    HSHVECMAT(I,I+1,LA,U,C,V,A);HSHVECMAT(I,I+1,LB,U,C,V,B)
end HSH2COL;
        eop