如图所示，F(x)就可以用erf（sqrt（2）*x）来表示了！这样我的最后程序。应该是：xs=(0:0.01:1);% 自变量采样序列
ys=0*xs; % 初始化因变量序列
F=@(t)(1/2*erf(sqrt(2)*t)); % 略去被积函数具体内容
for k=1:length(xs)
    x=xs(k);
    ys(k)=46160*x+quad(F,0,1000*x); % 我省略了所有系数，自己补上
end
plot(xs,ys); % 出图

不知道对否？还有我想把刻度弄小点，怎么改啊？0.1太大了！