F098 JINRLIB

ОБЪЕДИНЕННЫЙ ИНСТИТУТ ЯДЕРНЫХ ИССЛЕДОВАНИЙ

Библиотека программ JINRLIB

SOLIPSYST - решение системы линейных уравнений методом Жордана

F098

Автор: А.П.Сапожников

Язык: Фортран

Функция SOLIPSYST решает систему линейных уравнений A * X = B общего вида методом Жордана с выбором главного элемента в столбце.

Структура:

Тип:	-	FUNCTION
Имена входа для пользователя:	-	SOLIPSYST

Обращение:

DETERM = SOLIPSYST (A,N,M)

На входе:

A	-	(REAL*8) матрица размерностью (N,N+M+1), содержащая в первых N столбцах коэффициенты А системы, в следующих М столбцах - М комплектов правых частей В, последний столбец - рабочий;
N	-	(INTEGER) размерность системы (2,3,...);
M	-	(INTEGER) число различных вариантов правых частей (0,1,2,...).

На выходе:

Первые N столбцов матрицы А содержат матрицу, обратную исходной,
следующие М столбцов - М комплектов векторов Х, решений исходной системы.
При М=0 производится только обращение матрицы.

Результат, выдаваемый функцией, равен детерминанту системы или 0, если исходная матрица сингулярна.

Пример:

       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION A(3,6)
       DATA A/1.D0,4.D0,0.D0, 2.D0,5.D0,7.D0, 3.D0,6.D0,8.D0,
     *     6.D0,15.D0,15.D0, 12.D0,30.D0,30.D0, 0.D0,0.D0,0.D0/
       WRITE(*,*) ' DETERM=',SOLIPSYST(A,3,2)
       DO J=4,5
         WRITE(*,*) (A(I,J),I=1,3)
       ENDDO
       . . .

Результат:

       DETERM=      18.000000000000000
          1.000000000000000        1.000000000000000        1.000000000000000
          2.000000000000000        2.000000000000000        2.000000000000000