function B = rowcombine(A,i,j,k)
    % ROWCOMBINE
    %           Makes a new matrix B by replacing row i  of A with 
    %            (row i of A) + k * (row j of A)
    B = A;
    B(i,:) = A(i,:) + k * A(j,:);