source/maxlik/src_FPy.org/fibonacci.f

   1 C FILE: FIB1.F
   2       SUBROUTINE FIB(A,N)
   3 C
   4 C     CALCULATE FIRST N FIBONACCI NUMBERS
   5 C
   6       INTEGER N
   7       REAL*8 A(N)
   8       DO I=1,N
   9          IF (I.EQ.1) THEN
  10             A(I) = 0.0D0
  11          ELSEIF (I.EQ.2) THEN
  12             A(I) = 1.0D0
  13          ELSE
  14             A(I) = A(I-1) + A(I-2)
  15          ENDIF
  16       ENDDO
  17       END
  18 C END FILE FIB1.F