在 Matlab 中使用四边形时出错

Question

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

??? 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

与pi/8的解析值一致。