在Matlab中使用quad时出错

Question

我在matlab中有一个匿名函数：f=@(x) x/(1+x^4)。我想用quad将它从0集成到1，但是由于某种原因，mpower抱怨它不能接受x^4，因为“输入必须是标量和方阵。”如果我用x*x*x*x替换x^4，我会得到"Error using ==> mtime内部矩阵尺寸必须一致“。下面是完整的错误输出。

??? Error using ==> mpower
Inputs must be a scalar and a square matrix.

Error in ==> intapprox>@(x)x/(1+x^4) at 4
f=@(x) x/(1+x^4);

Error in ==> quad at 77
y = f(x, varargin{:});

Error in ==> intapprox at 27
area = quad(f,0,1);

Answer 1

使用f=@(x) x./(1+x.^4);，即./和.^，因为你不想做矩阵除法和幂运算。然后我就得到了

>> area = quad(f,0,1)
area =

    0.3927

这与π/8的分析值是一致的。