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

frad · 20.01.2013, 22:37

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

Gruvi · 20.01.2013, 22:51

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

Цитата:

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 · 20.01.2013, 23:06

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

Gruvi · 20.01.2013, 23:53

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

frad · 21.01.2013, 09:44

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

20.01.2013, 22:37	#1 (permalink)
frad Новичок Регистрация: 20.01.2013 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Определитель матрицы, Fortran Помогите пожалуйста. Нужен код на фортране для нахождения определителя матрицы любого размера (чтение из текстового файла) по формуле Лапласа. Спасибо!

20.01.2013, 22:37
Helpmaster Member Регистрация: 08.03.2016 Сообщений: 0	Вам не помешает прочесть схожие по содержанию темы Матрицы С++ Пожалуйста, помогите решить простую задачу по информатике в FORTRAN про массив Не работает Сompaq visual fortran 6.6 на windows 7 х32

20.01.2013, 23:06	#3 (permalink)
frad Новичок Регистрация: 20.01.2013 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Спасибо, но это не то! Мне нужно по формуле (разложению) Лапласа (в этом то вся соль), а это методом Гаусса.

20.01.2013, 23:53	#4 (permalink)
Gruvi VIP user Регистрация: 10.03.2011 Сообщений: 765 Записей в дневнике: 1 Сказал(а) спасибо: 10 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 3453	Вот тебе Теорема Лапласа =) смотри и учись. А по факту это - работа с двумерным массивом.

21.01.2013, 09:44	#5 (permalink)
frad Новичок Регистрация: 20.01.2013 Сообщений: 3 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Спасибо конечно, но эту теорему я знаю еще с 1 курса, а вот как ее осуществить - может какой блок-схемой, хотя бы поможете (

Ads
Яндекс Member Регистрация: 31.10.2006 Сообщений: 40200 Записей в дневнике: 0 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 55070