as.f90

来自「Visual Fortran常用数值算法集」· F90 代码 · 共 13 行

FUNCTION as(N,B)
DIMENSION A(15)
A(1)=EXP(-B)/B
IF (N==1) THEN
  AS=A(N)
ELSE
  DB=1./B
  DO I=2,N
    A(I)=A(1)+DB*(I-1)*A(I-1)
  END DO
  AS=A(N)
END IF
END FUNCTION as