Текст подпрограммы и версий
amazr_p.zip , amaze_p.zip 		Тексты тестовых примеров
tamazr_p.zip , tamaze_p.zip

Подпрограмма:  AMAZR (модуль AMAZR_p)

Назначение

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

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

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

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

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

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

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

Unit tamazr_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMAZR_p;

function tamazr: String; 

implementation

function tamazr: String;
var
N,_i :Integer;
C :Array [0..15] of Real;
const
A :Array [0..15] 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,1.0 );
B :Array [0..9] of Real = ( 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0 );
begin
Result := '';
N := 4;
AMAZR(A,N,B,C);
Result := Result + #$0D#$0A;
for _i:=0 to 15 do
 begin
  Result := Result + Format('%20.16f ',[C[_i]]);
  if ( ((_i+1) mod 2)=0 )
   then Result := Result + #$0D#$0A;
 end;
Result := Result + #$0D#$0A;
UtRes('tamazr',Result);  { вывод результатов в файл tamazr.res }
exit;
end;

end.


Результаты:

                | 2.0   2.0   2.0   2.0 |
                | 2.0   2.0   2.0   2.0 |
       C  =  | 2.0   2.0   2.0   2.0 |
                | 2.0   2.0   2.0   2.0 |