我正在运行 octave，我一直在尝试使用 ode45、ode54、ode23 等来整合方程`
当要积分的时间间隔增加时，执行时间非线性增加。从 t=0 积分到 t=T 需要 1.2 秒。10T 需要 12 秒，没问题！但是对于 100T，它需要 500 秒，而不是预期的 120 秒。这是为什么呢？我可以绕过它吗？
我必须将 1.2 秒的时间间隔积分大约 1000 倍。1200秒就好了。但4小时后我没有任何结果。
编辑：我在实际方程式中犯了几个错误，现在它是正确的。脚本如下所示：
U=1e-12; # Unit factor to make computation precise and fast.
B = 0.7 * 6.2 / 2 * 1e26*U^2;
w = 6*pi/800 * 1e17*U;
OPT = odeset ("InitialStep",1e-16/U,"MaxStep", 4e-16/U, ...
"RelTol", 1e-1, "AbsTol", 1e-1);
tSpan=[0 1e-12/U]; #150 ps => 15 deg from 0 deg.
init=[0 0];
tic
[t, g]=ode23( @(t,G)[G(2); B*cos(G(1))*sin(w*t)],tSpan,init,OPT);
t = t*U;
g=g(:,1);
toc
plot(t,g)