|
COMPLEX FUNCTION FILON(CFUN,CEX,ARG,N,FAC)
COMPLEX AI,CSUM,CFUN(1),CEX(1),DFUN,DEX,DARG,DPROD
DIMENSION ARG(1)
AI=CMPLX(0.,1E0)
CSUM=0E0
DO 10 I=1,N-1
DARG=-AI/(ARG(I+1)-ARG(I))
DFUN=CFUN(I+1)-CFUN(I)
DEX=CEX(I+1)-CEX(I)
DPROD=CFUN(I+1)*CEX(I+1)-CFUN(I)*CEX(I)
CSUM=CSUM+DARG*(DPROD-DARG*DFUN*DEX)
10 CONTINUE
FILON=FAC*CSUM
RETURN
END
请问上述子函数中,变量CFUN,CEX以及ARG的数组大小为N,为什么声明的时候用CFUN(1),CEX(1)和ARG(1)呢?
|
|