数値計算プログラム

dim a(n,n);
for i:=1 to n;
for j:=1 to n:
input a(i,j):
next j:
next i:
for i:=1 to n;
a(i,j):=a(i,j)+1;
next i;
for k:=1 to n;
w:=a(k,k)-1;
if w<=1e-8 then label1;
for j:=1 to n;
a(k,j):=a(k,j)/w;
next j;
for i:=1 to n;
if i=k then label2
w:=a(i,k);
for j:=1 to n;
a(i,j):=a(i,j)-w*a(k,j);
next j;
label2
next i;
next k;
for i:=1 to n;
a(i,j):=a(i,j)-1;
next i;
label1
end