Определитель матрицы, Fortran

[frad]

Определитель матрицы, Fortran

Помогите пожалуйста. Нужен код на фортране для нахождения определителя матрицы любого размера (чтение из текстового файла) по формуле Лапласа.
Спасибо!

 

[Gruvi]

Вот тебе программа. (и прикрепил тебе фаил еще)

REAL FUNCTION FindDet(matrix, n)
IMPLICIT NONE
REAL, DIMENSION(n,n) :: matrix
INTEGER, INTENT(IN) :: n
REAL :: m, temp
INTEGER :: i, j, k, l
LOGICAL :: DetExists = .TRUE.
l = 1
!Convert to upper triangular form
DO k = 1, n-1
IF (matrix(k,k) == 0) THEN
DetExists = .FALSE.
DO i = k+1, n
IF (matrix(i,k) /= 0) THEN
DO j = 1, n
temp = matrix(i,j)
matrix(i,j)= matrix(k,j)
matrix(k,j) = temp
END DO
DetExists = .TRUE.
l=-l
EXIT
ENDIF
END DO
IF (DetExists .EQV. .FALSE.) THEN
FindDet = 0
return
END IF
ENDIF
DO j = k+1, n
m = matrix(j,k)/matrix(k,k)
DO i = k+1, n
matrix(j,i) = matrix(j,i) - m*matrix(k,i)
END DO
END DO
END DO

!Calculate determinant by finding product of diagonal elements
FindDet = l
DO i = 1, n
FindDet = FindDet * matrix(i,i)
END DO

END FUNCTION FindDet

Это подпрограмма-функция, которая на входе получает квадратную матицу размерности n и считает ее дискриминант.

 

[frad]

Спасибо, но это не то! Мне нужно по формуле (разложению) Лапласа (в этом то вся соль), а это методом Гаусса.

 

[Gruvi]

Вот тебе Теорема Лапласа =) смотри и учись. А по факту это - работа с двумерным массивом.

 

[frad]

Спасибо конечно, но эту теорему я знаю еще с 1 курса, а вот как ее осуществить - может какой блок-схемой, хотя бы поможете (