fredin.for

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

      FUNCTION fredin(x,n,a,b,t,f,w,g,ak)
      INTEGER n
      REAL fredin,a,b,x,f(n),t(n),w(n),g,ak
      EXTERNAL ak,g
CU    USES ak,g
      INTEGER i
      REAL sum
      sum=0.
      do 11 i=1,n
        sum=sum+ak(x,t(i))*w(i)*f(i)
11    continue
      fredin=g(x)+sum
      return
      END