Расшифровка кода

[drNabla]

Здравствуйте.
Есть код из Делфи, но ни кто не знает по какому алгоритму эта программа работает, а алгоритм этот знать очень хочется, а человек который писал программу и изначально делал алгоритм уже давно расслабляется в америке... так как из моего окружения в делфи ни кто не разбирается то последняя надежда на вашу помощь.
Заранее спасибо за любую помощь.
Пример кода:

procedure faddy(ind:byte);
label 1;
var i,j,k,sup: byte;
k1,l,sgn: integer;
a,b,c:array[1 .. 25,1 .. 25] of real;
p,bnk:extended;
begin
i:=1;
case ind of
0: begin dim:= 4; sup:=tz; end;
{1,2: begin dim:=4; sup:=tz div 2; end;}
1,2: begin dim:=4; sup:=8; end;
end;
for k:=1 to sup do begin
k1:=k+1; a[k,k]:=1.0;
if k = sup then goto 1;
for l:=k1 to sup do begin
a[k,l]:=mss[i]; a[l,k]:=mss[i];
i:=i+1; end;
end;
1: sgn:=1; vol:=0.0; bnk:=1.0;
for i:=1 to sup do begin for j:=1 to sup do begin
b[i,j]:=a[i,j];end;end;
for k:=1 to dim-1 do begin
p:=0.0; bnk:=(bnk*(sup-k+1))/k;
for i:=1 to sup do begin
p:=p+b[i,i]; end;
p:=p/k;
for i:=1 to sup do begin
b[i,i]:=b[i,i]-p; end;
for i:=1 to sup do begin
for j:=1 to sup do begin
c[i,j]:=0.0;
for l:=1 to sup do begin
c[i,j]:=c[i,j]+a[i,l]*b[l,j]; end;
end;
end;
for j:=1 to sup do begin
for i:=1 to sup do begin
b[i,j]:=c[i,j];
end;
end;
sgn:=-sgn; end;
bnk:=bnk*(sup-dim+1);
for i:=1 to sup do begin
vs[i]:=sgn*b[i,i];
vol:=vol+vs[i]; end;
for i:=1 to sup do begin
vs[i]:=100*vs[i]/vol; end;
vol:=vol/bnk;
end;

[wm_leviathan]

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