在Octave

Question

在Octave 4.0.2中我定义了一个函数S，如下所示：

S = @(x) (Y(k)+((Y(k+1)-Y(k))/h(k)-(2*M(k)+M(k+1))*h(k)/6)*(x-X(k))+M(k)*(k-X(k))^2/2+(M(k+1)-M(k))*(x-X(k))^3/(6*h(k)))

当我调用它来评估区间[X（k），X（k + 1）]中的数字时，我得到了我期望的结果，但是当我尝试用命令绘制它时：

fplot(S, [X(k), X(k+1)]); hold on;

我收到错误＆＃34;错误：对于A ^ b，A必须是方阵。使用。^表示元素功率。＆＃34;

发生了什么事？

Answer 1

好吧，我发现当将函数传递给fplot时，它必须能够获取输入向量并返回输出向量。