倍频程中的弧长计算

Question

我正在尝试计算Octave中3D曲线的弧长。我找到了一个简短的脚本（https://www.mathworks.com/help/matlab/math/integration-to-find-arc-length.html）并将其修改为以下格式，并得到了“运算符*：不合格的参数（op1为33x1，op2为33x1）”错误。我迷路了。

代码：

# calculating the length of a 3D curve
t=0:0.1:0.25*pi;
r3=81.109;
m=5;
n=20;
f=@(t) sqrt(r3^2+(t*n.^t+m*t*n.^t*log(n))^2)
len=integral(f,0,0.25*pi)