code 34609; procedure HSH3ROW3(L,U,UX,J,A1,A2,A3,A,B,X); value L,U,UX,J,A1,A2,A3;integer L,J,U,UX;real A1,A2,A3;array A,B,X; if A2^=0 or A3^=0 then begin real R,T,C;array V[J:J+2];integer K; 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(L,U,J,J+2,C,V,A);HSHVECTAM(L,U,J,J+2,C,V,B); HSHVECTAM(L,UX,J,J+2,C,V,X) end HSH3ROW3; eop