as.f90

Visual Fortran常用数值算法集

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