Метод квадратного корня

[kitty_girl]

Разбираю чужой код. суть задачи - решение слау методом квадратного корня для ленточных матриц. в коде существует процедура,которая,насколько я поняла, умножает матрицу саму на себя для избежания выхода в комплексную область при извлечении корней. но матрица перемножается каким-то странным образом. может кто знает,что это за метод? сама процедура ниже

MatrixB = array[1..maxn, 1..maxl] of real;

Procedure TFormMain.MultiB (n,l: integer; Var B:MatrixB);
var i,j,k,u:integer;
temp: MatrixB;
begin
for i:=1 to maxn do
for j:=1 to maxl do
temp[i,j] := 0;

for i:=1 to n do
for j:=1 to l do
temp[i,1] := temp[i,1] + B[i,j]*B[i,j];
for i:=1 to n do
for j:=2 to l do
begin
u := 1;
for k:=j to l do
begin
temp[i,j] := temp[i,j] + B[i,k]*B[j+i-1,u];
inc(u);
end;
end;
for i:=1 to n do
for j:=1 to l do
B[i,j] := temp[i,j]
end;

[Abstraction]

Что-то странное. Первый столбец - суммы квадратов строк, последующие - суммы попарных произведений, но... Из интереса скармливаю матрицу 3х3, считаю. ~ - неопределённые значения, B[j+i-1,u] лезет за пределы массива.

  1 0 3        10 0 0        10  3 12
B 3 1 2  temp1 14 0 0  temp2 14 14  ~
  4 5 1        42 0 0        42  ~  ~

Оно правда работает?

[kitty_girl]

да))
дело в том,что фактический размер матрицы - maxN*maxL,но ненулевые элементы расположены в области N*L,т.е. выхода за границу нет,там ноль стоит