如何求这个数值积分，老是出错啊，不知道为什么？

[summer_xueer]

先建了一个M文件
function y= myfun(x)
y=exp(20./(5+1./sin(x).^2))/(5+sin(x).^2);
然后在命令窗口键入
Q=quad(@myfun,-pi/2,pi/2)
就会提示出错：??? Index exceeds matrix dimensions.
这是为什么啊？

[cfchenchen]

y=exp(20./(5+1./sin(x).^2))./(5+sin(x).^2);
你在/(5+sin(x).^2)前少一个点.上面是加上之后的.结果得:8.6554;
还有一个问题,就是在-pi/2与pi/2处,没有意义的,所以我们可以让它们分别加上一个与减去一个非常小的数,这样结果不变,也不会出现这样的情况了,比如说:可以从-pi/2+0.00000001到pi/2-0.00000001处积分,这样就行了,我也是新手噢!!

[summer_xueer]

首先非常感谢楼上的给予帮助，但是由于本人matlab学得太烂，按照你说的去再重新改过来，好像还是得不到答案，最后还是NaN

[cfchenchen]

我把那个东西完整地发给你看一下.
M文件中函数:
function y=myfun(x)
y=exp(20./(5+1./sin(x).^2))./(5+sin(x).^2);
然后在命令窗口键入
Q=quad(@myfun,-pi/2-0.000001,pi/2-0.000001)


结果如下:>> Q=quad(@myfun,-pi/2-0.000001,pi/2-0.000001)

Q =

8.6554
你再看一下,那儿出错了!

[summer_xueer]

按照你说的，我已经得出了结果，确实是需要注意添加那个.，现在我还有一个问题要请教，当积分表达式当中有另外一个变量（除了积分变量之外），而这个积分又如何求出来呢？是不是应该用一个for循环呢？可是我这么做的话会提示出错

[cfchenchen]

按你这样说的,好像不是求数值积分,好像是求符号积分的.对符号积分应该用int积分函数的!
例如:
syms x y
f=sqrt(exp(x)+1)*y;
F=int(f,x)
结果得:
F =

y*(2*(exp(x)+1)^(1/2)-2*atanh((exp(x)+1)^(1/2)))
这个是不定积分,如果要求定积分的话,可以这样改一下int函数
int(f,x,0,2)意思是对x积分是从0到2的.

另外对数值积而言,在表达式只能有一个变量了,只这样才能算出结果啊!
对于以上说的是符号积分就不一样了,它可以有多个变量的.