为什么我这个程序出来y值全是NAN？

Question

function dx= fun4(t,x) %建立函数 g=10; m=0.03; r=sqrt(x(1)^2+x(3)^2); r1=-x(1)/r; r2=-x(3)/r; dx=[x(2);g*m*r1/r^2;x(4);g*m*r2/r^2]; end x0=[0;0;0;0]; %调用函数 [t,y]=ode45(@fun4,[0,20],x0) plot(y(:,1),y(:,3))

Answer 1

初始值为全0的时候， r=sqrt(x(1)^2+x(3)^2)的取值为0了，r1=-x(1)/r; r2=-x(3)/r均为0/0，即为nan，nan与任何数运算返回值均为nan。解决办法是修改初始条件为非0值。

x0=[1;-1;2;-2]; %调用函数 
[t,y]=ode45(@fun4,[0,20],x0) 
plot(y(:,1),y(:,3))