c=randi(10,3); b=[2; 5; 4]; %c=[2 1 3; 4 5 -7; -1 2 2]; SPRAWDZENIE %b=[4 ;7; -3]; SPRAWDZENIE AA=gaussian(c,b); z=length(b)+1; bb=AA(:,z); xxx=backward(AA,bb); function Ag = gaussian(A,b) Ag = [A b]; n = size(Ag,1); for k=1:n-1 for i = k+1:n l = Ag(i,k) / Ag(k,k); for j = k:n+1 Ag(i,j) = Ag(i,j) - l * Ag(k,j); end end end end function x = backward(U,c) n = size(U-1,1); x = zeros(n,1); for i = n:-1:1 s = 0; for j = i+1:n s = s + U(i,j)*x(j); end x(i) = (c(i) - s) / U(i,i); end end