c-------------------------------------------------------------      
      subroutine eigtransf(n,ndim,a,b,d,c)
      implicit none
      integer n,ndim
      double precision a(ndim,n),b(ndim,n),c(ndim,n),d(ndim)
      integer i,j,k
      do i=1,n
        do j=1,n
          c(i,j)=0.0d0
          do k=1,n
            c(i,j)=c(i,j)+a(i,k)*b(k,j)*d(k)
          enddo
        enddo
      enddo
      return
      end