Текст подпрограммы и версий
amm5r_p.zip , amm5e_p.zip 		Тексты тестовых примеров
tamm5r_p.zip , tamm5e_p.zip

Подпрограмма:  AMM5R (модуль AMM5R_p)

Назначение

Умножение симметричных матриц, заданных в компактной форме.

Математическое описание

AMM5R вычисляет произведение симметричных матриц А и B порядка N, заданных в компактной форме. Результатом этого произведения является квадратная матрица C порядка N.

Использование 

procedure AMM5R(var A :Array of Real; var B :Array of Real;
                N :Integer; var C :Array of Real); 

   Параметры
A, B - вещественные векторы длины N * (N + 1)/2, в которых задаются симметричные матрицы А и B в компактной форме;
N - заданный порядок матриц А и B (тип: целый);
C - вещественный двумерный массив размера N*N, в котором запоминается матрица произведения. 
   Версии
AMM5E - умножение с расширенной (Extended) точностью симметричных матриц, заданных в компактной форме. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
  В подпрограмме АММ5R параметры А, B и C должны иметь тип Extended.

Пример использования 

Unit tamm5r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMM5R_p;

function tamm5r: String; 

implementation

function tamm5r: String;
var
N,_i :Integer;
C :Array [0..24] of Real;
const
A :Array [0..14] of Real = ( 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,
1.0,1.0,1.0 );
B :Array [0..14] of Real = ( 2.0,2.0,2.0,2.0,2.0,2.0,2.0,2.0,2.0,2.0,2.0,2.0,
2.0,2.0,2.0 );
begin
Result := '';
N := 5;
AMM5R(A,B,N,C);
Result := Result + #$0D#$0A;
for _i:=0 to 24 do
 begin
  Result := Result + Format('%20.16f ',[C[_i]]);
  if ( ((_i+1) mod 3)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('tamm5r',Result);  { вывод результатов в файл tamm5r.res }
exit;
end;

end.


Результаты:

                 |  10   10   10   10   10  |
                 |  10   10   10   10   10  |
      C  =    |  10   10   10   10   10  |
                 |  10   10   10   10   10  |
                 |  10   10   10   10   10  |