Пример рекурсии в фортране

! Д. К. Фаддеев, И. С. Соминский.
! Сборник задач по высшей алгебре. 1977.
! № 279. Вычислить определитель.
! Ответ: 900.
 
    integer, parameter :: n=4
    real*8 :: a(n,n) = &
! Обратите внимание, что матрица в фортране записывается по столбцам.
    (/ 1, -2,  3,  4, &
       2,  1, -4,  3, &
       3, -4, -1, -2, &
       4,  3,  2, -1  /)
    integer :: ia(n) = (/ (i, i=1,n) /)
    real*8 det
    print *, det(a, n, ia, n)
    end
    
    recursive real*8 function det(a, n, ia, k)
!    Функция вычисляет минор матрицы a размера n x n, стоящий на пересечении
!    k последних строк и k столбцов, номера которых перечислены в массиве ia.
    integer k,n
    real*8, dimension (n,n) :: a
    integer, dimension(k) :: ia
    integer i, j, sign
    integer, dimension(k-1) :: iar
    real*8 :: s
    if(k==1) then
        det=a(n,ia(1))
        return
    end if
    if(k==n) ia = (/ (i, i=1,n) /)
    s=0.0
    sign = 1
    do i = 1,k
        iar=(/ (ia(j), j=1,i-1), (ia(j), j=i+1,k) /)
        s = s + sign * a(n-k+1,ia(i)) * det(a, n, iar, k-1)
        sign = -sign
    end do
    det = s
    end function det