我试图在matlab中使用ode45和ode23s来确定我的确切问题的收敛顺序。
为此,请使用
options = odeset('InitialStep',dt,'RelTol',Inf,'AbsTol',Inf,'MaxStep',dt,'Stats','on');
[~,y] = ode23s(@funode,[time(1):dt:time(end)],init,options);
控制步长。
这对我在ode23s上测试过的问题很有效,但有些原因不适用于ode45。我是幸运的是它在ode23s上工作了吗?我究竟做错了什么?或者为什么不在ode45上工作?