vander.for

来自「Numerical Recipes一书中例子的源码所用到的函数集,William」· FOR 代码 · 共 36 行

      SUBROUTINE vander(x,w,q,n)
      INTEGER n,NMAX
      DOUBLE PRECISION q(n),w(n),x(n)
      PARAMETER (NMAX=100)
      INTEGER i,j,k
      DOUBLE PRECISION b,s,t,xx,c(NMAX)
      if(n.eq.1)then
        w(1)=q(1)
      else
        do 11 i=1,n
          c(i)=0.d0
11      continue
        c(n)=-x(1)
        do 13 i=2,n
          xx=-x(i)
          do 12 j=n+1-i,n-1
            c(j)=c(j)+xx*c(j+1)
12        continue
          c(n)=c(n)+xx
13      continue
        do 15 i=1,n
          xx=x(i)
          t=1.d0
          b=1.d0
          s=q(n)
          do 14 k=n,2,-1
            b=c(k)+xx*b
            s=s+q(k-1)*b
            t=xx*t+b
14        continue
          w(i)=s/t
15      continue
      endif
      return
      END