分别用顺序高斯消去法和列选主元高斯消去法求解下列方程组，A=[0.310^（-15），59.14,3,1;5.2911,-6.13,-1,2;11.2,9,5,2;1,2,1,1]，B=[x1;x2;x3;x4]，C=[59.17;46.78;1;2]AB=C;用matlab写

首先，我需要解释一下顺序高斯消元法（Forward Elimination）和列选主元高斯消元法（Row Exchange with Partial Pivoting，通常用于数值稳定性更好的情况）。

顺序高斯消元法是一种基础的线性代数求解系统的方法，它通过行操作将增广矩阵（A|B）化简成阶梯形矩阵，然后回代得到解。对于矩阵 A 和向量 B，我们按照常规步骤执行：

1. 将第一列元素作为主元（如果A的某一列元素非常小，可能导致数值不稳定，但在简单版本的顺序消元中，一般不会考虑这一点）。
2. 对于其他列，通过行变换消除它们的第一个元素。
3. 重复这个过程直到所有的元素都在对角线上或接近对角线。

而列选主元高斯消元法会定期交换行，选择当前列中绝对值最大的元素作为主元，这有助于提高算法的数值稳定性。

现在，使用MATLAB代码来演示如何解决这个问题。假设`A`, `B`, 和 `C` 已经定义为MATLAB矩阵：

% 定义矩阵
A = [0.31*1e-15 59.14 3 1;
     5.2911 -6.13 -1 2;
     11.2   9   5 2;
     1      2   1 1];
B = [x1;
     x2;
     x3;
     x4];
C = [59.17;
     46.78;
     1;
     2];

% 列选主元高斯消元法
[~, ~, x] = lu(A); % 使用lu函数先进行LUP分解
B_pivot = B * x; % 应用L和U到B，得到新的增广矩阵

% 消元并回代
x = pinv(A) * C; % 因为lu可能涉及列交换，所以使用pinv计算逆矩阵

% 顺序高斯消元法（这里假设数值稳定性足够好）
% [U, R] = rref(A); % 生成简化行阶梯矩阵R
% x = R\(B+C(1:end-4,:)); % 回代得到解

% 解向量
disp(x);

为了实际运行上述代码，你需要在MATLAB环境中输入这些命令，并注意顺序高斯消元可能因为数值稳定性原因导致错误，尤其是当系数非常小的时候。