按你的程序是有问题，如开平方出现虚数后，就无法再算下去。我改了一下，看能不能达到你的预想：
g=0.189^2/0.148-0.1;
x=0:pi/100:g;
m=(-0.189-(0.189^2-0.148*x).^0.5)*((-0.0074)*(-1));
n=(-0.189+(0.189^2-0.148*x).^0.5)*((-0.0074)*(-1));
ms=sort(m);
ns=sort(n);
xlen=length(x);
for k=1:xlen
c=linspace(m(k),n(k),xlen);
y=(((2*pi)^0.5)*6)^(-1)*exp(-((c-25).^2)*(72^(-1)));
z=trapz(c,y,2)
figure
plot(x,y)
end
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ号：24596024