数値計算プログラム

input n,e.j;
l1=0;
dim a(n,n+1),b(n),c(n);
for i:=1 to n;
c(i):=o;
for j:=1 to n+1;
input a(i,j);
next j;
next i;

for i:=1 to n;
label l8;
w:=a(i,i);
if w:=0 then l1;
for j:=1 to n+1;
a(i,j):=a(i,j)/w;
next j;
next i;

label l5;
for i:=1 to n;
b(i):=0;
for j:=1 to n;
if i=j then l2;
b(i):=b(i)-c(j)*a(i,j);
label l2;
next j;
b(i):=b(i)+a(i,j);
next i;
l1:=l1+1;

if l1=l then l3;
next i;
for i:=1 to n;
if abs(b(i)-c(i))>e then l4;
next i;

label l3;
for i:=1 to n;
print b(i);
next i;
end;

label l4;
for j:=1 to n;
c(j):=b(j);
next j;
goto l5;

label l1;
if i=n then l6;
for k:=i+1 to n;
if a(h,h)<>0 then i,h:goto l7;
next h;:goto l6;

label l7;
for l:=1 to n;
v:=a(i,l);
a(i,l):=a(h,l);
a(h,l):=v;
next l;
goto l8;

label l7;
print "funou";
end;