回复: 请高手帮我解一个数值积分，谢谢：）
仅作讨论，
x=[0:0.001:1];
y1=x.^(19081/10000).*besselk(16719/5000,1/562949953421312*7177179097769788095511012573184^(1/2).*x.^(1/2));
y2=(-pi^(1/2).*erf(10*500/707.*x)+pi^(1/2));
y3=y1.*y2
y=7.3810*y3;
m=trapz(x,y)
结果
m =
NaN
我用的是梯形积分，你可以画一下y=f(x)的图形，但积分后为何得到NaN，我就不知道了，或许梯形积分有什么限制我不知道吧
|