Текст подпрограммы и версий
ama2r_p.zip , ama2e_p.zip 		Тексты тестовых примеров
tama2r_p.zip , tama2e_p.zip

Подпрограмма:  AMA2R (модуль AMA2R_p)

Назначение

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

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

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

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

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

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

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

Unit TAMA2R_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMA2R_p;

function TAMA2R: String; 

implementation

function TAMA2R: String;
var
N,NC,_i :Integer;
B :Array [0..7] of Real;
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 );
begin
Result := '';
{ прототип оператора DАТА на FORTRANе }
B[0] := 0.0;   
B[1] := 1.0;   
B[2] := 1.0;   
B[3] := 1.0;   
B[4] := 1.0;   
B[5] := 1.0;   
B[6] := 1.0;   
B[7] := 1.0;   

N := 4;
NC := 1;
AMA2R(A,N,B,NC,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('TAMA2R',Result);  { вывод результатов в файл TAMA2R.res }
exit;
end;

end.

Результаты:

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