function x=Gauss(A,b);
n=size(A); v=[1:n]';
for k=1:n-1
    [i j]=PosMax(A(k:n,k:n));
    i=i+k-1;
    j=j+k-1;
    [A b]=Swap_linie(A,k,i,b);
    [A v]=Swap_col(A,k,j,v);
    T=MatGauss(A(:,k),k);
    A=T*A
    b=T*b;
end
x=zeros(n,1);
y=Triug(A,b);
for i=1:n
    x(v(i))=y(i);
end