Текст подпрограммы и версий
ama4r_p.zip , ama4e_p.zip 		Тексты тестовых примеров
tama4r_p.zip , tama4e_p.zip

Подпрограмма:  AMA4R (модуль AMA4R_p)

Назначение

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

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

AMA4R вычисляет сумму ленточной матрицы  А порядка  N, заданной в компактной форме, и ленточной симметричной матрицы  B того же порядка, заданной также в компактной форме. Результатом этого сложения является ленточная матрица  C порядка  N с  max (NLА, NС) нижними и  max (NUА, NС) верхними ко - диагоналями. Матрица  C запоминается в компактной форме.

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

procedure AMA4R(var A :Array of Real; N :Integer; NLA :Integer;
                NUA :Integer; var B :Array of Real; NC :Integer;
                var C :Array of Real); 

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

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

Unit tama4r_p;
interface
uses
SysUtils, Math, { Delphi }
LStruct, Lfunc, UtRes_p, AMA4R_p;

function tama4r: String; 

implementation

function tama4r: String;
var
N,NLA,NUA,NC,_i :Integer;
C :Array [0..19] of Real;
const
A :Array [0..19] of Real = ( 0.0,0.0,1.0,1.0,1.0,0.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,0.0 );
B :Array [0..9] of Real = ( 0.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0 );
begin
Result := '';
N := 5;
NLA := 2;
NUA := 1;
NC := 1;
AMA4R(A,N,NLA,NUA,B,NC,C);
Result := Result + #$0D#$0A;
for _i:=0 to 19 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('tama4r',Result);  { вывод результатов в файл tama4r.res }
exit;
end;

end.


Результаты:

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