SOLIPSYST                Библиотека "JINRLIB"          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