subroutine copy(a,b,n)
      implicit none
      integer n,i
      real*8 a(n),b(n)
      do i=1,n
       b(i)=a(i)
      enddo
      return
      end