%NOT FINISHED YET!

A=[3 -0.1 -0.2;0.1 7 -0.3;0.3 -0.2 10];
B=[7.85 -19.3 71.4]';

%initial guess
x=[1 1 1]';
n=length(B);
tol=0.0000001;
error_tol=1;
%x(1)=(B(1)-A(1,2)*x(2)-A(1,3)*x(3))/A(1,1);
while error_tol>tol 
for ii=1:n
xold=x;
	Ax=0;
	for jj=1:n
		if jj==ii
			%do nothing
		else
			Ax=Ax+A(ii,jj)*x(jj);
		end
	end
x(ii)=(B(ii)-Ax)/A(ii,ii);
end
x;
errorx=abs((x-xold)./x)*100;
error_tol=max(errorx);
end
x