A=[3 2 5;7 1 6;4 9 8];
b=[10;20;30];
B=[A b]
disp(A);
disp(b);
disp(B);
for (j=1:2)
sem=B(j,j);
B(j,:)=B(j,:)/B(j,j);
for (i=j+1:3)
B(i,:)=B(i,:)-B(i,j)*B(j,:);
end
end
disp(B);
z=B(3,4)/B(3,3);
disp('Nilai Z ')
disp(z);
y=(B(2,4)-B(2,3)*z)/B(2,2);
disp('Nilai Y ')
disp(y);
x=(B(1,4)-B(1,2)*y-B(1,3)*z)/B(1,1);
disp('Nilai X ')
disp(x);