История изменений
Исправление 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