c****************************************************
c update_rij(n)
c update r(i,j) into common/matrix2/ r(i,j)
      subroutine update_rij(n)
   
       implicit double precision(a-h, o-z)
       parameter(maxa=500,maxn=maxa*3)   
       integer n
       double precision rr,rr2
       common/matrix/x(maxa,3),bangle(maxa),rz(maxa)
       common/matrix2/ r(maxa,maxa)
 
       do i=1,n-1
          do j=i+1,n
             rr=0.0
             do k=1,3
             rr=(x(i,k)-x(j,k))**2+rr
             enddo
             rr2=dsqrt(rr)
             r(i,j)=rr2
             r(j,i)=rr2
           enddo
        enddo

        return 
        end