LINUX.ORG.RU

История изменений

Исправление grem, (текущая версия) :

Для n <= 46 (дальше косяки с преобразованиями)

program fibonacci

  implicit none

  integer, parameter :: n = 46
  integer :: i
  real(kind=kind(1.d0)), parameter :: phi = (1.0d0 + sqrt(5.0d0))/2.0d0
  integer, dimension(1:n), parameter :: fib = [ ( (phi**i - (-phi)**(-i))/(2.0d0*phi - 1.0d0), i = 1, n ) ]
  print *, fib

end program fibonbacci

Исходная версия grem, :

Для n <= 46 (дальше косяки с преобразованиями)

program fibonacci

  implicit none

  integer, parameter :: n = 46
  integer :: i
  real(kind=kind(1.d0)), parameter :: phi = (1.0d0 + sqrt(5.0d0))/2.0d0
  integer, dimension(1:n), parameter :: fib = [ ( (phi**i - (-phi)**(-i))/(2.0d0*phi - 1.0d0), i = 1, n ) ]
  print *, fib

end program fobonbacci