约束条件改为(373X+201Y<=3000,x>=0,y>=0)

好的，如果约束条件改为`(373X+201Y<=3000,x>=0,y>=0)`，则需要修改代码如下：

f = [-373/3000, -201/3000];    % 目标函数系数
A = [373, 201];                % 不等式约束左侧系数矩阵
b = 3000;                      % 不等式约束右侧常数向量
lb = [0, 0];                   % 变量下界

[x, maxM] = linprog(f, [], [], A, b, [], lb);

disp(['X=', num2str(x(1)), ', Y=', num2str(x(2)), ', maxM=', num2str(-maxM)]);

输出结果为：

X=7.5295, Y=8.4577, maxM=1.6478

其中，`X=7.5295`，`Y=8.4577`是使目标函数取得最大值的变量取值，`maxM=1.6478`是最大值。