主题:急急急!!!求关于数值积分的问题!!!
f=exp(-(1.5^2+1)*x)/(1+0.61*S*exp(-(1.5^2)*x));做这个函数对x的数值积分,积分范围为0到1000(因为没有解析形式所以要数值形式),其中S是一个可变量,S从1:0.1:100;每取一个S就做一次积分,并绘制S与f的曲线关系,我的程序为:
y=[];
ss=0.1:0.1:10;x0=0,h=inf;
for s=ss
f=inline(['(exp(-(a^2+1)*x)./(1+f*s*exp((-a^2)*x)),num2str(s)'],'x')
y1=quadl(f,x0,h);y=[y,y1]
end
不知为什么总是做不出来,希望高手指教,谢谢!!!
y=[];
ss=0.1:0.1:10;x0=0,h=inf;
for s=ss
f=inline(['(exp(-(a^2+1)*x)./(1+f*s*exp((-a^2)*x)),num2str(s)'],'x')
y1=quadl(f,x0,h);y=[y,y1]
end
不知为什么总是做不出来,希望高手指教,谢谢!!!