function B = gauss(A,row)
% function gauss
% Creates zeros under the leading entry in row of matrix A 

epsilon = 10^(-14);

[n,m]=size(A);

Anonzero = abs(A)>epsilon*ones(n,m);

pivotcolumn=find(Anonzero(row,:),1,'first');

for elimrow=(row+1):1:n
   A(elimrow,:)=A(elimrow,:)-A(elimrow,pivotcolumn)/A(row,pivotcolumn)*A(row,:);
end

B=A;


 
