ОБЪЕДИНЕННЫЙ   ИНСТИТУТ   ЯДЕРНЫХ   ИССЛЕДОВАНИЙ
lit
  Библиотека   программ   JINRLIB eng

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


home up e-mail