GNU Octave 中的定积分计算存在问题

任务是这样的：在 GNU Octave 中计算函数的值（菲涅耳变换）以及理解如何在 Octave 中计算积分的主要问题。请不要提供 C++、Python 和其他编程语言和程序，我已经知道这些选项，它们不适合我。

我附上代码：

clear all;
lambda = 0.0005;
T = 1;
f = 100;
z = 150;
k = 2*pi/lambda;
alpha = 1.5;
T = 1;

x = linspace (0, 2*pi, 200);
u = exp(-i*x.^2/2*f + i*alpha*x.^3);

for j = 1:200
 x_u(j) = x(j)-u(j);
endfor;

integ1 = exp(-i*x.^2/2*f + i*alpha*x.^3 + i*k*(x_u).^2);


integ2 = quad(integ1, -T, T);
G = -i/(lambda*z)*exp(i*k*z)*integ1;

因此，第一个参数倒数第二行的错误，Octave 要求将其设为字符串或函数。我真的不明白如何正确输入该功能，以便他停止对我发誓。