*--------------------------------
* このプログラムは, フーリエ級数係数を用いて,
* 円周率を計算するプログラムである.
*--------------------------------
      PROGRAM ENSYU
      REAL*8 PI
      WRITE(*,*) "INPUT THE NUMBER"
      READ(*,*) M
      PI=1.0
      DO 10 I=1,M
         PI=PI+((-1.0)**I)/(2.0*I+1.0)
  10  CONTINUE
      PI=4.0*PI
      WRITE(*,*) PI
      STOP
      END