function rmat = naiveGE(A)
  numRows = size(A)(1);
  numCols = size(A)(2);
  ndims = numRows;
  rmat = A;

  for currIdx = 1:(ndims - 1)
    pivot = rmat(currIdx, currIdx);
    for ii = (currIdx+1):ndims
      rmat(ii, currIdx+1:numCols) = rmat(ii, currIdx+1:numCols) - ...
             rmat(ii,currIdx) / pivot * rmat(currIdx, currIdx+1:numCols);
      rmat(ii, currIdx) = 0;
    end
  end
end