SUBROUTINE DMXTXF(M,N,A,LDA,nb,B,LDB)
      IMPLICIT double precision (A-H,O-Z)
      DIMENSION A(LDA,N),B(LDB,N)
      DO 20 J=1,N
      DO 20 I=j,N
      AA=0.0D0
      DO 10 L=1,M
   10 AA=AA+A(L,I)*A(L,J)
   20 B(j,i)=AA
      call dfull1(n,b,ldb)
      RETURN
      END