已知4.2 x1+ 3.3 x2+ 3.3 x3<=26，其中x1，x2，x3都是非负整数. matlab 编程求出方程组所有的解,用矩阵M存储.M的每一行表示一个解，第1、2、3列分别存储x1,x2,x3的值. 并且矩阵的元素分别按照第1列从小到大排列，如果第1列元素相同，则按照第⒉列递增排列，如果第2列相同，则按照第3列递增排列.

可以使用MATLAB的intlinprog函数来求解该问题，代码如下：

c = [0 0 0]; % 目标函数系数
A = [4.2 3.3 3.3]; % 约束矩阵
b = 26; % 约束向量
lb = [0 0 0]; % 变量下界
ub = [Inf Inf Inf]; % 变量上界
intcon = [1 2 3]; % 整数变量
options = optimoptions('intlinprog','Display','off'); % 设置intlinprog函数的选项
[x,fval] = intlinprog(c,intcon,A,b,[],[],lb,ub,options); % 求解整数线性规划问题
num_solutions = fval + 1; % 解的数量

M = zeros(num_solutions,3); % 初始化解矩阵
M(1,:) = x; % 第一个解
for i = 2:num_solutions
    % 添加新的约束条件，要求新解的目标函数值大于等于当前解的目标函数值加1
    A = [A; c];
    b = [b; fval+1];
    [x,fval] = intlinprog(c,intcon,A,b,[],[],lb,ub,options); % 求解整数线性规划问题
    M(i,:) = x; % 添加新解到解矩阵中
end

% 按照要求排序解矩阵
[~,idx] = sortrows(M,[1 2 3]);
M = M(idx,:);

解释一下代码：

1. 定义目标函数系数c、约束矩阵A、约束向量b、变量下界lb、变量上界ub和整数变量intcon。
2. 使用intlinprog函数求解整数线性规划问题，得到一个最优解x和目标函数值fval。
3. 初始化解矩阵M，将第一个解x存储到M的第一行。
4. 循环求解新的解，添加新的约束条件，要求新解的目标函数值大于等于当前解的目标函数值加1。将新的解添加到解矩阵M中。
5. 使用sortrows函数按照要求排序解矩阵M。

注意，由于整数线性规划问题的求解是NP难问题，可能会存在计算时间较长的情况。如果您的问题规模较大，建议使用专业的求解器来求解。