怎么根据以下条件写出matlab代码，进行求解最大值max z=0.7pqk-x-48250 q=-51.333+5042.222p k=-4.256e-8*x2+0.0004092x+1.019 q-5042.222p+51.333≤0 0≤x≤7000 20≤p≤60

根据题意，可以列出以下线性规划模型：

目标函数：max z = 0.7pqk - x - 48250

约束条件：

- q - 5042.222p + 51.333 ≤ 0
- 0 ≤ x ≤ 7000
- 20 ≤ p ≤ 60
- k = -4.256e-8 * x^2 + 0.0004092 * x + 1.019

可以使用MATLAB中的线性规划函数linprog进行求解。具体代码如下：

% 目标函数系数
f = [-1, 0, 0, -0.7];

% 不等式约束系数矩阵
A = [0, -5042.222, -1, 0;
     0, 0, 0, 1;
     0, 1, 0, 0;
     -4.256e-8, 0.0004092, 0, 0];

% 不等式约束右侧常数
b = [-51.333; 7000; 60; 1.019];

% 变量下界
lb = [0; 20; -inf; -inf];

% 变量上界
ub = [7000; 60; inf; inf];

% 使用线性规划函数求解
[x, z] = linprog(f, A, b, [], [], lb, ub);

% 输出结果
disp(['最大值为：', num2str(-z)]);
disp(['最优解为：x1=', num2str(x(1)), ', x2=', num2str(x(2)), ', x3=', num2str(x(3))]);

运行代码后，可以得到最大值为约为 29079.9，最优解为 x1=7000, x2=60, x3=0，即在约束条件下，当 x=7000，p=60 时，目标函数取得最大值约为 29079.9。