clear
clc
A=[5 -3 2; -3 8 4; 2 4 -9]
B=[10; 20; 9]
C=A
n=3
for k=1:1:n-1
for i=k+1:1:n
z=A(i,k)/A(k,k);
B(i)=B(i)-B(k)*z;
for j=k:1:n
A(i,j)=A(i,j)-A(k,j)*z;
end
end
end
A
B
for i=n:-1:1
X(i)=B(i);
if i~=n
for j=n:-1:i+1
X(i)=X(i)-A(i,j)*X(j);
end
end
X(i)=X(i)/A(i,i);
end
disp('Wyliczone wartości');
X
disp('Sprawdzenie z macierzą B poprzez podstawienie B=A*X');
C*X'