Функции Фортрана, возвращающие неожиданные типы и значения

Я работаю над проектом, в котором необходимо реализовать несколько числовых методов на Фортране. Для этого мне нужно написать несколько рекурсивных функций. Вот мой код.

!     
! File:   main.F95
!

RECURSIVE FUNCTION integrate(n) RESULT(rv)
    IMPLICIT NONE
    DOUBLE PRECISION :: rv
    INTEGER, INTENT(IN) :: n
    DOUBLE PRECISION, PARAMETER :: minusone = -1.0
    IF (n == 1) THEN
        rv = 10 !exp(minusone)
        RETURN
    ELSE
        rv = 1 - (n * integrate(n - 1))
        RETURN
    END IF
END FUNCTION integrate

RECURSIVE FUNCTION factorial(n) RESULT(res)
    INTEGER res, n
    IF (n.EQ. 0) THEN
        res = 1
    ELSE
        res = n * factorial(n - 1)
    END IF
END

PROGRAM main
    DOUBLE PRECISION :: rv1
    PRINT *, factorial(5)
    PRINT *, integrate(2)

    !READ *, rv1

END PROGRAM main

Для этой программы вывод:

         NaN
       1

Если я изменю порядок операторов печати (, строки 30 и 31 ), вывод будет:

         1
-19.000000

Вывод должен быть (для исходного порядка оператора печати ):

  120  
  -19 

. Я взял факториальную функцию со страницы Википедии Особенности языка Fortran 95 .

  • Компилятор :gfortran 4.5.3 с Cygwin
  • IDE :Netbeans 7.0.1
  • Платформа :Windows 7

6
задан francescalus 25 October 2019 в 09:25
поделиться