Question

同一微分方程组 不同的方法求解 结果为什么不一样

这两种方法求得同一个微分方程组的解不一样，请您看一下原因在哪，谢谢
1.
代码：
R=600;r=100;L=100;b=12.3;lam1=0.171;
vdp10001=@ (x,y)[-lam1*sin(y(1))*tan(y(1))/r;tan(y(1))/r];
[a,b] = ode45(vdp10001,[-L 0],[0.99*pi/2 0]);

2.
代码：
R=600;r=100;L=100;b=12.3;lam1=0.171;options=odeset('RelTol',1e-10,'AbsTol',1e-10);
m函数文件“wf1”:
function dy=wf1(x,y,options,R,r,lam1)
dy=[-lam1*sin(y(1))*tan(y(1))/r;tan(y(1))/r];
end
[m,n]=ode45('wf1',[-L 0],[pi*0.999/2 0],options,R,r,lam1);
说明一下：不是[0.99*pi/2 0]与[pi*0.999/2 0]的原因 我试过的 可以解决的请回答 解决不了的请不要说一些无关的 谢谢

Answer 1

求解 结果是一样的。可能没有清除上次程序执行变量，而导致变量窜用。

代码：

运行结果

追问

您好 问题很明显啊 运行结果 第一张曲线点相对稀疏，第二幅曲线点密集 我核对了数值解 发现有一点偏差 但不知道为啥