SUBROUTINE XYT (N1,N2,N3,L,M,N,X,Y,Z)
      DOUBLE PRECISION   X(N1,M)
      DOUBLE PRECISION   Y(N2,M)
      DOUBLE PRECISION   Z(N3,N)
      DOUBLE PRECISION   TEMP
C
C  DESCRIPTION:
C
C  FORM MATRIX PRODUCT Z = X * Y-TRANSPOSE
C
C  ON INPUT:
C  X IS L X M MATRIX IN N1 X M ARRAY
C  Y IS N X M MATRIX IN N2 X M ARRAY
C
C  ON OUTPUT:
C  Z IS L X N MATRIX IN N3 X N ARRAY
C
C  WRITTEN BY:  J. DOUGLAS BIRDWELL
C               DEPT. OF ELEC. ENGR.
C               FERRIS 311
C               UNIVERSITY OF TENNESSEE
C               KNOXVILLE, TN  37996-2100
C               (615) 974-5468
C
      DO 1000 J = 1, N
         DO 1000 I = 1, L
            TEMP = 0.D0
            DO 500 K = 1, M
500            TEMP = TEMP + X(I,K) * Y(J,K)
1000        Z(I,J) = TEMP
      RETURN
      END