c*******************************************************
c force_electro_init(econst,debye,ddna)
c
      subroutine force_electro_init(econst,debye,ddna)

      implicit double precision(a-h, o-z)
      parameter(ds=0.001,ndneigh=10/ds)
      parameter(maxa=500,maxn=maxa*3,ncmax=5)
       common/fe/neigh(maxa,maxa),fe(ndneigh)

      i0=dsqrt(ddna)/ds

      do i=1,ndneigh
         s=i*ds
         yy=debye*s
         fe(i)=econst*exp(-yy)*(yy+1)/(s*s*s)
c         write(111,*) fe(i)
      enddo
c      write(*,*) 'initialize the fe'
c      pause

      
      return 
      end