自己贴出答案吧！希望能帮到大家！
clear all; close all
F = @(x) 1/2*erf(x/sqrt(2));
xs = [0:0.05:1].'; % 变量x的采样序列
ys = [10000:5000:100000].'; % 变量y的采样序列
[XS, YS] = meshgrid(xs, ys);
ZS = 0*XS;
for g = 1:length(xs)
    x = xs(g);
    for  k = 1:length(ys)
        y = ys(k);
        xy = x*y;
        ZS(k, g) = 1.104*xy + 0.72*quad(F,0,xy/40);
    end
end
mesh(XS,YS,ZS);

不错！虽然程序可以改进的空间还很大。

那前辈能帮忙改改吗？呵呵

Matlab擅长的还是矩阵运算，for语句能不用就别用吧=。=

仅仅个人意见

那前辈对我的那个问题有什么高见啊？