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
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}