distcal.f

Multi Dimensional Scaling

      subroutine distcal(rval,dist,notus1,nd1)
c
      implicit none
c
      integer notus1,nd1
c
      integer iotus0,iotus1,id
c
      real dist(0:notus1,0:notus1),rval(0:nd1,0:notus1)
c
      do 100 iotus0=0,notus1
      do 100 iotus1=0,notus1
        dist(iotus0,iotus1)=0
        do 110 id=0,nd1
          dist(iotus0,iotus1)=dist(iotus0,iotus1)+
     &      (rval(id,iotus0)-rval(id,iotus1))**2
 110    continue
c
        dist(iotus0,iotus1)=dist(iotus0,iotus1)**0.5
c---
cc        write(*,*) iotus0,iotus1,dist(iotus0,iotus1)
cc        read(*,*)
c---
 100    continue
c
      return
c
      end