Fortran 95 source code to compute Napier's constant with FMLIB.

napier.f95

!****************************************************
! ネイピア数（自然対数の底 e）の計算
!
! date          name            version
! 2018.11.17    mk-mode.com     1.00 新規作成
!
! Copyright(C) 2018 mk-mode.com All Rights Reserved.
!****************************************************
!
program napier
  use FMZM
  implicit none
  integer, parameter  :: KETA = 1000  ! 小数点以下の桁数
  integer, parameter  :: N    =  451  ! 計算項数
  type(IM)            :: k            ! 階乗部分
  type(FM)            :: np           ! ネイビア数
  character(16)       :: fmt_s        ! 結果出力フォーマット
  character(KETA + 2) :: s            ! 結果文字列
  integer             :: i            ! ループインデックス

  call FM_SET(KETA + 2)
  k  = TO_IM('1')
  np = TO_FM('1')
  do i = 1, N
    k = k * i
    np = np + 1.0d0 / k
  end do
  !call FMPRINT(np)  ! <= FMLIB 標準の書式で出力
  write (fmt_s, '("F", I0, ".", I0)') KETA + 2, KETA
  call FM_FORM(fmt_s, np, s)
  print '(A)', s

  stop
end program napier