公式4：Y=18.464-0.0858*z1+24.08*0.2+0.5*z2;
公式5：Q=0.3*z1*z2
我只能编写一个公式的响应曲面及等值线图，两个公式如何编程，谢谢大家！很急
clear;clc
z1=50:10:160;z2=0.5:0.25:3;
[z1,z2]=meshgrid(z1,z2);
Y=18.464-0.0858*z1+24.08*0.2+0.5*z2;
surf(z1,z2,Y);
figure;
[c h]=contour(z1,z2,Y);
clabel(c);