clear all;
clc;
a=[27 6 -1;6 15 2 ; 1 1 54];
x=[0;0;0];
b=[85;72;110];
n=input ('n=');
for k=1:10
for i=1:n
s=0;
for j=1:n
if (i~=j)
s=s+x(j)*a(i,j);
end
end
x(i)=(b(i)-s)/a(i,i);
end
disp(x');
end