SUBROUTINE CROSS(A,B,AXB) DIMENSION A(3),B(3),AXB(3) DO 1 I=1,3 J=I+1 IF(J.EQ.4) J=1 K=J+1 IF(K.EQ.4) K=1 1 AXB(I)=A(J)*B(K) - B(J)*A(K) RETURN END