Gauss-elimináció kiíratások nélkül

% Gauss-elimináció kiíratások nélkül. A test_script.m számára

function x = gauss_test(A,b)

n=length(A);
A=[A,b];
for k=1:n-1
  for i=k+1:n
    for j=k+1:n+1
      if ( A(k,k) == 0 )
    [maximum, ind] = max(abs(A(k:n,k)));
    if ( maximum == 0 )
      return;
    end
    akt = A(ind,1:n+1);
    A(ind,1:n+1) = A(k,1:n+1);
    A(k,1:n+1) = akt;
      end
        A(i,j) = A(i,j) -  ( (A(i,k) / A(k,k)) * A(k,j) );
    end
    A(i,k)=0;
  end
end
x=zeros(n,1);
x(n,1) = A(n,n+1) / A(n,n);
temp=0;
for k=n-1:-1:1
  for j=k+1:n
    temp = temp + A(k,j) * x(j,1);
  end
    x(k,1) = (A(k,n+1) - temp) / A(k,k);
    temp=0;
end