c-------------------------------------------------------------      
      subroutine eigtransf1(n,ndim,ab,d,c)
      implicit none
      integer n,ndim
      double precision ab(ndim,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)+ab(k,j,i)*d(k)
          enddo
        enddo
      enddo
      return
      end