code 34610; procedure HSH3ROW2(LA,LB,U,J,A1,A2,A3,A,B); value LA,LB,U,J,A1,A2,A3;integer LA,LB,U,J;real A1,A2,A3;array A,B; if A2^=0 or A3^=0 then begin real R,T,C;array V[J:J+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[J+2]:=1;V[J+1]:=A2/T;V[J]:=A3/T; HSHVECTAM(LA,U,J,J+2,C,V,A);HSHVECTAM(LB,U,J,J+2,C,V,B) end HSH3ROW2; eop