对每行加个注释，解释更加详细一点

当然可以，下面是对整个代码添加了详细的行级注释的版本：

% 函数GaussXQLineMain接收两个输入参数，矩阵A和向量b，返回解向量x和简化版的系数矩阵XA
function [x,XA]=GaussXQLineMain(A,b)

    % 获取矩阵A的尺寸，n表示行数
    N = size(A);
    n = N(1);

    % 初始化一个变量用于保存当前行的主元索引
    index = 0;

    % 主循环遍历矩阵的第一列到最后一列
    for i=1:(n-1) % 从第一行开始到倒数第二行
        % 寻找当前行的最大绝对值元素作为主元
        main_elemant = max(abs(A(1:n,i))); 
        % 内部循环遍历当前列，找到最大绝对值元素对应的行索引
        for k=i:n
            if(abs(A(k,i))==main_elemant)
                % 找到主元所在行索引，跳出内部循环
                index = k; 
                break;
            end
        end
        
        % 临时存储当前行的数据，准备交换
        temp = A(i,1:n);
        
        % 将主元行与索引对应的行交换
        A(i,1:n) = A(index,1:n);
        A(index,1:n) = temp;
        
        % 交换向量b中的对应元素
        bb = b(index); 
        b(index)=b(i);
        % 将bb赋值给b(i)
        b(i) = bb; 

        % 检查主元是否为0，若为0则报错
        if(A(i,i)==0)
            disp('对角元素为0！');
            return;
        end

        % 消元过程，处理当前行以下的非主元行
        for j=(i+1):n
            % 线性消元，l=A(j,i)，m=A(i,i)
            l = A(j,i);
            m = A(i,i);
            % 将j行减去l乘以i行的倍数
            A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m;
            
            % 更新对应b的值
            b(j)=b(j)-l*b(i)/m; 
        end
    end
    
    % 回代阶段：从最后一行开始，逐步计算每个未知数x的值
    n = length(b); 
    x = zeros(n, 1); % 初始化解向量

    for i = n:-1:1 % 从最后一个元素到第一个元素
        % 计算第i个未知数x(i)的值
        x(i) = b(i) / A(i, i);
        
        % 更新剩余行的b，排除已经计算过的项
        b(1:i-1) = b(1:i-1) - A(1:i-1, i) * x(i);
    end 

    % 返回解向量x和简化后的系数矩阵XA（实际上就是原始矩阵A）
    XA = A;

end

这里对每行进行了注释，解释了代码的主要步骤和目的。