crank.for

Numerical Recipes一书中例子的源码所用到的函数集,William H. Press 和 Saul A. Teukolsky 所著

      SUBROUTINE crank(n,w,s)
      INTEGER n
      REAL s,w(n)
      INTEGER j,ji,jt
      REAL rank,t
      s=0.
      j=1
1     if(j.lt.n)then
        if(w(j+1).ne.w(j))then
          w(j)=j
          j=j+1
        else
          do 11 jt=j+1,n
            if(w(jt).ne.w(j))goto 2
11        continue
          jt=n+1
2         rank=0.5*(j+jt-1)
          do 12 ji=j,jt-1
            w(ji)=rank
12        continue
          t=jt-j
          s=s+t**3-t
          j=jt
        endif
      goto 1
      endif
      if(j.eq.n)w(n)=n
      return
      END