c init_linear.f
c ***********************************************************************
c 
      subroutine init_linear(n)
c Given out an initial configuration: a straight line 

      implicit double precision(a-h,o-z)
      parameter(maxa=500,maxn=maxa*3,pi2=2.d0*3.1415926535d0)
      integer n
      common/matrix/x(maxa,3),bangle(maxa),rz(maxa)
        common/matrix2/ r(maxa,maxa)
      common/param/temp,bendc,r0,beadr,ss1,interval,dtime,gamma,
     +amass,am(maxn),coef,dtime2
      	common/bf_coords/f(maxa,3),v(maxa,3),u(maxa,3)
        common/ddx/d_x(maxa,3)
  
      do j=1,3
	x(1,j)=0.d0
      enddo
	
      do i=2, n
	 x(i,1)=x(i-1,1)+r0
         x(i,2)=0.d0
	 x(i,3)=0.d0
      enddo

      call update_bond(n)
      call update_rij(n)

      return
      end