Текст подпрограммы и версий
ama3r_p.zip , ama3e_p.zip , ama3c_p.zip 		Тексты тестовых примеров
tama3r_p.zip , tama3e_p.zip , tama3c_p.zip

Подпрограмма:  AMA3R (модуль AMA3R_p)

Назначение

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

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

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

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

procedure AMA3R(var A :Array of Real; N :Integer;
                var B :Array of Real; NLB :Integer; NUB :Integer;
                var C :Array of Real); 

   Параметры
A - вещественный двумерный массив размера N на N, в котором задается матрица  A;
N - заданный порядок матриц  А и  B (тип: целый);
B - вещественный двумерный массив размера N на (NLВ + NUВ + 1), в котором в компактной форме задается ленточная матрица  B;
      NLB -       NUB   заданное число нижних и верхних ко - диагоналей матрицы  B соответственно (тип: целый);
C - вещественный двумерный массив порядка  N, в котором запоминается результирующая матрица сложения. 
   Версии
AMA3E - сложение с расширенной (Extended) точностью квадратной матрицы с ленточной матрицей, заданной в компактной форме.
AMA3C - сложение комплексной квадратной матрицы с ленточной комплексной матрицей, заданной в компактной форме. 
   Вызываемые подпрограммы:  нет. 
   Замечания по использованию
 

В подпрограмме АМА3E параметры А, B и C должны иметь тип Extended.

В подпрограмме AMA3C параметры А, B и C должны иметь тип Complex.

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

Unit TAMA3R_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMA3R_p;

function TAMA3R: String; 

implementation

function TAMA3R: String;
var
N,NLB,NUB,_i :Integer;
B :Array [0..15] 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] := 0.0;
B[2] := 1.0;
B[3] := 1.0;
B[4] := 0.0;
B[5] := 1.0;
B[6] := 1.0;
B[7] := 1.0;
B[8] := 1.0;
B[9] := 1.0;
B[10] := 1.0;
B[11] := 1.0;
B[12] := 1.0;
B[13] := 1.0;
B[14] := 1.0;
B[15] := 0.0;

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

end.

Результаты:

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