input n;
dim a(n,n+1),x(n);
input e;
for i:=1 to n;
for j:=1 to n;
input a(i,j);
next j;
next i;
for i:=1 to n;
input a(i,i+1);
next i;
n1:=n+1;
n2:=n-1;
for k:=1 to n2;
a1:=a(k,k);
for j:=k to n1;
a(j,j):=a(k,j)/a1;
next j;
for i:=k+1 to n;
a2:=a(i,k);
for j:=k to n1;
a(i,j):=a(i,j)-a2*a(k,j);
next j;
next i;
next k;
if abs(a(n,n1))<=e then l1;
a3:=a(n,n);
a(n,n)=a(n,n)/a3;
a(n,n1):=a(n,n1)/a3;
x(n):=a(n,n1);
for i:=1 to n2;
k:=n-i;
x(k):=a(k,n1);
for j:=k to n2;
x(k):=x(k)-a(k,j+1)*x(j+1);
next j;
next i;
for i:=1 to n;
print x(i);
next i;
end;
label l1;
print "error";
end;